a. Using the second-order Taylor method with the indicated step size h, calculate the solution estimates at x1, x2, and x3.
b. Write a MATLAB script file that uses Euler ODE and Taylor_2_ ODE (see Problem 12) to find the approximate values produced by Euler’s and second-order Taylor methods and returns a table that includes these values, as well as the exact values and the global percent relative error for both methods, at all mesh points in the given interval.
Write a user-defined function with function call y = Taylor_2_ ODE(f,fp,x,y0) that solves a first-order IVP in the form y′ = f(x,y), y(x0) = y0 using the second-order Taylor method. The input argument fp represents the first derivative of the function f, while the others are as in Euler ODE. Execute the function to solve