These are routines from the public domain slatec library,
https://www.netlib.org/slatec

We include only routines that we use, not the entire library.
