f.f90

这是一个可以用来计算洛仑兹方程解的程序

subroutine F(t0,x0,y0,z0,n,D)
real*8 theta,r,b
theta=10
r=24.74
b=8/3
		 if(n.EQ.1.0) then
		    D=theta*(y0-x0)
		 end if
		 if(n.EQ.2.0) then
		    D=x0*(r-z0)-y0
		 end if
		 if(n.EQ.3.0) then
		    D=x0*y0-b*z0
		 end if
end