Write a script using Matlab:

CP(T) = R (a+bT +cT^2 ) where R = 8.314 J mol1 K-1.

(Species) (a) (b) (c)

Methane ( 1.702 ) ( 9.081*10^-3 ) ( -2.164*10^-6 )

Ethanol ( 3.518 ) ( 20.001*10^-3 ) ( -6.002*10^-6 )

formaldehyde ( 4.191 ) ( 7.022*10^-3 ) ( -1.877*10^-6 )

Write a function start with:

function Cp = getCp(species, T)

to calculate the heat capacity Cp for a given species and temperature T (in kelvin). The input species should be a string, the input T should be a scalar (a single number), and the output Cp should be a scalar (a single number).

The result should be :

>> getCp('methane', 298)

ans = 35.051