% Draw a torus
[phi,theta] = meshgrid(linspace(0,2*pi,100));
x =  (cos(phi) + 3) .* cos(theta);
y =  (cos(phi) + 3) .* sin(theta);
z = sin(phi);
c = sin(3*theta); 
 
surf(x,y,z,c)
shading interp
title('this is not a doughnut');
