plot(pq2(1,:),pq2(2,:),'*r')
[px,py]=bezi(pq2);
t=linspace(0,1);
bx=polyval(px,t);
by=polyval(py,t);
hold
plot(bx,by)
hold
