% Programme dipole.m
% Calcul du potentiel d un dipole
%
% Licence de Physique, UL5A
% Universite Paul Sabatier
% Licence de Physique UL5, TD numeriques, 2003-2004
% Version 08/02/2004 LK
%----------------------------------------
clear all;
clg
nom='??'

% constantes
%------------------------

C = 299792458;                  % Vitesse de la lumi�re
us4pe0 = C*C*1E-7               % 1 sur (4 PI epsilon0) = C^2 10^-7

% Donnees du dipole
%------------------------

q=5*1.6E-19                     % Charge en Coulombs
q1 = q                          % Charge 1
q2 =-q                          % Charge 2 opposee a q1

a=??                          % Distance entre les charges (en m)
borne = ??                    % limite du champ calcule

% Initialisation de la grille
%------------------------------

ny=??                          % Nombre de points en y
ymin= -borne
ymax= borne
ypas=(ymax-ymin)/(ny-1)         % Calcul du pas en y

nz=??                         % Nombre de points en z
zmin= -borne
zmax= borne
zpas=(zmax-zmin)/(nz-1)         % Calcul du pas en z

% Calcul du potentiel
%-------------------------

for ly=1:ny
	for lz=1:nz
		y(ly)=ymin+(ly-1)*ypas;
		z(lz)=zmin+(lz-1)*zpas;

		r1=??
		r2=??
		r= ??

		V  (lz,ly) = ??
		Vap(lz,ly) = ??
	end
end

%Partie graphique
%--------------------

v=[-8,-4, -2, -1,-0.5, -0.2, 0, 0.2, 0.5, 1, 2, 4, 8]     %Isocontours

hold on
Cex = contour (y,z,V,v);
Cap = contour (y,z,Vap,v);

title=(['Potentiel : a=', a, nom])
xlabel('y')
ylabel('z')
clabel(Cex,'manual')
clabel(Cap,'manual')
hold off

