This chapter describes interface functions for
GNUPLOT ox server ox_sm1_gnuplot
.
These interface functions are defined in the file gnuplot
.
The file `gnuplot' is at
`$(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
The function
gnuplot_heat
(dt,step) demonstrates our gnuplot interface.
It numerically solves the heat equation
The segment [0,1] is devided into Heat_N
segments.
If the celebrated Courant-Friedrichs-Levi number
dt*Heat_N
*Heat_N
is less than or equal to 0.5,
then the explicit scheme is numerically stable.
One can observe the instability by changing CFL number.
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.