la fonction gaussienne

bonjour

je veux écrire un programme en utilisant les fichiers pour tracer la fonction gaussienne (c la solution initiale de mon problème ) quand je plote la fonction je trouve pas la bonne allure Pièce jointe 116191par contre c quand j'écris la formule de la fonction f(x) = exp((-(x+l)**2)/(2.D0))/sqrt((2.D0)*pi) je trouve la bonne allure Pièce jointe 116190 je vois pas c quoi le problème.

Si quelqu'un pourrait m'aider je serai vraiment reconnaissante

merci d'avance!

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28

program solution implicit none   real, dimension(:),allocatable :: x integer :: i,n real::l,xdeb,xfin,dx,pi   open(unit=10, file='parametre') read(10,*) n,xdeb,xfin close(10)   allocate(x(n))   l=xfin-xdeb dx = l*1.D0/(n+1) pi = 4.D0*datan(1.D0) do i=1,n x(i) = (1.D0/sqrt(2.D0*pi))*exp((-((i*dx)-(l/2.D0))*((i*dx)-(l/2.D0)))/2.D0) end do   open(unit=10, file='solution_ini.dat') do i = 1,n write(10,*) x(i) end do close(10)   end program solution

---

et pour les parametres j'ai pris 100 0.D0 1.D0

Salut.
C'est parce que les abscisses sont différents. Je retrouve la même courbe avec :

Code:

---

n=100 ; xdeb=-10.d0 ; xfin = 10.d0

---

Vérifie si l'argument de l'exponentiel est correct (facteur pi manquant ?)

salut __dardanos__

j'ai mis pi = 4.D0*datan(1.D0) (ligne 16 du prog) mais pour quoi ça marche pas avec xdeb = 0 et xfin = 1?!!!

On peut réécrire la boucle :

Code:

---

1 2 3 4

do i = 1, n absc = i * dx - l / 2.D0 x(i) = 1.D0 / sqrt(2.D0*pi) * exp(-.5d0 * absc**2) end do

---

La forme générale d'une gaussienne étant :

Code:

---

1.D0 / sqrt(2.D0*pi) / sigma * exp(-.5d0 * (absc-mu)**2 / sigma**2)

---

Tu as pris une gaussienne avec un écart type sigma de 1. Elle est centrée, mu=0 (elle est maximale en zéro). Rappelle toi que l'écart type est lié à la largeur à mi-hauteur.
Avec xdeb=0 et xfin=1, x varie de -0,49 à +0,49. Comme sigma=1, on a toujours |x|<sigma/2. D'où la courbe qui varie peu.
Avec xdeb=-10 et xfin=+10, x varie de -9,80 à +9,80. x prend des valeurs telles que |x|>>sigma. On retrouve la forme en cloche présentée dans les manuels.

merci __dardanos__ :ccool: