/usr/share/doc/libplplot11/examples/octave/x11c.m is in octave-plplot 5.9.9-2ubuntu2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 | ## Copyright (C) 1998, 1999, 2000 Joao Cardoso.
##
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by the
## Free Software Foundation; either version 2 of the License, or (at your
## option) any later version.
##
## This program 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
## General Public License for more details.
##
## This file is part of plplot_octave.
## It is based on the corresponding demo function of PLplot.
## Does a series of mesh plots for a given data set, with different
## viewing options in each plot.
1;
function cmap1_init()
i = [0; 1]; # left boundary, right boundary
h = [240; 0]; # blue -> green -> yellow -> red
l = [0.6; 0.6];
s = [0.8; 0.8];
plscmap1n(256);
plscmap1l(0, i, h, l, s, zeros(2,1));
endfunction
function ix11c
global DRAW_LINEXY MAG_COLOR BASE_CONT
XPTS = 35; ## Data points in x */
YPTS = 46; ## Datat points in y */
opt = [DRAW_LINEXY, DRAW_LINEXY ]';
alt = [33, 17]';
az = [24, 115]';
title = ["#frPLplot Example 11 - Alt=33, Az=24, Opt=3",
"#frPLplot Example 11 - Alt=17, Az=115, Opt=3"];
## Parse and process command line arguments */
## (void) plparseopts(&argc, argv, PL_PARSE_FULL);
## Initialize plplot */
plinit();
xx = 3.0*((0:XPTS-1) - fix(XPTS / 2)) / fix(XPTS / 2);
yy = 3.0*((0:YPTS-1) - fix(YPTS / 2)) / fix(YPTS / 2);
[x,y] = meshgrid (xx,yy);
z = 3 * (1-x).^2 .* exp(-(x.^2) - (y+1).^2) - ...
10 * (x/5 - x.^3 - y.^5) .* exp(-x.^2 - y.^2)- ...
1/3 * exp(-(x+1).^2 - y.^2);
nlevel = 10;
zmax = max(max(z));
zmin = min(min(z));
step = (zmax-zmin)/(nlevel+1);
clevel = linspace(zmin+step, zmax-step, nlevel)';
cmap1_init;
for k=1:2
for i=0:3
pladv(0);
plcol0(1);
plvpor(0, 1, 0, 0.9);
plwind(-1, 1, -1, 1.5);
plw3d(1, 1, 1.2, -3, 3, -3, 3, zmin, zmax, alt(k), az(k));
plbox3("bnstu", "x axis", 0, 0,
"bnstu", "y axis", 0, 0,
"bcdmnstuv", "z axis", 0, 4);
plcol0(2);
switch (i)
case 0
plmesh(xx', yy', z', opt(k));
case 1
plmesh(xx', yy', z', opt(k) + MAG_COLOR);
case 2
plot3d(xx', yy', z', opt(k) + MAG_COLOR, 1);
case 3
plmeshc(xx', yy', z', opt(k) + MAG_COLOR + BASE_CONT, clevel);
endswitch
plcol0(3);
plmtex("t", 1.0, 0.5, 0.5, deblank(title(k,:)));
endfor
endfor
plend1();
endfunction
ix11c
|