%% Escoamento em uma canto
%
% Essa solução para escoamento potêncial representa um escoamento que
% acomparanha um canto com ângulo entre duas superfícies sólidas. Essa
% solução é resultado de uma formulação complexa.
%
% Aqui a parte imaginária do número complexo é a função corrente e a parte
% real a função potêncial. O gráfico representa a função corrente, em preto
% as linhas de corrente e na cor a intensidade do valor da função.

A=1;
beta=3*pi/4;
xlist=-1:0.04:1;
ylist=0:0.04:2;
[x,y]=meshgrid(xlist,ylist);
z=x+1i*y;
w=A*z.^(pi/beta);
contourf(xlist,ylist,imag(w),0:0.125:2.5);