Oliver Kullmann <O.Kullmann <at> swansea.ac.uk> writes:
>
> Hello,
>
> I couldn't find information on whether the logarithmic integrals
>
> Li_m(x) = integral_0^x log(t)^(-m) dt
>
> for x >= 0 are available in R?
I saw your request only this weekend.
The first logarithmic integral can be computed using the exponential
integral Ei(x) per
li(x) = Ei(log(x))
and elliptic integrals are part of the 'gsl' package, so
library('gsl')
x <- seq(2, 10, by=0.5)
y <- expint_Ei(log(x))
y
See e.g. the Handbook of Mathematical Functions for how to reduce higher
logarithmic integrals.
Another possibility is to use the Web API of 'keisan', the calculation
library of Casio.
Regards
Hans Werner
> Best wishes
>
> Oliver
>
>
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
