plplot 0.0.0

data/examples/x25.rb

@@ -0,0 +1,114 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+
+#  Filling and clipping polygons.
+
+#--------------------------------------------------------------------------
+# Test program for filling polygons and proper clipping
+#--------------------------------------------------------------------------
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Initialize plplot
+
+plssub(3, 3)
+plinit
+
+xextreme = [
+  [-120.0,  120.0],
+  [-120.0,  120.0],
+  [-120.0,  120.0],
+  [ -80.0,   80.0],
+  [-220.0, -120.0],
+  [- 20.0,   20.0],
+  [ -20.0,   20.0],
+  [ -80.0,   80.0],
+  [  20.0,  120.0]
+]
+
+yextreme = [
+  [-120.0, 120.0],
+  [  20.0, 120.0],
+  [ -20.0, 120.0],
+  [ -20.0, 120.0],
+  [-120.0, 120.0],
+  [-120.0, 120.0],
+  [ -20.0,  20.0],
+  [ -80.0,  80.0],
+  [-120.0, 120.0],
+]
+
+
+2.times do |k|
+  4.times do |j|
+    x0 = []
+    y0 = []
+    if (j == 0)
+      # Polygon 1: a diamond
+      x0[0] = 0; y0[0] = -100
+      x0[1] = -100; y0[1] = 0
+      x0[2] = 0; y0[2] = 100
+      x0[3] = 100; y0[3] = 0
+    elsif (j == 1)
+      # Polygon 1: a diamond - reverse direction
+      x0[3] = 0; y0[3] = -100
+      x0[2] = -100; y0[2] = 0
+      x0[1] = 0; y0[1] = 100
+      x0[0] = 100; y0[0] = 0
+      npts  = 4
+    elsif (j == 2)
+      # Polygon 2: a square with punctures
+      x0[0] = -100; y0[0] = -100
+      x0[1] = -100; y0[1] = -80
+      x0[2] = 80; y0[2] = 0
+      x0[3] = -100; y0[3] = 80
+      x0[4] = -100; y0[4] = 100
+      x0[5] = -80; y0[5] = 100
+      x0[6] = 0; y0[6] = 80
+      x0[7] = 80; y0[7] = 100
+      x0[8] = 100; y0[8] = 100
+      x0[9] = 100; y0[9] = -100
+    elsif (j == 3)
+      # Polygon 2: a square with punctures - reversed direction
+      x0[9] = -100; y0[9] = -100
+      x0[8] = -100; y0[8] = -80
+      x0[7] = 80; y0[7] = 0
+      x0[6] = -100; y0[6] = 80
+      x0[5] = -100; y0[5] = 100
+      x0[4] = -80; y0[4] = 100
+      x0[3] = 0; y0[3] = 80
+      x0[2] = 80; y0[2] = 100
+      x0[1] = 100; y0[1] = 100
+      x0[0] = 100; y0[0] = -100
+    end
+
+    xextreme.each_index do |i|
+      pladv(0)
+      plvsta
+      plwind(xextreme[i][0], xextreme[i][1], yextreme[i][0], yextreme[i][1])
+
+      plcol0(2)
+      plbox("bc", 1.0, 0, "bcnv", 10.0, 0)
+      plcol0(1)
+      plpsty(0)
+      if (k == 0)
+        plfill(x0, y0)
+      else
+        plgradient(x0, y0, 45)
+      end
+      plcol0(2)
+      pllsty(1)
+      plline(x0, y0)
+    end
+  end
+end
+
+# Don't forget to call plend to finish off!
+
+plend

data/examples/x26.rb

