Please see Approximatrix Simply Fortran.