I converted the Fortran routines polint.for xpolint.for ratint.for and xratint.for to double precision, but then they do not work for me.
Any ideas?