@@ -0,0 +1,180 @@
+#!/usr/bin/env ruby
+# encoding: UTF-8
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include Math
+
+# Multi-lingual version of the first page of example 4.
+#
+# Copyright (C) 2006 Alan Irwin
+# Copyright (C) 2006 Andrew Ross
+#
+# Thanks to the following for providing translated strings for this example:
+# Valery Pipin (Russian)
+#
+# This file is part of PLplot.
+#
+# PLplot is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Library Public License as published
+# by the Free Software Foundation; either version 2 of the License, or
+# (at your option) any later version.
+#
+# PLplot is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU Library General Public License for more details.
+#
+# You should have received a copy of the GNU Library General Public License
+# along with PLplot; if not, write to the Free Software
+# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
+
+# This example designed just for devices (e.g., psttfc and the
+# cairo-related devices) that use the pango and fontconfig libraries. The
+# best choice of glyph is selected by fontconfig and automatically rendered
+# by pango in way that is sensitive to complex text layout (CTL) language
+# issues for each unicode character in this example. Of course, you must
+# have the appropriate TrueType fonts installed to have access to all the
+# required glyphs.
+#
+# Translation instructions: The strings to be translated are given by
+# x_label, y_label, alty_label, title_label, and line_label below.  The
+# encoding used must be UTF-8.
+#
+# The following strings to be translated involve some scientific/mathematical
+# jargon which is now discussed further to help translators.
+#
+# (1) dB is a decibel unit, see http://en.wikipedia.org/wiki/Decibel .
+# (2) degrees is an angular measure, see
+#    http://en.wikipedia.org/wiki/Degree_(angle) .
+# (3) low-pass filter is one that transmits (passes) low frequencies.
+# (4) pole is in the mathematical sense, see
+#    http://en.wikipedia.org/wiki/Pole_(complex_analysis) .  "Single Pole"
+#    means a particular mathematical transformation of the filter function has
+#    a single pole, see
+#    http://ccrma.stanford.edu/~jos/filters/Pole_Zero_Analysis_I.html .
+#    Furthermore, a single-pole filter must have an inverse square decline
+#    (or -20 db/decade). Since the filter plotted here does have that
+#    characteristic, it must by definition be a single-pole filter, see also
+#    http://www-k.ext.ti.com/SRVS/Data/ti/KnowledgeBases/analog/document/faqs/1p.htm
+# (5) decade represents a factor of 10, see
+#    http://en.wikipedia.org/wiki/Decade_(log_scale) .
+
+X_LABEL = [
+    "Frequency",
+    "Частота",
+]
+
+Y_LABEL = [
+    "Amplitude (dB)",
+    "Амплитуда (dB)",
+]
+
+ALTY_LABEL = [
+    "Phase shift (degrees)",
+    "Фазовый сдвиг (градусы)",
+]
+
+TITLE_LABEL = [
+    "Single Pole Low-Pass Filter",
+    "Однополюсный Низко-Частотный Фильтр",
+]
+
+LINE_LABEL = [
+    "-20 dB/decade",
+    "-20 dB/десяток",
+]
+
+# Makes transliteration from C example easier
+def pow(a,b)
+  a ** b
+end
+
+#--------------------------------------------------------------------------
+# plot1
+#
+# Log-linear plot.
+#--------------------------------------------------------------------------
+
+def plot1(type, x_label, y_label, alty_label, title_label, line_label)
+
+  pladv(0)
+
+  # Set up data for log plot
+
+  freql = NArray.float(101)
+  ampl  = NArray.float(101)
+  phase = NArray.float(101)
+
+  f0 = 1.0
+  101.times do |i|
+    freql[i] = -2.0 + i / 20.0
+    freq     = pow(10.0, freql[i])
+    ampl[i]  = 20.0 * log10(1.0 / sqrt(1.0 + pow(freq / f0, 2)))
+    phase[i] = -(180 / PI) * atan(freq / f0)
+  end
+
+  plvpor(0.15, 0.85, 0.1, 0.9)
+  plwind(-2.0, 3.0, -80.0, 0.0)
+
+  # Try different axis and labelling styles.
+
+  plcol0(1)
+  case (type)
+  when 0
+    plbox("bclnst", 0.0, 0, "bnstv", 0.0, 0)
+  when 1
+    plbox("bcfghlnst", 0.0, 0, "bcghnstv", 0.0, 0)
+  end
+
+  # Plot ampl vs freq
+
+  plcol0(2)
+  plline(freql, ampl)
+  plcol0(1)
+  plptex(1.6, -30.0, 1.0, -20.0, 0.5, line_label)
+
+  # Put labels on
+
+  plcol0(1)
+  plmtex("b", 3.2, 0.5, 0.5, x_label)
+  plmtex("t", 2.0, 0.5, 0.5, title_label)
+  plcol0(2)
+  plmtex("l", 5.0, 0.5, 0.5, y_label)
+
+  # For the gridless case, put phase vs freq on same plot
+
+  if (type == 0)
+    plcol0(1)
+    plwind(-2.0, 3.0, -100.0, 0.0)
+    plbox("", 0.0, 0, "cmstv", 30.0, 3)
+    plcol0(3)
+    plline(freql, phase)
+    plcol0(3)
+    plmtex("r", 5.0, 0.5, 0.5, alty_label)
+  end
+end
+
+#--------------------------------------------------------------------------
+# Illustration of logarithmic axes, and redefinition of window.
+#--------------------------------------------------------------------------
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Initialize plplot
+
+plinit()
+plfont(2)
+
+# Make log plots using two different styles.
+
+X_LABEL.each_index do |i|
+  plot1(i, X_LABEL[i], Y_LABEL[i], ALTY_LABEL[i], TITLE_LABEL[i],
+        LINE_LABEL[i])
+end
+
+plend

