I need help to write MATLAB codes of the following question. And also I need help with answering the questions in the text.

We want to estimate the integral:

with MATLAB program. The integral can be estimated analytically with at least ten decimal places with the value:

Use this value as a reference value. Begin estimating the integral over the entire interval with MATLAB's integral. Start with the standard tolerance. Then set a higher tolerance. Does that match the reference value? Then solve the integral with the trapezoidal rule. Start by dividing the interval into 20 parts. How much and why does the integral value differ from the reference value? Then halve the step length in the trapezoidal rule until:

How many halves did you need? How many function evaluations has the trapezoidal rule done? Was integral more effective? Why?

-2 e-()? dx. 0.0003 Ipef = 0.0039/T. Iref Irap