code: def Bezier(points, at):
      
      at = np.asarray(at)
      at_flat = at.ravel()
      N = len(points)
      curve = np.zeros((at_flat.shape[0], 2))
      for ii in range(N):
          curve += np.outer(Bernstein(N - 1, ii)(at_flat), points[ii])
      return curve.reshape(at.shape + (2,))
docs: Build Bézier curve from points.
      Deprecated. CatmulClark builds nicer splines