/usr/share/doc/libplplot11/examples/octave/p21.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 | ## Copyright (C) 2003 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.
function p21
global GRID_CSA GRID_DTLI GRID_NNI GRID_NNIDW GRID_NNLI GRID_NNAIDW
##ocmap = colormap;
colormap(bgr);
global pl_automatic_replot
t = pl_automatic_replot;
pl_automatic_replot = 0;
gx = 30;
gy = 40;
n = 600;
alg = ["Cubic Spline Approximation",
"Delaunay Linear Interpolation",
"Natural Neighbors Interpolation",
"KNN Inv. Distance Weighted",
"3NN Linear Interpolation",
"4NN Around Inv. Dist. Weighted"];
opt = [0, 0, -1, 50, 1.001, 0];
x = rand(n,1)*6-3;
y = rand(n,1)*6-3;
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);
xg = linspace (-2.5, 3, gx)';
yg = linspace (-3, 2.5, gy)';
multiplot(3,2);
j = 1;
for i = [GRID_CSA GRID_DTLI GRID_NNI GRID_NNIDW GRID_NNLI GRID_NNAIDW]
pladv(j++);
zg = griddata(x, y, z, xg, yg, i, opt(i));
## zg(isnan(zg)) = 0;
old_dofi = warning("query","Octave:fortran-indexing");
warning("off","Octave:fortran-indexing");
zg(isnan(zg)) = 0;
warning(old_dofi.state,"Octave:fortran-indexing");
title(alg(i,:));
meshc(xg, yg, zg');
endfor
pl_automatic_replot = t;
##colormap(ocmap);
endfunction
|