function jouk(s)

%  Disegna i profili di Joukowsky
%  Calcola e disegna la lunghezza della corda
%  il centro aerodinamico e l'angolo di portanza nulla 
%  by Nicola Morganti
%   
%  jouk(s)   
% 
%  s= Singolarità del profilo
% 
%  puo' essere anche complessa e in tal caso il profilo
%  non sara' piu' simmetrico rispetto all'asse reale


a=-(s+1)/2
th=[0:.001:2*pi];
z=exp(i.*th);
z1=z+((1+a).^2)./(z+a);
plot(z1,'r');
axis equal;
hold on;
c1=1+((1+a).^2)./(1+a)
c2=s+((1+a).^2)./(s+a)
plot(c1,'bo');
plot(c2,'bo');
lunghezza_corda=abs(c1-c2)
Angolo_di_portanza_nulla_in_gradi=atan((imag(c2)-imag(c1))/(real(c2)-real(c1)))*180/pi
line([real(c1),real(c2)],[imag(c1),imag(c2)]);
A=real(c2)+0.25*abs(c1-c2)*cos(atan(imag(c2)/(abs(c2)+1)));
B=imag(c1)/(real(c1)-1)*(A-1);
centro_aerodinamico=complex(A,B)
plot(centro_aerodinamico,'go')
hold off
figure(2)
axis equal;
x=[-8*pi/180:0.001:14*pi/180];
cl=8*pi/lunghezza_corda*sin(x-atan((imag(c2)-imag(c1))/(real(c2)-real(c1))));
plot(x,cl);
grid;