clear all;

A=400; U=248.61;

r=[0:0.1:2]; t=[-2*pi:0.1:2*pi];

f=U*cos(t); g=-U*sin(t);

[r,t]=meshgrid(r,t);

for n=0:200

    f=f+A.*(1./(2.^(n+1))).*((2.*r-2.*n.*(cos(t)+0.01.*sin(t)))./(r.^2+(1.+0.01.^2).*n.^2-2.*n.*r.*(cos(t)+0.01.*sin(t))));

end;

for n=0:200

    g=g+A.*(1./(2.^(n+1))).*((-2*n.*(0.01.*cos(t)-sin(t)))./(r.^2+(1+0.01.^2).*n.^2-2.*n.*r.*(cos(t)+0.01.*sin(t))));

end;    

contour(r,t,f,0:2:100);

hold on;

contour(r,t,g,0:2:100);