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