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: