Fonction et matrice

Version imprimable

Fonction et matrice

Resalut !
Desolé pour ces questions mais ça commence à m'énerver un peu des trucs que je ne connais pas encore.
J'ai une matrice A et je veux écrire une fonction qui me donne comme réponse une matrice aussi,
j'ai écrit la fonction Ah mais à l'execution ça me donne :

Code:

1
2
3
4
5
6
7
8
    bru1.f: Dans la fonction `ah':
bru1.f:52:
           function Ah(s,MM,phi)
                    1
bru1.f:61: (suite) :
                 Ah(j1,j2)=cmplx(0.,0.)
                 2
déclaration invalide ou référence au symboe « ah » à (2) [initialement vu à (1)]

Alors pouvez-vous me donner une indication de ce qui ce passe ?
Voilà le début de ma fonction

Code:

1
2
3
4
5
6
7
8
9
10
11
            function Ah(s,MM,phi)
        integer MM,j1,j2,l1,l2
        double precision h,pi,s
        complex phi(-MM,MM-1)
        parameter(ii=(0.,1.))
        parameter(pi=3.141592653d0)
        h=2.*real(s)/real(MM)     
        do j1=-MM,MM-1
           do j2=-MM,MM-1
              Ah(j1,j2)=cmplx(0.,0.)
........

merci

Pourquoi pas une subroutine? Ca serait plus propre à mon avis.

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
        subroutine jesaispas(Ah,s,MM,phi)
        implicit none
        integer ii
        complex Ah(-MM,MM-1)
        integer MM,j1,j2,l1,l2
        double precision h,pi,s
        complex phi(-MM,MM-1)
        parameter(ii=(0.,1.))
        parameter(pi=3.141592653d0)
        h=2.*real(s)/real(MM)     
        do j1=-MM,MM-1
           do j2=-MM,MM-1
              Ah(j1,j2)=cmplx(0.,0.)
........

21/04/2008, 16h10
feynman

t as raison, mais c est une petite partie d un tres long programme , dont il ne me donne pas des bons resultats pour cela j essay de modifier un peu les choses pour capter l erreur

Et si tu declare le type de retour de ta fonction.

Code:

1
2
3
4
5
6
7
8
9
10
11
 
function Ah(s,MM,phi)
implicit none
   complex :: Ah(-MM:MM-1,-MM:MM-1)
   ...
   do j1=-MM,MM-1
      do j2=-MM,MM-1
         Ah(j1,j2)=cmplx(0.,0.)
      enddo
   enddo
...

PS : c'est beaucoup mieux avec un implicit none

22/04/2008, 21h23
feynman

Ah merde! c est une erreur des betes.
desole et merci:lol:

quand meme j obtient toujours le meme probleme :

Code:

1
2
3
4
5
6
7
8
bru1.f: Dans la fonction `ah':
bru1.f:52:
           function Ah(s,MM,phi)
                    1
bru1.f:54: (suite) :
           complex Ah(-MM:MM-1,-MM:MM-1)
                   2
déclaration invalide ou référence au symboe « ah » à (2) [initialement vu à (1)]