Bezier Curve/Patch Subdivision

The function splits a k-th order Bezier curve or Bezier patch at user defined parameter(s).
Use:
CP = bezsubdiv(C, U) splits a k-th order Bezier curve with control points C at parameter U, and returns new control points CP.
[CP, CL, CR] = BEZSUBDIV(...) also returns k-by-3 matrices CL and CR for control points CL and CR for left and right segments.
CP = BEZSUBDIV(C,U,W) splits a Bezier patch with control points C at parameter U and W, and returns new control points CP.
[CP,CLU,CRU,CLW,CRW] = BEZSUBDIV(...) also returns ku-by-kw-3 matrices CLU, CRU, CLW and CRW for control points for top left, top right, bottom left and bottom right patches.