Ej1d

t=linspace(0,1);
Vx=[1 2 4 4.6];
Vy=[1 3 -1 1.5];
plot(Vx,Vy)
hold on
qx=[0];
qy=[0];
n=3;
for i=0:n
    qx=qx+Vx(i+1)*bernstein(n,i,t);
    qy=qy+Vy(i+1)*bernstein(n,i,t);
end
plot(qx,qy)
hold off