Discussion :

[cocanouar]

Bonjour,
J'ai une équation itérative et je suis débutant dans la programmation et je cherche un aide pour résoudre cette équation ( dans les pièces jointes):
avec

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
k0=1;
a=1;
k=0.5;
p=0.1

avec Re c'est la partie réal et Im c'est la partie imaginaire

et Voilà le petit programme que j'ai essayé :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
k0=1;
a=1;
k=0.5;
p=0.1;
for j=1:50,
    Re(j)=(1+p/(1+(pi^2/16)*(k0*a)^4*(j/k)^4))/((1+p/(1+(pi^2/16)*(k0*a)^4*(j/k)^4))^2+p^2*((j/k)+1/((pi/4)*(k0*a)^2*(j/k)^2))^2);
    Im(j)=-p*((j/k)+1/((pi/4)*(k0*a)^2*(j/k)^2))/((1+1/(1+(pi^2/16)*(k0*a)^4*(j/k)^4))^2+p^2*((j/k)+1/((pi/4)*(k0*a)^2*(j/k)^2))^2); 
end
plot(Re,'r');
 
pause;
print -djpeg figure 1
end

Merci d'avance

Bonne année