home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21fb.zip / octave / SCRIPTS.ZIP / scripts.fat / control / ax2dlim.m < prev    next >
Text File  |  1999-12-24  |  2KB  |  64 lines
  1. ## Copyright (C) 1998 Auburn University.  All Rights Reserved
  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 the 
  7. ## Free Software Foundation; either version 2, or (at your option) any 
  8. ## later version. 
  9. ## 
  10. ## Octave is distributed in the hope that it will be useful, but WITHOUT 
  11. ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 
  12. ## FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License 
  13. ## for more details.
  14. ## 
  15. ## You should have received a copy of the GNU General Public License 
  16. ## along with Octave; see the file COPYING.  If not, write to the Free 
  17. ## Software Foundation, 59 Temple Place, Suite 330, Boston, MA 02111 USA. 
  18.  
  19. ## -*- texinfo -*-
  20. ## @deftypefn{Function File } { @var{axvec} =} ax2dlim (@var{axdata})
  21. ##  determine axis limits for 2-d data(column vectors); leaves a 10% margin 
  22. ##  around the plots.
  23. ##  puts in margins of +/- 0.1 if data is one dimensional (or a single point)
  24. ## 
  25. ## @strong{Inputs}
  26. ##    @var{axdata} nx2 matrix of data [x,y]
  27. ## 
  28. ## @strong{Outputs}
  29. ##    @var{axvec} vector of axis limits appropriate for call to axis() function
  30. ## @end deftypefn
  31.  
  32. function axvec = ax2dlim (axdata)
  33.  
  34.   if(isempty(axdata))
  35.     axdata = 0;
  36.   endif
  37.  
  38.   ## compute axis limits
  39.   minv = min(axdata);
  40.   maxv = max(axdata);
  41.   delv = (maxv-minv)/2;      # breadth of the plot
  42.   midv = (minv + maxv)/2;    # midpoint of the plot
  43.   axmid = [midv(1), midv(1), midv(2), midv(2)];
  44.   axdel = [-0.1, 0.1,-0.1,0.1];   # default plot width (if less than 2-d data)
  45.   if(max(delv) == 0)
  46.     if(midv(1) != 0)
  47.       axdel(1:2) = [-0.1*midv(1),0.1*midv(1)];
  48.     endif
  49.     if(midv(2) != 0)
  50.       axdel(3:4) = [-0.1*midv(2),0.1*midv(2)];
  51.     endif
  52.   else
  53.     ## they're at least one-dimensional
  54.     if(delv(1) != 0)
  55.       axdel(1:2) = 1.1*[-delv(1),delv(1)];
  56.     endif
  57.     if(delv(2) != 0)
  58.       axdel(3:4) = 1.1*[-delv(2),delv(2)];
  59.     endif
  60.   endif
  61.   axvec = axmid + axdel; 
  62. endfunction
  63.  
  64.