Sets the gauss Laguerre roots and weights for 30 discretization points. Used for integration in this module. roots
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp), | intent(out), | dimension(no_roots) | :: | roots |
subroutine gauss_laguerre_30_roots(roots)
!! Sets the gauss Laguerre roots and weights for 30
!! discretization points. Used for integration in this module.
!! roots
!
! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
real(dp), dimension(no_roots), intent(out):: roots
roots= (/4.740718054080526184d-02,&
2.499239167531593919d-01,&
6.148334543927683749d-01,&
1.143195825666101451d+00,&
1.836454554622572344d+00,&
2.696521874557216147d+00,&
3.725814507779509288d+00,&
4.927293765849881879d+00,&
6.304515590965073635d+00,&
7.861693293370260349d+00,&
9.603775985479263255d+00,&
1.153654659795613924d+01,&
1.366674469306423489d+01,&
1.600222118898106771d+01,&
1.855213484014315029d+01,&
2.132720432178312819d+01,&
2.434003576453269346d+01,&
2.760555479678096091d+01,&
3.114158670111123683d+01,&
3.496965200824907072d+01,&
3.911608494906788991d+01,&
4.361365290848483056d+01,&
4.850398616380419980d+01,&
5.384138540650750571d+01,&
5.969912185923549686d+01,&
6.618061779443848991d+01,&
7.344123859555988076d+01,&
8.173681050672767867d+01,&
9.155646652253683726d+01,&
1.041575244310588886d+02/)
end subroutine gauss_laguerre_30_roots