Discussion :

[the_reward]

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 par 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 je vois pas c quoi le problème.

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

merci d'avance!

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
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