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