data/examples/x27.rb

@@ -0,0 +1,124 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include NMath
+
+# Drawing "spirograph" curves - epitrochoids, cycolids, roulettes
+#
+# Copyright (C) 2007  Arjen Markus
+#
+# This file is part of PLplot.
+#
+# PLplot is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Library Public License as published
+# by the Free Software Foundation; either version 2 of the License, or
+# (at your option) any later version.
+#
+# PLplot is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU Library General Public License for more details.
+#
+# You should have received a copy of the GNU Library General Public License
+# along with PLplot; if not, write to the Free Software
+# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
+
+def cycloid
+  # TODO
+end
+
+NPNT = 20000
+
+def spiro(params)
+
+    # Fill the coordinates
+
+    windings = params[3].to_i
+    steps    = NPNT / windings
+    dphi     = 8.0 * acos(-1.0) / steps
+
+    p01 = params[0] - params[1]
+    npnt = windings * steps + 1
+    phi  = NArray.float(npnt).indgen!.mul!(dphi)
+    phiw = NArray.float(npnt).indgen!.mul!(dphi * p01 / params[1])
+    xcoord = cos(phi).mul!(p01) + cos(phiw).mul!(params[2])
+    ycoord = sin(phi).mul!(p01) - sin(phiw).mul!(params[2])
+
+    xmin = xcoord.min
+    xmax = xcoord.max
+    ymin = ycoord.min
+    ymax = ycoord.max
+
+    if (xmax - xmin > ymax - ymin)
+        scale = xmax - xmin
+    else
+        scale = ymax - ymin
+    end
+    xmin = -0.65 * scale
+    xmax = 0.65 * scale
+    ymin = -0.65 * scale
+    ymax = 0.65 * scale
+
+    plwind(xmin, xmax, ymin, ymax)
+
+    plcol0(1)
+    plline(xcoord, ycoord)
+end
+
+#--------------------------------------------------------------------------
+# Generates two kinds of plots:
+#   - construction of a cycloid (animated)
+#   - series of epitrochoids and hypotrochoids
+#--------------------------------------------------------------------------
+
+# R, r, p, N
+params = [
+  [21.0,   7.0,  7.0,  3.0], # Deltoid
+  [21.0,   7.0, 10.0,  3.0],
+  [21.0,  -7.0, 10.0,  3.0],
+  [20.0,   3.0,  7.0, 20.0],
+  [20.0,   3.0, 10.0, 20.0],
+  [20.0,  -3.0, 10.0, 20.0],
+  [20.0,  13.0,  7.0, 20.0],
+  [20.0,  13.0, 20.0, 20.0],
+  [20.0, -13.0, 20.0, 20.0]
+]
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Initialize plplot
+
+plinit
+
+# Illustrate the construction of a cycloid
+
+cycloid
+
+# Loop over the various curves
+# First an overview, then all curves one by one
+
+plssub(3, 3) # Three by three window
+
+9.times do |i|
+  pladv(0)
+  plvpor(0.0, 1.0, 0.0, 1.0)
+  spiro(params[i])
+end
+
+pladv(0)
+plssub(1, 1) # One window per curve
+
+9.times do |i|
+  pladv(0)
+  plvpor(0.0, 1.0, 0.0, 1.0)
+  spiro(params[i])
+end
+
+# Don't forget to call plend to finish off!
+
+plend

data/examples/x28.rb

