?? rossler.f90
字號:
!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
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -