Gammatest

[lhs,rhs] =argn(0);
if rhs < 1 ;N=[10,20,50,100];end
x=linspace(0,1,100);
plot2d(cos(2*%pi*x)',sin(2*%pi*x)',1,"151"," ",[-2,-2,2,2]);
xtitle('Bezier curve approximating a circle')
icol=2
for n=N ;
t=sqrt(linspace(0,1,n));
p=[cos(2*%pi*t);sin(2*%pi*t)];
y=bezier(p,x);
plot2d(y(1,:)',y(2,:)',icol,"000");
icol=icol+1;
end;