function xFunction
% xFunction.m
% An Example of using Function
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/02/05 $

clear
clc
disp('Calculating the reaction rate of ammonia synthesis')
t = input('vector of temperatures (℃) = ')
T = t + 273.15;         % T: Temperatre, K
dH = ReactionHeat(T)

% Ploting the results
plot(T,dH,'o-')
xlabel('T (K)')
ylabel('(-ΔH_r)')
title(' Reaction rate VS temperature')

% ------------------------------------------------------------------
function dH = ReactionHeat(T)
% T: Temperatre, K
% dH: Reaction heat (-ΔHr), kJ/kmol
dH = (10000+0.219*T-2.845e-3*T.^2+0.9703e-6*T.^3)*4.184;