@@ -0,0 +1,369 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include Math
+
+# plmtex3, plptex3 demo.
+#
+# Copyright (C) 2007, 2008, 2009 Alan W. Irwin
+#
+# This file is part of PLplot.
+#
+# PLplot is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Library Public License as published
+# by the Free Software Foundation; either version 2 of the License, or
+# (at your option) any later version.
+#
+# PLplot is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU Library General Public License for more details.
+#
+# You should have received a copy of the GNU Library General Public License
+# along with PLplot; if not, write to the Free Software
+# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
+
+# Choose these values to correspond to tick marks. */
+XPTS        =  2
+YPTS        =  2
+NREVOLUTION = 16
+NROTATION   =  8
+NSHEAR      =  8
+
+#--------------------------------------------------------------------------
+# Demonstrates plotting text in 3D.
+#--------------------------------------------------------------------------
+
+
+xmin     = 0.0; xmax = 1.0; xmid = 0.5 * (xmax + xmin); xrange = xmax - xmin
+ymin     = 0.0; ymax = 1.0; ymid = 0.5 * (ymax + ymin); yrange = ymax - ymin
+zmin     = 0.0; zmax = 1.0; zmid = 0.5 * (zmax + zmin); zrange = zmax - zmin
+ysmin    = ymin + 0.1 * yrange
+ysmax    = ymax - 0.1 * yrange
+ysrange  = ysmax - ysmin
+dysrot   = ysrange / (NROTATION - 1)
+dysshear = ysrange / (NSHEAR - 1)
+zsmin    = zmin + 0.1 * zrange
+zsmax    = zmax - 0.1 * zrange
+zsrange  = zsmax - zsmin
+dzsrot   = zsrange / (NROTATION - 1)
+dzsshear = zsrange / (NSHEAR - 1)
+pstring   = "The future of our civilization depends on software freedom."
+
+# Allocate and define the minimal x, y, and z to insure 3D box
+x = NArray.float(XPTS).indgen!.mul!((xmax-xmin)/(XPTS-1)) + xmin
+y = NArray.float(YPTS).indgen!.mul!((ymax-ymin)/(YPTS-1)) + ymin
+z = NArray.float(YPTS, XPTS)
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+plinit
+
+# Page 1: Demonstrate inclination and shear capability pattern. */
+
+pladv(0)
+plvpor(-0.15, 1.15, -0.05, 1.05)
+plwind(-1.2, 1.2, -0.8, 1.5)
+plw3d(1.0, 1.0, 1.0, xmin, xmax, ymin, ymax, zmin, zmax,
+      20.0, 45.0)
+
+plcol0(2)
+plbox3("b", "", xmax - xmin, 0,
+       "b", "", ymax - ymin, 0,
+       "bcd", "", zmax - zmin, 0)
+
+# z = zmin
+plschr(0, 1)
+NREVOLUTION.times do |i|
+  omega         = 2 * PI * i / NREVOLUTION
+  sin_omega     = sin(omega)
+  cos_omega     = cos(omega)
+  x_inclination = 0.5 * xrange * cos_omega
+  y_inclination = 0.5 * yrange * sin_omega
+  z_inclination = 0.0
+  x_shear       = -0.5 * xrange * sin_omega
+  y_shear       = 0.5 * yrange * cos_omega
+  z_shear       = 0.0
+  plptex3(
+    xmid, ymid, zmin,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.0, "  revolution")
+end
+
+# x = xmax
+plschr(0, 1)
+NREVOLUTION.times do |i|
+  omega         = 2 * PI * i / NREVOLUTION
+  sin_omega     = sin(omega)
+  cos_omega     = cos(omega)
+  x_inclination = 0.0
+  y_inclination = -0.5 * yrange * cos_omega
+  z_inclination = 0.5 * zrange * sin_omega
+  x_shear       = 0.0
+  y_shear       = 0.5 * yrange * sin_omega
+  z_shear       = 0.5 * zrange * cos_omega
+  plptex3(
+    xmax, ymid, zmid,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.0, "  revolution")
+end
+
+# y = ymax
+plschr(0, 1)
+NREVOLUTION.times do |i|
+  omega         = 2 * PI * i / NREVOLUTION
+  sin_omega     = sin(omega)
+  cos_omega     = cos(omega)
+  x_inclination = 0.5 * xrange * cos_omega
+  y_inclination = 0.0
+  z_inclination = 0.5 * zrange * sin_omega
+  x_shear       = -0.5 * xrange * sin_omega
+  y_shear       = 0.0
+  z_shear       = 0.5 * zrange * cos_omega
+  plptex3(
+    xmid, ymax, zmid,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.0, "  revolution")
+end
+# Draw minimal 3D grid to finish defining the 3D box.
+plmesh(x, y, z, PL::DRAW_LINEXY)
+
+# Page 2: Demonstrate rotation of string around its axis.
+pladv(0)
+plvpor(-0.15, 1.15, -0.05, 1.05)
+plwind(-1.2, 1.2, -0.8, 1.5)
+plw3d(1.0, 1.0, 1.0, xmin, xmax, ymin, ymax, zmin, zmax,
+      20.0, 45.0)
+
+plcol0(2)
+plbox3("b", "", xmax - xmin, 0,
+       "b", "", ymax - ymin, 0,
+       "bcd", "", zmax - zmin, 0)
+
+# y = ymax
+plschr(0, 1)
+x_inclination = 1
+y_inclination = 0
+z_inclination = 0
+x_shear       = 0
+NROTATION.times do |i|
+  omega     = 2 * PI * i / NROTATION
+  sin_omega = sin(omega)
+  cos_omega = cos(omega)
+  y_shear   = 0.5 * yrange * sin_omega
+  z_shear   = 0.5 * zrange * cos_omega
+  zs        = zsmax - dzsrot * i
+  plptex3(
+    xmid, ymax, zs,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.5, "rotation for y = y#dmax#u")
+end
+
+# x = xmax
+plschr(0, 1)
+x_inclination = 0
+y_inclination = -1
+z_inclination = 0
+y_shear       = 0
+NROTATION.times do |i|
+  omega     = 2 * PI * i / NROTATION
+  sin_omega = sin(omega)
+  cos_omega = cos(omega)
+  x_shear   = 0.5 * xrange * sin_omega
+  z_shear   = 0.5 * zrange * cos_omega
+  zs        = zsmax - dzsrot * i
+  plptex3(
+    xmax, ymid, zs,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.5, "rotation for x = x#dmax#u")
+end
+
+# z = zmin
+plschr(0, 1)
+x_inclination = 1
+y_inclination = 0
+z_inclination = 0
+x_shear       = 0
+NROTATION.times do |i|
+  omega     = 2 * PI * i / NROTATION
+  sin_omega = sin(omega)
+  cos_omega = cos(omega)
+  y_shear   = 0.5 * yrange * cos_omega
+  z_shear   = 0.5 * zrange * sin_omega
+  ys        = ysmax - dysrot * i
+  plptex3(
+    xmid, ys, zmin,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.5, "rotation for z = z#dmin#u")
+end
+# Draw minimal 3D grid to finish defining the 3D box. */
+plmesh(x, y, z, PL::DRAW_LINEXY)
+
+# Page 3: Demonstrate shear of string along its axis.
+# Work around xcairo and pngcairo (but not pscairo) problems for
+# shear vector too close to axis of string. (N.B. no workaround
+# would be domega = 0.)
+domega = 0.05
+pladv(0)
+plvpor(-0.15, 1.15, -0.05, 1.05)
+plwind(-1.2, 1.2, -0.8, 1.5)
+plw3d(1.0, 1.0, 1.0, xmin, xmax, ymin, ymax, zmin, zmax,
+      20.0, 45.0)
+
+plcol0(2)
+plbox3("b", "", xmax - xmin, 0,
+       "b", "", ymax - ymin, 0,
+       "bcd", "", zmax - zmin, 0)
+
+# y = ymax
+plschr(0, 1)
+x_inclination = 1
+y_inclination = 0
+z_inclination = 0
+y_shear       = 0
+NSHEAR.times do |i|
+  omega     = domega + 2 * PI * i / NSHEAR
+  sin_omega = sin(omega)
+  cos_omega = cos(omega)
+  x_shear   = 0.5 * xrange * sin_omega
+  z_shear   = 0.5 * zrange * cos_omega
+  zs        = zsmax - dzsshear * i
+  plptex3(
+    xmid, ymax, zs,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.5, "shear for y = y#dmax#u")
+end
+
+# x = xmax
+plschr(0, 1)
+x_inclination = 0
+y_inclination = -1
+z_inclination = 0
+x_shear       = 0
+NSHEAR.times do |i|
+  omega     = domega + 2 * PI * i / NSHEAR
+  sin_omega = sin(omega)
+  cos_omega = cos(omega)
+  y_shear   = -0.5 * yrange * sin_omega
+  z_shear   = 0.5 * zrange * cos_omega
+  zs        = zsmax - dzsshear * i
+  plptex3(
+    xmax, ymid, zs,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.5, "shear for x = x#dmax#u")
+end
+
+# z = zmin
+plschr(0, 1)
+x_inclination = 1
+y_inclination = 0
+z_inclination = 0
+z_shear       = 0
+NSHEAR.times do |i|
+  omega     = domega + 2 * PI * i / NSHEAR
+  sin_omega = sin(omega)
+  cos_omega = cos(omega)
+  y_shear   = 0.5 * yrange * cos_omega
+  x_shear   = 0.5 * xrange * sin_omega
+  ys        = ysmax - dysshear * i
+  plptex3(
+    xmid, ys, zmin,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.5, "shear for z = z#dmin#u")
+end
+# Draw minimal 3D grid to finish defining the 3D box.
+plmesh(x, y, z, PL::DRAW_LINEXY)
+
+# Page 4: Demonstrate drawing a string on a 3D path.
+pladv(0)
+plvpor(-0.15, 1.15, -0.05, 1.05)
+plwind(-1.2, 1.2, -0.8, 1.5)
+plw3d(1.0, 1.0, 1.0, xmin, xmax, ymin, ymax, zmin, zmax,
+      40.0, -30.0)
+
+plcol0(2)
+plbox3("b", "", xmax - xmin, 0,
+       "b", "", ymax - ymin, 0,
+       "bcd", "", zmax - zmin, 0)
+
+plschr(0, 1.2)
+# domega controls the spacing between the various characters of the
+# string and also the maximum value of omega for the given number
+# of characters in *pstring.
+domega = 2 * PI / pstring.length
+omega  = 0
+# 3D function is a helix of the given radius and pitch
+radius = 0.5
+pitch  = 1 / (2 * PI)
+pstring.split('').each do |p1string|
+  sin_omega = sin(omega)
+  cos_omega = cos(omega)
+  xpos      = xmid + radius * sin_omega
+  ypos      = ymid - radius * cos_omega
+  zpos      = zmin + pitch * omega
+  # In general, the inclination is proportional to the derivative of
+  # the position wrt theta.
+  x_inclination = radius * cos_omega
+  y_inclination = radius * sin_omega
+  z_inclination = pitch
+  # The shear vector should be perpendicular to the 3D line with Z
+  # component maximized, but for low pitch a good approximation is
+  # a constant vector that is parallel to the Z axis.
+  x_shear   = 0
+  y_shear   = 0
+  z_shear   = 1
+  plptex3(
+    xpos, ypos, zpos,
+    x_inclination, y_inclination, z_inclination,
+    x_shear, y_shear, z_shear,
+    0.5, p1string)
+
+  omega += domega
+end
+# Draw minimal 3D grid to finish defining the 3D box.
+plmesh(x, y, z, PL::DRAW_LINEXY)
+
+# Page 5: Demonstrate plmtex3 axis labelling capability
+pladv(0)
+plvpor(-0.15, 1.15, -0.05, 1.05)
+plwind(-1.2, 1.2, -0.8, 1.5)
+plw3d(1.0, 1.0, 1.0, xmin, xmax, ymin, ymax, zmin, zmax,
+      20.0, 45.0)
+
+plcol0(2)
+plbox3("b", "", xmax - xmin, 0,
+       "b", "", ymax - ymin, 0,
+       "bcd", "", zmax - zmin, 0)
+
+plschr(0, 1)
+plmtex3("xp", 3.0, 0.5, 0.5, "Arbitrarily displaced")
+plmtex3("xp", 4.5, 0.5, 0.5, "primary X-axis label")
+plmtex3("xs", -2.5, 0.5, 0.5, "Arbitrarily displaced")
+plmtex3("xs", -1.0, 0.5, 0.5, "secondary X-axis label")
+plmtex3("yp", 3.0, 0.5, 0.5, "Arbitrarily displaced")
+plmtex3("yp", 4.5, 0.5, 0.5, "primary Y-axis label")
+plmtex3("ys", -2.5, 0.5, 0.5, "Arbitrarily displaced")
+plmtex3("ys", -1.0, 0.5, 0.5, "secondary Y-axis label")
+plmtex3("zp", 4.5, 0.5, 0.5, "Arbitrarily displaced")
+plmtex3("zp", 3.0, 0.5, 0.5, "primary Z-axis label")
+plmtex3("zs", -2.5, 0.5, 0.5, "Arbitrarily displaced")
+plmtex3("zs", -1.0, 0.5, 0.5, "secondary Z-axis label")
+# Draw minimal 3D grid to finish defining the 3D box.
+plmesh(x, y, z, PL::DRAW_LINEXY)
+
+# Clean up
+plend