read Words123:

read SCHUTZENBERGER:

K:=250:

for i from 2 to 4 do

print(` r is `, i):

lu:=SeqVR(i,K):

gu:=1+add(lu[j]*x^j,j=1..nops(lu)):

gu:=expand(gu^2):

lu:=[seq(coeff(gu,x,i),i=0..K)]:

ku:=numer(Empir(lu,x,P)):

if ku<>0 then

ku:=add(factor(coeff(ku,P,i))*P^i,i=0..degree(ku,P)):
ku:=subs(P=F^2,ku):
print(`The minimal equation is`):
print(ku):

print(`and in Maple input form`):

lprint(ku):
else
print(`Nothing found`):
fi:
od:


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


quit: