bernob Subroutine

public subroutine bernob(n, bn)

************80

! BERNOB computes the Bernoulli number Bn.

Licensing:

This routine is copyrighted by Shanjie Zhang and Jianming Jin.  However, 
they give permission to incorporate this routine into a user program 
provided that the copyright is acknowledged.

Modified:

11 July 2012

Author:

Shanjie Zhang, Jianming Jin

Reference:

Shanjie Zhang, Jianming Jin,
Computation of Special Functions,
Wiley, 1996,
ISBN: 0-471-11963-6,
LC: QA351.C45.

Parameters:

Input, integer ( kind = 4 ) N, the index.

Output, real ( kind = 8 ) BN, the value of the N-th Bernoulli number.

Arguments

Type IntentOptional Attributes Name
integer(kind=4) :: n
real(kind=8) :: bn(0:n)