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