{VERSION 3 0 "IBM INTEL NT" "3.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 
0 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }
{CSTYLE "" -1 256 "" 1 16 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 
257 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 0 
0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 
0 }{CSTYLE "" -1 260 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 
261 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 0 0 0 
0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 263 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 
0 }{CSTYLE "" -1 264 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 
265 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 266 "" 0 1 0 0 0 
0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 
0 }{CSTYLE "" -1 268 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 
269 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 270 "" 0 1 0 0 0 
0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 271 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 
0 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 
273 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 274 "" 0 1 0 0 0 
0 0 0 1 0 0 0 0 0 0 }{CSTYLE "" -1 275 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 
0 }{CSTYLE "" -1 276 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 }{PSTYLE "Normal
" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }0 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 
-1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }
{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }
1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "
" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 
-1 0 }}
{SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT 256 23 "Entropy and Temperature" 
}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 227 "This w
orksheet contains calculations that illustrate the statistical mechani
cs definition of entropy. Two coupled Einstein solids are considered f
ollowing the explanations by Thomas A Moore and Daniel V Schroeder (Am
. J. Phys. " }{TEXT 257 2 "65" }{TEXT -1 13 ", 26 (1997))." }}{PARA 0 
"" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 741 "In macroscopic t
hermodynamics entropy is defined as the quantity that increases by Q/T
 when a thermodynamic system receives the energy amount Q by heating w
hile at temperature T. This is not intuitive, and the additional expla
nations about entropy measuring order in the system are also ad-hoc. T
his definition of entropy is referred to as the Clausius definition. B
oltzmann is the originator of the statistical mechanics definition. In
 the Einstein solid the thermodynamic system is simplified by consider
ing a given number of oscillators at fixed locations. Each atom has th
ree degrees of freedom, and the oscillators can take up energy in quan
tized units. The quantum zero-point energy of the oscillators is irrel
evant for the discussion." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "
" 0 "" {TEXT -1 520 "The connection between the macroscopic and micros
copic descriptions is as follows. Macroscopically the system is descri
bed by a a few variables (For a tank of an atomic gas: Number of atoms
, Volume of the gas, and total internal energy U; other variables such
 as pressure and temperature follow from the ideal gas law). Microscop
ically there can be many states that are consistent with this macrosta
te. The number of consistent microstates is called the multiplicity. D
ifferent macrostates have different multiplicities." }}{PARA 0 "" 0 "
" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 435 "If a system is isolate
d in a given macrostate it is equally likely to be in any of the consi
stent microstates - they have all equal probability [fundamental assum
ption of stat mech]. The probability equals the inverse multiplicity f
or that macrostate. This assumption leads to the second law of thermod
ynamics, namely that macroscopic objects exhibit irreversible behaviou
r. To prove this point we consider two coupled Einstein solids." }}
{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 93 "First we \+
calculate the multiplicity for an Einstein solid (for a derivation see
: H B Callen, " }{TEXT 258 54 "Thermodynamics and an Introduction to T
hermostatistics" }{TEXT -1 10 ", p. 334)." }}{PARA 0 "" 0 "" {TEXT -1 
0 "" }}{PARA 0 "" 0 "" {TEXT -1 20 "We have a system of " }{TEXT 262 
1 "N" }{TEXT -1 53 " one-dimensional oscillators (the number of atoms \+
is " }{TEXT 263 1 "N" }{TEXT -1 45 "/3), each oscillator can be in a s
tate where " }{TEXT 264 1 "n" }{TEXT -1 29 " units of energy are store
d, " }{TEXT 259 1 "n" }{TEXT -1 50 "=0,1,2,...; the total energy of th
e system equals " }{TEXT 265 1 "q" }{TEXT -1 18 " units of energy (" }
{TEXT 266 1 "q" }{TEXT -1 53 " is also a non-negative integer). The un
it of energy " }{TEXT 260 1 "E" }{TEXT -1 4 " = h" }{TEXT 261 2 " f" }
{TEXT -1 116 " , where f is he natural frequency of the oscillator and
 h is Planck's constant. The multiplicity for a system with " }{TEXT 
268 1 "N" }{TEXT -1 42 " degrees of freedom and a total number of " }
{TEXT 267 1 "q" }{TEXT -1 85 " units of energy can be shown to be (at \+
least try to verify this by counting for low " }{TEXT 270 1 "N" }
{TEXT -1 2 ", " }{TEXT 269 1 "q" }{TEXT -1 2 "):" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 21 "restart; with(plots):" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 32 "Omega:=(N,q)->binomial(q+N-1,q);" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#>%&OmegaGR6$%\"NG%\"qG6\"6$%)operatorG%&arrowGF)
-%)binomialG6$,(9%\"\"\"9$F2!\"\"F2F1F)F)F)" }}}{EXCHG {PARA 0 "" 0 "
" {TEXT -1 140 "A system of 3 oscillators with 3 units of energy share
d can be in a state (300), (210), (111), (012), (030), ... a total of \+
10 combinations." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "Omega(3
,3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"#5" }}}{EXCHG {PARA 0 "" 0 "
" {TEXT -1 830 "The two coupled Einstein solids represent an interacti
ng system; we label the two by A and B. We can consider a macrostate w
here the total energy U is split between A and B as U = U_A + U_B. The
 reason for being able to specify U_A and U_B lies in the fact that th
e coupling between the two solids is weak: while within the solids the
 systems are fluctuating between the microstates consistent with the m
acrostates (e.g., U_A, N_A), the transfer of energy between A and B ha
ppens on a longer timescale and it is possible to know U_A and U_B. On
e calculates the multiplicities for A and B, and takes the product of \+
these to find the multiplicity for the coupled system. This is true si
nce the two solids are independent of each other, and any microstate i
n A can combine with any microstate in B (for given macrostates in A a
nd B)." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 
257 "What we need now is a counting scheme: for given N_A, N_B and tot
al number of excitation quanta q_t, we wish to find out how many micro
states there are for a given combination where q_A, and q_B units of e
nergy stored in A, B respectively (q_A + q_B = q_t)." }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 19 "N_A:=300; N_B:=200;" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%$N_AG\"$+$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$N_BG
\"$+#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "q_t:=20;" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#>%$q_tG\"#?" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT -1 114 "We work with lists in Maple (lists are ordered sets), an
d define ourselves a function to add an element to a list:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "ladd:=(L,e)->[op(L),e];" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#>%%laddGR6$%\"LG%\"eG6\"6$%)operatorG%&arrow
GF)7$-%#opG6#9$9%F)F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "
ladd([1,2,3],4);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7&\"\"\"\"\"#\"\"$
\"\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "O_A:=[]: O_B:=[]: \+
O_t:=[]: " }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "for q_A from 0 to q_t \+
do:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "q_B:=q_t-q_A;" }}{PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 20 "o_A:=Omega(N_A,q_A);" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 20 "o_B:=Omega(N_B,q_B);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 
0 13 "o_t:=o_A*o_B;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "O_A:=ladd(O_
A,o_A); O_B:=ladd(O_B,o_B); O_t:=ladd(O_t,o_t);" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 3 "od:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 115 "The lists
 can be printed. If q is large, however, the lists are very long, and \+
we have commented out the statement." }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 5 "#O_A;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "#O_
B;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "#O_t;" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "listplot(O_t);" }}{PARA 13 "" 1 "" 
{GLPLOT2D 631 225 225 {PLOTDATA 2 "6#-%'CURVESG6#777$$\"\"\"\"\"!$\"+R
04#3\"\"#>7$$\"\"#F*$\"+:;jkH\"#?7$$\"\"$F*$\"+I\"*p))Q\"#@7$$\"\"%F*$
\"+vL:ZK\"#A7$$\"\"&F*$\"+.(*)e$>\"#B7$$\"\"'F*$\"+&fD#f()FE7$$\"\"(F*
$\"+bd)47$\"#C7$$\"\")F*$\"+`$Rt'*)FP7$$\"\"*F*$\"+:*z,6#\"#D7$$\"#5F*
$\"+Lo,2TFen7$$\"#6F*$\"+-,]ZmFen7$$\"#7F*$\"+c7dj*)Fen7$$\"#8F*$\"+4(
p^+\"\"#E7$$\"#9F*$\"+YKJB$*Fen7$$\"#:F*$\"+RD+$3(Fen7$$\"#;F*$\"+EZjR
VFen7$$\"#<F*$\"+,L.%4#Fen7$$\"#=F*$\"+-$H)pwFP7$$F-F*$\"+9F01?FP7$$F3
F*$\"+yZzSLFE7$$F9F*$\"+OQGkEF?" 1 2 0 1 0 2 9 1 4 2 1.000000 
45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 141 "We want
 to look at the width of the distribution. We can write a procedure to
 generate the mean, and the deviation of the data from the mean." }}
{PARA 0 "" 0 "" {TEXT -1 108 "In the calculation of the average we sub
tract by one since the first value of the list corresponds to q_A=0." 
}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "MSD:=proc(L) local n,i,av
,dev,wt;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "n:=nops(L); wt:=add(L[i
],i=1..n):" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "av:=add((i)*L[i],i=1.
.n)/wt:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "dev:=evalf(sqrt(add((i-a
v)^2*L[i],i=1..n)/wt));" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "print(`a
verage: `,av-1,` deviation: `,dev);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 
4 "end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "MSD(O_t);" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6&%*average:~G\"#7%-~deviation:~G$\"+vu/
KA!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 237 "We calculated the expe
ctation value of the distance squared from the average value and took \+
the square root to obtain a distance. The expectation values for the a
verage are normalized by dividing the sums by the sum over all data po
ints." }}{PARA 0 "" 0 "" {TEXT -1 17 "N_A=300, N_B=200:" }}{PARA 0 "" 
0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 52 "We see from repeated
 calculations with increased q: " }}{PARA 0 "" 0 "" {TEXT -1 18 "q=100
