In matlab please help me

Bessel functions arise in a variety of engineering and science applications. These functions are not amenable to straightforward evaluations, so they are usually compiled in standard tables, for example: Use this data to estimate the 1st order Bessel function corresponding to x=2.13, namely J1(2.13) using: 1) A4th order interpolating polynomial generated using MATLAB polyfit. Calculate the true relative error (t) in this estimate. For the error calculation, determine the true value using the built-in function bessel j. See MATLAB help for syntax and more information about Bessel functions. The subscript 1 in J1(x) indicates the order of the Bessel function (the nu input you see in the help documentation. 2) A cubic spline with "not-a-knot" end conditions generated using MATLAB. Calculate the true relative error (t) in this estimate too. For the next 2 interpolations, formulate the interpolating polynomial by hand and calculate it in your m-file as an equation, or multiple equations if you want to break up the long equation for part (4) 3) A 1st order Lagrange interpolating polynomial. Calculate the true relative error (t) in this estimate. 4) A2nd order Lagrange interpolating polynomial Make sure to use the closest 3 data points. Calculate the true relative error (t) in this estimate. Your script should only output the estimate and t from parts (1) to (4) in 2-element row vectors