この節では GNUPLOT の ox サーバ ox_sm1_gnuplot
とのインタフェース関数を解説する.
これらの関数はファイル `gnuplot' で定義されている.
gnuplot は
`$(OpenXM_HOME)/lib/asir-contrib/' にある.
nobuki@yama:~$ asir
This is Risa/Asir, Version 20020802 (Kobe Distribution).
Copyright (C) 1994-2000, all rights reserved, FUJITSU LABORATORIES LIMITED.
Copyright 2000,2001, Risa/Asir committers, http://www.openxm.org/.
GC 6.1(alpha5) copyright 2001, H-J. Boehm, A. J. Demers, Xerox, SGI, HP.
PARI 2.2.1(alpha), copyright (C) 2000,
C. Batut, K. Belabas, D. Bernardi, H. Cohen and M. Olivier.
OpenXM/Risa/Asir-Contrib(20020804), Copyright 2000-2002, OpenXM.org
help("keyword"); ox_help(0); ox_help("keyword"); ox_grep("keyword");
for help messages (unix version only).
[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
Go to the first, previous, next, last section, table of contents.