next up previous contents
Next: F. MathematicaによるSIRモデルのシミュレーション Up: E..4 Runge-Kutta法のプログラム Previous: E..4.1 Julia 言語によるプログラム例 rk1.jl,

E..4.2 Julia 言語によるプログラム例 lotka-vortella.jl

多次元の例として、やや唐突だが、Lotka-Volterra の方程式

(25)   $\displaystyle \frac{\D x}{\D t}=ax-bxy,$
(26)   $\displaystyle \frac{\D y}{\D t}=-cy+dxy$

を解いてみよう。


\begin{breakitembox}[l] {Julia言語によるプログラム例 \texttt{lotka-vo... ...footnotesize \verbatimfile{ode_for_ms/lotka-volterra-gr.jl} } \end{breakitembox} 
% julia
julia> includej("lotka-volterra-gr.jl")
julia> gr()
julia> testlv()
julia> testlv(1,100)


桂田 祐史