tpo élec

% 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='DOMENECH FOURCAULT';

% -------------- 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=  0.5*10^(-10);          % Decalage entre les pseudo sph�res
borne = 2*10^(-9);           % 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=500;
ypas=(ymax-ymin)/(ny-1);

zmin= -borne;
zmax=  borne;
nz=500;
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)= -q*us4pe0/(r1);
    elseif r1 <= r0
        V1(lz, ly)= -q*us4pe0*(r1^2)/(r0^3);
    end

    r2=sqrt(y(ly).^2 + (z(lz)+a/2).^2);
    if r2 > r0
        V2(lz, ly)= q*us4pe0/(r2);
    elseif r2 <= r0
        V2(lz, ly)= q*us4pe0*(r2^2)/(r0^3);
    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)= q*us4pe0*(1/r1 - 1/r2);
        elseif r<= r0
        Vap(lz,ly)= q*us4pe0*a*z(lz)/(r0^3);
    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,'--');
xlabel('y')
ylabel('z')
title=[nom]
clabel(C, 'manual');
hold off