Discussion :

[feynman]

Bonjour,
j ai définit une fonction numérique qui vant 1 au dehors de [-3/4,3/4] et une fonction numérique dans cette interval , j ai ecrit le program, mais j ai obtenu que ma fonction vaut toujours 1!!!!!!!!!!!!!!!
Alors pouvez m aider et dire ou est ce truc génant?

voila mon code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
           program ess
         integer i,j
         parameter(N=16)
         real gama(N-1),h
         h=1/N
         do i=1,N-1
            gama(i)=gamma(i*h)
         enddo
         do j=1,N-1
            print*, j,gama(j)
         enddo
         end
 
          real function gamma(x)
          real x
          if(x.GT.-3/4.AND.x.LT.3/4) then
              gamma =(10*((x**2-0.75**2)**4)*cos(2*pi*x)+1)**2
          else
              gamma =1.
          end if
          return
          end

meercci

[DaZumba]

3/4 est une division entiere qui vaut 0. Ta condition ne risque pas d'etre verifiee... 3.0/4.0 par contre vaut 0.75.

[feynman]

Bah j ai remplacé 3/4 par 3.0/4.0 et j ai obtrenu :


1 4.00451803
2 4.00451803
3 4.00451803
4 4.00451803
5 4.00451803
6 4.00451803
7 4.00451803
8 4.00451803
9 4.00451803
10 4.00451803
11 4.00451803
12 4.00451803
13 4.00451803
14 4.00451803
15 4.00451803
ce qui est toujours bizare!

[feynman]

Sayez c est résolu , tout le truc était ce 1.0/n au lieu de 1/N :@

[feynman]

Mais quand meme je ne suis pas convaincu des résultats,
en comparant la fonction avec scilab je n obtient pas la meme chose
alors pourquoi?

[qwilldune]

Bonjour,

As tu résolu ton problème?

Loin de moi la prétention d'avoir LA solution. Cela fait une décennie que je n'ai pas fait de FORTRAN et je m'y remet pour le accelerer des codes Matlab. Ton post m'a intrigué et c'est typiquement un bon exercice pour se replonger dans le sujet.

J'ai d'abord joué sous Matlab pour avoir les bons résultats. J'ai ensuite revu le source pour obtenir les mêmes résultats.

Je pense que le soucis vennait de ce qui est en argument de "gamma". (h est toujours égal à zéro). me trompe je?

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
program ess
 
integer i
parameter(n=16)
real gama(n-1),h
 
h=1/real(n)
do i=1,n-1
      	gama(i)=gamma(i*h)
       	print*, i*h,gama(i)
enddo
end
 
function gamma(x)
real x,pi
pi=3.14159265359
 
if(x.GT.-0.75.AND.x.LT.0.75) then
       	gamma =(10*((x**2-0.75**2)**4)*cos(2*pi*x)+1)**2
else
      	gamma =1.
end if
return
end