rossler.f90

来自「因开发环境中没有fortran,故选C++代替」· F90 代码 · 共 64 行

!This program is to study(i) the properity(i) of rossler model
!The reference paper is   .....

program rossler

!!---------------------------------------------------------!
implicit none 
integer  N
integer, parameter:: NN=7
real x, y, z, dx, dy, dz, dt
Real a, b, c, x_c(NN),c1,c2
integer  i,j,k
! The values of parameters about the sy(i)s
character aa*1,  bb*1, cc*2
!---------------------------------------------------------!
a=0.2; b=0.2; N=1E6;  dt=0.01; !c=2.6 !,3.5, 4.1, 4.18, 4.21, 4.6, 5.7
x_c(1:7) =(/2.6, 3.5, 4.1, 4.18, 4.21, 4.6, 5.7/)  
!------------------------------------!
do j=1, 40000

  x=0.10
  y=0.30
  z=0.2
end do

do k=1,NN 
  	 c=x_c(k)
	 
     !c=c*100.0
	 c1=(c*100)/100
	 c2=c*100.0-100*int(c)

	write(aa, '(i1.1)')  int(k)	
    write(bb, '(i1.1)')  int(c1)
	write(cc, '(i2.2)')  int(c2*1.01)
	
   open (unit=1,file='Data\rossler_c'//aa//'='//bb//'.'//cc//'.dat',status='unknown') ! 

  do 	i=1,40000  	
        dx=dt*(  -y-z   )
        dy=dt*(  x+a*y  )
        dz=dt*(  b+(x-c)*z     )
        x=x+dx
        y=y+dy
        z=z+dz
      if(	 MOD( i, 40 )==0   ) then	  
	!	print*,   'i=', i, 'c=',  c, 'c2=', c2

         write (1,'(1x, 3f12.5)')   x, y, z
      end if
   
   	
    

   end do 
   
  close (1)
  

    
  end do

end program rossler