ruby-pgplot 0.1.5

data/README

@@ -0,0 +1,100 @@
+      Ruby/PGPLOT ver 0.1.4 by Masahiro TANAKA  2012-01-04
+
+= Description
+
+  * Ruby interface to PGPLOT.
+    It provides PGPLOT functions as a Ruby module.
+
+= Requirement
+
+  * PGPLOT ver 5.2.x
+     (http://astro.caltech.edu/~tjp/pgplot/)
+  * Libraries which PGPLOT drivers use:
+    ** X11 library
+    ** PNG library
+    ** GrWin Graphics Library for MS-Windows.
+     (http://spdg1.sci.shizuoka.ac.jp/grwinlib/english/)
+  * Ruby ver 1.8.x/1.9.x
+  * NArray ver 0.5.x/0.6.x
+     (http://www.ruby-lang.org/en/raa-list.rhtml?name=NArray)
+
+= Installation
+
+  * install PGPLOT library
+
+    - FORTRAN77 compiler required for source-compile.
+    - C interface library (libcpgplot.a) is also required.
+
+  * install Ruby, NArray
+
+  * install Ruby/PGPLOT
+
+    - extract distribution.
+
+    - compile and install by:
+
+      ruby extconf.rb [options (see below)]
+      make
+      make site-install
+      (or make install)
+
+    - Options for `ruby extconf.rb':
+       --with-pgplot-include=path : path to cpgplot.h
+       --with-pgplot-lib=path	  : path to PGPLOT libraries.
+       --with-x11-dir=path	  : path to X11
+       --with-grwin : for using GrWin (cygwin/mingw).
+       --with-sunws : PGPLOT is compiled with Sun WorkShop f77.
+		      gfortran is assumed without this option.
+
+      example:
+       ruby extconf.rb --with-x11-dir=/usr/X11R6 \
+	--with-pgplot-include=/usr/local/pgplot \
+	--with-pgplot-lib=/usr/local/pgplot
+
+= Usage
+
+  * Environment variables for PGPLOT (csh):
+
+      setenv PGPLOT_DIR  /usr/local/lib/pgplot/ # pgxwin_server, rgb.txt
+      setenv PGPLOT_FONT ${PGPLOT_DIR}/grfont.dat
+      setenv PGPLOT_DEV  /xwin
+
+  * To load the Ruby/PGPLOT extension;
+
+      require "pgplot"
+
+  * Module name is `Pgplot'.  If you want to omit the module name;
+
+      include Pgplot
+
+  * Pgplot module method names are lower case of corresponding
+    PGPLOT functions.  To start plotting:
+
+      pgbeg('/xwin')
+
+    See PGPLOT manual for function usage.
+    Arguments are not always same as the FORTRAN PGPLOT.
+    "FuncUsage" file includes a list of functions with
+    different arguments.
+
+= Platforms tested
+
+  * ruby 1.9.2p180 (2011-02-18 revision 30909) [x86_64-linux]
+  * ruby 1.8.7 (2011-02-18 patchlevel 334) [x86_64-linux]
+  * NArray 0.6.0.1
+  * gcc/gfortran version 4.5.1 20100924 (Red Hat 4.5.1-4) (GCC)
+
+= License
+
+  This program is free software.
+  You can distribute/modify this program
+  under the same terms as Ruby itself.
+  NO WARRANTY.
+
+= Acknowledgment
+
+  M.T. thanks to Tim Pearson for developing the useful PGPLOT library
+  and giving me the permission to distribute the Ruby version of pgdemos.
+  M.T. also thanks to the developers of the Ruby language
+  for providing the wonderful language.
+

data/README.ja

@@ -0,0 +1,91 @@
+      Ruby/PGPLOT ver 0.1.4 by Masahiro TANAKA  2012-01-04
+
+= 概要
+
+  * PGPLOTをRubyから使うための拡張ライブラリ。
+
+= 必要なもの
+
+  * PGPLOT ver 5.2.x
+     (http://astro.caltech.edu/~tjp/pgplot/)
+  * PGPLOTのドライバが使うライブラリ (必要に応じて):
+    ** X11ライブラリ
+    ** PNGライブラリ
+    ** GrWinグラフィクスライブラリ for MS-Windows (Cygwinでのみ動作確認)
+     (http://spdg1.sci.shizuoka.ac.jp/grwinlib/)
+  * Ruby ver 1.8.x/1.9.x
+  * NArray ver 0.5.x/0.6.x
+     (http://www.ruby-lang.org/en/raa-list.rhtml?name=NArray)
+
+= インストール方法
+
+  * PGPLOT をインストール
+
+    - FORTRAN77が必要です (gfortranでもOK)。
+    - C言語インタフェースライブラリ(cpgplot)を作成する必要があります。
+
+  * Ruby, NArray をインストール
+
+  * Ruby/PGPLOT をインストール
+
+    - ソースを展開
+
+    - コンパイル・インストール:
+      ruby extconf.rb [オプション(下記参照)]
+      make
+      make site-install
+      (または make install)
+
+    - ruby extconf.rb のオプション:
+       --with-pgplot-include=path : cpgplot.hのディレクトリ。
+       --with-pgplot-lib=path	  : PGPLOTライブラリのディレクトリ。
+       --with-x11-dir=path	  : X11のpath
+       --with-grwin : for using GrWin (cygwin/mingw).
+       --with-sunws : PGPLOTをSun WorkShopのf77を使って
+		      コンパイルした場合。指定しない時はg77。
+     例:
+       ruby extconf.rb --with-x11-dir=/usr/X11R6 \
+	--with-pgplot-include=/usr/local/pgplot \
+	--with-pgplot-lib=/usr/local/pgplot
+
+= 使用方法
+
+  * PGPLOTの環境変数設定(csh)
+
+      setenv PGPLOT_DIR  /usr/local/lib/pgplot/ # pgxwin_server, rgb.txt
+      setenv PGPLOT_FONT ${PGPLOT_DIR}/grfont.dat
+      setenv PGPLOT_DEV  /xwin
+
+  * 拡張ライブラリをロードするには:
+
+      require "pgplot"
+
+  * モジュール名は`Pgplot'です。モジュール名を省略したいときは:
+
+      include Pgplot
+
+  * Pgplotモジュールのメソッドは、PGPLOTの関数と同じ名前(小文字)です。
+    プロットを開始するには:
+
+      pgbeg('/xwin')
+
+    詳細は PGPLOTのマニュアルを参照して下さい。
+    引数は必ずしもFORTRAN版と同じではありません。
+    引数が違う関数のリストは、"FuncUsage"というファイルにあります。
+
+= 動作確認
+
+  * ruby 1.9.2p180 (2011-02-18 revision 30909) [x86_64-linux]
+  * ruby 1.8.7 (2011-02-18 patchlevel 334) [x86_64-linux]
+  * NArray 0.6.0.1
+  * gcc/gfortran version 4.5.1 20100924 (Red Hat 4.5.1-4) (GCC)
+
+= 配布条件
+
+  Ruby本体と同じです。
+  無保証です。
+
+= 謝辞
+
+  PGPLOTを開発し、Ruby版pgdemoの配布許可を下さった Tim Pearson氏
+  に感謝します。Ruby開発に携わった皆様に感謝します。

data/cogen.rb

@@ -0,0 +1,157 @@
+# cogen.rb : Code generation script for Ruby/PGPLOT
+#
+#   Copyright (c) 2000,2001 Masahiro TANAKA <masa@ir.isas.ac.jp>
+#
+#   This program is free software.
+#   You can distribute/modify this program
+#   under the same terms as Ruby itself.
+#   NO WARRANTY.
+
+$pgfuncs = %w(
+
+pgend::
+pgbbuf::
+pgebuf::
+pgpage::
+pgpap:1,1:
+pgupdt::
+pgpanl:1,1:
+pgclos::
+
+pgbox:2,1,0,2,1,0:
+pgtbox:2,1,0,2,1,0:
+pgvsiz:1,1,1,1:
+pgvstd::
+pgwnad:1,1,1,1:
+pgsubp:0,0:
+
+pgwedg:2,1,1,1,1,2:
+
+# Draw Funcs
+pgdraw:1,1:
+pgmove:1,1:
+pgrect:1,1,1,1:
+pgarro:1,1,1,1:
+pgcirc:1,1,1:
+pgpt1:1,1,0:
+pgerr1:0,1,1,1,1:
+pglab:2,2,2:
+pgptxt:1,1,1,1,2:
+pgtext:1,1,2:
+pgmtxt:2,1,1,1,2:
+pgetxt::
+
+pgiden::
+pgldev::
+pgsave::
+pgunsa::
+pgeras::
+
+# Set Funcs
+pgsch:1:
+pgscf:0:
+pgsci:0:
+pgsfs:0:
+pgsls:0:
+pgslw:0:
+pgsclp:0:
+pgsitf:0:
+pgslct:0:
+pgstbg:0:
+pgscr:0,1,1,1:
+pgshls:0,1,1,1:
+pgsah:0,1,1:
+pgscrl:1,1:
+pgscir:0,0:
+pgscrn:0,2:0
+pgshs:1,1,1:
+pgsvp:1,1,1,1:
+pgswin:1,1,1,1:
+
+# Query Funcs
+pgqch::1
+pgqcf::0
+pgqci::0
+pgqfs::0
+pgqls::0
+pgqlw::0
+pgqclp::0
+pgqid::0
+pgqitf::0
+pgqndt::0
+pgqtbg::0
+pgqcr:0:1,1,1
+pgqvp:0:1,1,1,1
+pgqwin::1,1,1,1
+pgqcol::0,0
+pgqcir::0,0
+pgqpos::1,1
+pgqvsz:0:1,1,1,1
+
+).grep(/:.*:/).collect{|i| i.split(":",3)}
+
+
+def pgfuncgen(name, inp, out)
+  inp = inp.split(",").collect{|i| i.to_i}
+  out = out.split(",").collect{|i| i.to_i}
+  ninp = inp.size
+  nout = out.size
+  # int->0, float->1
+  val2 = ["NUM2INT","NUM2DBL","StringValuePtr"]
+  type = ["int","float",nil]
+  conv = ["INT2NUM","rb_float_new",nil]
+  # Initialize Array
+  prot = ["VALUE obj"]
+  pass = []
+  vars = []
+  retn = []
+  inp.each_with_index { |i,x|
+    prot << "VALUE arg#{x}"
+    pass << "#{val2[i]}(arg#{x})"
+  }
+  out.each_with_index { |i,x|
+    vars << "#{type[i]} var#{x};"
+    pass << "&var#{x}"
+    retn << "#{conv[i]}(var#{x})"
+  }
+  if nout==0 then
+    retn = "Qtrue";
+  elsif nout>1 then
+    retn = "rb_ary_new3(#{nout},"+retn.join(",")+")"
+  else
+    retn = retn.join("")
+  end
+
+  vars = vars.join("")
+  prot = prot.join(",")
+  pass = pass.join(",")
+
+  return "
+static VALUE
+  rb_pgplot_#{name}(#{prot})
+{
+  #{vars}
+  c#{name}(#{pass});
+  return #{retn};
+}
+"
+end
+
+
+def cogen_pgplot
+  fin  = open("rb_pgplot.c.in","r")
+  fout = open("rb_pgplot.c","w")
+  while l = fin.gets
+    if    /--- auto-generated funcs will be placed here ---/ =~ l
+      $pgfuncs.each{|x| fout.print pgfuncgen(*x)}
+    elsif /--- auto-generated defs will be placed here ---/  =~ l
+      $pgfuncs.each{|x|
+	n = x[1].split(",").size
+	fout.print "  rb_define_module_function(mPgplot,\"#{x[0]}\",rb_pgplot_#{x[0]},#{n});\n"}
+    else
+      fout.print
+    end
+  end
+end
+
+cogen_pgplot

data/demo/pgdemo1.rb

@@ -0,0 +1,414 @@
+require 'narray'
+require 'pgplot'
+include Pgplot
+include NMath
+PI    = Math::PI
+TWOPI = PI*2
+BLACK,WHITE,RED,GREEN,BLUE,CYAN,MAGENT,YELLOW = (0..7).to_a
+FULL,DASH,DOTDSH,DOTTED,FANCY = (1..5).to_a
+NORMAL,ROMAN,ITALIC,SCRIPT = (1..4).to_a
+SOLID,HOLLOW = (1..2).to_a
+
+# ====== Utility function ======
+
+def indgen arg
+  if arg.kind_of?(Range)
+    return NArray.sfloat(arg.size).indgen!(arg.first)
+  elsif arg.kind_of?(Numeric)
+    return NArray.sfloat(arg).indgen!
+  else
+    raise ArgumentError, "invalid argument"
+  end
+end
+
+def randomn n=1
+  rr = NArray.sfloat(n)
+  xx = NArray.sfloat(n)
+  idx= NArray.int(n).indgen!
+  i  = 0
+  while i<n
+    x = NArray.sfloat(n*4/3).random!(1) * 2 - 1
+    y = NArray.sfloat(n*4/3).random!(1) * 2 - 1
+    r = x**2 + y**2
+    idx = (r<1).where
+    siz = [n-i,idx.size-1].min
+    rr[i] = r[idx[0...siz]]
+    xx[i] = x[idx[0...siz]]
+    #printf "i=%d,siz=%d,idx.size=%d\n",i,siz,idx.size
+    i += siz
+  end
+  return xx * sqrt(-2*log(rr)/rr)
+end
+
+# ====== Demo function ======
+
+def pgex0
+  %w(
+    version
+    state
+    user
+    now
+    device
+    file
+    type
+    dev/type
+    hardcopy
+    terminal
+    cursor
+  ).each{|i| print " ",i," = ",pgqinf(i),"\n"}
+
+  x1,x2,y1,y2 = pgqvsz(1)
+  x = x2-x1
+  y = y2-y1
+  printf " Plot dimensions (x,y; inches): %9.2f, %9.2f
+                          (mm): %9.2f, %9.2f\n",x, y, x*25.4, y*25.4
+end
+
+def pgex1
+  pgenv 0, 10, 0, 20, 0, 1
+  pglab '(x)', '(y)', 'PGPLOT Example 1:  y = x\u2'
+
+  xs = [1.0,2.0,3.0,4.0,5.0]
+  ys = [1.0,4.0,9.0,16.0,25.0]
+  pgpt xs,ys,9
+
+  n = 60
+  xr = indgen(n)*0.1
+  yr = xr**2
+  pgline xr, yr
+end
+
+def pgex2
+  pgenv -2.0,10.0,-0.4,1.2, 0,1
+  pglab '(x)', 'sin(x)/x', 'PGPLOT Example 2:  Sinc Function'
+
+  xr = (indgen(100)-20.5)/6.0
+  yr = sin(xr)/xr
+  pgline xr, yr
+end
+
+def pgex3
+  pgenv 0.0,720.0,-2.0,2.0, 0,-2
+  pgsave
+  pgsci 14
+  pgbox 'G',30.0, 0,'G', 0.2, 0
+  pgsci 5
+  pgbox 'ABCTSN',90.0, 3,'ABCTSNV', 0.0,0
+  pgsci 3
+  pglab 'x (degrees)','f(x)','PGPLOT Example 3'
+
+  xr = indgen(360) * 2
+  arg= xr/180*PI
+  yr = sin(arg) + cos(arg*2)*0.5 + cos(arg*1.5+PI/3)*0.5
+
+  pgsci 6
+  pgsls 2
+  pgslw 3
+  pgline xr, yr
+  pgunsa
+end
+
+def pgex4
+  seed = -5678921
+  data = randomn(1000)
+  pgsave
+  pghist data, 31, -3.1..3.1, 0
+  data = randomn(200)*0.5 + 1 
+  pgsci 15
+  pghist data, 31, -3.1..3.1, 3
+  pgsci 0
+  pghist data, 31, -3.1..3.1, 1
+  pgsci 1
+  pgbox 'BST', 0.0, 0, ' ', 0.0, 0
+
+  pglab 'Variate', ' ','PGPLOT Example 4:  Histograms (Gaussian)'
+
+  x = indgen(620)*0.01 - 3.1
+  y = exp(-(x**2)*0.5) * (0.2*1000/sqrt(2*PI))
+
+  pgline x,y
+  pgunsa
+end
+
+def pgex5
+  np=15
+  freq = NArray[ 26.0, 38.0, 80.0, 160.0, 178.0, 318.0, 365.0, 408.0,
+                 750.0, 1400.0, 2695.0, 2700.0, 5000.0, 10695.0, 14900.0 ]
+  flux = NArray[ 38.0, 66.4, 89.0, 69.8, 55.9, 37.4, 46.8, 42.4, 27.0,
+		 15.8, 9.09, 9.17, 5.35, 2.56, 1.73 ]
+  err  = NArray[ 6.0, 6.0, 13.0, 9.1, 2.9, 1.4, 2.7, 3.0, 0.34, 0.8,
+		 0.2, 0.46, 0.15, 0.08, 0.01 ]
+  pgsave
+  pgsci CYAN
+  pgenv -2.0,2.0,-0.5,2.5,1,30
+  pglab 'Frequency, \gn (GHz)',
+        'Flux Density, S\d\gn\u (Jy)',
+        'PGPLOT Example 5:  Log-Log plot'
+  x = indgen(100)*0.03 + 1.3
+  xp = x-3
+  yp = -x*1.15 - exp(-x)*7.72 + 5.18
+  pgsci RED
+  pgline xp,yp
+
+  xp = log10(freq)-3.0
+  yp = log10(flux)
+  pgsci GREEN
+  pgpt xp,yp,17
+
+  yhi = log10(flux+2*err)
+  ylo = log10(flux-2*err)
+  pgerry xp,ylo,yhi
+  pgunsa
+end
+
+def pgex6
+  twopi = PI*2
+  npol = 6
+
+  n1 = [ 0, 3, 4, 5, 5, 6, 8 ]  # index-0 is dummy
+  n2 = [ 0, 1, 1, 1, 2, 1, 3 ]
+  lab = [ '', 'Fill style 1 (solid)',
+	  'Fill style 2 (outline)',
+          'Fill style 3 (hatched)',
+          'Fill style 4 (cross-hatched)']
+  pgbbuf
+  pgsave
+  pgpage
+  pgsvp 0.0, 1.0, 0.0, 1.0
+  pgwnad 0.0, 10.0, 0.0, 10.0
+  pgsci 1
+  pgmtxt 'T', -2.0, 0.5, 0.5, 
+    'PGPLOT fill area: routines PGPOLY, PGCIRC, PGRECT'
+
+  for k in 1..4
+    pgsci 1
+    y0 = 10.0 - 2.0*k
+    pgtext 0.2, y0+0.6, lab[k]
+    pgsfs k
+    for i in 1..npol
+      pgsci i
+      angle = indgen(n1[i])*(n2[i]*twopi/n1[i])
+      x = cos(angle)*0.5 + i
+      y = sin(angle)*0.5 + y0
+      pgpoly x,y
+    end
+    pgsci 7
+    pgcirc 7.0, y0, 0.5
+    pgsci 8
+    pgrect 7.8, 9.5, y0-0.5, y0+0.5
+  end
+  pgunsa
+  pgebuf
+end
+
+def pgex7
+  pgbbuf
+  pgsave
+  pgsci 1
+  pgenv 0.0,5.0, -0.3,0.6, 0,1
+  pglab '\fix', '\fiy', 'PGPLOT Example 7: scatter plot'
+
+  xs = NArray.sfloat(300).random!(1.0)*5
+  ys = randomn(300)*0.05 + xs*exp(-xs)
+
+  pgsci 3
+  pgpt xs[0..99],ys[0..99], 3
+  pgpt xs[100..199],ys[100..199], 17
+  pgpt xs[200..299],ys[200..299], 21
+
+  xr = indgen(101)*0.05
+  yr = xr*exp(-xr)
+  pgsci 2
+  pgline xr,yr
+
+  xp = xs[100]
+  yp = ys[100]
+  xsig = 0.2
+  ysig = 0.1
+  pgsci 5
+  pgsch 3.0
+  pgerr1 5, xp, yp, xsig, 1.0
+  pgerr1 6, xp, yp, ysig, 1.0
+  pgpt1 xp,yp,21
+  pgunsa
+  pgebuf
+end
+
+def pgex8
+  pgpage
+  pgbbuf
+  pgsave
+  pgsvp(0.1,0.6,0.1,0.6)
+  pgswin(0.0, 630.0, -2.0, 2.0)
+  pgsci(CYAN)
+  pgbox('ABCTS', 90.0, 3, 'ABCTSV', 0.0, 0)
+  pgsci(RED)
+  pgbox('N',90.0, 3, 'VN', 0.0, 0)
+  xr = indgen(360)*2
+  yr = sin(xr/57.29577951)
+  pgsci(MAGENT)
+  pgsls(DASH)
+  pgline(xr,yr)
+  pgswin(90.0, 720.0, -2.0, 2.0)
+  pgsci(YELLOW)
+  pgsls(DOTTED)
+  pgline(xr,yr)
+  pgsls(FULL)
+  pgsvp(0.45,0.85,0.45,0.85)
+  pgswin(0.0, 180.0, -2.0, 2.0)
+  pgsci(0)
+  pgrect(0.0, 180.0, -2.0, 2.0)
+  pgsci(BLUE)
+  pgbox('ABCTSM', 60.0, 3, 'VABCTSM', 1.0, 2)
+  pgsci(WHITE)
+  pgsls(DASH)
+  pgline(xr,yr)
+  pgunsa
+  pgebuf
+end
+
+def pgex9
+  pgbbuf
+  pgsave
+  pgsci(5)
+  #pgfunt(fx,fy,360,0.0,TWOPI,0)  # pgfunt is not implemented
+  t = indgen(361)/360*TWOPI
+  x = sin(t*5)
+  y = sin(t*4)
+  pgenv -1,1,-1,1
+  pgline x,y
+  pgsci(3)
+  pglab('x','y','PGPLOT Example 9:  routine PGFUNT')
+  pgunsa
+  pgebuf
+end
+
+def bessel_j0 arg
+
+  r = NArray.sfloat(arg.size)
+  x = arg.abs
+  idx1,idx2 = (x<=3).where2
+  if idx1.size>0 then
+    xo3 = x[idx1]/3.0
+    t   = xo3**2
+    r[idx1] = 1.0 +
+	t*(-2.2499997 +
+	t*( 1.2656208 +
+	t*(-0.3163866 +
+	t*( 0.0444479 +
+	t*(-0.0039444 +
+	t*( 0.0002100))))))
+  end
+  if idx2.size>0 then
+    xx = x[idx2]
+    t  = 3.0/xx
+    f0 = 0.79788456 +
+	t*(-0.00000077 + 
+	t*(-0.00552740 +
+	t*(-0.00009512 +
+	t*( 0.00137237 +
+	t*(-0.00072805 +
+	t*( 0.00014476))))))
+    theta0 = xx - 0.78539816 +
+      	t*(-0.04166397 +
+      	t*(-0.00003954 +
+      	t*( 0.00262573 +
+      	t*(-0.00054125 +
+      	t*(-0.00029333 +
+      	t*( 0.00013558))))))
+    r[idx2] = f0*cos(theta0)/sqrt(xx)
+  end
+  return r
+end
+
+def bessel_j1 arg
+  r = NArray.sfloat(arg.size)
+  x = arg.abs
+  idx1,idx2 = (x<=3).where2
+  if idx1.size>0 then
+    xo3 = x[idx1]/3.0
+    t = xo3**2
+    f = 0.5 + t*(-0.56249985 +
+	t*( 0.21093573 +
+	t*(-0.03954289 +
+	t*( 0.00443319 +
+	t*(-0.00031761 +
+	t*( 0.00001109))))))
+    r[idx1] = f * arg[idx1]
+  end
+  if idx2.size>0 then
+    xx = x[idx2]
+    t  = 3.0/xx
+    f1 = 0.79788456 +
+        t*( 0.00000156 +
+	t*( 0.01659667 + 
+        t*( 0.00017105 +
+        t*(-0.00249511 +
+        t*( 0.00113653 + 
+        t*(-0.00020033))))))
+    theta1 = xx - 2.35619449 + 
+        t*( 0.12499612 +
+        t*( 0.00005650 +
+        t*(-0.00637879 +
+        t*( 0.00074348 +
+        t*( 0.00079824 +
+        t*(-0.00029166))))))
+    r[idx2] = f1*cos(theta1)/sqrt(xx)
+  end
+  idx = (arg<0).where
+  #p idx
+  #p r[idx]
+  r[idx] = -r[idx] if idx.size>0
+  return r
+end
+
+def pgex10
+  pgbbuf
+  pgsave
+  pgsci(YELLOW)
+  # PGFUNX(PGBSJ0,500,0.0,10.0*PI,0)
+  x = indgen(500)/50*PI
+  y = bessel_j0(x)
+  pgenv 0,PI*10, y.min,y.max
+  pgline x,y
+  pgsci(RED)
+  pgsls(DASH)
+  # PGFUNX(PGBSJ1,500,0.0,10.0*PI,1)
+  pgline x, bessel_j1(x)
+  pgsci(GREEN)
+  pgsls(FULL)
+  pglab('\fix', '\fiy',
+        '\frPGPLOT Example 10: routine PGFUNX')
+  pgmtxt('T', -4.0, 0.5, 0.5,
+	 '\frBessel Functions')
+  pgarro(8.0, 0.7, 1.0, bessel_j0(NArray[1.0])[0])
+  pgarro(12.0, 0.5, 9.0, bessel_j1(NArray[9.0])[0])
+  pgstbg(GREEN)
+  pgsci(0)
+  pgptxt(8.0, 0.7, 0.0, 0.0, ' \fiy = J\d0\u(x)')
+  pgptxt(12.0, 0.5, 0.0, 0.0, ' \fiy = J\d1\u(x)')
+  pgunsa
+  pgebuf
+end
+
+# ====== Demo start ======
+
+raise "device not found" if pgopen<0
+
+pgex0
+pgex1
+pgex2
+pgex3
+pgsubp 2,1
+pgex4
+pgex5
+pgsubp 1,1
+pgex6
+pgex7
+pgex8
+pgex9
+pgex10
+
+pgclos
+exit
+