read Cfinite:

print(`The following verifies the Kasteleyn's formula (Eq. (15) in his classic paper), for various numeric z0 for all even n from 2 to 10`):

for z0 from 1 to 4 do
print(`When z0=`, z0):
for n from 2 to 10 by 2 do
print(`When n=,`, n, `the formula is`):
print(evalb(KastF1(n,z0)=DimerData(z0)[n])):
od:

od:

print(`This took`, time(), `seconds. `):



quit: