add a new function interpolant(mu,solution,M)