Pretty much working

Uses:
double precision variables
spline + 3 power laws to give number of lines (correct)
reproduces Ly alpha forest reliably
includes LLS ( log NHI > 17 )
robust to change in number of points for integration
shows no tendency to crash, although not tested