Legrangian interpolation using mathematica part 2