F=inline('(r.^5)* cos(5*theta) + (r.^3)*cos(3*theta)/2 +1','r','theta');
G=inline('(r.^5)* sin(5*theta) + (r.^3)*sin(3*theta)/2','r','theta');
r=0:2/100:2;
r=r'; 
theta=0:pi/100:pi;
FF=F(r,theta); 
GG=G(r,theta); 
mesh(theta,r,FF)
figure
mesh(theta,r,GG)
u=sqrt(FF.^2+GG.^2);
figure
mesh(theta,r,u)
%axis([0 pi 0 2 0 1]) 