2011-08-26から1日間の記事一覧

微分方程式の数値解析 #10

有名どころでLorenz方程式を解いてみる。 方程式のパラメータはLorenzが最初の論文で示したものを使用する。 program ode use runge_kutta_4 implicit none double precision, parameter :: a = 10, b = 28, c = 8d0/3 ! Lorenz's parameters double precisi…