function xFunction
clear; clc
disp('a')
t=input('t')
I=t+273.15;
dH=ReactionHeat(I)
plot(I,dH,'k-.')
xlabel('I(K)'),
ylabel('(-\DeltaH\_r)'),
title('Reaction rate VS tempature')

function dH=ReactionHeat(I)
dH=(10000+0.219*I-2.845e-3*I.^2+0.9703e-6*I.^-3)*4.184;