Help:=proc(): print(`FCmt(L)`): end: #inputs a list L, conjectured to be asmpt. to L[n]=C*mu^n*n^theta #tried to find appx. for C, mu, and theta #log(L[n])=log(C)+n *log(mu)+ theta*log(n) FCmt:=proc(L) local n,eq,var,logC, logmu, theta: n:=nops(L): var:={logC,logmu,theta}: eq:=evalf({log(L[n])-logC-n*logmu-theta*log(n), log(L[n-1])-logC-(n-1)*logmu-theta*log(n-1), log(L[n-2])-logC-(n-2)*logmu-theta*log(n-2)}): var:=solve(eq,var): logC:=subs(var,logC): logmu:=subs(var,logmu): theta:=subs(var,theta): [exp(logC),exp(logmu),theta]: end: