vx=[1,2,4,4.6]
vy=[1,3,-1,1.5]
plot(vx,vy)
hold on
t=linspace(0,1);
n= 3;
qx=[0];
qy=[0];
for i=0:n;
qx=qx+vx(i+1)*berstein(n,i,t);
qy=qy+vy(i+1)*berstein(n,i,t);
plot(qx,qy)
hold on
end
hold off
vx =
1.0000 2.0000 4.0000 4.6000
vy =
1.0000 3.0000 -1.0000 1.5000