h=0.1                  % paso de muestreo
 
u=1:h:2;               % valor u según el muestreo en el intervalo [1,2]
 
v=0:h:2*pi+h;            % valor v según el muestreo en el intervalo [0,2*pi]
 
[rho,theta]=meshgrid(u,v); % matrices de las componentes polares de la placa
 
 xx=rho.*cos(theta);        % parametrización
 
 yy=rho.*sin(theta);
 
subplot(1,2,1)
 
mesh(xx,yy,0*yy)
 
axis([-4,4,-4,4])

subplot(1,2,2)
 
fx=(sin((pi*theta)/2))./(30*rho);
 
fy=0*yy;
 
mesh(xx+fx,yy+fy,yy*0)
 
axis([-4,4,-4,4])
