この節では 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("plot sin(x**2);"); 0
関数
gnuplot_heat
(dt,step) はわれわれの GNUPLOT インタフェース関数のデモ
である.
この関数は熱伝導方程式
区間 [0,1] は Heat_N
個に分割される.
有名な Courant-Friedrichs-Levi 数
dt*Heat_N
*Heat_N
が 0.5 以下であれば,
陽的差分スキームは安定である.
CFL を変えることにより, 不安定性が生じるのを観察できる.
Heat_N=20; gnuplot_heat(0.001,30); (CFL number is 0.4) 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.