Title: Fortran subroutines for computing the optimal interpolation formula