fencha.f90

来自「Rossler分叉程序」· F90 代码 · 共 46 行

	program main

	use DFLIB

	implicit none

	integer Nmax,n,dummy
	real*8 c,x0,y0,z0,x1,y1,z1,deta_t,x_max,temp

	OPEN (UNIT=1, FILE ='date.dat')

	Nmax=10000000
	
	x_max=-1000

	deta_t=1.d-2
	
	do c=2.5,6.5,0.01
		x0=1.d0
    	y0=1.1d0
	    z0=1.2d0
	  do n=1,Nmax
 	    x1=x0-(y0+z0)*deta_t
	    y1=y0+(x0+(2.d-1)*y0)*deta_t
	    z1=z0+((2.d-1)+x0*z0-c*z0)*deta_t

	    x0=x1
	    y0=y1
	    z0=z1

	    if(n>100000)then
!		  write(*,*) x0
		  temp=(y0+z0)
		  if(dabs(temp)<1.d-4.and.x0>0)then
            write(1,*) c,x0
	        dummy=SETPIXELRGB(int2((70+100*c)),int2((600-40*x0)),130000)
		  endif 
	    endif
	  enddo
	  
	  x_max=-1000
	enddo
	
	end