Lagrange interpolation Matlab

The user should be able to enter: 1. A function named f. 2. A number of points (nodes) with their respective values. 3. A point x0 The output should be: 1. A Lagrange polynomial (function of x) that approximate the function with agreement in the points. 2. A figure representing both the curve of the function and that of the polynomial in the interval [xmin, xmax] where xmin (xmax) is the minimum (resp. maximum) nodes value. 3. An approximation of f(x0) by Lagrange polynomial. 4. The approximation absolute and relative errors. The code should be well documented explaining all the steps needed to achieve the final result.