*** Version 1.2 *** Auteur : Philippe Roy

  program chaos

  implicit none
  integer j,n,k,p
  parameter (n=1000)
  real*4 x2(n),r,dx,rmax,dr,f,i
  parameter (rmax=5)
  external f

  open(11,file='r.res')

  call initial(x2,dr,n)

  do 10 i=1,rmax,dr
    do 20 j=2,n,1
      p=j-1
      x2(j)=f(i,x2,n,p)
      if (j.gt.800) then
        write(11,*) i,x2(j)
      endif
    20 continue
  10 continue

  close(11)
  end

  subroutine initial(x,dr,n)

  implicit none
  integer n
  real*4 x(n),dr

 30 print *,''
  print *,'Quelle est la valeur initiale de x (entre 0 et 1) ?'
  read(*,*) x(1)

  if ((x(1).lt.0) .or. (x(1).gt.1)) then
    print *,'La valeur doit etre comprise entre 0 et 1'
    goto 30
  else
    goto 40
  endif

 40 print *,'Quelle est la valeur du pas d''itÅÈration dr ?'
  read(*,*) dr

  end

  real*4 function f(r,x,n,p)

  implicit none
  integer n,i,p
  real*4 r,x(n)

  end

Retour à la section chaos

Retour à la page principale 1