: 60, dev 5.4" }}{PARA 0 "" 0 "" {TEXT -1 19 "q=200: 120, dev 8.2" }}
{PARA 0 "" 0 "" {TEXT -1 20 "q=300: 180, dev 10.7" }}{PARA 0 "" 0 "" 
{TEXT -1 20 "q=400: 240, dev 13.1" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}
{PARA 0 "" 0 "" {TEXT -1 20 "q=800: 480, dev 22.3" }}{PARA 0 "" 0 "" 
{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 294 "For fixed number of degr
ees of freedom we have increased the amount of energy deposited. The c
ombination of energy quanta distributed between A and B follows some o
rder: q_A = 6/10 * q_t is the location with the largest multiplicity f
or the choice of N_A vs N_B. How does the width grow with q?" }}{PARA 
0 "" 0 "" {TEXT -1 74 "This has to be considered carefully. The number
 of states increases like q" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 
0 "" 0 "" {TEXT -1 18 " For 20/30 we get " }}{PARA 0 "" 0 "" {TEXT -1 
20 "q=800: 480. dev 56.6" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "v1:=56.6/480;" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%#v1G$\"+nm;z6!#5" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 14 "v2:=111.5/960;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#
v2G$\"+L$e9;\"!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "v2/v1;
" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+;L#)\\)*!#5" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 17 "evalf(1/sqrt(2));" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"+5y1rq!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
15 "v3:=166.3/1440;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#v3G$\"+66'[:
\"!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "v3/v1;" }}{PARA 11 
"" 1 "" {XPPMATH 20 "6#$\"+W^(Qz*!#5" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 17 "evalf(1/sqrt(3));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#
$\"+$p-Nx&!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "v4:=13.1/2
40;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#v4G$\"+LLLea!#6" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "v4/v1;" }}{PARA 11 "" 1 "" {XPPMATH 
20 "6#$\"+j_(*GY!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "eval
f(1/sqrt(4));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+++++]!#5" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "v8:=22.3/480;" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#>%#v8G$\"+LL$ek%!#6" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 6 "v8/v1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+F$H*RR
!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "evalf(1/sqrt(8));" }
}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+0R`NN!#5" }}}{EXCHG {PARA 0 "" 0 
"" {TEXT -1 162 "The relative width is not decreasing as 1/sqrt(q), bu
t we are also not in the limit of large q, yet. If N was larger than q
, the rel. width would go as 1/sqrt(N)." }}{PARA 0 "" 0 "" {TEXT -1 0 
"" }}{PARA 0 "" 0 "" {TEXT -1 694 "For many degrees of freedom (N appr
oaching 10^24, and implying also a large number of energy units deposi
ted) the distribution becomes very narrowly peaked. Thus, there is an \+
almost uniquely defined macrostate with a very high multiplicity compa
red to the others. The second law of thermodynamics (that states that \+
entropy increases) can be cast into the form that the system evolves a
lways towards the macrostate with the largest multiplicity. It happens
 since in the course of fluctuations the system is bound to evolve int
o the state with the largest probability. Given that all microstates a
re assumed to be equally probable, it must evolve towards the macrosta
te with largest multiplicity." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 305 "Entropy is defined now as the log
arithm of the multiplicity. This ensures the additivity of entropies i
n weakly coupled systems (since their combined multiplicity was given \+
as a product of the individual ones). In addition, the entropy definit
ion contains Boltzmann's constant as a proportionality factor." }}
{PARA 0 "" 0 "" {TEXT -1 15 "S = k ln(Omega)" }}{PARA 0 "" 0 "" {TEXT 
-1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 116 "For plotting purposes we chose
 k=1 (a unit system). We can calculate the entropies for the subsystem
s and the total." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 37 "P1:=lis
tplot(map(log,O_A),color=red):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 38 "P2:=listplot(map(log,O_B),color=blue):" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 39 "P3:=listplot(map(log,O_t),color=green):" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "display(P1,P2,P3);" }}{PARA 
13 "" 1 "" {GLPLOT2D 684 262 262 {PLOTDATA 2 "6%-%'CURVESG6$777$$\"\"
\"\"\"!F*7$$\"\"#F*$\"+vCy.d!\"*7$$\"\"$F*$\"+cXxr5!\")7$$\"\"%F*$\"+H
g&H`\"F67$$\"\"&F*$\"+t)*pl>F67$$\"\"'F*$\"+_)ekP#F67$$\"\"(F*$\"+$39$
pFF67$$\"\")F*$\"+y:3ZJF67$$\"\"*F*$\"+*>A=^$F67$$\"#5F*$\"+>(4^'QF67$
$\"#6F*$\"+Q`=3UF67$$\"#7F*$\"+SI0UXF67$$\"#8F*$\"+m;an[F67$$\"#9F*$\"
+\\qM&=&F67$$\"#:F*$\"+N;1'\\&F67$$\"#;F*$\"+9f>+eF67$$\"#<F*$\"+1V>)4
'F67$$\"#=F*$\"+$>Z/R'F67$$\"#>F*$\"+%>+tn'F67$$\"#?F*$\"+N91fpF67$$\"
#@F*$\"+<t+OsF6-%'COLOURG6&%$RGBG$\"*++++\"F6F*F*-F$6$777$F($\"+Yx7bkF
67$F,$\"++Qz:iF67$F2$\"+#>)yrfF67$F8$\"+Kc$Gs&F67$F=$\"+E\"H'oaF67$FB$
\"+&>C)3_F67$FG$\"+9;.V\\F67$FL$\"+I\"33n%F67$FQ$\"+RWk\"R%F67$FV$\"+!
H\\\\5%F67$Fen$\"+k!G+\"QF67$Fjn$\"+\\J01NF67$F_o$\"+)z@?>$F67$Fdo$\"+
tSpmGF67$Fio$\"+ru\\GDF67$F^p$\"+?CPv@F67$Fcp$\"+7U]/=F67$Fhp$\"+]I\"=
T\"F67$F]q$\"+%4v%3**F07$Fbq$\"+ntJ)H&F07$FgqF*-F\\r6&F^rF*F*F_r-F$6$7
7Fdr7$F,$\"+[?<'y'F67$F2$\"+[FcVqF67$F8$\"+h;zbsF67$F=$\"+***GVV(F67$F
B$\"+ZIG&e(F67$FG$\"+(pXBr(F67$FL$\"+3(*)y\"yF67$FQ$\"+PmY.zF67$FV$\"+
5!f+(zF67$Fen$\"+-M@=!)F67$Fjn$\"+*=1\"[!)F67$F_o$\"+lMcf!)F67$Fdo$\"+
B6/_!)F67$Fio$\"+2\"fX-)F67$F^p$\"+M$ob(zF67$Fcp$\"+=&)p-zF67$Fhp$\"+V
-E-yF67$F]q$\"+/x9owF67$Fbq$\"+rJ*))[(F6Ffq-F\\r6&F^rF*F_rF*" 1 2 0 1 
0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" 
{TEXT -1 112 "The graph shows that at equilibrium (which happens at q_
A=60) the total entropy is stationary and has a maximum." }}{PARA 0 "
" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 208 "We understand now
 that the evolution of the macrosystem can be described as consistent \+
with the 2nd law of thermodynamics in the form that the isolated syste
m evolves towards a state with the largest entropy." }}{PARA 0 "" 0 "
" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 99 "Note that the total ent
ropy S_t is fairly flat in the vicinity of the maximum as a function o
f q_A." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT 271 30 "How can we define temperature?" }}{PARA 0 "" 0 "" {TEXT 
-1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 206 "One can adopt a definition bas
ed on energy flow and thermal equilibrium: two objects in thermal cont
act are at the same temperature if they are in thermal equilibrium, in
 which case there is no spontaneous " }{TEXT 272 3 "net" }{TEXT -1 
532 " flow between them. If there is a spontaneous net flow from A to \+
B, then A loses energy and is at a higher temperature, and B gains ene
rgy and is at a lower temperature. Now to relate temperature to entrop
y we have to find in the above figure at the equilibrium point two qua
ntitites that become equal. These are the slopes of the curves for S_A
(q_A) and S_B(q_A), which are equal in magnitude and opposite in sign.
 S_B(q_B) has the same slope as S_A(q_A), i.e., the sign is also equal
. Thus, temperature has to be related to dS/dq." }}{PARA 0 "" 0 "" 
{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 163 "On dimensional grounds (
given the choice of constant k in the definition of S) it is the inver
se temperature that equals the rate of change of entropy with energy.
" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 824 "The \+
above figure is used by Moore and Schroeder to justify this. Assume th
at you are on the right of the equilibrium point q_A=60 units, i.e., s
olid A has more energy than what it would have at equilibrium. In this
 region the curve for S_A (red) has a positive slope that is levelling
 off. The (blue) curve has a negative slope that is increasing in magn
itude. What does this mean? If a small amount of energy (e.g. one unit
 dq=1) were to pass from A to B in this region (a movement to the left
 in the graph, or a move towards equilibrium), the increase in S_B wou
ld be bigger than the decrease in S_A. The total entropy would increas
e (green curve). The second law of TD tells us that this process will \+
happen spontaneously (there are many more microstates in this region, \+
i.e., a fluctuation will find this regime easily)." }}{PARA 0 "" 0 "" 
{TEXT -1 253 "The steeper the slope in the entropy versus energy curve
 the more the system wants to obey the second law and increase its ene
rgy (B), while the a shallower entropy-energy curve for the other syst
em (A) means that it doesn't mind to give up some energy." }}{PARA 0 "
" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 374 "To summarize: to \+
the right of the equilibrium point A has more energy than at equilibri
um and energy flows from A to B to reach equilibrium. A thus has a hig
her temperature than B there. Its entropy-energy curve has a smaller s
lope, while the slope for |S_B(q_A)|  is larger. Thus, the inverse of \+
the temperature is proportional to the rate of change of entropy with \+
energy:" }}{PARA 0 "" 0 "" {TEXT -1 11 "1/T = dS/dU" }}{PARA 0 "" 0 "
" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 171 "When the volume of the
 system is allowed to change (as is true for gases) one has to take pa
rtial derivatives, and one cannot simply rearrange the above eq. into \+
dS=dU/T !" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 
256 "For the Einstein solid we can use the above equation to calculate
 the temperature as T = dU/dS by using neighbouring values in the entr
opy as the energy is changed by one unit up and down. A central differ
ence formula is used, the spacing in U equals dq=2." }}{PARA 0 "" 0 "
" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 81 "We take a single solid \+
with 50 degrees of freedom and vary q between zero and 100" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "N:=50;" }}{PARA 11 "" 1 "" {XPPMATH 
20 "6#>%\"NG\"#]" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "O1:=[]:" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 23 "for q from 0 to 100 do:" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 15 "o1:=Omega(N,q);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 16 
"O1:=ladd(O1,o1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "S1:=map(evalf,map(log,O1)):" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT -1 37 "Now we can calculate the temperature:" }}}{EXCHG {PARA 0 
"> " 0 "" {MPLTEXT 1 0 7 "T1:=[]:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 22 "for q from 2 to 99 do:" }}{PARA 0 "> " 0 "" {MPLTEXT 
1 0 24 "t1:=2/(S1[q+1]-S1[q-1]);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 
"T1:=ladd(T1,t1); od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "li
stplot(T1);" }}{PARA 13 "" 1 "" {GLPLOT2D 695 227 227 {PLOTDATA 2 "6#-
%'CURVESG6#7^q7$$\"\"\"\"\"!$\"+-%Gpz#!#57$$\"\"#F*$\"+QgO$G$F-7$$\"\"
$F*$\"+3-vyOF-7$$\"\"%F*$\"+o#z$HSF-7$$\"\"&F*$\"+/OW_VF-7$$\"\"'F*$\"
+q3ccYF-7$$\"\"(F*$\"+O*en%\\F-7$$\"\")F*$\"+y-EE_F-7$$\"\"*F*$\"+_!os
\\&F-7$$\"#5F*$\"+q:OhdF-7$$\"#6F*$\"+][r>gF-7$$\"#7F*$\"+]hAtiF-7$$\"
#8F*$\"+O*)fAlF-7$$\"#9F*$\"+7]RonF-7$$\"#:F*$\"+w/26qF-7$$\"#;F*$\"+q
2+^sF-7$$\"#<F*$\"+3))\\)[(F-7$$\"#=F*$\"+q!GQs(F-7$$\"#>F*$\"+/E@dzF-
7$$\"#?F*$\"+#[W))=)F-7$$\"#@F*$\"+E&*))=%)F-7$$\"#AF*$\"+A>\\Z')F-7$$
\"#BF*$\"+9zxu))F-7$$\"#CF*$\"+![e35*F-7$$\"#DF*$\"+3<$eK*F-7$$\"#EF*$
\"+5Zy\\&*F-7$$\"#FF*$\"+k_zs(*F-7$$\"#GF*$\"+=I$\\***F-7$$\"#HF*$\"+^
gi@5!\"*7$$\"#IF*$\"+AMoV5Fdt7$$\"#JF*$\"+=0nl5Fdt7$$\"#KF*$\"+x>f(3\"
Fdt7$$\"#LF*$\"+1?X46Fdt7$$\"#MF*$\"+nWDJ6Fdt7$$\"#NF*$\"+wG+`6Fdt7$$
\"#OF*$\"+Y/qu6Fdt7$$\"#PF*$\"+q,N'>\"Fdt7$$\"#QF*$\"+eZ&z@\"Fdt7$$\"#
RF*$\"+=n^R7Fdt7$$\"#SF*$\"+0%Q5E\"Fdt7$$\"#TF*$\"+f>_#G\"Fdt7$$\"#UF*
$\"+x$pRI\"Fdt7$$\"#VF*$\"+BDQD8Fdt7$$\"#WF*$\"+?JwY8Fdt7$$\"#XF*$\"+w
F6o8Fdt7$$\"#YF*$\"+(*HV*Q\"Fdt7$$\"#ZF*$\"+\">D2T\"Fdt7$$\"#[F*$\"+k1
*>V\"Fdt7$$\"#\\F*$\"+c1B`9Fdt7$$\"#]F*$\"+3jWu9Fdt7$$\"#^F*$\"+*pQc\\
\"Fdt7$$\"#_F*$\"+f)3o^\"Fdt7$$\"#`F*$\"+^x&z`\"Fdt7$$\"#aF*$\"+fi3f:F
dt7$$\"#bF*$\"+\\_>!e\"Fdt7$$\"#cF*$\"+QbG,;Fdt7$$\"#dF*$\"+VyNA;Fdt7$
$\"#eF*$\"+3HTV;Fdt7$$\"#fF*$\"+@9Xk;Fdt7$$\"#gF*$\"++SZ&o\"Fdt7$$\"#h
F*$\"+k7[1<Fdt7$$\"#iF*$\"+-QZF<Fdt7$$\"#jF*$\"+[@X[<Fdt7$$\"#kF*$\"+>
oTp<Fdt7$$\"#lF*$\"+[$o.z\"Fdt7$$\"#mF*$\"+]rI6=Fdt7$$\"#nF*$\"+,PBK=F
dt7$$\"#oF*$\"+b%[J&=Fdt7$$\"#pF*$\"+$y^S(=Fdt7$$\"#qF*$\"+\"4W\\*=Fdt
7$$\"#rF*$\"+Dd#e\">Fdt7$$\"#sF*$\"+iqpO>Fdt7$$\"#tF*$\"+\\%ev&>Fdt7$$
\"#uF*$\"+k,Ty>Fdt7$$\"#vF*$\"+zDD**>Fdt7$$\"#wF*$\"+ef3??Fdt7$$\"#xF*
$\"+i0\"4/#Fdt7$$\"#yF*$\"+;nsh?Fdt7$$\"#zF*$\"+QY`#3#Fdt7$$\"#!)F*$\"
+yXL.@Fdt7$$\"#\")F*$\"+7o7C@Fdt7$$\"##)F*$\"+K:\"\\9#Fdt7$$\"#$)F*$\"
+e*)ol@Fdt7$$\"#%)F*$\"+_$fk=#Fdt7$$\"#&)F*$\"+oGA2AFdt7$$\"#')F*$\"+%
pzzA#Fdt7$$\"#()F*$\"+_+t[AFdt7$$\"#))F*$\"+OTZpAFdt7$$\"#*)F*$\"+_?@!
H#Fdt7$$\"#!*F*$\"+QS%4J#Fdt7$$\"#\"*F*$\"+E-nJBFdt7$$\"##*F*$\"+s2R_B
Fdt7$$\"#$*F*$\"+[e5tBFdt7$$\"#%*F*$\"+[b\"QR#Fdt7$$\"#&*F*$\"+u+_9CFd
t7$$\"#'*F*$\"+Q&>_V#Fdt7$$\"#(*F*$\"+SS\"fX#Fdt7$$\"#)*F*$\"+UPgwCFdt
" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 7 "S1[13];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6
#$\"+&))3'=G!\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "T1[13];
" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+O*)fAl!#5" }}}{EXCHG {PARA 0 "
" 0 "" {TEXT -1 90 "We are out by one step when comparing the entropy \+
and temperature lists! (cf fig. 5 in MS)" }}{PARA 0 "" 0 "" {TEXT -1 
164 "The graph shows how the temperature increases with the energy amo
unt stored in the solid. Traditionally one measures the heat capacity \+
as a function of temperature." }}{PARA 0 "" 0 "" {TEXT -1 216 "The hea
t capacity is the rate of change of the amount of energy with temperat
ure. We can calculate it also by differencing. One normalizes it usual
ly by the number of oscillators (to take out the size of the system).
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "C1:=[]:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "for q from 2 to 97 do:" }}{PARA 0 "
> " 0 "" {MPLTEXT 1 0 26 "c1:=2/(T1[q+1]-T1[q-1])/N;" }}{PARA 0 "> " 
0 "" {MPLTEXT 1 0 20 "C1:=ladd(C1,c1); od:" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT -1 100 "The simple listplot graphs the specific heat capacity (a
t fixed volume) as a function of energy (q)." }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 13 "listplot(C1);" }}{PARA 13 "" 1 "" {GLPLOT2D 711 
243 243 {PLOTDATA 2 "6#-%'CURVESG6#7\\q7$$\"\"\"\"\"!$\"+KQ1OX!#57$$\"
\"#F*$\"+/_$=O&F-7$$\"\"$F*$\"+W!>u$fF-7$$\"\"%F*$\"+/\"QxP'F-7$$\"\"&
F*$\"+/QVInF-7$$\"\"'F*$\"+WlC@qF-7$$\"\"(F*$\"+)Q,gE(F-7$$\"\")F*$\"+
C-AvuF-7$$\"\"*F*$\"+s>GcwF-7$$\"#5F*$\"+#*oc9yF-7$$\"#6F*$\"+s#>T&zF-
7$$\"#7F*$\"+sB0y!)F-7$$\"#8F*$\"+3(3))=)F-7$$\"#9F*$\"+S(Q$)G)F-7$$\"
#:F*$\"+S3Ay$)F-7$$\"#;F*$\"+7%[(f%)F-7$$\"#<F*$\"+cC*R`)F-7$$\"#=F*$
\"+Sp%=g)F-7$$\"#>F*$\"+ou1k')F-7$$\"#?F*$\"+;eH@()F-7$$\"#@F*$\"+'RxS
x)F-7$$\"#AF*$\"+%y#)G#))F-7$$\"#BF*$\"+O=7o))F-7$$\"#CF*$\"+#*)\\,\"*
)F-7$$\"#DF*$\"+'\\v#\\*)F-7$$\"#EF*$\"++$pd)*)F-7$$\"#FF*$\"+_D()>!*F
-7$$\"#GF*$\"+cqz^!*F-7$$\"#HF*$\"+o4t\"3*F-7$$\"#IF*$\"+O?%)4\"*F-7$$
\"#JF*$\"+G:GO\"*F-7$$\"#KF*$\"++;=h\"*F-7$$\"#LF*$\"+C:m%=*F-7$$\"#MF
*$\"+?f$o?*F-7$$\"#NF*$\"+?jzF#*F-7$$\"#OF*$\"+GHjZ#*F-7$$\"#PF*$\"+SZ
Vm#*F-7$$\"#QF*$\"+'4lUG*F-7$$\"#RF*$\"+;>>,$*F-7$$\"#SF*$\"+OEG<$*F-7
$$\"#TF*$\"+W$)eK$*F-7$$\"#UF*$\"+o3;Z$*F-7$$\"#VF*$\"+%=\\5O*F-7$$\"#
WF*$\"+cYHu$*F-7$$\"#XF*$\"+3d$pQ*F-7$$\"#YF*$\"+k9,*R*F-7$$\"#ZF*$\"+
?Yb5%*F-7$$\"#[F*$\"+_xf@%*F-7$$\"#\\F*$\"+WI<K%*F-7$$\"#]F*$\"++?IU%*
F-7$$\"#^F*$\"+sv+_%*F-7$$\"#_F*$\"+g(>8Y*F-7$$\"#`F*$\"+?zDq%*F-7$$\"
#aF*$\"+O\\$)y%*F-7$$\"#bF*$\"+;B3([*F-7$$\"#cF*$\"+gZ,&\\*F-7$$\"#dF*
$\"+'>PE]*F-7$$\"#eF*$\"+%ew*4&*F-7$$\"#fF*$\"+3s/<&*F-7$$\"#gF*$\"+?W
&Q_*F-7$$\"#hF*$\"+sjTI&*F-7$$\"#iF*$\"+#R[n`*F-7$$\"#jF*$\"+!y\\Ga*F-
7$$\"#kF*$\"+%)Gu[&*F-7$$\"#lF*$\"+?*RWb*F-7$$\"#mF*$\"+o7$*f&*F-7$$\"
#nF*$\"+C-Cl&*F-7$$\"#oF*$\"+crPq&*F-7$$\"#pF*$\"+OsMv&*F-7$$\"#qF*$\"
+%ec,e*F-7$$\"#rF*$\"+![/[e*F-7$$\"#sF*$\"+3KJ*e*F-7$$\"#tF*$\"+;mn$f*
F-7$$\"#uF*$\"+Cv*yf*F-7$$\"#vF*$\"+K6+-'*F-7$$\"#wF*$\"+sI(fg*F-7$$\"
#xF*$\"+KH#)4'*F-7$$\"#yF*$\"+Wxc8'*F-7$$\"#zF*$\"+[Q><'*F-7$$\"#!)F*$
\"++Wr?'*F-7$$\"#\")F*$\"+C@9C'*F-7$$\"##)F*$\"++1YF'*F-7$$\"#$)F*$\"+
/\\oI'*F-7$$\"#%)F*$\"+;?$Qj*F-7$$\"#&)F*$\"+#*f)oj*F-7$$\"#')F*$\"+g_
%)R'*F-7$$\"#()F*$\"+3mtU'*F-7$$\"#))F*$\"+!o\\bk*F-7$$\"#*)F*$\"+)Gy#
['*F-7$$\"#!*F*$\"+w8%4l*F-7$$\"#\"*F*$\"+;&HNl*F-7$$\"##*F*$\"+Sq0c'*
F-7$$\"#$*F*$\"+go^e'*F-7$$\"#%*F*$\"+gY!4m*F-7$$\"#&*F*$\"+cTCj'*F-7$
$\"#'*F*$\"+)eBbm*F-" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 
45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 19 "To plot C(T) we do
:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "plot([[T1[j],C1[j+1]] \+
$j=1..95],view=[0..3,0..1]);" }}{PARA 13 "" 1 "" {GLPLOT2D 706 283 
283 {PLOTDATA 2 "6%-%'CURVESG6$7[q7$$\"1+++-%Gpz#!#;$\"1*****R?N=O&F*7
$$\"1+++QgO$G$F*$\"1+++W!>u$fF*7$$\"1+++3-vyOF*$\"1+++/\"QxP'F*7$$\"1+
++o#z$HSF*$\"1+++/QVInF*7$$\"1+++/OW_VF*$\"1+++WlC@qF*7$$\"1+++q3ccYF*
$\"1+++)Q,gE(F*7$$\"1+++O*en%\\F*$\"1+++C-AvuF*7$$\"1+++y-EE_F*$\"1+++
s>GcwF*7$$\"1*****>0os\\&F*$\"1+++#*oc9yF*7$$\"1+++q:OhdF*$\"1+++s#>T&
zF*7$$\"1+++][r>gF*$\"1+++sB0y!)F*7$$\"1+++]hAtiF*$\"1+++3(3))=)F*7$$
\"1+++O*)fAlF*$\"1+++S(Q$)G)F*7$$\"1+++7]RonF*$\"1+++S3Ay$)F*7$$\"1+++
w/26qF*$\"1+++7%[(f%)F*7$$\"1+++q2+^sF*$\"1+++cC*R`)F*7$$\"1+++3))\\)[
(F*$\"1+++Sp%=g)F*7$$\"1******p!GQs(F*$\"1+++ou1k')F*7$$\"1+++/E@dzF*$
\"1+++;eH@()F*7$$\"1*****>[W))=)F*$\"1+++'RxSx)F*7$$\"1+++E&*))=%)F*$
\"1+++%y#)G#))F*7$$\"1+++A>\\Z')F*$\"1+++O=7o))F*7$$\"1+++9zxu))F*$\"1
+++#*)\\,\"*)F*7$$\"1+++![e35*F*$\"1+++'\\v#\\*)F*7$$\"1+++3<$eK*F*$\"
1,+++$pd)*)F*7$$\"1+++5Zy\\&*F*$\"1+++_D()>!*F*7$$\"1+++k_zs(*F*$\"1,+
+cqz^!*F*7$$\"1+++=I$\\***F*$\"1+++o4t\"3*F*7$$\"1+++^gi@5!#:$\"1+++O?
%)4\"*F*7$$\"1+++AMoV5Fat$\"1+++G:GO\"*F*7$$\"1+++=0nl5Fat$\"1++++;=h
\"*F*7$$\"1+++x>f(3\"Fat$\"1+++C:m%=*F*7$$\"1+++1?X46Fat$\"1,++?f$o?*F
*7$$\"1+++nWDJ6Fat$\"1+++?jzF#*F*7$$\"1+++wG+`6Fat$\"1+++GHjZ#*F*7$$\"
1+++Y/qu6Fat$\"1+++SZVm#*F*7$$\"1+++q,N'>\"Fat$\"1*****f4lUG*F*7$$\"1+
++eZ&z@\"Fat$\"1+++;>>,$*F*7$$\"1+++=n^R7Fat$\"1,++OEG<$*F*7$$\"1+++0%
Q5E\"Fat$\"1+++W$)eK$*F*7$$\"1+++f>_#G\"Fat$\"1+++o3;Z$*F*7$$\"1+++x$p
RI\"Fat$\"1+++%=\\5O*F*7$$\"1+++BDQD8Fat$\"1+++cYHu$*F*7$$\"1+++?JwY8F
at$\"1+++3d$pQ*F*7$$\"1+++wF6o8Fat$\"1+++k9,*R*F*7$$\"1+++(*HV*Q\"Fat$
\"1+++?Yb5%*F*7$$\"1+++\">D2T\"Fat$\"1+++_xf@%*F*7$$\"1+++k1*>V\"Fat$
\"1+++WI<K%*F*7$$\"1+++c1B`9Fat$\"1++++?IU%*F*7$$\"1+++3jWu9Fat$\"1+++
sv+_%*F*7$$\"1+++*pQc\\\"Fat$\"1+++g(>8Y*F*7$$\"1+++f)3o^\"Fat$\"1+++?
zDq%*F*7$$\"1+++^x&z`\"Fat$\"1+++O\\$)y%*F*7$$\"1+++fi3f:Fat$\"1+++;B3
([*F*7$$\"1+++\\_>!e\"Fat$\"1+++gZ,&\\*F*7$$\"1+++QbG,;Fat$\"1+++'>PE]
*F*7$$\"1+++VyNA;Fat$\"1*****Rew*4&*F*7$$\"1+++3HTV;Fat$\"1+++3s/<&*F*
7$$\"1+++@9Xk;Fat$\"1+++?W&Q_*F*7$$\"1++++SZ&o\"Fat$\"1+++sjTI&*F*7$$
\"1+++k7[1<Fat$\"1+++#R[n`*F*7$$\"1+++-QZF<Fat$\"1+++!y\\Ga*F*7$$\"1++
+[@X[<Fat$\"1+++%)Gu[&*F*7$$\"1+++>oTp<Fat$\"1+++?*RWb*F*7$$\"1+++[$o.
z\"Fat$\"1+++o7$*f&*F*7$$\"1+++]rI6=Fat$\"1+++C-Cl&*F*7$$\"1+++,PBK=Fa
t$\"1*****f:x.d*F*7$$\"1+++b%[J&=Fat$\"1+++OsMv&*F*7$$\"1+++$y^S(=Fat$
\"1+++%ec,e*F*7$$\"1+++\"4W\\*=Fat$\"1+++![/[e*F*7$$\"1+++Dd#e\">Fat$
\"1+++3KJ*e*F*7$$\"1+++iqpO>Fat$\"1+++;mn$f*F*7$$\"1+++\\%ev&>Fat$\"1+
++Cv*yf*F*7$$\"1+++k,Ty>Fat$\"1+++K6+-'*F*7$$\"1+++zDD**>Fat$\"1+++sI(
fg*F*7$$\"1+++ef3??Fat$\"1*****>$H#)4'*F*7$$\"1+++i0\"4/#Fat$\"1+++Wxc
8'*F*7$$\"1+++;nsh?Fat$\"1+++[Q><'*F*7$$\"1+++QY`#3#Fat$\"1++++Wr?'*F*
7$$\"1+++yXL.@Fat$\"1+++C@9C'*F*7$$\"1+++7o7C@Fat$\"1++++1YF'*F*7$$\"1
+++K:\"\\9#Fat$\"1+++/\\oI'*F*7$$\"1+++e*)ol@Fat$\"1+++;?$Qj*F*7$$\"1+
++_$fk=#Fat$\"1,++#*f)oj*F*7$$\"1+++oGA2AFat$\"1+++g_%)R'*F*7$$\"1+++%
pzzA#Fat$\"1+++3mtU'*F*7$$\"1+++_+t[AFat$\"1+++!o\\bk*F*7$$\"1+++OTZpA
Fat$\"1+++)Gy#['*F*7$$\"1+++_?@!H#Fat$\"1+++w8%4l*F*7$$\"1+++QS%4J#Fat
$\"1+++;&HNl*F*7$$\"1+++E-nJBFat$\"1+++Sq0c'*F*7$$\"1+++s2R_BFat$\"1++
+go^e'*F*7$$\"1+++[e5tBFat$\"1+++gY!4m*F*7$$\"1+++[b\"QR#Fat$\"1+++cTC
j'*F*7$$\"1+++u+_9CFat$\"1+++)eBbm*F*-%'COLOURG6&%$RGBG$\"#5!\"\"\"\"!
Feil-%+AXESLABELSG6$%!GFiil-%%VIEWG6$;Feil$\"\"$Feil;Feil$\"\"\"Feil" 
1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 
"" 0 "" {TEXT -1 260 "Note that dimensionful constants such as k, and \+
the energy quantum E have been set to unity. The temperature unit is a
ctually E/k. The heat capacity is constant above E/k=1, and drops to z
ero below. At high temperatures the specific heat capacity approaches \+
k." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 75 "Can
 we investigate the low-temperature behaviour? We need many oscillator
s." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 21 "N:=1000000; q_m:=200;" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%\"NG\"(+++\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$q_
mG\"$+#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "O1:=[]:" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "for q from 0 to q_m do:" }}
{PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "o1:=Omega(N,q);" }}{PARA 0 "> " 0 "
" {MPLTEXT 1 0 16 "O1:=ladd(O1,o1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 
3 "od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "S1:=map(evalf,map
(log,O1)):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 37 "Now we can calculat
e the temperature:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "T1:=[]
:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "for q from 2 to q_m-1 \+
do:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "t1:=2/(S1[q+1]-S1[q-1]);" }}
{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "T1:=ladd(T1,t1); od:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "listplot(T1);" }}{PARA 13 "" 1 "" 
{GLPLOT2D 695 227 227 {PLOTDATA 2 "6#-%'CURVESG6#7bw7$$\"\"\"\"\"!$\"+
E1\\Cu!#67$$\"\"#F*$\"+_$e,u(F-7$$\"\"$F*$\"+wM^`zF-7$$\"\"%F*$\"+[IV=
\")F-7$$\"\"&F*$\"+ssGa#)F-7$$\"\"'F*$\"++b_q$)F-7$$\"\"(F*$\"+]i`s%)F
-7$$\"\")F*$\"+#z2Pc)F-7$$\"\"*F*$\"+m)>jk)F-7$$\"#5F*$\"+g])>s)F-7$$
\"#6F*$\"+y)))=z)F-7$$\"#7F*$\"+)*)Hp&))F-7$$\"#8F*$\"++m!y\"*)F-7$$\"
#9F*$\"+'Gt](*)F-7$$\"#:F*$\"+_x<H!*F-7$$\"#;F*$\"+or[!3*F-7$$\"#<F*$
\"+SlIH\"*F-7$$\"#=F*$\"+%G#*e<*F-7$$\"#>F*$\"+u?Y?#*F-7$$\"#?F*$\"+#Q
-KE*F-7$$\"#@F*$\"++UF/$*F-7$$\"#AF*$\"+kv\"QM*F-7$$\"#BF*$\"+5_&>Q*F-
7$$\"#CF*$\"+g^z=%*F-7$$\"#DF*$\"+IIVa%*F-7$$\"#EF*$\"+?S&*)[*F-7$$\"#
FF*$\"+-VVA&*F-7$$\"#GF*$\"+YB%\\b*F-7$$\"#HF*$\"+y(Rle*F-7$$\"#IF*$\"
+5CG<'*F-7$$\"#JF*$\"+34AZ'*F-7$$\"#KF*$\"+W9Sw'*F-7$$\"#LF*$\"+#Gm[q*
F-7$$\"#MF*$\"+sRlK(*F-7$$\"#NF*$\"+[+!)f(*F-7$$\"#OF*$\"++sL'y*F-7$$
\"#PF*$\"+UcH7)*F-7$$\"#QF*$\"+ELqP)*F-7$$\"#RF*$\"+ogei)*F-7$$\"#SF*$
\"+Kz'p))*F-7$$\"#TF*$\"+O8(3\"**F-7$$\"#UF*$\"+yqJM**F-7$$\"#VF*$\"+)
pCt&**F-7$$\"#WF*$\"+KC\"*z**F-7$$\"#XF*$\"+M(4-+\"!#57$$\"#YF*$\"+T&*
Q-5Fdy7$$\"#ZF*$\"+t@`/5Fdy7$$\"#[F*$\"+Y!Rm+\"Fdy7$$\"#\\F*$\"+&\\6(3
5Fdy7$$\"#]F*$\"+v2v55Fdy7$$\"#^F*$\"+y!eF,\"Fdy7$$\"#_F*$\"+MXt95Fdy7
$$\"#`F*$\"+#>\"o;5Fdy7$$\"#aF*$\"+k!*f=5Fdy7$$\"#bF*$\"+5\"*[?5Fdy7$$
\"#cF*$\"+RANA5Fdy7$$\"#dF*$\"+8$*=C5Fdy7$$\"#eF*$\"+Z6+E5Fdy7$$\"#fF*
$\"+>&)yF5Fdy7$$\"#gF*$\"+y@bH5Fdy7$$\"#hF*$\"+LGHJ5Fdy7$$\"#iF*$\"+j6
,L5Fdy7$$\"#jF*$\"+4yqM5Fdy7$$\"#kF*$\"+*Q$QO5Fdy7$$\"#lF*$\"++&Q!Q5Fd
y7$$\"#mF*$\"+/PnR5Fdy7$$\"#nF*$\"+N&*GT5Fdy7$$\"#oF*$\"+Bl)G/\"Fdy7$$
\"#pF*$\"+k^YW5Fdy7$$\"#qF*$\"+Df-Y5Fdy7$$\"#rF*$\"+s#pv/\"Fdy7$$\"#sF
*$\"+_c4\\5Fdy7$$\"#tF*$\"+'[010\"Fdy7$$\"#uF*$\"+n\"*4_5Fdy7$$\"#vF*$
\"+)4xN0\"Fdy7$$\"#wF*$\"+j'R]0\"Fdy7$$\"#xF*$\"+7s[c5Fdy7$$\"#yF*$\"+
*4?z0\"Fdy7$$\"#zF*$\"+f'Q$f5Fdy7$$\"#!)F*$\"+GKug5Fdy7$$\"#\")F*$\"+=
T8i5Fdy7$$\"##)F*$\"+G;^j5Fdy7$$\"#$)F*$\"+Xg([1\"Fdy7$$\"#%)F*$\"+gwA
m5Fdy7$$\"#&)F*$\"+ancn5Fdy7$$\"#')F*$\"+uN*)o5Fdy7$$\"#()F*$\"+*Q3-2
\"Fdy7$$\"#))F*$\"+Z9^r5Fdy7$$\"#*)F*$\"+!)H!G2\"Fdy7$$\"#!*F*$\"+JK3u
5Fdy7$$\"#\"*F*$\"+<CNv5Fdy7$$\"##*F*$\"+b2hw5Fdy7$$\"#$*F*$\"+r%ey2\"
Fdy7$$\"#%*F*$\"+bd4z5Fdy7$$\"#&*F*$\"+0GK!3\"Fdy7$$\"#'*F*$\"+N)R:3\"
Fdy7$$\"#(*F*$\"+8qu#3\"Fdy7$$\"#)*F*$\"+0X%R3\"Fdy7$$\"#**F*$\"+bD8&3
\"Fdy7$$\"$+\"F*$\"+)G6j3\"Fdy7$$\"$,\"F*$\"+=3[(3\"Fdy7$$\"$-\"F*$\"+
J9k)3\"Fdy7$$\"$.\"F*$\"+#=$z*3\"Fdy7$$\"$/\"F*$\"+Xi$44\"Fdy7$$\"$0\"
F*$\"+'zq?4\"Fdy7$$\"$1\"F*$\"+8q>$4\"Fdy7$$\"$2\"F*$\"+d\\J%4\"Fdy7$$
\"$3\"F*$\"+5[U&4\"Fdy7$$\"$4\"F*$\"+(pEl4\"Fdy7$$\"$5\"F*$\"+U2i(4\"F
dy7$$\"$6\"F*$\"+Nrq)4\"Fdy7$$\"$7\"F*$\"+Vfy*4\"Fdy7$$\"$8\"F*$\"+)Hd
35\"Fdy7$$\"$9\"F*$\"+67#>5\"Fdy7$$\"$:\"F*$\"+'*z(H5\"Fdy7$$\"$;\"F*$
\"+qw-/6Fdy7$$\"$<\"F*$\"+1.206Fdy7$$\"$=\"F*$\"+/h516Fdy7$$\"$>\"F*$
\"+\"3Nr5\"Fdy7$$\"$?\"F*$\"+zt:36Fdy7$$\"$@\"F*$\"+RJ<46Fdy7$$\"$A\"F
*$\"+#Q#=56Fdy7$$\"$B\"F*$\"+!>&=66Fdy7$$\"$C\"F*$\"+5<=76Fdy7$$\"$D\"
F*$\"+D?<86Fdy7$$\"$E\"F*$\"+'GcT6\"Fdy7$$\"$F\"F*$\"+aW8:6Fdy7$$\"$G
\"F*$\"+\"o1h6\"Fdy7$$\"$H\"F*$\"+<J2<6Fdy7$$\"$I\"F*$\"+FP.=6Fdy7$$\"
$J\"F*$\"+/'))*=6Fdy7$$\"$K\"F*$\"+-z$*>6Fdy7$$\"$L\"F*$\"+];)37\"Fdy7
$$\"$M\"F*$\"+z)>=7\"Fdy7$$\"$N\"F*$\"+5GvA6Fdy7$$\"$O\"F*$\"+].oB6Fdy
7$$\"$P\"F*$\"+%f-Y7\"Fdy7$$\"$Q\"F*$\"+l(>b7\"Fdy7$$\"$R\"F*$\"+4<VE6
Fdy7$$\"$S\"F*$\"+]'Qt7\"Fdy7$$\"$T\"F*$\"+C1CG6Fdy7$$\"$U\"F*$\"+0w8H
6Fdy7$$\"$V\"F*$\"+>)H+8\"Fdy7$$\"$W\"F*$\"+yr\"48\"Fdy7$$\"$X\"F*$\"+
#y*zJ6Fdy7$$\"$Y\"F*$\"+OxnK6Fdy7$$\"$Z\"F*$\"+x5bL6Fdy7$$\"$[\"F*$\"+
Z)>W8\"Fdy7$$\"$\\\"F*$\"+%3%GN6Fdy7$$\"$]\"F*$\"+JQ9O6Fdy7$$\"$^\"F*$
\"+$>**p8\"Fdy7$$\"$_\"F*$\"+7-&y8\"Fdy7$$\"$`\"F*$\"+IppQ6Fdy7$$\"$a
\"F*$\"+b%R&R6Fdy7$$\"$b\"F*$\"+,xPS6Fdy7$$\"$c\"F*$\"+T=@T6Fdy7$$\"$d
\"F*$\"+b=/U6Fdy7$$\"$e\"F*$\"+)ynG9\"Fdy7$$\"$f\"F*$\"+\\(*oV6Fdy7$$
\"$g\"F*$\"+>x]W6Fdy7$$\"$h\"F*$\"+5=KX6Fdy7$$\"$i\"F*$\"+P>8Y6Fdy7$$
\"$j\"F*$\"+x#Qp9\"Fdy7$$\"$k\"F*$\"+z3uZ6Fdy7$$\"$l\"F*$\"+e'R&[6Fdy7
$$\"$m\"F*$\"+%zM$\\6Fdy7$$\"$n\"F*$\"+.i7]6Fdy7$$\"$o\"F*$\"+NR\"4:\"
Fdy7$$\"$p\"F*$\"+q\")p^6Fdy7$$\"$q\"F*$\"+E)yC:\"Fdy7$$\"$r\"F*$\"+`f
D`6Fdy7$$\"$s\"F*$\"++'HS:\"Fdy7$$\"$t\"F*$\"+>)*za6Fdy7$$\"$u\"F*$\"+
$fmb:\"Fdy7$$\"$v\"F*$\"+S+Lc6Fdy7$$\"$w\"F*$\"+X,4d6Fdy7$$\"$x\"F*$\"
+$*o%y:\"Fdy7$$\"$y\"F*$\"+//ge6Fdy7$$\"$z\"F*$\"+H2Nf6Fdy7$$\"$!=F*$
\"+(y(4g6Fdy7$$\"$\"=F*$\"+J;%3;\"Fdy7$$\"$#=F*$\"+\\Ceh6Fdy7$$\"$$=F*
$\"+$4?B;\"Fdy7$$\"$%=F*$\"+;Y0j6Fdy7$$\"$&=F*$\"+Thyj6Fdy7$$\"$'=F*$
\"+@Z^k6Fdy7$$\"$(=F*$\"+4-Cl6Fdy7$$\"$)=F*$\"+eE'f;\"Fdy7$$\"$*=F*$\"
+HBom6Fdy7$$\"$!>F*$\"+t!*Rn6Fdy7$$\"$\">F*$\"+yG6o6Fdy7$$\"$#>F*$\"+*
zB)o6Fdy7$$\"$$>F*$\"+h>`p6Fdy7$$\"$%>F*$\"+`tBq6Fdy7$$\"$&>F*$\"+&*)R
4<\"Fdy7$$\"$'>F*$\"+6(R;<\"Fdy7$$\"$(>F*$\"+goLs6Fdy7$$\"$)>F*$\"+h7.
t6Fdy" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 41 "We have clearly reached low temperatures.
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "C1:=[]:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "for q from 2 to q_m-3 do:" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 26 "c1:=2/(T1[q+1]-T1[q-1])/N;" }}{PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 20 "C1:=ladd(C1,c1); od:" }}}{EXCHG {PARA 0 "" 0 "
" {TEXT -1 100 "The simple listplot graphs the specific heat capacity \+
(at fixed volume) as a function of energy (q)." }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 13 "listplot(C1);" }}{PARA 13 "" 1 "" {GLPLOT2D 
711 243 243 {PLOTDATA 2 "6#-%'CURVESG6#7`w7$$\"\"\"\"\"!$\"+O]b!y$!#87
$$\"\"#F*$\"+iI;(G&F-7$$\"\"$F*$\"+Qa^\\mF-7$$\"\"%F*$\"+9tfLzF-7$$\"
\"&F*$\"+!*e%Q;*F-7$$\"\"'F*$\"+y:HN5!#77$$\"\"(F*$\"+xk&3:\"FG7$$\"\"
)F*$\"+L`gj7FG7$$\"\"*F*$\"+%y=RP\"FG7$$\"#5F*$\"+tq3#[\"FG7$$\"#6F*$
\"+t(Q$)e\"FG7$$\"#7F*$\"+t\"eGp\"FG7$$\"#8F*$\"+urz&z\"FG7$$\"#9F*$\"
+?LG(*=FG7$$\"#:F*$\"+`vU(*>FG7$$\"#;F*$\"+SOK'4#FG7$$\"#<F*$\"++A0%>#
FG7$$\"#=F*$\"+!y&o!H#FG7$$\"#>F*$\"+W&*G'Q#FG7$$\"#?F*$\"+mA#4[#FG7$$
\"#@F*$\"+==juDFG7$$\"#AF*$\"+#\\jum#FG7$$\"#BF*$\"+'*[YfFFG7$$\"#CF*$
\"+eMn]GFG7$$\"#DF*$\"+a97THFG7$$\"#EF*$\"+A\"R3.$FG7$$\"#FF*$\"+y#e)>
JFG7$$\"#GF*$\"+K(4#3KFG7$$\"#HF*$\"+Gv\"fH$FG7$$\"#IF*$\"+5_+$Q$FG7$$
\"#JF*$\"+c\"*[pMFG7$$\"#KF*$\"+odRbNFG7$$\"#LF*$\"+C![2k$FG7$$\"#MF*$
\"+#Qebs$FG7$$\"#NF*$\"+YN%)4QFG7$$\"#OF*$\"+%e;O*QFG7$$\"#PF*$\"+'4+p
(RFG7$$\"#QF*$\"+C(4(fSFG7$$\"#RF*$\"+=,0UTFG7$$\"#SF*$\"+=5%RA%FG7$$
\"#TF*$\"+SMR0VFG7$$\"#UF*$\"+E=T'Q%FG7$$\"#VF*$\"+%e8qY%FG7$$\"#WF*$
\"+)p8sa%FG7$$\"#XF*$\"+#H<qi%FG7$$\"#YF*$\"+90V1ZFG7$$\"#ZF*$\"+-_Y&y
%FG7$$\"#[F*$\"+)zLT'[FG7$$\"#\\F*$\"+i:WU\\FG7$$\"#]F*$\"+9SQ?]FG7$$
\"#^F*$\"+M5*z4&FG7$$\"#_F*$\"+yNFv^FG7$$\"#`F*$\"+SR@__FG7$$\"#aF*$\"
+)f@)G`FG7$$\"#bF*$\"+a)3^S&FG7$$\"#cF*$\"+y=4\"[&FG7$$\"#dF*$\"+M6ycb
FG7$$\"#eF*$\"+3r;KcFG7$$\"#fF*$\"+?6D2dFG7$$\"#gF*$\"+)*H/#y&FG7$$\"#
hF*$\"+KlbceFG7$$\"#iF*$\"+/1!3$fFG7$$\"#jF*$\"+u*fZ+'FG7$$\"#kF*$\"+w
`WygFG7$$\"#lF*$\"+9q)=:'FG7$$\"#mF*$\"+)yi]A'FG7$$\"#nF*$\"+3\\'zH'FG
7$$\"#oF*$\"+9+kqjFG7$$\"#pF*$\"+!zuIW'FG7$$\"#qF*$\"+[tB:lFG7$$\"#rF*
$\"+s4:(e'FG7$$\"#sF*$\"+'pk)emFG7$$\"#tF*$\"+gPNInFG7$$\"#uF*$\"+7Od,
oFG7$$\"#vF*$\"+)*HesoFG7$$\"#wF*$\"+#H$RVpFG7$$\"#xF*$\"+EF)R,(FG7$$
\"#yF*$\"+#)zL%3(FG7$$\"#zF*$\"+W,YarFG7$$\"#!)F*$\"+S&*QCsFG7$$\"#\")
F*$\"+;-9%H(FG7$$\"##)F*$\"+#\\$ojtFG7$$\"#$)F*$\"+O^)HV(FG7$$\"#%)F*$
\"+su6-vFG7$$\"#&)F*$\"+w'o5d(FG7$$\"#')F*$\"+7+zRwFG7$$\"#()F*$\"+7SM
3xFG7$$\"#))F*$\"+kGqwxFG7$$\"#*)F*$\"+y#p[%yFG7$$\"#!*F*$\"+IB)G\"zFG
7$$\"#\"*F*$\"++Sn!)zFG7$$\"##*F*$\"+Q(*G[!)FG7$$\"#$*F*$\"+1Iy:\")FG7
$$\"#%*F*$\"+)>QI=)FG7$$\"#&*F*$\"+_r6]#)FG7$$\"#'*F*$\"+Y&\\rJ)FG7$$
\"#(*F*$\"+w*QQQ)FG7$$\"#)*F*$\"+_))G]%)FG7$$\"#**F*$\"+!eLp^)FG7$$\"$
+\"F*$\"+5U;$e)FG7$$\"$,\"F*$\"+7t:\\')FG7$$\"$-\"F*$\"+)G.`r)FG7$$\"$
.\"F*$\"+)=>6y)FG7$$\"$/\"F*$\"+g7bY))FG7$$\"$0\"F*$\"+il,7*)FG7$$\"$1
\"F*$\"+['yu(*)FG7$$\"$2\"F*$\"++UlU!*FG7$$\"$3\"F*$\"+[,v2\"*FG7$$\"$
4\"F*$\"+=!pC<*FG7$$\"$5\"F*$\"+=x,P#*FG7$$\"$6\"F*$\"+)Q?;I*FG7$$\"$7
\"F*$\"+]d]m$*FG7$$\"$8\"F*$\"+)>d3V*FG7$$\"$9\"F*$\"+Eog%\\*FG7$$\"$:
\"F*$\"+')H\")e&*FG7$$\"$;\"F*$\"+W$HEi*FG7$$\"$<\"F*$\"+IaG'o*FG7$$\"
$=\"F*$\"+;N.](*FG7$$\"$>\"F*$\"+5ME8)*FG7$$\"$?\"F*$\"+y&Gl()*FG7$$\"
$@\"F*$\"+#*e5S**FG7$$\"$A\"F*$\"+8hL+5!#67$$\"$B\"F*$\"+;hi15Ffam7$$
\"$C\"F*$\"+wc(G,\"Ffam7$$\"$D\"F*$\"+^^9>5Ffam7$$\"$E\"F*$\"+5LVD5Ffa
m7$$\"$F\"F*$\"+:)Q;.\"Ffam7$$\"$G\"F*$\"+Wz&y.\"Ffam7$$\"$H\"F*$\"+&*
37W5Ffam7$$\"$I\"F*$\"+JOK]5Ffam7$$\"$J\"F*$\"+j$)\\c5Ffam7$$\"$K\"F*$
\"+WAri5Ffam7$$\"$L\"F*$\"+$ob)o5Ffam7$$\"$M\"F*$\"++m*\\2\"Ffam7$$\"$
N\"F*$\"+>t?\"3\"Ffam7$$\"$O\"F*$\"+8AI(3\"Ffam7$$\"$P\"F*$\"+i\\U$4\"
Ffam7$$\"$Q\"F*$\"+RZd*4\"Ffam7$$\"$R\"F*$\"+Y^j06Ffam7$$\"$S\"F*$\"+[
lv66Ffam7$$\"$T\"F*$\"+V6#y6\"Ffam7$$\"$U\"F*$\"+2^'Q7\"Ffam7$$\"$V\"F
*$\"+Cr'*H6Ffam7$$\"$W\"F*$\"+sM+O6Ffam7$$\"$X\"F*$\"+a?,U6Ffam7$$\"$Y
\"F*$\"+tC.[6Ffam7$$\"$Z\"F*$\"+GM1a6Ffam7$$\"$[\"F*$\"+uN5g6Ffam7$$\"
$\\\"F*$\"+Wm5m6Ffam7$$\"$]\"F*$\"+&=r?<\"Ffam7$$\"$^\"F*$\"+m(Q!y6Ffa
m7$$\"$_\"F*$\"+%GiR=\"Ffam7$$\"$`\"F*$\"+)4I**=\"Ffam7$$\"$a\"F*$\"+j
h*e>\"Ffam7$$\"$b\"F*$\"+_=\"=?\"Ffam7$$\"$c\"F*$\"+B\"ox?\"Ffam7$$\"$
d\"F*$\"+(HqO@\"Ffam7$$\"$e\"F*$\"+1Nc>7Ffam7$$\"$f\"F*$\"+#RWaA\"Ffam
7$$\"$g\"F*$\"++=OJ7Ffam7$$\"$h\"F*$\"+6iEP7Ffam7$$\"$i\"F*$\"+))G0V7F
fam7$$\"$j\"F*$\"+C?#*[7Ffam7$$\"$k\"F*$\"+b?xa7Ffam7$$\"$l\"F*$\"+l3g
g7Ffam7$$\"$m\"F*$\"+]8^m7Ffam7$$\"$n\"F*$\"+R8Hs7Ffam7$$\"$o\"F*$\"+'
GW!y7Ffam7$$\"$p\"F*$\"+ne(QG\"Ffam7$$\"$q\"F*$\"+t&y'*G\"Ffam7$$\"$r
\"F*$\"+;#\\aH\"Ffam7$$\"$s\"F*$\"+,@C,8Ffam7$$\"$t\"F*$\"+`2+28Ffam7$
$\"$u\"F*$\"+]=s78Ffam7$$\"$v\"F*$\"+H$>&=8Ffam7$$\"$w\"F*$\"+nhFC8Ffa
m7$$\"$x\"F*$\"+`A$*H8Ffam7$$\"$y\"F*$\"+@OmN8Ffam7$$\"$z\"F*$\"+&zp9M
\"Ffam7$$\"$!=F*$\"+hz5Z8Ffam7$$\"$\"=F*$\"+Arv_8Ffam7$$\"$#=F*$\"+C:a
e8Ffam7$$\"$$=F*$\"+8=@k8Ffam7$$\"$%=F*$\"+.Xwp8Ffam7$$\"$&=F*$\"+#\\^
aP\"Ffam7$$\"$'=F*$\"+nbF\"Q\"Ffam7$$\"$(=F*$\"+Ur%oQ\"Ffam7$$\"$)=F*$
\"+t_N#R\"Ffam7$$\"$*=F*$\"+fY1)R\"Ffam7$$\"$!>F*$\"+\"pyPS\"Ffam7$$\"
$\">F*$\"+M3O49Ffam7$$\"$#>F*$\"+KM([T\"Ffam7$$\"$$>F*$\"+K6_?9Ffam7$$
\"$%>F*$\"+_$phU\"Ffam7$$\"$&>F*$\"+'4v;V\"Ffam7$$\"$'>F*$\"+OiCP9Ffam
" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 
0 "" 0 "" {TEXT -1 19 "To plot C(T) we do:" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 58 "plot([[T1[j],C1[j+1]] $j=1..q_m-5],view=[0..0.2,0..
0.02]);" }}{PARA 13 "" 1 "" {GLPLOT2D 706 283 283 {PLOTDATA 2 "6%-%'CU
RVESG6$7_w7$$\"1*****fi!\\Cu!#<$\"1+++iI;(G&!#>7$$\"1+++_$e,u(F*$\"1++
+Qa^\\mF-7$$\"1+++wM^`zF*$\"1+++9tfLzF-7$$\"1+++[IV=\")F*$\"1+++!*e%Q;
*F-7$$\"1*****>F(Ga#)F*$\"1+++y:HN5!#=7$$\"1,+++b_q$)F*$\"1+++xk&3:\"F
B7$$\"1+++]i`s%)F*$\"1+++L`gj7FB7$$\"1+++#z2Pc)F*$\"1+++%y=RP\"FB7$$\"
1+++m)>jk)F*$\"1+++tq3#[\"FB7$$\"1+++g])>s)F*$\"1+++t(Q$)e\"FB7$$\"1++
+y)))=z)F*$\"1+++t\"eGp\"FB7$$\"1+++)*)Hp&))F*$\"1+++urz&z\"FB7$$\"1++
++m!y\"*)F*$\"1+++?LG(*=FB7$$\"1+++'Gt](*)F*$\"1+++`vU(*>FB7$$\"1,++_x
<H!*F*$\"1+++SOK'4#FB7$$\"1+++or[!3*F*$\"1++++A0%>#FB7$$\"1+++SlIH\"*F
*$\"1+++!y&o!H#FB7$$\"1,++%G#*e<*F*$\"1+++W&*G'Q#FB7$$\"1,++u?Y?#*F*$
\"1+++mA#4[#FB7$$\"1,++#Q-KE*F*$\"1+++==juDFB7$$\"1++++UF/$*F*$\"1+++#
\\jum#FB7$$\"1+++kv\"QM*F*$\"1+++'*[YfFFB7$$\"1******4_&>Q*F*$\"1+++eM
n]GFB7$$\"1+++g^z=%*F*$\"1+++a97THFB7$$\"1,++IIVa%*F*$\"1+++A\"R3.$FB7
$$\"1+++?S&*)[*F*$\"1+++y#e)>JFB7$$\"1+++-VVA&*F*$\"1+++K(4#3KFB7$$\"1
+++YB%\\b*F*$\"1+++Gv\"fH$FB7$$\"1+++y(Rle*F*$\"1+++5_+$Q$FB7$$\"1+++5
CG<'*F*$\"1+++c\"*[pMFB7$$\"1+++34AZ'*F*$\"1+++odRbNFB7$$\"1+++W9Sw'*F
*$\"1+++C![2k$FB7$$\"1,++#Gm[q*F*$\"1+++#Qebs$FB7$$\"1*****>(RlK(*F*$
\"1+++YN%)4QFB7$$\"1+++[+!)f(*F*$\"1+++%e;O*QFB7$$\"1++++sL'y*F*$\"1++
+'4+p(RFB7$$\"1+++UcH7)*F*$\"1+++C(4(fSFB7$$\"1+++ELqP)*F*$\"1+++=,0UT
FB7$$\"1+++ogei)*F*$\"1+++=5%RA%FB7$$\"1+++Kz'p))*F*$\"1+++SMR0VFB7$$
\"1+++O8(3\"**F*$\"1+++E=T'Q%FB7$$\"1+++yqJM**F*$\"1+++%e8qY%FB7$$\"1+
++)pCt&**F*$\"1+++)p8sa%FB7$$\"1+++KC\"*z**F*$\"1+++#H<qi%FB7$$\"1+++M
(4-+\"!#;$\"1+++90V1ZFB7$$\"1+++T&*Q-5Fcy$\"1+++-_Y&y%FB7$$\"1+++t@`/5
Fcy$\"1+++)zLT'[FB7$$\"1+++Y!Rm+\"Fcy$\"1+++i:WU\\FB7$$\"1+++&\\6(35Fc
y$\"1+++9SQ?]FB7$$\"1+++v2v55Fcy$\"1+++M5*z4&FB7$$\"1+++y!eF,\"Fcy$\"1
+++yNFv^FB7$$\"1+++MXt95Fcy$\"1+++SR@__FB7$$\"1+++#>\"o;5Fcy$\"1+++)f@
)G`FB7$$\"1+++k!*f=5Fcy$\"1+++a)3^S&FB7$$\"1+++5\"*[?5Fcy$\"1+++y=4\"[
&FB7$$\"1+++RANA5Fcy$\"1+++M6ycbFB7$$\"1+++8$*=C5Fcy$\"1+++3r;KcFB7$$
\"1+++Z6+E5Fcy$\"1+++?6D2dFB7$$\"1+++>&)yF5Fcy$\"1+++)*H/#y&FB7$$\"1++
+y@bH5Fcy$\"1+++KlbceFB7$$\"1+++LGHJ5Fcy$\"1+++/1!3$fFB7$$\"1+++j6,L5F
cy$\"1+++u*fZ+'FB7$$\"1+++4yqM5Fcy$\"1+++w`WygFB7$$\"1+++*Q$QO5Fcy$\"1
+++9q)=:'FB7$$\"1++++&Q!Q5Fcy$\"1+++)yi]A'FB7$$\"1+++/PnR5Fcy$\"1+++3
\\'zH'FB7$$\"1+++N&*GT5Fcy$\"1+++9+kqjFB7$$\"1+++Bl)G/\"Fcy$\"1+++!zuI
W'FB7$$\"1+++k^YW5Fcy$\"1+++[tB:lFB7$$\"1+++Df-Y5Fcy$\"1+++s4:(e'FB7$$
\"1+++s#pv/\"Fcy$\"1+++'pk)emFB7$$\"1+++_c4\\5Fcy$\"1+++gPNInFB7$$\"1+
++'[010\"Fcy$\"1+++7Od,oFB7$$\"1+++n\"*4_5Fcy$\"1+++)*HesoFB7$$\"1+++)
4xN0\"Fcy$\"1+++#H$RVpFB7$$\"1+++j'R]0\"Fcy$\"1+++EF)R,(FB7$$\"1+++7s[
c5Fcy$\"1+++#)zL%3(FB7$$\"1+++*4?z0\"Fcy$\"1+++W,YarFB7$$\"1+++f'Q$f5F
cy$\"1+++S&*QCsFB7$$\"1+++GKug5Fcy$\"1+++;-9%H(FB7$$\"1+++=T8i5Fcy$\"1
+++#\\$ojtFB7$$\"1+++G;^j5Fcy$\"1+++O^)HV(FB7$$\"1+++Xg([1\"Fcy$\"1+++
su6-vFB7$$\"1+++gwAm5Fcy$\"1+++w'o5d(FB7$$\"1+++ancn5Fcy$\"1+++7+zRwFB
7$$\"1+++uN*)o5Fcy$\"1+++7SM3xFB7$$\"1+++*Q3-2\"Fcy$\"1+++kGqwxFB7$$\"
1+++Z9^r5Fcy$\"1,++y#p[%yFB7$$\"1+++!)H!G2\"Fcy$\"1,++IB)G\"zFB7$$\"1+
++JK3u5Fcy$\"1++++Sn!)zFB7$$\"1+++<CNv5Fcy$\"1+++Q(*G[!)FB7$$\"1+++b2h
w5Fcy$\"1,++1Iy:\")FB7$$\"1+++r%ey2\"Fcy$\"1+++)>QI=)FB7$$\"1+++bd4z5F
cy$\"1,++_r6]#)FB7$$\"1+++0GK!3\"Fcy$\"1+++Y&\\rJ)FB7$$\"1+++N)R:3\"Fc
y$\"1,++w*QQQ)FB7$$\"1+++8qu#3\"Fcy$\"1+++_))G]%)FB7$$\"1+++0X%R3\"Fcy
$\"1+++!eLp^)FB7$$\"1+++bD8&3\"Fcy$\"1,++5U;$e)FB7$$\"1+++)G6j3\"Fcy$
\"1+++7t:\\')FB7$$\"1+++=3[(3\"Fcy$\"1*****zG.`r)FB7$$\"1+++J9k)3\"Fcy
$\"1+++)=>6y)FB7$$\"1+++#=$z*3\"Fcy$\"1+++g7bY))FB7$$\"1+++Xi$44\"Fcy$
\"1+++il,7*)FB7$$\"1+++'zq?4\"Fcy$\"1+++['yu(*)FB7$$\"1+++8q>$4\"Fcy$
\"1++++UlU!*FB7$$\"1+++d\\J%4\"Fcy$\"1,++[,v2\"*FB7$$\"1+++5[U&4\"Fcy$
\"1+++=!pC<*FB7$$\"1+++(pEl4\"Fcy$\"1+++=x,P#*FB7$$\"1+++U2i(4\"Fcy$\"
1+++)Q?;I*FB7$$\"1+++Nrq)4\"Fcy$\"1+++]d]m$*FB7$$\"1+++Vfy*4\"Fcy$\"1*
****z>d3V*FB7$$\"1+++)Hd35\"Fcy$\"1+++Eog%\\*FB7$$\"1+++67#>5\"Fcy$\"1
+++')H\")e&*FB7$$\"1+++'*z(H5\"Fcy$\"1*****RMHEi*FB7$$\"1+++qw-/6Fcy$
\"1+++IaG'o*FB7$$\"1+++1.206Fcy$\"1+++;N.](*FB7$$\"1+++/h516Fcy$\"1+++
5ME8)*FB7$$\"1+++\"3Nr5\"Fcy$\"1*****zdGl()*FB7$$\"1+++zt:36Fcy$\"1***
**>*e5S**FB7$$\"1+++RJ<46Fcy$\"1+++8hL+5F*7$$\"1+++#Q#=56Fcy$\"1+++;hi
15F*7$$\"1+++!>&=66Fcy$\"1+++wc(G,\"F*7$$\"1+++5<=76Fcy$\"1+++^^9>5F*7
$$\"1+++D?<86Fcy$\"1+++5LVD5F*7$$\"1+++'GcT6\"Fcy$\"1+++:)Q;.\"F*7$$\"
1+++aW8:6Fcy$\"1+++Wz&y.\"F*7$$\"1+++\"o1h6\"Fcy$\"1+++&*37W5F*7$$\"1+
++<J2<6Fcy$\"1+++JOK]5F*7$$\"1+++FP.=6Fcy$\"1+++j$)\\c5F*7$$\"1+++/'))
*=6Fcy$\"1+++WAri5F*7$$\"1+++-z$*>6Fcy$\"1+++$ob)o5F*7$$\"1+++];)37\"F
cy$\"1++++m*\\2\"F*7$$\"1+++z)>=7\"Fcy$\"1+++>t?\"3\"F*7$$\"1+++5GvA6F
cy$\"1+++8AI(3\"F*7$$\"1+++].oB6Fcy$\"1+++i\\U$4\"F*7$$\"1+++%f-Y7\"Fc
y$\"1+++RZd*4\"F*7$$\"1+++l(>b7\"Fcy$\"1+++Y^j06F*7$$\"1+++4<VE6Fcy$\"
1+++[lv66F*7$$\"1+++]'Qt7\"Fcy$\"1+++V6#y6\"F*7$$\"1+++C1CG6Fcy$\"1+++
2^'Q7\"F*7$$\"1+++0w8H6Fcy$\"1+++Cr'*H6F*7$$\"1+++>)H+8\"Fcy$\"1+++sM+
O6F*7$$\"1+++yr\"48\"Fcy$\"1+++a?,U6F*7$$\"1+++#y*zJ6Fcy$\"1+++tC.[6F*
7$$\"1+++OxnK6Fcy$\"1+++GM1a6F*7$$\"1+++x5bL6Fcy$\"1+++uN5g6F*7$$\"1++
+Z)>W8\"Fcy$\"1+++Wm5m6F*7$$\"1+++%3%GN6Fcy$\"1+++&=r?<\"F*7$$\"1+++JQ
9O6Fcy$\"1+++m(Q!y6F*7$$\"1+++$>**p8\"Fcy$\"1+++%GiR=\"F*7$$\"1+++7-&y
8\"Fcy$\"1+++)4I**=\"F*7$$\"1+++IppQ6Fcy$\"1+++jh*e>\"F*7$$\"1+++b%R&R
6Fcy$\"1+++_=\"=?\"F*7$$\"1+++,xPS6Fcy$\"1+++B\"ox?\"F*7$$\"1+++T=@T6F
cy$\"1+++(HqO@\"F*7$$\"1+++b=/U6Fcy$\"1+++1Nc>7F*7$$\"1+++)ynG9\"Fcy$
\"1+++#RWaA\"F*7$$\"1+++\\(*oV6Fcy$\"1++++=OJ7F*7$$\"1+++>x]W6Fcy$\"1+
++6iEP7F*7$$\"1+++5=KX6Fcy$\"1+++))G0V7F*7$$\"1+++P>8Y6Fcy$\"1+++C?#*[
7F*7$$\"1+++x#Qp9\"Fcy$\"1+++b?xa7F*7$$\"1+++z3uZ6Fcy$\"1+++l3gg7F*7$$
\"1+++e'R&[6Fcy$\"1+++]8^m7F*7$$\"1+++%zM$\\6Fcy$\"1+++R8Hs7F*7$$\"1++
+.i7]6Fcy$\"1+++'GW!y7F*7$$\"1+++NR\"4:\"Fcy$\"1+++ne(QG\"F*7$$\"1+++q
\")p^6Fcy$\"1+++t&y'*G\"F*7$$\"1+++E)yC:\"Fcy$\"1+++;#\\aH\"F*7$$\"1++
+`fD`6Fcy$\"1+++,@C,8F*7$$\"1++++'HS:\"Fcy$\"1+++`2+28F*7$$\"1+++>)*za
6Fcy$\"1+++]=s78F*7$$\"1+++$fmb:\"Fcy$\"1+++H$>&=8F*7$$\"1+++S+Lc6Fcy$
\"1+++nhFC8F*7$$\"1+++X,4d6Fcy$\"1+++`A$*H8F*7$$\"1+++$*o%y:\"Fcy$\"1+
++@OmN8F*7$$\"1+++//ge6Fcy$\"1+++&zp9M\"F*7$$\"1+++H2Nf6Fcy$\"1+++hz5Z
8F*7$$\"1+++(y(4g6Fcy$\"1+++Arv_8F*7$$\"1+++J;%3;\"Fcy$\"1+++C:ae8F*7$
$\"1+++\\Ceh6Fcy$\"1+++8=@k8F*7$$\"1+++$4?B;\"Fcy$\"1+++.Xwp8F*7$$\"1+
++;Y0j6Fcy$\"1+++#\\^aP\"F*7$$\"1+++Thyj6Fcy$\"1+++nbF\"Q\"F*7$$\"1+++
@Z^k6Fcy$\"1+++Ur%oQ\"F*7$$\"1+++4-Cl6Fcy$\"1+++t_N#R\"F*7$$\"1+++eE'f
;\"Fcy$\"1+++fY1)R\"F*7$$\"1+++HBom6Fcy$\"1+++\"pyPS\"F*7$$\"1+++t!*Rn
6Fcy$\"1+++M3O49F*7$$\"1+++yG6o6Fcy$\"1+++KM([T\"F*7$$\"1+++*zB)o6Fcy$
\"1+++K6_?9F*7$$\"1+++h>`p6Fcy$\"1+++_$phU\"F*7$$\"1+++`tBq6Fcy$\"1+++
'4v;V\"F*7$$\"1+++&*)R4<\"Fcy$\"1+++OiCP9F*-%'COLOURG6&%$RGBG$\"#5!\"
\"\"\"!F[in-%+AXESLABELSG6$%!GF_in-%%VIEWG6$;F[in$\"\"#Fjhn;F[in$Fein!
\"#" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 95 "The curve for the specific heat C(T) turn
s around at low T to appraoch zero with a small slope!" }}{PARA 0 "" 
0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 274 26 "The Boltzma
nn distribution" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" 
{TEXT -1 374 "An important quantity in statistical mechanics is the Bo
ltzmann factor: for a system at a given temperature it provides the pr
obability to find the system in a microstate of energy E. The system i
s at temperature T while in thermal contact with a heat bath (a large \+
system). This relative probability is given by exp(-E/kT), and is deri
ved in texts on statistical mechanics." }}{PARA 0 "" 0 "" {TEXT -1 0 "
" }}{PARA 0 "" 0 "" {TEXT -1 302 "The results from the previous sectio
n can be used to demonstrate this behaviour. We can consider a coupled
 system of a N=1000 oscillator Einstein solid coupled to a single osci
llator. We vary the temperature by changing the amount of energy units
 shared between the single oscillator and the heat bath." }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "restart; with(plots):" }}{PARA 0 ">
 " 0 "" {MPLTEXT 1 0 56 "Omega:=(N,q)->binomial(q+N-1,q); ladd:=(L,e)-
>[op(L),e];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&OmegaGR6$%\"NG%\"qG6
\"6$%)operatorG%&arrowGF)-%)binomialG6$,(9%\"\"\"9$F2!\"\"F2F1F)F)F)" 
}}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%laddGR6$%\"LG%\"eG6\"6$%)operator
G%&arrowGF)7$-%#opG6#9$9%F)F)F)" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 
73 "The 'solid' A is the single oscillator, while B represents the hea
t bath." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "N_A:=1; N_B:=100
0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$N_AG\"\"\"" }}{PARA 11 "" 1 "
" {XPPMATH 20 "6#>%$N_BG\"%+5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 9 "q_t:=100;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$q_tG\"$+\"" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "O_A:=[]: O_B:=[]: O_t:=[]:" 
}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "for q_A from 0 to q_t do:
" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "q_B:=q_t-q_A;" }}{PARA 0 "> " 
0 "" {MPLTEXT 1 0 20 "o_A:=Omega(N_A,q_A);" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 20 "o_B:=Omega(N_B,q_B);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 
0 13 "o_t:=o_A*o_B;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "O_A:=ladd(O_
A,o_A); O_B:=ladd(O_B,o_B); O_t:=ladd(O_t,o_t);" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 3 "od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "O_t:
=map(evalf@log,O_t):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "P1:
=listplot(O_t,view=[0..30,250..340],color=blue): display(P1);" }}
{PARA 13 "" 1 "" {GLPLOT2D 627 297 297 {PLOTDATA 2 "6$-%'CURVESG6$7aq7
$$\"\"\"\"\"!$\"+voH=L!\"(7$$\"\"#F*$\"+<qK%H$F-7$$\"\"$F*$\"+fdEqKF-7
$$\"\"%F*$\"+*37hC$F-7$$\"\"&F*$\"+v\\'=A$F-7$$\"\"'F*$\"+gL_(>$F-7$$
\"\"(F*$\"+ph3tJF-7$$\"\")F*$\"+-Bb[JF-7$$\"\"*F*$\"+N1#R7$F-7$$\"#5F*
$\"+@+>*4$F-7$$\"#6F*$\"+(GfV2$F-7$$\"#7F*$\"+KsU\\IF-7$$\"#8F*$\"+KER
CIF-7$$\"#9F*$\"+JUD**HF-7$$\"#:F*$\"+Z2,uHF-7$$\"#;F*$\"+n3m[HF-7$$\"
#<F*$\"+[K?BHF-7$$\"#=F*$\"+9lj(*GF-7$$\"#>F*$\"+b#f>(GF-7$$\"#?F*$\"+
H+<YGF-7$$\"#@F*$\"+etE?GF-7$$\"#AF*$\"+C(\\Uz#F-7$$\"#BF*$\"+vb6oFF-7
$$\"#CF*$\"+<L'=u#F-7$$\"#DF*$\"+98\\:FF-7$$\"#EF*$\"+))y**)o#F-7$$\"#
FF*$\"+;8QiEF-7$$\"#GF*$\"+I)Rcj#F-7$$\"#HF*$\"+6;x3EF-7$$\"#IF*$\"+\"
zu<e#F-7$$\"#JF*$\"+^ukaDF-7$$\"#KF*$\"+9wQFDF-7$$\"#LF*$\"+\\K***\\#F
-7$$\"#MF*$\"+jAYsCF-7$$\"#NF*$\"+.DzWCF-7$$\"#OF*$\"+]<)pT#F-7$$\"#PF
*$\"+;x-*Q#F-7$$\"#QF*$\"+V!G4O#F-7$$\"#RF*$\"+)H!oKBF-7$$\"#SF*$\"+r>
G/BF-7$$\"#TF*$\"+o/tvAF-7$$\"#UF*$\"+5J-ZAF-7$$\"#VF*$\"+Gr:=AF-7$$\"
#WF*$\"+e'H\"*=#F-7$$\"#XF*$\"+Ox$*f@F-7$$\"#YF*$\"+$Hy08#F-7$$\"#ZF*$
\"+Z\"[55#F-7$$\"#[F*$\"+.SMr?F-7$$\"#\\F*$\"+QCYT?F-7$$\"#]F*$\"+,**R
6?F-7$$\"#^F*$\"+.F:\")>F-7$$\"#_F*$\"+2qr]>F-7$$\"#`F*$\"+A))3?>F-7$$
\"#aF*$\"+!*RE*)=F-7$$\"#bF*$\"+z\"Q#e=F-7$$\"#cF*$\"+ro+F=F-7$$\"#dF*
$\"+\\`c&z\"F-7$$\"#eF*$\"+!o3Rw\"F-7$$\"#fF*$\"+4<.K<F-7$$\"#gF*$\"+N
!H**p\"F-7$$\"#hF*$\"+'*\\fn;F-7$$\"#iF*$\"+]O-N;F-7$$\"#jF*$\"+c(3Ag
\"F-7$$\"#kF*$\"+YP9p:F-7$$\"#lF*$\"++<#e`\"F-7$$\"#mF*$\"+B`B-:F-7$$
\"#nF*$\"+.pPo9F-7$$\"#oF*$\"+$GQUV\"F-7$$\"#pF*$\"+?3\")*R\"F-7$$\"#q
F*$\"+M`3l8F-7$$\"#rF*$\"+l?0I8F-7$$\"#sF*$\"+61q%H\"F-7$$\"#tF*$\"+l)
>!f7F-7$$\"#uF*$\"+Qz*HA\"F-7$$\"#vF*$\"+x?i'=\"F-7$$\"#wF*$\"+j&y)\\6
F-7$$\"#xF*$\"+.Ev76F-7$$\"#yF*$\"+%>G_2\"F-7$$\"#zF*$\"+ozGP5F-7$$\"#
!)F*$\"+[,8*)**!\")7$$\"#\")F*$\"+,m#3g*Fcdl7$$\"##)F*$\"+D@u2#*Fcdl7$
$\"#$)F*$\"+/li4))Fcdl7$$\"#%)F*$\"+SC?1%)Fcdl7$$\"#&)F*$\"+64;(*zFcdl
7$$\"#')F*$\"+%Rb@e(Fcdl7$$\"#()F*$\"+'f%zgrFcdl7$$\"#))F*$\"+zJjKnFcd
l7$$\"#*)F*$\"+H(frH'Fcdl7$$\"#!*F*$\"+s3y`eFcdl7$$\"#\"*F*$\"+R)*z,aF
cdl7$$\"##*F*$\"+YoQS\\Fcdl7$$\"#$*F*$\"+fp`oWFcdl7$$\"#%*F*$\"+C!3])R
Fcdl7$$\"#&*F*$\"+/`A)[$Fcdl7$$\"#'*F*$\"+pp7wHFcdl7$$\"#(*F*$\"+Ig*eW
#Fcdl7$$\"#)*F*$\"+(Q]M*=Fcdl7$$\"#**F*$\"+)GOBJ\"Fcdl7$$\"$+\"F*$\"+z
_v2p!\"*7$$\"$,\"F*F*-%'COLOURG6&%$RGBGF*F*$\"*++++\"Fcdl-%%VIEWG6$;F*
Fet;$\"$]#F*$\"$S$F*" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 
45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 491 "The plot is shift
ed by one unit on the horizontal axis, i.e., the point at (1,335) corr
esponds to q_A=0. The graph shows the logarithm of the probability for
 the microstate (its multiplicity) as a function of the number of unit
s of energy. The function falls off exponentially which implies a line
ar fall-off of the logarithm. Thus, we have demonstrated the Boltzmann
 factor variation with E. [At the high end of the curve near q_A=q_t o
ne can see a deviation due to the finite number q_t]. " }}{PARA 0 "" 
0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 307 "Next we observe the
 change with temperature. We estimate the temparature of the heat bath
 at the centre of the graph shown above (it varies a little due to the
 finiteness of the calculation, i.e., the removal of 1-30 units of ene
rgy out of a bath of 1000 oscillators that share 100 units is not negl
igible)." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "Temp1:=1/(O_t[1
4]-O_t[15]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&Temp1G$\"+@$=9'R!#5
" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 44 "Now we quadruple the number o
f energy units:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "q_t:=400;
" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$q_tG\"$+%" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 26 "O_A:=[]: O_B:=[]: O_t:=[]:" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 25 "for q_A from 0 to q_t do:" }}{PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 13 "q_B:=q_t-q_A;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 
0 20 "o_A:=Omega(N_A,q_A);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "o_B:=
Omega(N_B,q_B);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "o_t:=o_A*o_B;" }
}{PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "O_A:=ladd(O_A,o_A); O_B:=ladd(O_B,
o_B); O_t:=ladd(O_t,o_t);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od:" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "O_t:=map(evalf@log,O_t):" }
}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 65 "P2:=listplot(O_t,view=[0..
30,790..840],color=green): display(P2);" }}{PARA 13 "" 1 "" {GLPLOT2D 
627 297 297 {PLOTDATA 2 "6$-%'CURVESG6$7]dl7$$\"\"\"\"\"!$\"+]G%\\L)!
\"(7$$\"\"#F*$\"+lBUA$)F-7$$\"\"$F*$\"++S))4$)F-7$$\"\"%F*$\"+(pFtH)F-
7$$\"\"&F*$\"+(R`ZG)F-7$$\"\"'F*$\"+U5;s#)F-7$$\"\"(F*$\"+u0bf#)F-7$$
\"\")F*$\"+L>#pC)F-7$$\"\"*F*$\"+h]FM#)F-7$$\"#5F*$\"+(*)4;A)F-7$$\"#6
F*$\"+#QE*3#)F-7$$\"#7F*$\"+cWA'>)F-7$$\"#8F*$\"+eS]$=)F-7$$\"#9F*$\"+
G^wq\")F-7$$\"#:F*$\"+.w+e\")F-7$$\"#;F*$\"+A9BX\")F-7$$\"#<F*$\"+ClVK
\")F-7$$\"#=F*$\"+YGi>\")F-7$$\"#>F*$\"+E.z1\")F-7$$\"#?F*$\"+,*QR4)F-
7$$\"#@F*$\"+2&o53)F-7$$\"#AF*$\"+\"3z\"o!)F-7$$\"#BF*$\"+e0Fb!)F-7$$
\"#CF*$\"+vGMU!)F-7$$\"#DF*$\"+nfRH!)F-7$$\"#EF*$\"+o(Hk,)F-7$$\"#FF*$
\"+9UW.!)F-7$$\"#GF*$\"+Q#R/*zF-7$$\"#HF*$\"+uZTxzF-7$$\"#IF*$\"+b2Pkz
F-7$$\"#JF*$\"+;rI^zF-7$$\"#KF*$\"+)yB#QzF-7$$\"#LF*$\"+.27DzF-7$$\"#M
F*$\"+%z(*>\"zF-7$$\"#NF*$\"+#*\\&))*yF-7$$\"#OF*$\"+GAp&)yF-7$$\"#PF*
$\"+M%4D(yF-7$$\"#QF*$\"+QlIfyF-7$$\"#RF*$\"+rM3YyF-7$$\"#SF*$\"+k,%G$
yF-7$$\"#TF*$\"+Vld>yF-7$$\"#UF*$\"+RDH1yF-7$$\"#VF*$\"+!3))Hz(F-7$$\"
#WF*$\"+%4j'zxF-7$$\"#XF*$\"+2vJmxF-7$$\"#YF*$\"+Z7&Hv(F-7$$\"#ZF*$\"+
SUcRxF-7$$\"#[F*$\"+8k:ExF-7$$\"#\\F*$\"+\"pFFr(F-7$$\"#]F*$\"+**zF*p(
F-7$$\"#^F*$\"+js!eo(F-7$$\"#_F*$\"+1aJswF-7$$\"#`F*$\"+_B!)ewF-7$$\"#
aF*$\"+D!o_k(F-7$$\"#bF*$\"+[BrJwF-7$$\"#cF*$\"+V_8=wF-7$$\"#dF*$\"+Km
`/wF-7$$\"#eF*$\"+Ok\"4f(F-7$$\"#fF*$\"+xXFxvF-7$$\"#gF*$\"+w4hjvF-7$$
\"#hF*$\"+_b#*\\vF-7$$\"#iF*$\"+D#=i`(F-7$$\"#jF*$\"+8*)[AvF-7$$\"#kF*
$\"+Ovt3vF-7$$\"#lF*$\"+7S'\\\\(F-7$$\"#mF*$\"+e#o6[(F-7$$\"#nF*$\"+\"
>]tY(F-7$$\"#oF*$\"+G(4NX(F-7$$\"#pF*$\"+%yY'RuF-7$$\"#qF*$\"+w7wDuF-7
$$\"#rF*$\"+<J&=T(F-7$$\"#sF*$\"+BA#zR(F-7$$\"#tF*$\"+2&oRQ(F-7$$\"#uF
*$\"+$)=**ptF-7$$\"#vF*$\"+jA*fN(F-7$$\"#wF*$\"+g&p>M(F-7$$\"#xF*$\"+%
oBzK(F-7$$\"#yF*$\"+[X&QJ(F-7$$\"#zF*$\"+h?w*H(F-7$$\"#!)F*$\"+Lhk&G(F
-7$$\"#\")F*$\"+um]rsF-7$$\"##)F*$\"+$fVtD(F-7$$\"#$)F*$\"+'zcJC(F-7$$
\"#%)F*$\"+$>Y*GsF-7$$\"#&)F*$\"+!p6Z@(F-7$$\"#')F*$\"+$>`/?(F-7$$\"#(
)F*$\"+21<'=(F-7$$\"#))F*$\"+QQ'=<(F-7$$\"#*)F*$\"+!zKv:(F-7$$\"#!*F*$
\"+nt<VrF-7$$\"#\"*F*$\"+suzGrF-7$$\"##*F*$\"+3IR9rF-7$$\"#$*F*$\"+vQ'
**4(F-7$$\"#%*F*$\"+w*4b3(F-7$$\"#&*F*$\"+57.rqF-7$$\"#'*F*$\"+yu_cqF-
7$$\"#(*F*$\"+y')*>/(F-7$$\"#)*F*$\"+5ZWFqF-7$$\"#**F*$\"+pa'G,(F-7$$
\"$+\"F*$\"+b3E)*pF-7$$\"$,\"F*$\"+i2j$)pF-7$$\"$-\"F*$\"+(3v*opF-7$$
\"$.\"F*$\"+BPHapF-7$$\"$/\"F*$\"+mleRpF-7$$\"$0\"F*$\"+4N&[#pF-7$$\"$
1\"F*$\"+WW45pF-7$$\"$2\"F*$\"+j#4`*oF-7$$\"$3\"F*$\"+cy\\!)oF-7$$\"$4
\"F*$\"+;,mloF-7$$\"$5\"F*$\"+Ifz]oF-7$$\"$6\"F*$\"+(=0f$oF-7$$\"$7\"F
*$\"+xx)4#oF-7$$\"$8\"F*$\"+%eVg!oF-7$$\"$9\"F*$\"+(\\s5z'F-7$$\"$:\"F
*$\"++W2wnF-7$$\"$;\"F*$\"+y\"\\5w'F-7$$\"$<\"F*$\"+:n*fu'F-7$$\"$=\"F
*$\"+%*o\"4t'F-7$$\"$>\"F*$\"+(f4er'F-7$$\"$?\"F*$\"+/Zn+nF-7$$\"$@\"F
*$\"+)47bo'F-7$$\"$A\"F*$\"+c;KqmF-7$$\"$B\"F*$\"+dK5bmF-7$$\"$C\"F*$
\"+!yc)RmF-7$$\"$D\"F*$\"++@eCmF-7$$\"$E\"F*$\"+%4z#4mF-7$$\"$F\"F*$\"
+Pw%Rf'F-7$$\"$G\"F*$\"+-weylF-7$$\"$H\"F*$\"+i))>jlF-7$$\"$I\"F*$\"+!
H\"yZlF-7$$\"$J\"F*$\"+dZLKlF-7$$\"$K\"F*$\"+K\"fo^'F-7$$\"$L\"F*$\"+%
Ga8]'F-7$$\"$M\"F*$\"+#3?e['F-7$$\"$N\"F*$\"+$Rc-Z'F-7$$\"$O\"F*$\"+#3
jYX'F-7$$\"$P\"F*$\"+9+/RkF-7$$\"$Q\"F*$\"+aqQBkF-7$$\"$R\"F*$\"+kSq2k
F-7$$\"$S\"F*$\"+14*>R'F-7$$\"$T\"F*$\"+SuCwjF-7$$\"$U\"F*$\"+ENZgjF-7
$$\"$V\"F*$\"+A!pYM'F-7$$\"$W\"F*$\"+'yL)GjF-7$$\"$X\"F*$\"+tw'HJ'F-7$
$\"$Y\"F*$\"+R02(H'F-7$$\"$Z\"F*$\"+PA9\"G'F-7$$\"$[\"F*$\"+?E=liF-7$$
\"$\\\"F*$\"+S:>\\iF-7$$\"$]\"F*$\"+X)oJB'F-7$$\"$^\"F*$\"+(Q9r@'F-7$$
\"$_\"F*$\"+6!G5?'F-7$$\"$`\"F*$\"+k&4\\='F-7$$\"$a\"F*$\"+#*)e(ohF-7$
$\"$b\"F*$\"+Qed_hF-7$$\"$c\"F*$\"+Y-OOhF-7$$\"$d\"F*$\"+b>6?hF-7$$\"$
e\"F*$\"+13$Q5'F-7$$\"$f\"F*$\"+Qm^(3'F-7$$\"$g\"F*$\"+)Gp62'F-7$$\"$h
\"F*$\"+!f)yagF-7$$\"$i\"F*$\"+\"Qu$QgF-7$$\"$j\"F*$\"+#\\E>-'F-7$$\"$
k\"F*$\"+bZW0gF-7$$\"$l\"F*$\"++!H*))fF-7$$\"$m\"F*$\"+b!zB(fF-7$$\"$n
\"F*$\"+\\ZzbfF-7$$\"$o\"F*$\"+0f<RfF-7$$\"$p\"F*$\"+[B_AfF-7$$\"$q\"F
*$\"++R$e!fF-7$$\"$r\"F*$\"+$Q5\"*)eF-7$$\"$s\"F*$\"+:;NseF-7$$\"$t\"F
*$\"+9ubbeF-7$$\"$u\"F*$\"+'fF(QeF-7$$\"$v\"F*$\"+v>'=#eF-7$$\"$w\"F*$
\"+j.'\\!eF-7$$\"$x\"F*$\"+tD-)y&F-7$$\"$y\"F*$\"+7%[5x&F-7$$\"$z\"F*$
\"+)oPSv&F-7$$\"$!=F*$\"+2-*pt&F-7$$\"$\"=F*$\"+td!*>dF-7$$\"$#=F*$\"+
(=%y-dF-7$$\"$$=F*$\"+]_i&o&F-7$$\"$%=F*$\"+g(G%ocF-7$$\"$&=F*$\"+8X>^
cF-7$$\"$'=F*$\"+.B#Rj&F-7$$\"$(=F*$\"+C>h;cF-7$$\"$)=F*$\"+lJE*f&F-7$
$\"$*=F*$\"+:e(=e&F-7$$\"$!>F*$\"+g'\\Wc&F-7$$\"$\">F*$\"+&[%)pa&F-7$$
\"$#>F*$\"+s+[HbF-7$$\"$$>F*$\"++i$>^&F-7$$\"$%>F*$\"+\\EN%\\&F-7$$\"$
&>F*$\"+$>HnZ&F-7$$\"$'>F*$\"+1c1faF-7$$\"$(>F*$\"+g;OTaF-7$$\"$)>F*$
\"+BrhBaF-7$$\"$*>F*$\"+h<$eS&F-7$$\"$+#F*$\"+S`+)Q&F-7$$\"$,#F*$\"+@w
8q`F-7$$\"$-#F*$\"+j$GAN&F-7$$\"$.#F*$\"+BtFM`F-7$$\"$/#F*$\"+dUG;`F-7
$$\"$0#F*$\"+:*[#)H&F-7$$\"$1#F*$\"+Z5<!G&F-7$$\"$2#F*$\"+*R]?E&F-7$$
\"$3#F*$\"+<n)QC&F-7$$\"$4#F*$\"+S(zcA&F-7$$\"$5#F*$\"+3#Hu?&F-7$$\"$6
#F*$\"+c[8*=&F-7$$\"$7#F*$\"+=kzq^F-7$$\"$8#F*$\"+BOT_^F-7$$\"$9#F*$\"
+*>')R8&F-7$$\"$:#F*$\"+pQ^:^F-7$$\"$;#F*$\"+cj*p4&F-7$$\"$<#F*$\"+wLV
y]F-7$$\"$=#F*$\"+XY#)f]F-7$$\"$>#F*$\"+u)p6/&F-7$$\"$?#F*$\"+t(oC-&F-
7$$\"$@#F*$\"+Y5s.]F-7$$\"$A#F*$\"+'RE\\)\\F-7$$\"$B#F*$\"+?X3m\\F-7$$
\"$C#F*$\"+9^>Z\\F-7$$\"$D#F*$\"+qyDG\\F-7$$\"$E#F*$\"+wCF4\\F-7$$\"$F
#F*$\"+;'Q-*[F-7$$\"$G#F*$\"+qf:r[F-7$$\"$H#F*$\"+<U-_[F-7$$\"$I#F*$\"
+GI%G$[F-7$$\"$J#F*$\"+t?h8[F-7$$\"$K#F*$\"+=5L%z%F-7$$\"$L#F*$\"+B&**
\\x%F-7$$\"$M#F*$\"+ZshbZF-7$$\"$N#F*$\"+TQ=OZF-7$$\"$O#F*$\"+b*)p;ZF-
7$$\"$P#F*$\"+MA;(p%F-7$$\"$Q#F*$\"+<LdxYF-7$$\"$R#F*$\"+R=$zl%F-7$$\"
$S#F*$\"+LuBQYF-7$$\"$T#F*$\"+C(*[=YF-7$$\"$U#F*$\"+L$)o)f%F-7$$\"$V#F
*$\"+zG$)yXF-7$$\"$W#F*$\"+sH#*eXF-7$$\"$X#F*$\"+>#e*QXF-7$$\"$Y#F*$\"
+B#Q*=XF-7$$\"$Z#F*$\"+!ei))\\%F-7$$\"$[#F*$\"+\")3tyWF-7$$\"$\\#F*$\"
+7FaeWF-7$$\"$]#F*$\"+`wHQWF-7$$\"$^#F*$\"+z_*zT%F-7$$\"$_#F*$\"+f^j(R
%F-7$$\"$`#F*$\"+co@xVF-7$$\"$a#F*$\"+G*RnN%F-7$$\"$b#F*$\"+DR?OVF-7$$
\"$c#F*$\"+#R3cJ%F-7$$\"$d#F*$\"+oG&\\H%F-7$$\"$e#F*$\"+%)oBuUF-7$$\"$
f#F*$\"+m*fMD%F-7$$\"$g#F*$\"+L;iKUF-7$$\"$h#F*$\"+(R@<@%F-7$$\"$i#F*$
\"+h(e2>%F-7$$\"$j#F*$\"+DKtpTF-7$$\"$k#F*$\"+xUk[TF-7$$\"$l#F*$\"++9
\\FTF-7$$\"$m#F*$\"+pSF1TF-7$$\"$n#F*$\"+_<*\\3%F-7$$\"$o#F*$\"+3RkjSF
-7$$\"$p#F*$\"+')*HA/%F-7$$\"$q#F*$\"+J%\\2-%F-7$$\"$r#F*$\"+v;?**RF-7
$$\"$s#F*$\"+WhexRF-7$$\"$t#F*$\"+aA!f&RF-7$$\"$u#F*$\"+6%\\T$RF-7$$\"
$v#F*$\"+9qK7RF-7$$\"$w#F*$\"+]WV!*QF-7$$\"$x#F*$\"+(4r%oQF-7$$\"$y#F*
$\"+BjVYQF-7$$\"$z#F*$\"+'[HV#QF-7$$\"$!GF*$\"+K*\\@!QF-7$$\"$\"GF*$\"
+)*p*)zPF-7$$\"$#GF*$\"+3+ddPF-7$$\"$$GF*$\"+w#o^t$F-7$$\"$%GF*$\"+06p
7PF-7$$\"$&GF*$\"+$yP,p$F-7$$\"$'GF*$\"+)e2vm$F-7$$\"$(GF*$\"+&y*zWOF-
7$$\"$)GF*$\"+EO,AOF-7$$\"$*GF*$\"+\\$[\"*f$F-7$$\"$!HF*$\"+!=.id$F-7$
$\"$\"HF*$\"+Ht<`NF-7$$\"$#HF*$\"+$**p+`$F-7$$\"$$HF*$\"+`.)o]$F-7$$\"
$%HF*$\"+wvg$[$F-7$$\"$&HF*$\"+73DgMF-7$$\"$'HF*$\"+(>4oV$F-7$$\"$(HF*
$\"+[=G8MF-7$$\"$)HF*$\"+nym*Q$F-7$$\"$*HF*$\"+Pj'fO$F-7$$\"$+$F*$\"+C
j<ULF-7$$\"$,$F*$\"+voH=LF-7$$\"$-$F*$\"+<qK%H$F-7$$\"$.$F*$\"+fdEqKF-
7$$\"$/$F*$\"+*37hC$F-7$$\"$0$F*$\"+v\\'=A$F-7$$\"$1$F*$\"+gL_(>$F-7$$
\"$2$F*$\"+ph3tJF-7$$\"$3$F*$\"+-Bb[JF-7$$\"$4$F*$\"+N1#R7$F-7$$\"$5$F
*$\"+@+>*4$F-7$$\"$6$F*$\"+(GfV2$F-7$$\"$7$F*$\"+KsU\\IF-7$$\"$8$F*$\"
+KERCIF-7$$\"$9$F*$\"+JUD**HF-7$$\"$:$F*$\"+Z2,uHF-7$$\"$;$F*$\"+n3m[H
F-7$$\"$<$F*$\"+[K?BHF-7$$\"$=$F*$\"+9lj(*GF-7$$\"$>$F*$\"+b#f>(GF-7$$
\"$?$F*$\"+H+<YGF-7$$\"$@$F*$\"+etE?GF-7$$\"$A$F*$\"+C(\\Uz#F-7$$\"$B$
F*$\"+vb6oFF-7$$\"$C$F*$\"+<L'=u#F-7$$\"$D$F*$\"+98\\:FF-7$$\"$E$F*$\"
+))y**)o#F-7$$\"$F$F*$\"+;8QiEF-7$$\"$G$F*$\"+I)Rcj#F-7$$\"$H$F*$\"+6;
x3EF-7$$\"$I$F*$\"+\"zu<e#F-7$$\"$J$F*$\"+^ukaDF-7$$\"$K$F*$\"+9wQFDF-
7$$\"$L$F*$\"+\\K***\\#F-7$$\"$M$F*$\"+jAYsCF-7$$\"$N$F*$\"+.DzWCF-7$$
\"$O$F*$\"+]<)pT#F-7$$\"$P$F*$\"+;x-*Q#F-7$$\"$Q$F*$\"+V!G4O#F-7$$\"$R
$F*$\"+)H!oKBF-7$$\"$S$F*$\"+r>G/BF-7$$\"$T$F*$\"+o/tvAF-7$$\"$U$F*$\"
+5J-ZAF-7$$\"$V$F*$\"+Gr:=AF-7$$\"$W$F*$\"+e'H\"*=#F-7$$\"$X$F*$\"+Ox$
*f@F-7$$\"$Y$F*$\"+$Hy08#F-7$$\"$Z$F*$\"+Z\"[55#F-7$$\"$[$F*$\"+.SMr?F
-7$$\"$\\$F*$\"+QCYT?F-7$$\"$]$F*$\"+,**R6?F-7$$\"$^$F*$\"+.F:\")>F-7$
$\"$_$F*$\"+2qr]>F-7$$\"$`$F*$\"+A))3?>F-7$$\"$a$F*$\"+!*RE*)=F-7$$\"$
b$F*$\"+z\"Q#e=F-7$$\"$c$F*$\"+ro+F=F-7$$\"$d$F*$\"+\\`c&z\"F-7$$\"$e$
F*$\"+!o3Rw\"F-7$$\"$f$F*$\"+4<.K<F-7$$\"$g$F*$\"+N!H**p\"F-7$$\"$h$F*
$\"+'*\\fn;F-7$$\"$i$F*$\"+]O-N;F-7$$\"$j$F*$\"+c(3Ag\"F-7$$\"$k$F*$\"
+YP9p:F-7$$\"$l$F*$\"++<#e`\"F-7$$\"$m$F*$\"+B`B-:F-7$$\"$n$F*$\"+.pPo
9F-7$$\"$o$F*$\"+$GQUV\"F-7$$\"$p$F*$\"+?3\")*R\"F-7$$\"$q$F*$\"+M`3l8
F-7$$\"$r$F*$\"+l?0I8F-7$$\"$s$F*$\"+61q%H\"F-7$$\"$t$F*$\"+l)>!f7F-7$
$\"$u$F*$\"+Qz*HA\"F-7$$\"$v$F*$\"+x?i'=\"F-7$$\"$w$F*$\"+j&y)\\6F-7$$
\"$x$F*$\"+.Ev76F-7$$\"$y$F*$\"+%>G_2\"F-7$$\"$z$F*$\"+ozGP5F-7$$\"$!Q
F*$\"+[,8*)**!\")7$$\"$\"QF*$\"+,m#3g*F_br7$$\"$#QF*$\"+D@u2#*F_br7$$
\"$$QF*$\"+/li4))F_br7$$\"$%QF*$\"+SC?1%)F_br7$$\"$&QF*$\"+64;(*zF_br7
$$\"$'QF*$\"+%Rb@e(F_br7$$\"$(QF*$\"+'f%zgrF_br7$$\"$)QF*$\"+zJjKnF_br
7$$\"$*QF*$\"+H(frH'F_br7$$\"$!RF*$\"+s3y`eF_br7$$\"$\"RF*$\"+R)*z,aF_
br7$$\"$#RF*$\"+YoQS\\F_br7$$\"$$RF*$\"+fp`oWF_br7$$\"$%RF*$\"+C!3])RF
_br7$$\"$&RF*$\"+/`A)[$F_br7$$\"$'RF*$\"+pp7wHF_br7$$\"$(RF*$\"+Ig*eW#
F_br7$$\"$)RF*$\"+(Q]M*=F_br7$$\"$*RF*$\"+)GOBJ\"F_br7$$\"$+%F*$\"+z_v
2p!\"*7$$\"$,%F*F*-%'COLOURG6&%$RGBGF*$\"*++++\"F_brF*-%%VIEWG6$;F*Fet
;$\"$!zF*$\"$S)F*" 1 2 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 
}}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 31 "Let us compare the two slopes:
 " }}{PARA 0 "" 0 "" {TEXT -1 57 "In case 1 (q_t=100) :  70 units vert
