% Programme sphere : calcule
% le champ electrique cree par la matiere
% polarisee.

% Universite Paul Sabatier
% Licence de Physique UL5, TD numeriques, 2003-2004
% Version 08/02/2004 LK

% -------------- Initialisation ------------------
clear all; 
clf reset

nom='??'

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

r0= 0.5E-9             % Rayon de la sphere
a=  ??            % Decalage entre les pseudo sph�res
borne = ??           % limite (en m) de la zone de calcul

q = 5 * 1.6E-19        % Charge
q1=  q
q2= -q

% ------ Definition de la grille de travail ------
ymin= -borne
ymax=  borne
ny=??
ypas=(ymax-ymin)/(ny-1)

zmin= -borne
zmax=  borne
nz=??
zpas=(zmax-zmin)/(nz-1)

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

    % --------- calcul "spheres decalees" --------
	r1=sqrt(y(ly).^2 + (z(lz)-a/2).^2);
	if r1 > r0
		V1(lz,ly)= ??
	elseif r1 <= r0 
		V1(lz, ly)= ??
	end

	r2=sqrt(y(ly).^2 + (z(lz)+a/2).^2);
	if r2 > r0  
		V2(lz, ly)= ??
	elseif r2 <= r0
		V2(lz, ly)= ??
	end
	V(lz, ly) = V1(lz,ly) + V2(lz,ly);
  
    % --------- calcul "approche dipole" ---------
	r = sqrt(y(ly)^2+z(lz)^2);
	if r > r0
		Vap(lz,ly)= ??
		elseif r<= r0
		Vap(lz,ly)= ??
	end
end
end

% --------------- Partie graphique  --------------
hold on
mini= min(V)*100;
maxi= max(V)*100;
pas= abs((maxi-mini)/40);

v = mini:pas:maxi;

C=contour(y,z,V,v);
%C=contour(y,z,Vap,v);
title([nom])
xlabel('y')
ylabel('z')
clabel(C, 'manual');
hold off

