octave_packages/symbolic-1.1.0/splot.m

   1 ## Copyright (C) 2002 Ben Sapp <bsapp@nua.lampf.lanl.gov>
   2 ##
   3 ## This program is free software; you can redistribute it and/or modify
   4 ## it under the terms of the GNU General Public License as published by
   5 ## the Free Software Foundation; either version 2 of the License, or
   6 ## (at your option) any later version.
   7 ##
   8 ## This program is distributed in the hope that it will be useful,
   9 ## but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  11 ## GNU General Public License for more details.
  12 ##
  13 ## You should have received a copy of the GNU General Public License
  14 ## along with this program; If not, see <http://www.gnu.org/licenses/>.
  15
  16 ## -*- texinfo -*-
  17 ## @deftypefn {Function File} splot(@var{f},@var{x},@var{range})
  18 ## Plot a symbolic function f(x) over range.
  19 ## @end deftypefn
  20
  21 function splot(expression,symbol,range)
  22   ## we should be a little smarter about this
  23   t = linspace(min(range),max(range),400);
  24   x = zeros(size(t));
  25   j = 1;
  26   for i = t
  27     x(j) = to_double(subs(expression,symbol,vpa(t(j))));
  28     j++;
  29   endfor
  30
  31   plot(t,x);
  32 endfunction