for a=0:112; 
   for b=0:60; 
      u=a/2; 
      w=b/2; 
      v=w/50-0.3; 
      X(a+1,b+1)=cos(u)+v*cos(u/2)*cos(u); 
      Y(a+1,b+1)=sin(u)+v*cos(u/2)*sin(u); 
      Z(a+1,b+1)=v*sin(u/2); 
   end 
end

surf(X,Y,Z) 
view(55, 50) 
shading interp 
light 
lighting phong 
material dull

set(gcf, 'color', 'w') 
colormap hsv