Legrangian interpolation using mathematica part 1