rho=linspace(1,5,30);                             %Definimos rho y theta
th=linspace(0,2*pi,30);

[U,V]=meshgrid(rho,th);                           %Creamos la malla

hold on
X=U.*cos(V);                                      %Parametrización de la superficie
Y=U.*sin(V);

f=@(rho,th)(rho+(1./rho)).*cos(th) +sqrt(2).*log(rho);  %Definimos la función potencial 
Z=f(U,V);                                         %Aplicamos la función potencial
surf(X,Y,Z);                                      %Dibujamos la función
plot(1*cos(th),1*sin(th),'k','lineWidth',3);      %Representamos el obstáculo
 
view(2);  
axis([-5,5,-5,5]);
colorbar;
title ('Función potencial');
xlabel ('EJE X');
ylabel ('EJE Y');
axis equal
hold off