この節では GNUPLOT の ox サーバ ox_sm1_gnuplot
とのインタフェース関数を解説する.
これらの関数はファイル `gnuplot' で定義されている.
gnuplot は
`$(OpenXM_HOME)/lib/asir-contrib/' にある.
[nobuki@yama doc]$ asir
This is Asir, Version 991116.
Copyright (C) FUJITSU LABORATORIES LIMITED.
3 March 1994. All rights reserved.
[155] load("xm");
1
xm Version 19991109. ox_help(0); ox_help("keyword"); for help message
[255] gnuplot_start();
0
[257] gnuplot("plot sin(x**2);");
0
関数
heat(dt,step) はわれわれの GNUPLOT インタフェース関数のデモ
である.
この関数は熱伝導方程式
区間 [0,1] は Heat_N 個に分割される.
有名な Courant-Friedrichs-Levi 数
dt*Heat_N*Heat_N が 0.5 以下であれば,
陽的差分スキームは安定である.
CFL を変えることにより, 不安定性が生じるのを観察できる.
Heat_N=20; heat(0.001,30); (CFL number is 0.4) Heat_N=20; heat(0.003,30); (CFL > 0.5 unstable)
Author of GNUPLOT: Thomas Williams, Colin Kelley
Go to the first, previous, next, last section, table of contents.