ical vs 30 horizontal" }}{PARA 0 "" 0 "" {TEXT -1 57 "In case 2 (q_t=4
00) :  35 units vertical vs 30 horizontal" }}{PARA 0 "" 0 "" {TEXT -1 
0 "" }}{PARA 0 "" 0 "" {TEXT -1 279 "We quadrupled the number of energ
y units in the heat bath made up of 1000 oscillators and the slope in \+
the log(Omega) plot halved. If Omega is proportional to exp(-E/kT) as \+
claimed we need to show that the temperature increases like the square
 root of the number of energy units." }}{PARA 0 "" 0 "" {TEXT -1 0 "" 
}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 27 "Temp2:=1/(O_t[14]-O_t[15]);" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%&Temp2G$\"+'*4^Qy!#5" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT -1 275 "The temperature of the heat bath is approximately twice \+
of what it was in the first case. Thus, we have demonstrated that the \+
probability distribution as calculated from the multiplicity associate
d with the heat bath follows the dependence given by the Boltzmann fac
tor exp(-" }{TEXT 275 1 "E" }{TEXT -1 2 "/k" }{TEXT 276 1 "T" }{TEXT 
-1 2 ")." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT 273 43 "Thermal behaviour of a two-state paramagnet" }}{PARA 0 "
" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 539 "Moore and Schroed
er discuss this problem as a an example that leads to further implicat
ions of the theoretical definition of temperature given above. In conv
entional thermodynamics temperature is the measure of the amount of ki
netic motion (linear motion of atoms in a gas, vibrational motion of d
iatomic or other molecules, vibrations of stationary atoms in the Eins
tein solid, etc.). In the case of a magnet the degrees of freedom are \+
contained in the spin orientation of individual atoms with respect to \+
an externally defined direction." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}
{PARA 0 "" 0 "" {TEXT -1 484 "The energy of an atom with an outer elec
tron with unpaired spin in an outer magnetic field depends on the orie
ntation of the electron spin, which gives rise to a magnetic moment (i
f the electron is in an s state,or in a state with zero magnetic quant
um number there is no contribution from orbital angular momentum to th
e magnetic moment). The magnetic moment vector is opposite to the spin
 vector due to the negative charge of the electron (see e.g., Wolfson \+
and Pasachoff, p. 1098)." }}{PARA 0 "" 0 "" {TEXT -1 962 "The energy o
f a magnetic dipole in an external field is given as U = - mu . B , an
d therefore the state with a magnetic dipole oriented in the same dire
ction as the external field has a lower energy than the state that has
 a counteraligned magnetic dipole moment. The magnetic dipole moment i
s defined such that the torque produced by the magnetic field on it ac
ts to align the dipole with the magnetic field vector. (The evidence f
or this is obtained from tracing the magnetic field of a big solenoid \+
or bar magnet by placing a compass at various locations: near the Nort
h and South poles it is evident how the compass needle (its dipole) al
igns with the field; at the sides of the 'big' magnet it apparently co
unteraligns: this is a result of the inhomogeneity of the field: the f
ield lines outside the bar magnet run in the opposite direction to for
m closed loops - this causes two permanent magnets side-by-side to sna
p together in a counteraligned fashion)." }}{PARA 0 "" 0 "" {TEXT -1 
0 "" }}{PARA 0 "" 0 "" {TEXT -1 309 "For our problem this complication
 does not exist. The external field is homogeneous. The total potentia
l energy of a system of N dipoles depends on how many are aligned (poi
nting up) N_u and how many are counteraligned (pointing down) N_d, wit
h N_u + N_d = N . If we chose units such that mu=1 and B=1 we have" }}
{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 66 "The energ
ies of the individual up and down states for the dipoles:" }}{PARA 0 "
" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 24 "E_u = -mu B, E_d =
 mu B." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 59 
"U = N_u E_u + N_d E_d = mu B (N_d - N_u) = mu B (N - 2 N_u)" }}{PARA 
0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 45 "The magnetizati
on of the ensemble is given as" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}
{PARA 0 "" 0 "" {TEXT -1 30 "M = mu (N_u - N_d) = - U / B ." }}{PARA 
0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 147 "These relation
s imply that the macrostate is specified by the total number N, the to
tal energy U, or equivalently by specifying N and N_u (or N_d)." }}
{PARA 0 "" 0 "" {TEXT -1 107 "The multiplicity of such a macrostate is
 given by the number of possibilities to chose N_u out of N states." }
}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "restart; with(plots): ladd:=
(L,e)->[op(L),e];" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "Omega:=(N,N_u)
->binomial(N,N_u);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%laddGR6$%\"LG
%\"eG6\"6$%)operatorG%&arrowGF)7$-%#opG6#9$9%F)F)F)" }}{PARA 11 "" 1 "
" {XPPMATH 20 "6#>%&OmegaG%)binomialG" }}}{EXCHG {PARA 0 "" 0 "" 
{TEXT -1 19 "The number of spins" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 7 "N:=100;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"NG\"$+
\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 264 "The energy range is from -
100 to 100 (or -N to N) given that the spins can all be aligned up to \+
counteraligned. First we calculate the energies, multiplicities and en
tropy: given that a single spin can only be aligne or counteraligned t
he loop runs over the spins:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 24 "U1:=-[]: O1:=[]: S1:=[]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 29 "for N_u from N to 0 by -1 do:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 
0 14 "u1:=(N-2*N_u);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "o1:=Omega(N
,N_u);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "s1:=log(o1);" }}{PARA 0 "
> " 0 "" {MPLTEXT 1 0 50 "U1:=ladd(U1,u1); O1:=ladd(O1,o1); S1:=ladd(S
1,s1);" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od:" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 3 "U1;" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#7aq!$+\"
!#)*!#'*!#%*!##*!#!*!#))!#')!#%)!##)!#!)!#y!#w!#u!#s!#q!#o!#m!#k!#i!#g
!#e!#c!#a!#_!#]!#[!#Y!#W!#U!#S!#Q!#O!#M!#K!#I!#G!#E!#C!#A!#?!#=!#;!#9!
#7!#5!\")!\"'!\"%!\"#\"\"!\"\"#\"\"%\"\"'\"\")\"#5\"#7\"#9\"#;\"#=\"#?
\"#A\"#C\"#E\"#G\"#I\"#K\"#M\"#O\"#Q\"#S\"#U\"#W\"#Y\"#[\"#]\"#_\"#a\"
#c\"#e\"#g\"#i\"#k\"#m\"#o\"#q\"#s\"#u\"#w\"#y\"#!)\"##)\"#%)\"#')\"#)
)\"#!*\"##*\"#%*\"#'*\"#)*\"$+\"" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 32 "plot([[U1[j],S1[j]] $j=1..N+1]);" }}{PARA 13 "" 1 "" 
{GLPLOT2D 743 227 227 {PLOTDATA 2 "6%-%'CURVESG6$7aq7$$!$+\"\"\"!F*7$$
!#)*F*$\"1#4))f=q^g%!#:7$$!#'*F*$\"1NFcbG92&)F07$$!#%*F*$\"1?lb/)\\$*>
\"!#97$$!##*F*$\"1p[Hm9>=:F;7$$!#!*F*$\"1V#)>%\\#o8=F;7$$!#))F*$\"1\"
\\NkB%*)*3#F;7$$!#')F*$\"1gpv*pK'\\BF;7$$!#%)F*$\"1-V!\\\\[\\f#F;7$$!#
#)F*$\"1%ev[*[SFGF;7$$!#!)F*$\"1lyAOBB[IF;7$$!#yF*$\"1a5)fxB%eKF;7$$!#
wF*$\"1oa(zu'zeMF;7$$!#uF*$\"1Nrn$\\N+l$F;7$$!#sF*$\"1o5es0sKQF;7$$!#q
F*$\"1)>'4#G]t+%F;7$$!#oF*$\"1^7_NlNuTF;7$$!#mF*$\"1h***4)o6MVF;7$$!#k
F*$\"1/+*fwjp[%F;7$$!#iF*$\"1&y*z#z\">LYF;7$$!#gF*$\"1I;\"4[jIx%F;7$$!
#eF*$\"1wmg+R\"o!\\F;7$$!#cF*$\"1Yv^SWlM]F;7$$!#aF*$\"1\"f$f,fxc^F;7$$
!#_F*$\"1kTug5Nt_F;7$$!#]F*$\"1xfG7o`%Q&F;7$$!#[F*$\"1gu$)pfZ!\\&F;7$$
!#YF*$\"1Xub#z)H\"f&F;7$$!#WF*$\"1jZl&GCro&F;7$$!#UF*$\"1@xb981ydF;7$$
!#SF*$\"1Pc4k&4U'eF;7$$!#QF*$\"1f?&yOgc%fF;7$$!#OF*$\"17=.Gu\\AgF;7$$!
#MF*$\"1vFSUuz%4'F;7$$!#KF*$\"1b-)=lIE;'F;7$$!#IF*$\"1cR$*>81EiF;7$$!#
GF*$\"14z2`\"[^G'F;7$$!#EF*$\"1a%\\,nW*RjF;7$$!#CF*$\"1of\"o_*\\!R'F;7
$$!#AF*$\"18vq+o&oV'F;7$$!#?F*$\"1^MrTi0zkF;7$$!#=F*$\"1I_E\"\\Lr^'F;7
$$!#;F*$\"1lu#QF?6b'F;7$$F;F*$\"1]FJjX/\"e'F;7$$!#7F*$\"1!Q/nsIpg'F;7$
$!#5F*$\"1j3!ok*zGmF;7$$!\")F*$\"1*>vc#)omk'F;7$$!\"'F*$\"1@1;qCbgmF;7
$$!\"%F*$\"1X]Ug:YqmF;7$$!\"#F*$\"1C@Z-RSwmF;7$F*$\"1V<?lTQymF;7$$\"\"
#F*Fiz7$$\"\"%F*Fdz7$$\"\"'F*F_z7$$\"\")F*Fjy7$$\"#5F*Fey7$$\"#7F*F`y7
$$\"#9F*F[y7$$\"#;F*Fgx7$$\"#=F*Fbx7$$\"#?F*F]x7$$\"#AF*Fhw7$$\"#CF*Fc
w7$$\"#EF*F^w7$$\"#GF*Fiv7$$\"#IF*Fdv7$$\"#KF*F_v7$$\"#MF*Fju7$$\"#OF*
Feu7$$\"#QF*F`u7$$\"#SF*F[u7$$\"#UF*Fft7$$\"#WF*Fat7$$\"#YF*F\\t7$$\"#
[F*Fgs7$$\"#]F*Fbs7$$\"#_F*F]s7$$\"#aF*Fhr7$$\"#cF*Fcr7$$\"#eF*F^r7$$
\"#gF*Fiq7$$\"#iF*Fdq7$$\"#kF*F_q7$$\"#mF*Fjp7$$\"#oF*Fep7$$\"#qF*F`p7
$$\"#sF*F[p7$$\"#uF*Ffo7$$\"#wF*Fao7$$\"#yF*F\\o7$$\"#!)F*Fgn7$$\"##)F
*FX7$$\"#%)F*FS7$$\"#')F*FN7$$\"#))F*FI7$$\"#!*F*FD7$$\"##*F*F?7$$\"#%
*F*F97$$\"#'*F*F47$$\"#)*F*F.7$$\"$+\"F*F*-%'COLOURG6&%$RGBG$F\\\\l!\"
\"F*F*-%+AXESLABELSG6$%!GF]el-%%VIEWG6$%(DEFAULTGFael" 1 2 0 1 0 2 9 
1 4 2 1.000000 45.000000 45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT 
-1 211 "The magnetization is trivial: it is simply related to the ener
gy, but we are interested in the temperature and the specific heat cap
acity. The temperature is given as the inverse of the slope of the abo
ve graph." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "S1[3];" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#-%#lnG6#\"%]\\" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 7 "T1:=[]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 45 "for q from 2 to N do: dif:=(S1[q+1]-S1[q-1]);" }}{PARA 0 "> " 
0 "" {MPLTEXT 1 0 51 "if dif<>0 then t1:=evalf((U1[q+1]-U1[q-1])/dif) \+
fi;" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "T1:=ladd(T1,t1): od:" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "#listplot(T1);" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "T1[10];" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"+7fL!G*!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 112 "j:='j': plot([[T1[j],-U1[j]/N] $j=1..N-2],style=point,view=[-10
..10,-1..1],title=`Magnetization per site M(T)`);" }}{PARA 13 "" 1 "" 
{GLPLOT2D 713 296 296 {PLOTDATA 2 "6'-%'CURVESG6$7^q7$$\"1+++_6$>q%!#;
$\"\"\"\"\"!7$$\"1+++/\\%RT&F*$\"1+++++++)*F*7$$\"1+++OUr#*fF*$\"1++++
+++'*F*7$$\"1+++'**H6^'F*$\"1*************R*F*7$$\"1+++!oUm*pF*$\"1+++
++++#*F*7$$\"1+++3*zLY(F*$\"1+++++++!*F*7$$\"1+++K6%*>zF*$\"1+++++++))
F*7$$\"1+++_3>s$)F*$\"1+++++++')F*7$$\"1******>G\\C))F*$\"1+++++++%)F*
7$$\"1+++7fL!G*F*$\"1+++++++#)F*7$$\"1+++OcoU(*F*$\"1+++++++!)F*7$$\"1
+++&*)>9-\"!#:$\"1+++++++yF*7$$\"1+++(fO(p5F]o$\"1+++++++wF*7$$\"1+++%
3h%>6F]o$\"1+++++++uF*7$$\"1+++Cp$3<\"F]o$\"1+++++++sF*7$$\"1+++?`6C7F
]o$\"1+++++++qF*7$$\"1+++O4cz7F]o$\"1+++++++oF*7$$\"1+++`vXP8F]o$\"1++
+++++mF*7$$\"1+++TR6)R\"F]o$\"1+++++++kF*7$$\"1+++z,(=Y\"F]o$\"1++++++
+iF*7$$\"1+++G[5H:F]o$\"1+++++++gF*7$$\"1+++IKC+;F]o$\"1+++++++eF*7$$
\"1+++Stwv;F]o$\"1,++++++cF*7$$\"1+++;wAc<F]o$\"1+++++++aF*7$$\"1+++Pv
DU=F]o$\"1+++++++_F*7$$\"1+++S;fM>F]o$\"1+++++++]F*7$$\"1+++.\")3M?F]o
$\"1+++++++[F*7$$\"1+++GtvT@F]o$\"1+++++++YF*7$$\"1+++C()zeAF]o$\"1+++
++++WF*7$$\"1+++I$[mQ#F]o$\"1+++++++UF*7$$\"1+++#\\Tq_#F]o$\"1+++++++S
F*7$$\"1+++Qj4#o#F]o$\"1+++++++QF*7$$\"1+++&eEW&GF]o$\"1+++++++OF*7$$
\"1+++NoHZIF]o$\"1+++++++MF*7$$\"1+++!pR[E$F]o$\"1+++++++KF*7$$\"1+++V
fO7NF]o$\"1+++++++IF*7$$\"1+++l!=oz$F]o$\"1+++++++GF*7$$\"1+++c-XFTF]o
$\"1+++++++EF*7$$\"1+++!Qzo^%F]o$\"1+++++++CF*7$$\"1+++![mF)\\F]o$\"1+
++++++AF*7$$\"1+++/#>1b&F]o$\"1+++++++?F*7$$\"1+++gkpeiF]o$\"1+++++++=
F*7$$\"1+++oh6nrF]o$\"1+++++++;F*7$$\"1+++S:2w$)F]o$\"1+++++++9F*7$$\"
1+++y1f15!#9$\"1+++++++7F*7$$\"1+++80tf7Fcy$\"1+++++++5F*7$$\"1+++Db=
\"o\"Fcy$Fin!#<7$$\"1+++A8VBDFcy$F\\rF_z7$$\"1+++;W%)[]Fcy$F^uF_z7$Fez
$F`xF_z7$$!1+++;W%)[]FcyF-7$$!1+++A8VBDFcy$!1+++++++?F_z7$$!1+++Db=\"o
\"Fcy$!1+++++++SF_z7$$!1+++80tf7Fcy$!1+++++++gF_z7$$!1+++y1f15Fcy$!1++
+++++!)F_z7$$!1+++S:2w$)F]o$!1+++++++5F*7$$!1+++oh6nrF]o$!1+++++++7F*7
$$!1+++gkpeiF]o$!1+++++++9F*7$$!1+++/#>1b&F]o$!1+++++++;F*7$$!1+++![mF
)\\F]o$!1+++++++=F*7$$!1+++!Qzo^%F]o$Fa[lF*7$$!1+++c-XFTF]o$!1+++++++A
F*7$$!1+++l!=oz$F]o$!1+++++++CF*7$$!1+++VfO7NF]o$!1+++++++EF*7$$!1+++!
pR[E$F]o$!1+++++++GF*7$$!1+++NoHZIF]o$!1+++++++IF*7$$!1+++&eEW&GF]o$!1
+++++++KF*7$$!1+++Qj4#o#F]o$!1+++++++MF*7$$!1+++#\\Tq_#F]o$!1+++++++OF
*7$$!1+++I$[mQ#F]o$!1+++++++QF*7$$!1+++C()zeAF]o$Ff[lF*7$$!1+++GtvT@F]
o$!1+++++++UF*7$$!1+++.\")3M?F]o$!1+++++++WF*7$$!1+++S;fM>F]o$!1++++++
+YF*7$$!1+++PvDU=F]o$!1+++++++[F*7$$!1+++;wAc<F]o$!1+++++++]F*7$$!1+++
Stwv;F]o$!1+++++++_F*7$$!1+++IKC+;F]o$!1+++++++aF*7$$!1+++G[5H:F]o$!1,
++++++cF*7$$!1+++z,(=Y\"F]o$!1+++++++eF*7$$!1+++TR6)R\"F]o$F[\\lF*7$$!
1+++`vXP8F]o$!1+++++++iF*7$$!1+++O4cz7F]o$!1+++++++kF*7$$!1+++?`6C7F]o
$!1+++++++mF*7$$!1+++Cp$3<\"F]o$!1+++++++oF*7$$!1+++%3h%>6F]o$!1++++++
+qF*7$$!1+++(fO(p5F]o$!1+++++++sF*7$$!1+++&*)>9-\"F]o$!1+++++++uF*7$$!
1+++OcoU(*F*$!1+++++++wF*7$$!1+++7fL!G*F*$!1+++++++yF*7$$!1******>G\\C
))F*$F`\\lF*7$$!1+++_3>s$)F*$!1+++++++#)F*7$$!1+++K6%*>zF*$!1+++++++%)
F*7$$!1+++3*zLY(F*$!1+++++++')F*7$$!1+++!oUm*pF*$!1+++++++))F*7$$!1+++
'**H6^'F*$!1+++++++!*F*7$$!1+++OUr#*fF*$!1+++++++#*F*7$$!1+++/\\%RT&F*
$!1*************R*F*-%'COLOURG6&%$RGBG$\"#5!\"\"F-F--%&TITLEG6#%<Magne
tization~per~site~M(T)G-%+AXESLABELSG6$%!GFbjl-%&STYLEG6#%&POINTG-%%VI
EWG6$;$!#5F-$FiilF-;$FjilF-F+" 1 5 0 1 0 2 9 1 4 2 1.000000 45.000000 
45.000000 0 }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 124 "When all spins ar
e aligned the temperature T is at its lowest. At high temperatures the
 average magnetization drops to zero." }}{PARA 0 "" 0 "" {TEXT -1 131 
