Avis sur la façon de programmer une équation

Bonjour,
J'ai besoin de votre avis sur la façon de programmer l'équation en attachement.

J'ai écrit le code mais je doute quand même de sa validité vu que je n'obtiens pas la bonne réponse.
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 do f=0,TT call matbar(2*MM,mu) call matwise(2*MM,TR,mu,D) call dft2d(2*MM,D,Dhat) call dft2d(2*MM,Gh,Ghat) call matwise(2*MM,Ghat,Dhat,DD) call idft2d(2*MM,DD,DDinvhat) do i=0,(2*MM-1) do j=0,(2*MM-1) mu(i,j)=III(i,j)-h**2*DDinvhat(i,j) enddo enddo enddo
Avec matbar qui est une subroutine qui transforme une matrice en son conjugué complexe, matwise pour le produit de Frobinus, dft2d pour la transformée de Fourier discrète, idft2d pour la transformée de Fourier discrète inverse.

Alors merci de me donner votre avis de la programmation de cette formule et si c'est possible vos idées ou vos suggestions.

Merci de nouveau