itjyb Subroutine

public subroutine itjyb(x, tj, ty)

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

! ITJYB computes integrals of Bessel functions J0(t) and Y0(t).

Discussion:

This procedure integrates Bessel functions J0(t) and Y0(t)
with respect to t from 0 to x.

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:

25 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, real ( kind = 8 ) X, the upper limit of the integral.

Output, real ( kind = 8 ) TJ, TY, the integrals of J0(t) and Y0(t) 
from 0 to x.

Arguments

Type IntentOptional Attributes Name
real(kind=8) :: x
real(kind=8) :: tj
real(kind=8) :: ty