 Rmax=6; hmax=2; n=4; hr=0.25*hmax/Rmax;
 r=([0:1:Rmax]);
 t=([0:pi/20:2*pi]);
 [th,rr]=meshgrid(t,r);
 [xx,yy]=pol2cart(th,rr);

 %   zz=hmax+hr*rr.*cos(n*th);   % Fisketuren
 %   zz=hmax/2+hr*rr.*cos(n*th)+xx*hr;
 for fi=0:-pi/12:-2*pi
   %zz=hmax+hr*rr.*cos(n*th + fi);   % Fisketuren
   zz=hmax/2+hr*rr.*cos(n*th+fi).*(cos(th+fi)+1)/2;+cos(th+fi)*3;
   meshz(xx,yy,zz);
   axis([ -Rmax Rmax -Rmax Rmax 0 hmax*1.5]);
   title('Virvelvindens golv')   
   pause (1)  % (1)
 end  
       

 view(az,el)
 for az=10:10:360;
   view(az,el);
   pause(0.5)
 end
 
 
