/usr/share/freemat/toolbox/graph/subplot.m is in freemat-data 4.0-5.
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 103 104 105 106 107 108 109 110 | % SUBPLOT SUBPLOT Subplot Function
%
% Usage
%
% This function divides the current figure into a 2-dimensional
% grid, each of which can contain a plot of some kind. The function
% has a number of syntaxes. The first version
%
% subplot(row,col,num)
%
% which either activates subplot number num, or
% sets up a subplot grid of size row x col, and then
% activates num. You can also set up subplots that cover multiple
% grid elements
%
% subplot(row,col,[vec])
%
% where vec is a set of indexes covered by the new subplot.
% Finally, as a shortcut, you can specify a string with three
% components
%
% subplot('mnp')
%
% or using the alternate notation
%
% subplot mnp
%
% where m is the number of rows, n is the number of columns
% and p is the index. You can also specify the location of the
% subplot explicitly using the syntax
%
% subplot('position',[left bottom width height])
%
%
% Copyright (c) 2002-2006 Samit Basu
% Licensed under the GPL
function hout = subplot(varargin)
m = 1; n = 1; p = 1;
if ((nargin == 2) && (isa(varargin{1},'char')))
if (strcmp(lower(varargin{1}),'position'))
h = axes('outerposition',varargin{2});
if (nargout > 0), hout = h; end;
return;
end
end
if (nargin == 1)
if (isa(varargin{1},'char'))
str = varargin{1};
if (length(str) ~= 3), error('Subplot with a string argument requires a string of the type mnp'); end
m = str(1)-'0';
n = str(2)-'0';
p = str(3)-'0';
else
axes(varargin{1}(1));
h = varargin{1}(1);
end
else
if (nargin >= 1), m = round(varargin{1}(1)); end
if (nargin >= 2), n = round(varargin{2}(1)); end
if (nargin >= 3), p = round(varargin{3}); p = p(:); end
end
row = m+1-round(idiv(p-1,n)+1);
col = round(mod(p-1,n)+1);
width = 1.0/n;
height = 1.0/m;
left = (col-1)*width;
bottom = (row-1)*height;
position = [left,bottom,left+width,bottom+height];
% Get the envelope of the axis
if (length(p) > 1)
position = [min(position(:,1:2)),max(position(:,3:4))];
end
position = [position(1),position(2),position(3)-position(1),...
position(4)-position(2)];
fig = gcf;
children = get(fig,'children');
found = 0;
deletelist = [];
for (i=1:length(children))
if (ishandle(children(i),'axes'))
outerpos = get(children(i),'outerposition');
if (all(outerpos==position))
axes(children(i));
found = 1;
elseif (intersects(outerpos,position))
deletelist = [deletelist,i];
end
end
end
children(deletelist) = [];
set(fig,'children',children);
if (~found)
h = axes('outerposition',position);
end
set(fig,'nextplot','add');
if (nargout > 0) hout = h; end;
function b = intersects(rect1,rect2)
nleft = max(rect1(1),rect2(1));
nbottom = max(rect1(2),rect2(2));
nright = min(rect1(1)+rect1(3),rect2(1)+rect2(3));
ntop = min(rect1(2)+rect1(4),rect2(2)+rect2(4));
if ((nright <= nleft) | (ntop <= nbottom))
b = 0;
else
area = (nright-nleft)*(ntop-nbottom);
b = area > .01;
end
|