x = linspace(-10,10);
y = linspace(-5,5);
[X,Y] = meshgrid(x,y);
q=1;
d=2;
% fonction de courant uniforme + source ŕ l'origine
RA =Y+q.*atan2(Y,X) ;
% fonction de courant uniforme + source en x=-d + puits en x=d
RA2=Y+q.*(atan2(Y,(X+d))-atan2(Y,(X-d)))
% 
e=1
X0=2;
rc=2;
rc2=rc*rc;
rc3=1
% fonction de courant uniforme + dipôle ŕ l'origine
CY=-Y.*(1-rc2./(X.*X+Y.*Y))
figure
contour(X,Y,RA,20)
figure(2)
contour(X,Y,RA2,20)
% fonction de courant uniforme + dipôle ŕ l'origine + dipôle en X0
RA4=-Y.*(1-rc2./(X.*X+Y.*Y)-rc3./((X-X0).*(X-X0)+Y.*Y))
figure(3)
contour(X,Y,RA4,40)
figure(4)
contour(X,Y,CY,50)