function ktrefftz(c,alpha)

% Calcolo di un profilo di Karman - Trefftz data
% la singolarità del punto c e l'angolo al bordo 
% di uscita alpha
% Calcola inoltre il centro aerodinamico
% la lunghezza della corda
% e l'angolo di portanza nulla
%
% Autore Nicola Morganti


k=(360-alpha)/180;
tetha=[0:0.001:2*pi];
Z=exp(i*tetha);
z=k*(c-1)./(((Z-1)./(Z-c)).^k-1);
plot(z,'r')
axis equal;
hold on;
c1=k*(c-1)./(((1-1)./(1-c)).^k-1);
lunghezza_corda=abs(c1)
Angolo_di_portanza_nulla_in_gradi=angle(c1)*180/pi
line([0,real(c1)],[0,imag(c1)]);
centro_aerodinamico=complex(0.25*abs(c1)*cos(angle(c1)),0.25*abs(c1)*sin(angle(c1)))
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-angle(c1));
plot(x,cl);
grid;