Not sure exactly what you want.
syms b ; % 1/T
syms N ;
Z = 2*(2*cosh(b))^N ;
F = log(Z)/N;
E_norm = -diff(F,b);
E_norm = simplify(E_norm) ;
fplot(E_norm,[ .1, 6]);
hold on
fplot(1/b,[.1,6]);
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!