solution in MATLAB

The equation CP=a+bT+cT2+dT3 is an empirical polynomial that describes the behavior of the heat capacity CP as a function of temperature in kelvins. The change in enthalpy (a measure of energy) as a gas is heated from T1 to T2 is the integral of this equation with respect to T : h=T1T2CPdT Find the change in enthalpy of oxygen gas as it is heated from 300 to 1000 K, using the MATLAB quadrature function integral. The values of a,b,c, and d for oxygen are as follows: a=25.48b=1.520102c=0.7155105d=1.312109