!Checking Logistic Map !Give the value of rmu,x0 and n here !where rmu is parameter value, x0 is initial value and n is number of iteration !Output is stored in fort.60 !After running this file run distrib.f using gfortran to get invariant density PROGRAM DISRIBUTION implicit real*8(a-h,o-z) pi=(4d0)*datan(1d0) rmu=4d0 x0=0.2222d0 n=50000 DO i1=1,n WRITE(60,*)i1,rmu,x0 !Output is stored in fort.60 WRITE(65,*)x0 !Only iterated points are stored here x0new=rmu*x0*(1d0-x0) x0=x0new END DO END