i=1:20;
a=cos(i/20*2*pi);
a=a'*a;
mesh(a)
pause
mattosurf(a,'test',[0 1 0 1])
!topo test
