t=linspace(-2,2,100);
[x,y]=meshgrid(t);
z=x.^2+y.^2;
surf2stl('taso.stl',x,y,z)