function map=joukowsky(a,b,c)

% Maps circle of center "a", radius "b" in z plane
% under joukowsky mapping z1= z + c^2/z
%
% "a" can be complex
%
% examples:
%      a=0.1*exp(3*i*pi/4); b=1+abs(a); joukowsky(a,b,1)
%      a=0.1*exp(3*i*pi/4); b=abs(1-a); joukowsky(a,b,1)

% Fabian Waleffe 1998/12/03 for Math 321

t=[0:0.02:1]*2*pi;
z=a+b*exp(i*t); z1=z+c^2./z;
x=real(z);y=imag(z); x1=real(z1);y1=imag(z1);

plot(x,y,x1,y1,cos(t),sin(t),':'), axis image
