Go to the first, previous, next, last section, table of contents.


GNUPLOT 函数

この節では 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.