E..3.5 Julia言語によるプログラム例 euler1-gr.jl

上で紹介した gnuplot を使うという方法は、分りやすいが、 今一つ使いにくく感じる。 Plots というパッケージをフロントエンド、 GR をバックエンドにする、というのを勧めている人がいるので、 それでやってみよう。

(動作がイマイチ安定していないような気がする。 バージョンによって挙動が違ったり。)

% julia
julia> include("euler1-gr.jl")
f (generic function with 1 method)

julia> gr()
julia> testeuler()
(出力される数値は省略)
julia> testeuler(1,100)
(出力される数値は省略)
julia> control+D
%

最初の testeuler() を実行すると、 結果が出るまで少し待たされる。 これはコンパイルをしているためで、 それ以降の testeuler() 呼び出しはすぐに実行される。

図: $\frac{\D x}{\D t}=x$, $x(0)=1$ を $[0,1]$ で解く