octave-nkf: scripts/specfun/ellipke.m annotate

annotate scripts/specfun/ellipke.m @ 16597:49832f60282e

if available, use ioctl to give up controlling terminal * octave-gui.cc (dissociate_terminal): If sys/ioctl.h is available and TIOCNOTTY is defined, use ioctl to give up controlling terminal on Unixy systems.

author	Ben Abbott <bpabbott@mac.com>
date	Wed, 01 May 2013 02:02:21 -0400
parents	a3fdd6041e64
children	e2de3c8882be

rev	line source
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	1 ## Copyright (C) 2001 David Billinghurst <David.Billinghurst@riotinto.com>
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	2 ## Copyright (C) 2001 Paul Kienzle <pkienzle@users.sf.net>
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	3 ## Copyright (C) 2003 Jaakko Ruohio
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	4 ##
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	5 ## This file is part of Octave.
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	6 ##
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	7 ## Octave is free software; you can redistribute it and/or modify it
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	8 ## under the terms of the GNU General Public License as published by
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	9 ## the Free Software Foundation; either version 3 of the License, or (at
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	10 ## your option) any later version.
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	11 ##
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	12 ## Octave is distributed in the hope that it will be useful, but
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	13 ## WITHOUT ANY WARRANTY; without even the implied warranty of
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	14 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	15 ## General Public License for more details.
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	16 ##
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	17 ## You should have received a copy of the GNU General Public License
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	18 ## along with Octave; see the file COPYING. If not, see
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	19 ## <http://www.gnu.org/licenses/>.
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	20
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	21 ## -- texinfo --
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	22 ## @deftypefn {Function File} {} ellipke (@var{m})
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	23 ## @deftypefnx {Function File} {} ellipke (@var{m}, @var{tol})
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	24 ## @deftypefnx {Function File} {[@var{k}, @var{e}] =} ellipke (@dots{})
16587 a3fdd6041e64 Fix ellipj, ellipke, and expint docstrings Mike Miller <mtmiller@ieee.org> parents: 16586 diff changeset	25 ## Compute complete elliptic integral of the first K(@var{m}) and second
a3fdd6041e64 Fix ellipj, ellipke, and expint docstrings Mike Miller <mtmiller@ieee.org> parents: 16586 diff changeset	26 ## E(@var{m}) kind.
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	27 ##
16587 a3fdd6041e64 Fix ellipj, ellipke, and expint docstrings Mike Miller <mtmiller@ieee.org> parents: 16586 diff changeset	28 ## @var{m} is either real array or scalar with 0 <= m <= 1.
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	29 ##
16587 a3fdd6041e64 Fix ellipj, ellipke, and expint docstrings Mike Miller <mtmiller@ieee.org> parents: 16586 diff changeset	30 ## @var{tol} is currently ignored (@sc{Matlab} uses this to allow faster,
a3fdd6041e64 Fix ellipj, ellipke, and expint docstrings Mike Miller <mtmiller@ieee.org> parents: 16586 diff changeset	31 ## less accurate approximation).
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	32 ##
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	33 ## Ref: Abramowitz, Milton and Stegun, Irene A. Handbook of Mathematical
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	34 ## Functions, Dover, 1965, Chapter 17.
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	35 ## @seealso{ellipj}
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	36 ## @end deftypefn
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	37
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	38 function [k, e] = ellipke (m)
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	39
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	40 if (nargin < 1 \|\| nargin > 2)
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	41 print_usage ();
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	42 endif
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	43
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	44 k = e = zeros (size (m));
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	45 m = m(:);
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	46 if any (!isreal (m))
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	47 error ("ellipke must have real m");
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	48 endif
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	49 if any (m > 1)
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	50 error ("ellipke must have m <= 1");
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	51 endif
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	52
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	53 Nmax = 16;
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	54 idx = find (m == 1);
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	55 if (!isempty (idx))
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	56 k(idx) = Inf;
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	57 e(idx) = 1;
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	58 endif
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	59
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	60 idx = find (m == -Inf);
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	61 if (!isempty (idx))
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	62 k(idx) = 0;
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	63 e(idx) = Inf;
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	64 endif
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	65
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	66 ## Arithmetic-Geometric Mean (AGM) algorithm
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	67 ## ( Abramowitz and Stegun, Section 17.6 )
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	68 idx = find (m != 1 & m != -Inf);
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	69 if (!isempty (idx))
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	70 idx_neg = find (m < 0 & m != -Inf);
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	71 mult_k = 1./sqrt (1 - m(idx_neg));
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	72 mult_e = sqrt (1 - m(idx_neg));
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	73 m(idx_neg) = -m(idx_neg)./(1 - m(idx_neg));
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	74 a = ones (length (idx), 1);
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	75 b = sqrt (1 - m(idx));
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	76 c = sqrt (m(idx));
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	77 f = 0.5;
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	78 sum = fc.c;
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	79 for n = 2:Nmax
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	80 t = (a + b)/2;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	81 c = (a - b)/2;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	82 b = sqrt (a.*b);
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	83 a = t;
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	84 f = f * 2;
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	85 sum = sum + fc.c;
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	86 if (all (c./a < eps)) break; endif
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	87 endfor
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	88 if (n >= Nmax) error ("ellipke: not enough workspace"); endif
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	89 k(idx) = 0.5*pi./a;
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	90 e(idx) = 0.5pi.(1 - sum)./a;
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	91 k(idx_neg) = mult_k.*k(idx_neg);
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	92 e(idx_neg) = mult_e.*e(idx_neg);
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	93 endif
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	94
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	95 endfunction
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	96
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	97 %% Test complete elliptic functions of first and second kind
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	98 %% against "exact" solution from Mathematica 3.0
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	99 %!test
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	100 %! m = [0.0; 0.01; 0.1; 0.5; 0.9; 0.99; 1.0 ];
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	101 %! [k,e] = ellipke (m);
16584 2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	102 %!
2f766ceeb03e Add ellipj, ellipke, and expint functions from Octave Forge Mike Miller <mtmiller@ieee.org> parents: diff changeset	103 %! # K(1.0) is really infinity - see below
16586 f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	104 %! k_exp = [1.5707963267948966192;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	105 %! 1.5747455615173559527;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	106 %! 1.6124413487202193982;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	107 %! 1.8540746773013719184;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	108 %! 2.5780921133481731882;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	109 %! 3.6956373629898746778;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	110 %! 0.0 ];
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	111 %! e_exp = [1.5707963267948966192;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	112 %! 1.5668619420216682912;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	113 %! 1.5307576368977632025;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	114 %! 1.3506438810476755025;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	115 %! 1.1047747327040733261;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	116 %! 1.0159935450252239356;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	117 %! 1.0 ];
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	118 %! if k(7)==Inf, k(7)=0; endif;
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	119 %! assert (k, k_exp, 8*eps);
f423873d3275 Style fixes for ellipj.cc, ellipke.m, and expint.m Mike Miller <mtmiller@ieee.org> parents: 16585 diff changeset	120 %! assert (e, e_exp, 8*eps);
16585 1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	121
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	122 %% Test against A&S Table 17.1
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	123 %!test
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	124 %! m = [0:5:50]'/100;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	125 %! k_exp = [1.570796326794897;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	126 %! 1.591003453790792;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	127 %! 1.612441348720219;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	128 %! 1.635256732264580;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	129 %! 1.659623598610528;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	130 %! 1.685750354812596;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	131 %! 1.713889448178791;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	132 %! 1.744350597225613;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	133 %! 1.777519371491253;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	134 %! 1.813883936816983;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	135 %! 1.854074677301372 ];
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	136 %! e_exp = [1.570796327;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	137 %! 1.550973352;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	138 %! 1.530757637;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	139 %! 1.510121831;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	140 %! 1.489035058;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	141 %! 1.467462209;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	142 %! 1.445363064;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	143 %! 1.422691133;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	144 %! 1.399392139;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	145 %! 1.375401972;
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	146 %! 1.350643881 ];
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	147 %! [k,e] = ellipke (m);
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	148 %! assert (k, k_exp, 1e-15);
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	149 %! assert (e, e_exp, 1e-8);
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	150
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	151 %% Test input validation
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	152 %!error ellipke ()
1a3bfb14b5da Add and fix tests for ellipj, ellipke, and expint Mike Miller <mtmiller@ieee.org> parents: 16584 diff changeset	153 %!error ellipke (1,2,3)

Mercurial > hg > octave-nkf

annotate scripts/specfun/ellipke.m @ 16597:49832f60282e