[x,y,z] = cylinder(1:10);
[nx, ny, nz]=surfnorm(x,y,z);
surfnorm(x,y,z)
save('C:\STIMATH\Output\nx.dat','nx','-ascii')
save('C:\STIMATH\Output\ny.dat','ny','-ascii')
save('C:\STIMATH\Output\nz.dat','nz','-ascii')
