3ProblemSet2Problem9 PDF
3ProblemSet2Problem9 PDF
\\base\root...\CBE353ProblemSet2Problem9.m
for i=1:length(T)
Cv(i)=R*((((E^2)/((K^2)*(T(i)^2)))*exp(-E/(K*T(i))))/(1+exp(-E/(K*T(i))))^2);
end
plot(T,Cv)%Plotting the values
xlabel('Temperature in Kelvin')
ylabel('Molar Heat Capacity')
%Calculating the maximum
f = @(x) R*((((E^2)/((K^2)*(x^2)))*exp(-E/(K*x)))/(1+exp(-E/(K*x)))^2);
max_x = fminbnd(@(x) -f(x),1,5000)
1 of 1