read `BiVariateMoms.txt`:

f:=1/(1-t-t^2*w):

print(`This for the curve for the limit of the log of the coefficient of`, t^n*w^(n*x) , ` divided by n, in the rational function`,f):
print(``):

plot(fmkNdata(f,t,w,15,100));

print(``):

t0:=time():
print(`This took`, t0, `seconds. `):
print(``):
print(`--------------------------------------------------`):
print(``):



f:=1/(1-(1+w)*t):

print(`This for the curve for the limit of the log of the coefficient of`, t^n*w^(n*x) , ` divided by n, in the rational function`,f):
print(``):

plot(fmkNdata(f,t,w,15,100));

print(``):

t1:=time():
print(`This took`, t1-t0, `seconds. `):
print(``):
print(`--------------------------------------------------`):
print(``):



f:=1/(1-t-t^3*w):

print(`This for the curve for the limit of the log of the coefficient of`, t^n*w^(n*x) , ` divided by n, in the rational function`,f):
print(``):

plot(fmkNdata(f,t,w,15,100));

print(``):

t2:=time():
print(`This took`, t2-t1, `seconds. `):
print(``):
print(`--------------------------------------------------`):
print(``):



f:=MDktw(2,t,w):
print(`This is the generating function for monomer-dimer tilings of an n by 4 rectangle`):
print(``):
print(`This for the curve for the limit of the log of the coefficient of`, t^n*w^(n*x) , ` divided by n, in the rational function`,f):
print(``):

plot(fmkNdata(f,t,w,10,30));

print(``):

t3:=time():
print(`This took`, t3-t2, `seconds. `):
print(``):
print(`--------------------------------------------------`):
print(``):



quit: