%APARTADO 7
 ro=linspace(1,3,50); 
 tt=linspace(0,pi/2,50);             
 [Ro,Tt]=meshgrid(ro,tt);
 x1=Ro.*cos(Tt);
 y2=Ro.*sin(Tt);
 div=((Ro-1)./4);
 mesh(x1,y2,div)
 shading flat
 title('Divergencia del campo')
 xlabel('Eje x')
 ylabel('Eje y')
 zlabel('Eje z')
 

 
  
 
 
 