Create Gauss-Legendre integration points and weights for one-dimensional integration on the interval [x1,x2].

Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|

real(kind=rk), | intent(in) | :: | x1 | lower limit of integration interval |
||

real(kind=rk), | intent(in) | :: | x2 | upper limit of integration interval |
||

real(kind=rk), | intent(out) | :: | x(:) | The coordinates of the gauss points on the interval [x1,x2]. The array has the length nIntP. |
||

real(kind=rk), | intent(out) | :: | w(:) | The quadrature weights. The array has the length nIntP. |
||

integer, | intent(in) | :: | nIntP | The number of integration points. |

Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|

real(kind=rk), | private | :: | z1 | some working variables |
|||

real(kind=rk), | private | :: | z | some working variables |
|||

real(kind=rk), | private | :: | xm | some working variables |
|||

real(kind=rk), | private | :: | xl | some working variables |
|||

real(kind=rk), | private | :: | pp | some working variables |
|||

real(kind=rk), | private | :: | p3 | some working variables |
|||

real(kind=rk), | private | :: | p2 | some working variables |
|||

real(kind=rk), | private | :: | p1 | some working variables |
|||

real(kind=rk), | private | :: | EPS | the relative precision of the points |
|||

integer, | private | :: | m | ||||

integer, | private | :: | i | ||||

integer, | private | :: | j |