From: <car...@us...> - 2011-11-20 21:08:46
|
Revision: 9151 http://octave.svn.sourceforge.net/octave/?rev=9151&view=rev Author: carandraug Date: 2011-11-20 21:08:39 +0000 (Sun, 20 Nov 2011) Log Message: ----------- specfun: help text for texinfo, update GP version Modified Paths: -------------- trunk/octave-forge/main/specfun/inst/Ci.m trunk/octave-forge/main/specfun/inst/Si.m trunk/octave-forge/main/specfun/inst/cosint.m trunk/octave-forge/main/specfun/inst/dirac.m trunk/octave-forge/main/specfun/inst/ellipj.m trunk/octave-forge/main/specfun/inst/ellipke.m trunk/octave-forge/main/specfun/inst/erfcinv.m trunk/octave-forge/main/specfun/inst/expint.m trunk/octave-forge/main/specfun/inst/expint_E1.m trunk/octave-forge/main/specfun/inst/expint_Ei.m trunk/octave-forge/main/specfun/inst/heaviside.m trunk/octave-forge/main/specfun/inst/laguerre.m trunk/octave-forge/main/specfun/inst/lambertw.m trunk/octave-forge/main/specfun/inst/psi.m trunk/octave-forge/main/specfun/inst/sinint.m trunk/octave-forge/main/specfun/inst/zeta.m Modified: trunk/octave-forge/main/specfun/inst/Ci.m =================================================================== --- trunk/octave-forge/main/specfun/inst/Ci.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/Ci.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,15 +13,18 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## CI compute the cosine integral function define by: -## +## -*- texinfo -*- +## @deftypefn {Function File} {@var{y} =} Ci (@var{z}) +## Compute the cosine integral function defined by: +## @verbatim ## Inf ## / ## Ci(x) = | cos(t)/t dt ## / ## x -## -##See also : cosint, Si, sinint, expint, expint_Ei. +## @end verbatim +## @seealso{cosint, Si, sinint, expint, expint_Ei} +## @end deftypefn function y = Ci(z) if (nargin != 1) @@ -33,4 +36,4 @@ y(real(z) == 0 & imag(z) <0) = 0.5*(expint_Ei(i.*y(real(z) == 0 & imag(z) <0))+expint_Ei(-i.*y(real(z) == 0 & imag(z) <0)))-i*pi./2; y(real(z)>=0) = -0.5.*(expint_E1(i.*y(real(z)>=0) )+expint_E1(-i.*y(real(z)>=0) )); y(real(z)<0) = -0.5.*(expint_E1(-i.*y(real(z)<0))+expint_E1(i.*y(real(z)<0)))+i*pi; -endfunction; +endfunction Modified: trunk/octave-forge/main/specfun/inst/Si.m =================================================================== --- trunk/octave-forge/main/specfun/inst/Si.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/Si.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,14 +13,17 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -##SI compute the sine integral function define by: -## +## -*- texinfo -*- +## @deftypefn {Function File} {@var{y} =} Si (@var{x}) +## Compute the sine integral defined by: +## @verbatim ## x ## / ## Si(x) = | sin(t)/t dt ## / ## 0 -## +## @end verbatim +## @end deftypefn function y = Si(x) if (nargin != 1) Modified: trunk/octave-forge/main/specfun/inst/cosint.m =================================================================== --- trunk/octave-forge/main/specfun/inst/cosint.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/cosint.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,19 +13,22 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## COSINT compute the cosine integral function define by: -## +## -*- texinfo -*- +## @deftypefn {Function File} {@var{y} =} cosint (@var{z}) +## Compute the cosine integral function defined by: +## @verbatim ## Inf ## / ## cosint(x) = | cos(t)/t dt ## / ## x -## -##See also : Ci, Si, sinint, expint, expint_Ei. +## @end verbatim +## @seealso{Ci, Si, sinint, expint, expint_Ei} +## @end deftypefn function y = cosint(z) if (nargin != 1) print_usage; endif y = Ci(z); -endfunction; +endfunction Modified: trunk/octave-forge/main/specfun/inst/dirac.m =================================================================== --- trunk/octave-forge/main/specfun/inst/dirac.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/dirac.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -14,7 +14,7 @@ ## along with this program; If not, see <http://www.gnu.org/licenses/>. ## -*- texinfo -*- -## @deftypefn {Function File} {} dirac(@var{x}) +## @deftypefn {Function File} {@var{y} =} dirac(@var{x}) ## Compute the dirac delta function. ## @seealso{heaviside} ## @end deftypefn @@ -25,4 +25,4 @@ endif y = double(x == 0.); y(x == 0) = Inf; -endfunction; +endfunction Modified: trunk/octave-forge/main/specfun/inst/ellipj.m =================================================================== --- trunk/octave-forge/main/specfun/inst/ellipj.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/ellipj.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,26 +13,29 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## usage: [sn,cn,dn] = ellipj(u,m[,tol]) +## -*- texinfo -*- +## @deftypefn {Function File} {[@var{sn}, @var{cn}, @var{dn}] =} ellipj (@var{u}, @var{m}[,@var{tol}]) +## Compute the Jacobi elliptic functions @var{sn}(@var{u}|@var{m}), +## @var{cn}(@var{u}|@var{m}) and @var{dn}(@var{u}|@var{m}) for complex argument +## @var{u} and parameter 0 <= @var{m} <= 1. ## -## Compute the Jacobi elliptic functions sn(u|m), cn(u|m) and dn(u|m) -## for complex argument u and parameter 0 <= m <= 1. +## @emph{warning}: the approximation blows up for abs(u)>20 near m=1. ## -## WARNING: the approximation blows up for abs(u)>20 near m=1. +## @var{tol} is accepted for compatibility, but ignored ## -## tol is accepted for compatibility, but ignored -## -## Ref: Abramowitz, Milton and Stegun, Irene A -## Handbook of Mathematical Functions, Dover, 1965 -## Chapter 16 (Sections 16.4, 16.13 and 16.15) -## -## Example +## @example +## @group ## m = linspace(0,1,200); u=linspace(-10,10,200); ## [U,M] = meshgrid(u,m); ## [sn, cn, dn] = ellipj(U,M); ## imagesc(sn); +## @end group +## @end example ## -## See also: ellipke +## Ref: Abramowitz, Milton and Stegun, Irene A. Handbook of Mathematical +## Functions, Dover, 1965, Chapter 16 (Sections 16.4, 16.13 and 16.15) +## @seealso{ellipke} +## @end deftypefn ## Author: David Billinghurst <Dav...@ri...> ## 2001-01-31 David Billinghurst <Dav...@ri...> Modified: trunk/octave-forge/main/specfun/inst/ellipke.m =================================================================== --- trunk/octave-forge/main/specfun/inst/ellipke.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/ellipke.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,22 +13,19 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## Compute: -## complete elliptic integral of first K(m) -## complete elliptic integral of second E(m) +## -*- texinfo -*- +## @deftypefn {Function File} {[@var{k}, @var{e}] =} ellipke (@var{m}[,@var{tol}]) +## Compute complete elliptic integral of first K(@var{m}) and second E(@var{m}). ## -## usage: [k,e] = ellipke(m[,tol]) +## @var{m} is either real array or scalar with 0 <= m <= 1 ## -## m is either real array or scalar with 0 <= m <= 1 -## -## tol Ignored. -## (Matlab uses this to allow faster, less accurate approximation) +## @var{tol} will be ignored (@sc{Matlab} uses this to allow faster, less +## accurate approximation) ## -## Ref: Abramowitz, Milton and Stegun, Irene A -## Handbook of Mathematical Functions, Dover, 1965 -## Chapter 17 -## -## See also: ellipj +## Ref: Abramowitz, Milton and Stegun, Irene A. Handbook of Mathematical +## Functions, Dover, 1965, Chapter 17. +## @seealso{ellipj} +## @end deftypefn ## Author: David Billinghurst <Dav...@ri...> ## Created: 31 January 2001 @@ -38,7 +35,6 @@ ## 2003-1-18 Jaakko Ruohio ## * extended for m < 0 - function [k,e] = ellipke( m ) if (nargin < 1 || nargin > 2) Modified: trunk/octave-forge/main/specfun/inst/erfcinv.m =================================================================== --- trunk/octave-forge/main/specfun/inst/erfcinv.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/erfcinv.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -24,4 +24,4 @@ print_usage; endif y = erfinv(1-x); -endfunction; +endfunction Modified: trunk/octave-forge/main/specfun/inst/expint.m =================================================================== --- trunk/octave-forge/main/specfun/inst/expint.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/expint.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,19 +13,22 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## EXPINT compute the exponential integral, -## +## -*- texinfo -*- +## @deftypefn {Function File} {@var{y} =} expint (@var{x}) +## Compute the exponential integral, +## @verbatim ## infinity ## / ## expint(x) = | exp(t)/t dt ## / ## x -## -## See also expint_Ei, expint_E1. +## @end verbatim +## @seealso{expint_E1, expint_Ei} +## @end deftypefn function y = expint(x) if (nargin != 1) print_usage; endif y = expint_E1(x); -endfunction; +endfunction Modified: trunk/octave-forge/main/specfun/inst/expint_E1.m =================================================================== --- trunk/octave-forge/main/specfun/inst/expint_E1.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/expint_E1.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,15 +13,18 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## EXPINT_E1 compute the exponential integral, -## +## -*- texinfo -*- +## @deftypefn {Function File} {@var{y} =} expint_E1 (@var{x}) +## Compute the exponential integral, +## @verbatim ## infinity ## / ## expint(x) = | exp(t)/t dt ## / ## x -## -## See also expint_Ei, expint. +## @end verbatim +## @seealso{expint, expint_Ei} +## @end deftypefn function y = expint_E1(x) if (nargin != 1) @@ -32,4 +35,4 @@ y(imag(x) < 0 & imag(x) != 0) = -expint_Ei(-y(imag(x) < 0 & imag(x) != 0)) +i.*pi; y(real(x) >= 0 & imag(x)==0) = -expint_Ei(-y(real(x) >= 0 & imag(x)==0)); y(real(x) < 0 & imag(x)==0) = -expint_Ei(-y(real(x) < 0 & imag(x)==0)) -i.*pi; -endfunction; +endfunction Modified: trunk/octave-forge/main/specfun/inst/expint_Ei.m =================================================================== --- trunk/octave-forge/main/specfun/inst/expint_Ei.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/expint_Ei.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,15 +13,18 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## EXPINT_EI compute the exponential integral, -## +## -*- texinfo -*- +## @deftypefn {Function File} {@var{y} =} expint_Ei (@var{x}) +## Compute the exponential integral, +## @verbatim ## infinity ## / ## expint_Ei(x) = - | exp(t)/t dt ## / ## -x -## -## See also expint, expint_E1. +## @end verbatim +## @seealso{expint, expint_E1} +## @end deftypefn function y = expint_Ei(x) if (nargin != 1) @@ -48,13 +51,13 @@ y(t) = conj(expint_Ei(conj(x(t)))); endif; ## Serie Expansion - else + else for k = 1:100; y(t) = y(t) + x(t).^k./(k.*factorial(k)); endfor y(t) = 0.577215664901532860606512090082402431 + log(x(t)) + y(t); endif endif - endif; + endif endfor -endfunction; +endfunction Modified: trunk/octave-forge/main/specfun/inst/heaviside.m =================================================================== --- trunk/octave-forge/main/specfun/inst/heaviside.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/heaviside.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, Modified: trunk/octave-forge/main/specfun/inst/laguerre.m =================================================================== --- trunk/octave-forge/main/specfun/inst/laguerre.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/laguerre.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, Modified: trunk/octave-forge/main/specfun/inst/lambertw.m =================================================================== --- trunk/octave-forge/main/specfun/inst/lambertw.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/lambertw.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -1,36 +1,9 @@ -%% usage: lambertw(z) or lambertw(n,z) -%% -%% Compute the Lambert W function of z. This function satisfies -%% W(z).*exp(W(z)) = z, and can thus be used to express solutions -%% of transcendental equations involving exponentials or logarithms. -%% -%% n must be integer, and specifies the branch of W to be computed; -%% W(z) is a shorthand for W(0,z), the principal branch. Branches -%% 0 and -1 are the only ones that can take on non-complex values. -%% -%% If either n or z are non-scalar, the function is mapped to each -%% element; both may be non-scalar provided their dimensions agree. -%% -%% This implementation should return values within 2.5*eps of its -%% counterpart in Maple V, release 3 or later. Please report any -%% discrepancies to the author, Nici Schraudolph <sch...@in...>. -%% -%% For further details, see: -%% -%% Corless, Gonnet, Hare, Jeffrey, and Knuth (1996), "On the Lambert -%% W Function", Advances in Computational Mathematics 5(4):329-359. - -%% Author: Nicol N. Schraudolph <sch...@in...> -%% Version: 1.0 -%% Created: 07 Aug 1998 -%% Keywords: Lambert W Omega special transcendental function - %% Copyright (C) 1998 by Nicol N. Schraudolph %% %% This program is free software; you can redistribute and/or %% modify it under the terms of the GNU General Public %% License as published by the Free Software Foundation; -%% either version 2, or (at your option) any later version. +%% either version 3, or (at your option) any later version. %% %% This program is distributed in the hope that it will be %% useful, but WITHOUT ANY WARRANTY; without even the implied @@ -42,6 +15,36 @@ %% License along with this software; see the file COPYING. If not, %% see <http://www.gnu.org/licenses/>. +## -*- texinfo -*- +## @deftypefn {Function File} {@var{x} = } lambertw (@var{z}) +## @deftypefnx {Function File} {@var{x} = } lambertw (@var{z}, @var{n}) +## Compute the Lambert W function of @var{z}. +## +## This function satisfies W(z).*exp(W(z)) = z, and can thus be used to express +## solutions of transcendental equations involving exponentials or logarithms. +## +## @var{n} must be integer, and specifies the branch of W to be computed; +## W(z) is a shorthand for W(0,z), the principal branch. Branches +## 0 and -1 are the only ones that can take on non-complex values. +## +## If either @var{n} or @var{z} are non-scalar, the function is mapped to each +## element; both may be non-scalar provided their dimensions agree. +## +## This implementation should return values within 2.5*eps of its +## counterpart in Maple V, release 3 or later. Please report any +## discrepancies to the author, Nici Schraudolph <sch...@in...>. +## +## For further details, see: +## +## Corless, Gonnet, Hare, Jeffrey, and Knuth (1996), "On the Lambert +## W Function", Advances in Computational Mathematics 5(4):329-359. +## @end deftypefn + +%% Author: Nicol N. Schraudolph <sch...@in...> +%% Version: 1.0 +%% Created: 07 Aug 1998 +%% Keywords: Lambert W Omega special transcendental function + function w = lambertw(b,z) if (nargin == 1) z = b; Modified: trunk/octave-forge/main/specfun/inst/psi.m =================================================================== --- trunk/octave-forge/main/specfun/inst/psi.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/psi.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, Modified: trunk/octave-forge/main/specfun/inst/sinint.m =================================================================== --- trunk/octave-forge/main/specfun/inst/sinint.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/sinint.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -2,7 +2,7 @@ ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,9 +13,12 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -##SININT compute the sine integral function. -##See also: Si. +## -*- texinfo -*- +## @deftypefn {Function File} {@var{y} =} sinint (@var{x}) +## Compute the sine integral function. +## @seealso{Si} +## @end deftypefn -function y = sinint(x); - y = Si(x); -endfunction; +function y = sinint (x) + y = Si (x); +endfunction Modified: trunk/octave-forge/main/specfun/inst/zeta.m =================================================================== --- trunk/octave-forge/main/specfun/inst/zeta.m 2011-11-20 20:26:20 UTC (rev 9150) +++ trunk/octave-forge/main/specfun/inst/zeta.m 2011-11-20 21:08:39 UTC (rev 9151) @@ -1,8 +1,8 @@ -# Copyright (C) 2006 Sylvain Pelissier <syl...@gm...> +## Copyright (C) 2006 Sylvain Pelissier <syl...@gm...> ## ## This program is free software; you can redistribute it and/or modify ## it under the terms of the GNU General Public License as published by -## the Free Software Foundation; either version 2 of the License, or +## the Free Software Foundation; either version 3 of the License, or ## (at your option) any later version. ## ## This program is distributed in the hope that it will be useful, @@ -13,7 +13,11 @@ ## You should have received a copy of the GNU General Public License ## along with this program; If not, see <http://www.gnu.org/licenses/>. -## ZETA compute the Riemann's Zeta function. +## -*- texinfo -*- +## @deftypefn {Function File} {@var{z} =} zeta (@var{t}) +## Compute the Riemann's Zeta function. +## @seealso{Si} +## @end deftypefn function z = zeta(t) if (nargin != 1) @@ -39,4 +43,4 @@ z(j) = 2.^t(j).*pi.^(t(j)-1).*sin(pi.*t(j)./2).*gamma(1-t(j)).*zeta(1-t(j)); endif endfor -endfunction; +endfunction This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site. |