plot(pis3(1,:),pis3(2,:),'*r')
t=linspace(0,1,8);
sq0=spline(t,pis3);
tt=linspace(0,1);
ar=ppval(sq0,tt);
sx=ar(1,:);
sy=ar(2,:);
hold
plot(sx,sy)
hold
