Mercurial > octave-nkf
annotate scripts/plot/fplot.m @ 10549:95c3e38098bf
Untabify .m scripts
| author | Rik <code@nomad.inbox5.com> |
|---|---|
| date | Fri, 23 Apr 2010 11:28:50 -0700 |
| parents | b14fd5116c29 |
| children | be55736a0783 |
| rev | line source |
|---|---|
| 8920 | 1 ## Copyright (C) 2005, 2006, 2007, 2008, 2009 Paul Kienzle |
| 5820 | 2 ## |
| 3 ## This file is part of Octave. | |
| 4 ## | |
| 5 ## Octave is free software; you can redistribute it and/or modify it | |
| 6 ## under the terms of the GNU General Public License as published by | |
| 7016 | 7 ## the Free Software Foundation; either version 3 of the License, or (at |
| 8 ## your option) any later version. | |
| 5820 | 9 ## |
| 10 ## Octave is distributed in the hope that it will be useful, but | |
| 11 ## WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
| 13 ## General Public License for more details. | |
| 14 ## | |
| 15 ## You should have received a copy of the GNU General Public License | |
| 7016 | 16 ## along with Octave; see the file COPYING. If not, see |
| 17 ## <http://www.gnu.org/licenses/>. | |
| 5820 | 18 |
| 19 ## -*- texinfo -*- | |
| 20 ## @deftypefn {Function File} {} fplot (@var{fn}, @var{limits}) | |
| 6781 | 21 ## @deftypefnx {Function File} {} fplot (@var{fn}, @var{limits}, @var{tol}) |
| 5820 | 22 ## @deftypefnx {Function File} {} fplot (@var{fn}, @var{limits}, @var{n}) |
| 6895 | 23 ## @deftypefnx {Function File} {} fplot (@dots{}, @var{fmt}) |
| 6699 | 24 ## Plot a function @var{fn}, within the defined limits. @var{fn} |
| 5820 | 25 ## an be either a string, a function handle or an inline function. |
| 26 ## The limits of the plot are given by @var{limits} of the form | |
| 27 ## @code{[@var{xlo}, @var{xhi}]} or @code{[@var{xlo}, @var{xhi}, | |
|
9040
dbd0c77e575e
Cleanup documentation file plot.texi
Rik <rdrider0-list@yahoo.com>
parents:
8920
diff
changeset
|
28 ## @var{ylo}, @var{yhi}]}. @var{tol} is the default tolerance to use for the |
| 6781 | 29 ## plot, and if @var{tol} is an integer it is assumed that it defines the |
| 6895 | 30 ## number points to use in the plot. The @var{fmt} argument is passed |
| 31 ## to the plot command. | |
| 5820 | 32 ## |
| 33 ## @example | |
|
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
9040
diff
changeset
|
34 ## @group |
| 6699 | 35 ## fplot ("cos", [0, 2*pi]) |
| 36 ## fplot ("[cos(x), sin(x)]", [0, 2*pi]) | |
|
9051
1bf0ce0930be
Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents:
9040
diff
changeset
|
37 ## @end group |
| 5820 | 38 ## @end example |
| 6781 | 39 ## @seealso{plot} |
| 5820 | 40 ## @end deftypefn |
| 41 | |
| 7017 | 42 ## Author: Paul Kienzle <pkienzle@users.sf.net> |
| 43 | |
| 6781 | 44 function fplot (fn, limits, n, linespec) |
| 8847 | 45 if (nargin < 2 || nargin > 4) |
| 6046 | 46 print_usage (); |
| 5820 | 47 endif |
| 48 | |
|
10400
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
49 if (!isreal (limits) || (numel (limits) != 2 && numel (limits) != 4)) |
|
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
50 error ("fplot: second input argument must be a real vector with 2 or 4 elements"); |
|
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
51 endif |
|
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
52 |
| 5820 | 53 if (nargin < 3) |
| 6781 | 54 n = 0.002; |
| 5820 | 55 endif |
| 56 | |
| 6781 | 57 have_linespec = true; |
| 58 if (nargin < 4) | |
| 59 have_linespec = false; | |
| 60 endif | |
| 5820 | 61 |
| 6781 | 62 if (ischar (n)) |
| 63 have_linespec = true; | |
| 64 linespec = n; | |
| 65 n = 0.002; | |
| 66 endif | |
| 67 | |
| 6699 | 68 if (strcmp (typeinfo (fn), "inline function")) |
| 69 fn = vectorize (fn); | |
| 70 nam = formula (fn); | |
| 71 elseif (isa (fn, "function_handle")) | |
| 72 nam = func2str (fn); | |
| 5820 | 73 elseif (all (isalnum (fn))) |
| 6781 | 74 nam = fn; |
|
10400
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
75 elseif (ischar (fn)) |
|
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
76 fn = vectorize (inline (fn)); |
|
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
77 nam = formula (fn); |
| 5820 | 78 else |
|
10400
b14fd5116c29
Ensure that 'limits' is a 2 or 4 vector, and that 'fn' is a function
Soren Hauberg <hauberg@gmail.com>
parents:
9051
diff
changeset
|
79 error ("fplot: first input argument must be a function handle, inline function or string"); |
| 6781 | 80 endif |
| 81 | |
| 82 if (floor(n) != n) | |
| 83 tol = n; | |
| 7280 | 84 x0 = linspace (limits(1), limits(2), 5)'; |
| 6781 | 85 y0 = feval (fn, x0); |
| 86 err0 = Inf; | |
| 7280 | 87 n = 8; |
| 6781 | 88 x = linspace (limits(1), limits(2), n)'; |
| 89 y = feval (fn, x); | |
| 90 | |
| 91 while (n < 2 .^ 20) | |
| 92 y00 = interp1 (x0, y0, x, "linear"); | |
| 93 err = 0.5 * max (abs ((y00 - y) ./ (y00 + y))(:)); | |
| 94 if (err == err0 || 0.5 * max (abs ((y00 - y) ./ (y00 + y))(:)) < tol) | |
| 10549 | 95 break; |
| 6781 | 96 endif |
| 97 x0 = x; | |
| 98 y0 = y; | |
| 99 err0 = err; | |
| 100 n = 2 * (n - 1) + 1; | |
| 101 x = linspace (limits(1), limits(2), n)'; | |
| 102 y = feval (fn, x); | |
| 103 endwhile | |
| 104 else | |
| 105 x = linspace (limits(1), limits(2), n)'; | |
| 106 y = feval (fn, x); | |
| 5820 | 107 endif |
| 108 | |
| 6781 | 109 if (have_linespec) |
| 110 plot (x, y, linespec); | |
| 111 else | |
| 112 plot (x, y); | |
| 113 endif | |
|
8145
7ef5b1b4e029
fplot.m: call axis after plot
John W. Eaton <jwe@octave.org>
parents:
7280
diff
changeset
|
114 |
|
7ef5b1b4e029
fplot.m: call axis after plot
John W. Eaton <jwe@octave.org>
parents:
7280
diff
changeset
|
115 if (length (limits) > 2) |
|
7ef5b1b4e029
fplot.m: call axis after plot
John W. Eaton <jwe@octave.org>
parents:
7280
diff
changeset
|
116 axis (limits); |
|
7ef5b1b4e029
fplot.m: call axis after plot
John W. Eaton <jwe@octave.org>
parents:
7280
diff
changeset
|
117 endif |
|
7ef5b1b4e029
fplot.m: call axis after plot
John W. Eaton <jwe@octave.org>
parents:
7280
diff
changeset
|
118 |
| 8507 | 119 if (isvector (y)) |
| 6781 | 120 legend (nam); |
| 121 else | |
| 8507 | 122 for i = 1:columns (y) |
| 6781 | 123 nams{i} = sprintf ("%s(:,%i)", nam, i); |
| 124 endfor | |
| 125 legend (nams{:}); | |
| 126 endif | |
| 5820 | 127 endfunction |
| 7245 | 128 |
| 129 %!demo | |
| 130 %! fplot ("cos", [0, 2*pi]) | |
| 131 | |
| 132 %!demo | |
| 133 %! fplot ("[cos(x), sin(x)]", [0, 2*pi]) |