この節では GNUPLOT の ox サーバ ox_sm1_gnuplot
とのインタフェース関数を解説する.
これらの関数はファイル `gnuplot.rr' で定義されている.
gnuplot.rr は
`$(OpenXM_HOME)/lib/asir-contrib/' にある.
[255] gnuplot.start();
0
[257] gnuplot.gnuplot("plot sin(x**2);");
0
関数
gnuplot.heat(dt,step) はわれわれの GNUPLOT インタフェース関数のデモ
である.
この関数は熱伝導方程式
区間 [0,1] は Heat_N 個に分割される.
static 変数 Heat_N は
関数 gnuplot.set_heat_N で設定する.
有名な Courant-Friedrichs-Levi 数
dt*Heat_N*Heat_N が 0.5 以下であれば,
陽的差分スキームは安定である.
CFL を変えることにより, 不安定性が生じるのを観察できる.
gnuplot.set_heat_N(20); gnuplot.heat(0.001,30); (CFL number is 0.4) gnuplot.set_heat_N(20); gnuplot.heat(0.003,30); (CFL > 0.5 unstable)
Author of GNUPLOT: Thomas Williams, Colin Kelley.
http://www.gnuplot.info
Go to the first, previous, next, last section, table of contents.