1 /* specfunc/gsl_sf_zeta.h
3 * Copyright (C) 1996, 1997, 1998, 1999, 2000, 2004 Gerard Jungman
5 * This program is free software; you can redistribute it and/or modify
6 * it 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.
10 * This program 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.
15 * You should have received a copy of the GNU General Public License
16 * along with this program; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
20 /* Author: G. Jungman */
22 #ifndef __GSL_SF_ZETA_H__
23 #define __GSL_SF_ZETA_H__
25 #include <gsl/gsl_sf_result.h>
30 # define __BEGIN_DECLS extern "C" {
31 # define __END_DECLS }
33 # define __BEGIN_DECLS /* empty */
34 # define __END_DECLS /* empty */
40 /* Riemann Zeta Function
41 * zeta(n) = Sum[ k^(-n), {k,1,Infinity} ]
44 * exceptions: GSL_EDOM, GSL_EOVRFLW
46 int gsl_sf_zeta_int_e(const int n, gsl_sf_result * result);
47 double gsl_sf_zeta_int(const int n);
50 /* Riemann Zeta Function
51 * zeta(x) = Sum[ k^(-s), {k,1,Infinity} ], s != 1.0
54 * exceptions: GSL_EDOM, GSL_EOVRFLW
56 int gsl_sf_zeta_e(const double s, gsl_sf_result * result);
57 double gsl_sf_zeta(const double s);
60 /* Riemann Zeta Function minus 1
61 * useful for evaluating the fractional part
62 * of Riemann zeta for large argument
65 * exceptions: GSL_EDOM, GSL_EOVRFLW
67 int gsl_sf_zetam1_e(const double s, gsl_sf_result * result);
68 double gsl_sf_zetam1(const double s);
71 /* Riemann Zeta Function minus 1 for integer arg
72 * useful for evaluating the fractional part
73 * of Riemann zeta for large argument
76 * exceptions: GSL_EDOM, GSL_EOVRFLW
78 int gsl_sf_zetam1_int_e(const int s, gsl_sf_result * result);
79 double gsl_sf_zetam1_int(const int s);
82 /* Hurwitz Zeta Function
83 * zeta(s,q) = Sum[ (k+q)^(-s), {k,0,Infinity} ]
86 * exceptions: GSL_EDOM, GSL_EUNDRFLW, GSL_EOVRFLW
88 int gsl_sf_hzeta_e(const double s, const double q, gsl_sf_result * result);
89 double gsl_sf_hzeta(const double s, const double q);
93 * eta(n) = (1-2^(1-n)) zeta(n)
95 * exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
97 int gsl_sf_eta_int_e(int n, gsl_sf_result * result);
98 double gsl_sf_eta_int(const int n);
102 * eta(s) = (1-2^(1-s)) zeta(s)
104 * exceptions: GSL_EUNDRFLW, GSL_EOVRFLW
106 int gsl_sf_eta_e(const double s, gsl_sf_result * result);
107 double gsl_sf_eta(const double s);
112 #endif /* __GSL_SF_ZETA_H__ */