read bezier:
#
#	Compute a full trajectory with independent x, y parametric functions.
#
x := [0, 2, 1, 0.5, 4, 3, 2.5, 3.5]:
y := [0, 3, 4, 1, 2, 3, 1, 4]:

Bzy := BezierApprox(y,t):
Bzx := BezierApprox(x,t):

with(plots):
plotBz := plot({[Bzx,Bzy,t=0..7]}, colour=YELLOW):
plotP  := plot({pairs(x,y)},style=point):
plotL  := plot({[pairs(x,y)]},style=line, colour=CYAN):
display({plotBz,plotP,plotL}, title=`Degree 7 Bezier space curve and control polygon`);