"The definition of temperature allows for two signs, depending on the \+
sign of the change of rate of the energy U with the entropy S." }}
{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 26 "What is t
he heat capacity?" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "C1:=[]:" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 23 "for q from 2 to N-2 do:" }}{PARA 0 "> " 0 "" 
{MPLTEXT 1 0 42 "c1:=(U1[q+1]-U1[q-1])/(T1[q+1]-T1[q-1])/N;" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 20 "C1:=ladd(C1,c1); od:" }}}{EXCHG {PARA 0 "
> " 0 "" {MPLTEXT 1 0 110 "j:='j': plot([[T1[j+1],C1[j]] $j=1..N-3],st
yle=point,view=[0..10,0..0.5],title=`Heat Capacity per site C(T)`);" }
}{PARA 13 "" 1 "" {GLPLOT2D 723 224 224 {PLOTDATA 2 "6'-%'CURVESG6$7]q
7$$\"1+++/\\%RT&!#;$\"1+++NS*))4$F*7$$\"1+++OUr#*fF*$\"1+++\"*HpXOF*7$
$\"1+++'**H6^'F*$\"1+++8xM%)RF*7$$\"1+++!oUm*pF*$\"1+++ozd+UF*7$$\"1++
+3*zLY(F*$\"1+++GOHKVF*7$$\"1+++K6%*>zF*$\"1+++[aN,WF*7$$\"1+++_3>s$)F
*$\"1+++z+3AWF*7$$\"1******>G\\C))F*$\"1+++yEe/WF*7$$\"1+++7fL!G*F*$\"
1+++kMQcVF*7$$\"1+++OcoU(*F*$\"1+++[LG$G%F*7$$\"1+++&*)>9-\"!#:$\"1+++
]T))*=%F*7$$\"1+++(fO(p5Fgn$\"1+++!)y\"*zSF*7$$\"1+++%3h%>6Fgn$\"1+++P
fYcRF*7$$\"1+++Cp$3<\"Fgn$\"1+++\")*4@#QF*7$$\"1+++?`6C7Fgn$\"1++++,/z
OF*7$$\"1+++O4cz7Fgn$\"1+++MY8HNF*7$$\"1+++`vXP8Fgn$\"1+++H#=SP$F*7$$
\"1+++TR6)R\"Fgn$\"1+++mx5:KF*7$$\"1+++z,(=Y\"Fgn$\"1+++Tvk`IF*7$$\"1+
++G[5H:Fgn$\"1+++*GO2*GF*7$$\"1+++IKC+;Fgn$\"1+++')*\\ts#F*7$$\"1+++St
wv;Fgn$\"1+++I#fVc#F*7$$\"1+++;wAc<Fgn$\"1+++dQa-CF*7$$\"1+++PvDU=Fgn$
\"1+++3ZgUAF*7$$\"1+++S;fM>Fgn$\"1+++MM<&3#F*7$$\"1+++.\")3M?Fgn$\"1++
+s;#3$>F*7$$\"1+++GtvT@Fgn$\"1+++ur1!y\"F*7$$\"1+++C()zeAFgn$\"1+++N(z
Lj\"F*7$$\"1+++I$[mQ#Fgn$\"1+++@l=\"\\\"F*7$$\"1+++#\\Tq_#Fgn$\"1+++6h
(QN\"F*7$$\"1+++Qj4#o#Fgn$\"1+++HF!=A\"F*7$$\"1+++&eEW&GFgn$\"1+++h()G
&4\"F*7$$\"1+++NoHZIFgn$\"1,++tuFY(*!#<7$$\"1+++!pR[E$Fgn$\"1*****>Gs3
g)Fhu7$$\"1+++VfO7NFgn$\"1*****Hj-\">vFhu7$$\"1+++l!=oz$Fgn$\"1+++$\\t
J]'Fhu7$$\"1+++c-XFTFgn$\"1+++'[#3bbFhu7$$\"1+++!Qzo^%Fgn$\"1+++3HjwYF
hu7$$\"1+++![mF)\\Fgn$\"1+++(eX%pQFhu7$$\"1+++/#>1b&Fgn$\"1+++o#o\\8$F
hu7$$\"1+++gkpeiFgn$\"1+++`l[uCFhu7$$\"1+++oh6nrFgn$\"1+++=:8*)=Fhu7$$
\"1+++S:2w$)Fgn$\"1+++*z&))z8Fhu7$$\"1+++y1f15!#9$\"1+++\"*H!fZ*!#=7$$
\"1+++80tf7F^y$\"1+++h\\[HfFay7$$\"1+++Db=\"o\"F^y$\"1+++'>1`;$Fay7$$
\"1+++A8VBDF^y$\"1+++`(ox=\"Fay7$$\"1+++;W%)[]F^y$\"1+++(G**Qe\"Fay7$F
bz$!1+++*\\-8'R!#>7$$!1+++;W%)[]F^y$!1+++<yU#G&Fiz7$$!1+++A8VBDF^yF_z7
$$!1+++Db=\"o\"F^yFjy7$$!1+++80tf7F^yFey7$$!1+++y1f15F^yF_y7$$!1+++S:2
w$)FgnFix7$$!1+++oh6nrFgnFdx7$$!1+++gkpeiFgnF_x7$$!1+++/#>1b&FgnFjw7$$
!1+++![mF)\\FgnFew7$$!1+++!Qzo^%FgnF`w7$$!1+++c-XFTFgnF[w7$$!1+++l!=oz
$FgnFfv7$$!1+++VfO7NFgnFav7$$!1+++!pR[E$FgnF\\v7$$!1+++NoHZIFgnFfu7$$!
1+++&eEW&GFgnFau7$$!1+++Qj4#o#FgnF\\u7$$!1+++#\\Tq_#FgnFgt7$$!1+++I$[m
Q#FgnFbt7$$!1+++C()zeAFgnF]t7$$!1+++GtvT@FgnFhs7$$!1+++.\")3M?FgnFcs7$
$!1+++S;fM>FgnF^s7$$!1+++PvDU=FgnFir7$$!1+++;wAc<FgnFdr7$$!1+++Stwv;Fg
nF_r7$$!1+++IKC+;FgnFjq7$$!1+++G[5H:FgnFeq7$$!1+++z,(=Y\"FgnF`q7$$!1++
+TR6)R\"FgnF[q7$$!1+++`vXP8FgnFfp7$$!1+++O4cz7FgnFap7$$!1+++?`6C7FgnF
\\p7$$!1+++Cp$3<\"FgnFgo7$$!1+++%3h%>6FgnFbo7$$!1+++(fO(p5FgnF]o7$$!1+
++&*)>9-\"FgnFhn7$$!1+++OcoU(*F*FX7$$!1+++7fL!G*F*FS7$$!1******>G\\C))
F*FN7$$!1+++_3>s$)F*FI7$$!1+++K6%*>zF*FD7$$!1+++3*zLY(F*F?7$$!1+++!oUm
*pF*F:7$$!1+++'**H6^'F*F57$$!1+++OUr#*fF*F07$$!1+++/\\%RT&F*F+-%'COLOU
RG6&%$RGBG$\"#5!\"\"\"\"!Fcdl-%&TITLEG6#%<Heat~Capacity~per~site~C(T)G
-%+AXESLABELSG6$%!GF[el-%&STYLEG6#%&POINTG-%%VIEWG6$;Fcdl$FadlFcdl;Fcd
l$\"\"&Fbdl" 1 5 0 1 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 }}}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "C1[1];" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"+NS*))4$!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 6 "T1[1];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+_6$>q%!#5" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 209 "The disatvantage of using lists i
nstead of tables: indexing is out of step if finite differencing is us
ed (loop starts at two instead of one; thus T[2] corresponds to S[1], \+
U[1], C[1] corresponds to T[2], etc." }}{PARA 0 "" 0 "" {TEXT -1 0 "" 
}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 377 "NOTE: WE HAVE SIMPLY MADE THE \+
QUANTITIES DIMENSIONLESS FOR PLOTTING PURPOSES by setting k=1, e=1 (en
ergy quantum), mu=1. To check the dimensions one should put the right \+
units back in. Where entropy is shown, it is really S/k; where the hea
t capacity is shown, it is C/k, where the temperature is shown, it is \+
kT/e for the Einstein solid, and kT/(mu B) for the paramagnet, etc." }
}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "122 1" 0 }
{VIEWOPTS 1 1 0 1 1 1803 }