Develop a MATLAB code to solve for cubic spline with parabolic run$-$out $y_0^{\text{'}\text{'}}=y_1^{\text{'}\text{'}}$;$y_3^{\text{'}\text{'}}=y_4^{\text{'}\text{'}},$ to fit the

data $(x_0,y_0)=(0,0),(x_1,y_1)=(1,1),(x_2,y_2)=(2,4),(x_3,y_3)=(3,1),(x_4,y_4)=(4,0).$ Submit the code, which

should include a graph of the piecewise cubic functions $($Use more grids for the interpolating function$).$