[windows251]

>> z=x.*exp(-x.^2-y.^2);
>> [u,v,w]=surfnorm(x,y,z);
>> surf(x,y,z); axis equal
>> mat2ans('curve.lgw',x,y,z);