pfys.f90

用fortran 语言编写的平方映射的程序

!This program is to study the bifurcation properity of pfys model
!The reference paper is  p63 .....Lutongxin
!! the equation is: x=miu*x*(1-x)
!! the results can be seen in p108,115	
program pfys
implicit none 
integer, parameter::  N=2e4
integer, parameter::  NN=60 !400 	  ！周期三窗口用NN=60；整个分岔用 NN=400
Real x,   miu

integer  i,   ii 
!character  aa*2, bb*2, cc*1
! The values of parameters about the sy(i)s	   

do ii=1, NN
   miu=3.82+0.001*ii	!观察周期三窗口(3.82--3.88)
  ! miu=0.0+0.01*ii		!观察整个分岔现象
   x=0.2	    

open (unit=1,file='Data\pfys_k1'//'.dat',status='unknown')	 


  
 do i=1,N
   
	    x=miu*x*(1-x)       
   
	 if(i>=500 .and. mod(i,5)==0) then	  
        write (1,'(1x, 4f12.5)')  miu, x
     end if

 end do  

 	

end do
   close (1)
end program pfys