Mercurial > octave-libgccjit
annotate scripts/plot/ellipsoid.m @ 8920:eb63fbe60fab
update copyright notices
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Sat, 07 Mar 2009 10:41:27 -0500 |
| parents | 967c15c5e265 |
| children | dbd0c77e575e |
| rev | line source |
|---|---|
| 8920 | 1 ## Copyright (C) 2007, 2008 Sylvain Pelissier |
| 7306 | 2 ## |
| 3 ## This file is part of Octave. | |
| 7305 | 4 ## |
| 7306 | 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 | |
| 7 ## the Free Software Foundation; either version 3 of the License, or (at | |
| 8 ## your option) any later version. | |
| 7305 | 9 ## |
| 7306 | 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. | |
| 7305 | 14 ## |
| 15 ## You should have received a copy of the GNU General Public License | |
| 7306 | 16 ## along with Octave; see the file COPYING. If not, see |
| 17 ## <http://www.gnu.org/licenses/>. | |
| 7305 | 18 |
| 19 ## -*- texinfo -*- | |
| 20 ## @deftypefn {Function File} {[@var{x}, @var{y}, @var{z}] =} ellipsoid (@var{xc},@var{yc}, @var{zc}, @var{xr}, @var{yr}, @var{zr}, @var{n}) | |
| 21 ## @deftypefnx {Function File} {} ellipsoid (@var{h}, @dots{}) | |
| 7306 | 22 ## Generate three matrices in @code{meshgrid} format that define an |
| 23 ## ellipsoid. Called with no return arguments, @code{ellipsoid} calls | |
| 24 ## directly @code{surf (@var{x}, @var{y}, @var{z})}. If an axes handle | |
| 25 ## is passed as the first argument, the the surface is plotted to this | |
| 26 ## set of axes. | |
| 7305 | 27 ## @seealso{sphere} |
| 28 ## @end deftypefn | |
| 29 | |
| 7306 | 30 ## Author: Sylvain Pelissier <sylvain.pelissier@gmail.com> |
| 31 | |
| 32 function [xx, yy, zz] = ellipsoid (varargin) | |
| 7305 | 33 |
| 34 [h, varargin, nargin] = __plt_get_axis_arg__ ((nargout > 0), "ellipsoid", | |
| 35 varargin{:}); | |
| 36 | |
| 37 if (nargin != 6 && nargin != 7) | |
| 38 print_usage (); | |
| 39 endif | |
| 40 | |
| 41 xc = varargin{1}; | |
| 42 yc = varargin{2}; | |
| 43 zc = varargin{3}; | |
| 44 xr = varargin{4}; | |
| 45 yr = varargin{5}; | |
| 46 zr = varargin{6}; | |
| 47 | |
|
8274
967c15c5e265
ellipsoid.m: fix default argument check
Soren Hauberg <hauberg@gmail.com>
parents:
7306
diff
changeset
|
48 if (nargin == 6) |
| 7305 | 49 n = 20; |
| 50 else | |
| 51 n = varargin{7}; | |
| 52 endif | |
| 53 | |
| 54 theta = linspace (0, 2 * pi, n + 1); | |
| 55 phi = linspace (-pi / 2, pi / 2, n + 1); | |
| 56 [theta, phi] = meshgrid (theta, phi); | |
| 57 | |
| 58 x = xr .* cos (phi) .* cos (theta) + xc; | |
| 59 y = yr .* cos (phi) .* sin (theta) + yc; | |
| 60 z = zr .* sin (phi) + zc; | |
| 61 | |
| 62 if (nargout > 0) | |
| 63 xx = x; | |
| 64 yy = y; | |
| 65 zz = z; | |
| 66 else | |
| 67 surf (h, x, y, z); | |
| 68 endif | |
| 69 | |
| 70 endfunction | |
| 71 | |
| 72 %!demo | |
| 73 %! ellipsoid (0, 0, 1, 2, 3, 4, 20); |