%% Perfil Parabólico
%
% Solução exata de escoamento potencial 2D para um arqueamento parabólico
% alpha deve ser em radianos
%
function[xy,gamma,Gamma,Cl,Cm0,dCp]=PerfilParabolico(U0,epsilon,alpha,c,N)

x  = linspace(0,c,N);
xy = [x' 4*epsilon*x'/c.*(1-x'/c)];

dCp = 4*alpha*sqrt((c-x)./x) + 32 *epsilon/c^2*sqrt(x*c - x.^2);
gamma = U0/2*dCp;

Gamma = U0*pi*c*(alpha+2*epsilon/c);

Cl = 2*pi*(alpha+2*epsilon/c);
Cm0 = -pi/2*(alpha+4*epsilon/c);