Aide sur l'équation de Schrödinger via MatLab

Fermé
Braindisturb Messages postés 1 Date d'inscription vendredi 10 mai 2013 Statut Membre Dernière intervention 10 mai 2013 - Modifié par Braindisturb le 10/05/2013 à 12:38
Bonjour,

J'ai un petit soucis, j'ai à comprendre un code qui serait apparemment lié à l'équation de Schrödinger :

D = 0.0001e-9; % On choisi notre Delta à cette valeur car on cherche un rayon de Bhor à 0.0529 nm
r = 0:D:0.5e-9; % Voici l'échelle des abscisses limité à 0,5 car après cette valeur de r les courbe s'éloigne du zéro
m = 9.1094e-31; % La masse d'un électron en kg
che = 1.60218e-19; % La charge d'un électron en Coulomb
E=(-13.619:0.0001:-13.6185)*che; % Calcul de E, les courbes d'energies nécéssaire à la condition de normalisation
epsilon = 8.85e-12; % La constante Epsilon
h = 6.6261e-34; % La constante de Plank
hbar = h/(2*pi); % Calcul de hbar
k = 0; % on fixe k=O du repere O,i,k pour une future utilisation dans un boucle for

Bon jusque la rien de compliquer, des déclaration de variable, c'est ensuite que je ne comprend pas.. ce qu'il se passe dans les boucles "for"

for j=1:1:length(E) % on trace les valeur de E
R(1,j)=0;
R(2,j)=1;
for i=2:1:length(r)-1 % on trace les valeur de R
Vr=-((che)^2)/(4*pi*epsilon*r(i));
R(i+1,j)=R(i,j)*(2+((2*m*D^2)/((hbar)^2))*(Vr-E(j)))-R(i-1,j);
end;
if (R(length(E),j)>0) %on prend la valeur qui coupe a 0
k=k+1;
end
figure(1)
plot(r,R)
title(['Les différentes courbes de E']);

end





si quelqu'un pouvait m'éclairer merci beaucoup. : )