[f] = mfeigensearch([54,4],[[5,3]]);
[g] = mfeigensearch([54,2],[[5,3]]);
[g1] = mfeigensearch([54,2],[[5,-3]]);
g2 = mffromell(S(0,1))[2];

primeList = primes(9)[3..9];

print("Hecke eigenvalues of f, g, g1 and g2:");
apply(p->print([p,mfcoef(f,p),mfcoef(g,p),mfcoef(g1,p),mfcoef(g2,p)]),primeList);