read FreeSAP:
 
print(`The Series Expansion for `):
print(` 2D (square-lattice) Self-Avoiding Polygons of Locally Bounded Width `):
K:=6:
L:=100:
for n from 1 to K do
lprint(`The first `, L, `coeffs. of the generating function (in the `):
lprint(` variable s), enumerating  2D (square-lattice) `):
lprint(`Self-Avoiding Polygons `):
lprint(`of width <=`, n , `equals `):
lu:=gfSeriesf(n,L):
print(lu):
print(``):
print(` The square roots of the `):
print(`the two last ratios (giving appx. to the "connective consts.")`):
print(` are: `):
print(evalf(sqrt(lu[L]/lu[L-1])),evalf(sqrt(lu[L-1]/lu[L-2])) ):
od:
quit: