[R] Discontinuities in a simple graph (machine precision?)

[Simon Cullen]

Thanks to Prof. Ripley, Peter Dalgaard and Uwe Ligges for pointing out  
that all I needed to do to improve the precision was replace
> inted <- integrate(integrand, a, Inf)
with
> inted <- integrate(integrand, a, Inf,rel.tol = .Machine$double.eps^0.5)

Thanks for the advice!

-- 
SC

Simon Cullen
Room 3030
Dept. Of Economics
Trinity College Dublin

Ph. (608)3477
Email cullens at tcd.ie