/*
  This code can be loaded, or copied and pasted, into Magma.
  It will load the data associated to the HMF, including
  the field, level, and Hecke and Atkin-Lehner eigenvalue data.
  At the *bottom* of the file, there is code to recreate the
  Hilbert modular form in Magma, by creating the HMF space
  and cutting out the corresponding Hecke irreducible subspace.
  From there, you can ask for more eigenvalues or modify as desired.
  It is commented out, as this computation may be lengthy.
*/

P<x> := PolynomialRing(Rationals());
g := P![25, 5, -11, -1, 1];
F<w> := NumberField(g);
ZF := Integers(F);

NN := ideal<ZF | {25, 5, 2/5*w^3 - 2/5*w^2 - 12/5*w + 1}>;

primesArray := [
[5, 5, 1/5*w^3 + 4/5*w^2 - 11/5*w - 7],
[5, 5, 1/5*w^3 - 6/5*w^2 - 1/5*w + 4],
[11, 11, w + 1],
[11, 11, 1/5*w^3 - 1/5*w^2 - 11/5*w + 2],
[16, 2, 2],
[19, 19, -1/5*w^3 + 1/5*w^2 + 1/5*w + 1],
[19, 19, 1/5*w^3 - 1/5*w^2 - 11/5*w],
[19, 19, -2/5*w^3 + 2/5*w^2 + 17/5*w],
[19, 19, w - 1],
[29, 29, w^2 - w - 8],
[29, 29, 2/5*w^3 - 2/5*w^2 - 7/5*w + 2],
[31, 31, 1/5*w^3 - 1/5*w^2 - 1/5*w + 2],
[31, 31, 2/5*w^3 - 2/5*w^2 - 17/5*w + 3],
[61, 61, -2/5*w^3 - 3/5*w^2 + 12/5*w + 5],
[61, 61, -4/5*w^3 - 6/5*w^2 + 39/5*w + 15],
[61, 61, 3/5*w^3 + 2/5*w^2 - 18/5*w - 3],
[61, 61, 7/5*w^3 + 3/5*w^2 - 62/5*w - 13],
[71, 71, 1/5*w^3 - 6/5*w^2 - 6/5*w + 4],
[71, 71, w^2 - 8],
[79, 79, 3/5*w^3 + 12/5*w^2 - 33/5*w - 22],
[79, 79, 1/5*w^3 + 4/5*w^2 - 11/5*w - 9],
[81, 3, -3],
[101, 101, 3/5*w^3 + 7/5*w^2 - 28/5*w - 15],
[101, 101, 1/5*w^3 + 9/5*w^2 - 16/5*w - 16],
[121, 11, 3/5*w^3 - 3/5*w^2 - 18/5*w + 1],
[131, 131, 1/5*w^3 + 4/5*w^2 - 16/5*w - 8],
[131, 131, w^2 - 2*w - 2],
[149, 149, 4/5*w^3 + 1/5*w^2 - 29/5*w - 5],
[149, 149, 4/5*w^3 - 9/5*w^2 - 19/5*w + 6],
[151, 151, 2/5*w^3 + 3/5*w^2 - 17/5*w - 3],
[151, 151, -2/5*w^3 + 7/5*w^2 + 7/5*w - 8],
[151, 151, 4/5*w^3 + 1/5*w^2 - 29/5*w - 2],
[151, 151, -4/5*w^3 + 9/5*w^2 + 19/5*w - 9],
[179, 179, -3/5*w^3 + 8/5*w^2 + 18/5*w - 6],
[179, 179, 2/5*w^3 + 3/5*w^2 - 12/5*w - 6],
[211, 211, 1/5*w^3 + 4/5*w^2 - 6/5*w - 10],
[211, 211, 2/5*w^3 - 7/5*w^2 - 12/5*w + 2],
[229, 229, 6/5*w^3 - 11/5*w^2 - 41/5*w + 11],
[229, 229, 2*w^2 - 2*w - 9],
[229, 229, -4/5*w^3 - 1/5*w^2 + 19/5*w + 2],
[229, 229, -3/5*w^3 + 3/5*w^2 + 28/5*w - 3],
[239, 239, 11/5*w^3 + 4/5*w^2 - 96/5*w - 19],
[239, 239, -6/5*w^3 + 1/5*w^2 + 51/5*w + 3],
[241, 241, 1/5*w^3 + 4/5*w^2 - 1/5*w - 6],
[241, 241, 4/5*w^3 - 4/5*w^2 - 29/5*w - 1],
[251, 251, -1/5*w^3 + 6/5*w^2 + 11/5*w - 12],
[251, 251, -1/5*w^3 + 11/5*w^2 + 1/5*w - 15],
[269, 269, 2/5*w^3 + 3/5*w^2 - 17/5*w - 2],
[269, 269, -2/5*w^3 + 7/5*w^2 + 7/5*w - 9],
[271, 271, 1/5*w^3 - 1/5*w^2 + 4/5*w],
[271, 271, 3/5*w^3 - 3/5*w^2 - 28/5*w + 2],
[281, 281, w^3 - 3*w^2 - 4*w + 13],
[281, 281, -3/5*w^3 + 3/5*w^2 + 28/5*w - 6],
[281, 281, -6/5*w^3 + 26/5*w^2 + 11/5*w - 16],
[281, 281, w^3 + w^2 - 8*w - 9],
[289, 17, 2/5*w^3 + 3/5*w^2 - 7/5*w - 5],
[289, 17, 4/5*w^3 - 9/5*w^2 - 29/5*w + 8],
[311, 311, -w^2 + 2*w + 11],
[311, 311, -1/5*w^3 + 6/5*w^2 + 1/5*w - 1],
[311, 311, 1/5*w^3 + 4/5*w^2 - 16/5*w + 1],
[311, 311, -1/5*w^3 - 4/5*w^2 + 11/5*w + 10],
[349, 349, -1/5*w^3 - 9/5*w^2 + 11/5*w + 16],
[349, 349, 2/5*w^3 - 12/5*w^2 - 7/5*w + 7],
[359, 359, -4/5*w^3 + 4/5*w^2 + 19/5*w - 1],
[359, 359, -w^3 + 3*w^2 + 5*w - 12],
[359, 359, w^3 - w^2 - 7*w + 2],
[359, 359, 4/5*w^3 + 6/5*w^2 - 29/5*w - 11],
[379, 379, 4/5*w^3 + 6/5*w^2 - 34/5*w - 9],
[379, 379, 3/5*w^3 - 3/5*w^2 - 28/5*w + 7],
[389, 389, w^2 - 3*w - 8],
[389, 389, -4/5*w^3 + 14/5*w^2 + 19/5*w - 16],
[401, 401, -3*w^2 + 2*w + 23],
[401, 401, 1/5*w^3 - 16/5*w^2 + 4/5*w + 11],
[409, 409, -w^3 + 2*w^2 + 3*w - 7],
[409, 409, 7/5*w^3 - 2/5*w^2 - 57/5*w - 2],
[419, 419, -1/5*w^3 + 1/5*w^2 + 16/5*w - 2],
[419, 419, 1/5*w^3 - 1/5*w^2 - 16/5*w],
[421, 421, 2/5*w^3 + 3/5*w^2 - 27/5*w - 3],
[421, 421, 1/5*w^3 + 4/5*w^2 + 4/5*w - 3],
[439, 439, -w^3 + 3*w^2 + 5*w - 16],
[439, 439, -4/5*w^3 - 6/5*w^2 + 29/5*w + 7],
[439, 439, -3/5*w^3 + 13/5*w^2 + 3/5*w - 11],
[439, 439, 4/5*w^3 + 6/5*w^2 - 39/5*w - 10],
[449, 449, 2/5*w^3 - 12/5*w^2 - 2/5*w + 5],
[449, 449, 4/5*w^3 - 24/5*w^2 - 4/5*w + 15],
[461, 461, -1/5*w^3 - 4/5*w^2 - 4/5*w + 4],
[461, 461, 2*w^2 - 2*w - 11],
[461, 461, 3/5*w^3 + 2/5*w^2 - 33/5*w - 4],
[479, 479, -w^3 + 3*w^2 + 6*w - 14],
[479, 479, -3/5*w^3 + 3/5*w^2 + 28/5*w + 3],
[499, 499, 2/5*w^3 + 8/5*w^2 - 17/5*w - 6],
[499, 499, -3/5*w^3 + 13/5*w^2 + 13/5*w - 17],
[509, 509, -w^3 + 2*w^2 + 4*w - 6],
[509, 509, 6/5*w^3 - 1/5*w^2 - 46/5*w - 4],
[521, 521, -2/5*w^3 + 2/5*w^2 + 7/5*w - 7],
[521, 521, -3/5*w^3 + 3/5*w^2 + 23/5*w - 8],
[529, 23, -7/5*w^3 + 17/5*w^2 + 42/5*w - 15],
[529, 23, 2/5*w^3 + 13/5*w^2 - 27/5*w - 21],
[541, 541, 1/5*w^3 - 11/5*w^2 - 1/5*w + 10],
[541, 541, 3/5*w^3 - 8/5*w^2 - 13/5*w + 13],
[541, 541, -2*w^2 + w + 13],
[541, 541, 3/5*w^3 + 2/5*w^2 - 23/5*w + 2],
[569, 569, w^3 - 2*w^2 - 5*w + 7],
[569, 569, -1/5*w^3 + 1/5*w^2 + 1/5*w - 6],
[571, 571, w^2 + w - 9],
[571, 571, -2/5*w^3 + 12/5*w^2 + 7/5*w - 12],
[571, 571, 1/5*w^3 + 9/5*w^2 - 11/5*w - 11],
[571, 571, -2/5*w^3 + 7/5*w^2 + 17/5*w - 4],
[601, 601, -3/5*w^3 + 13/5*w^2 + 18/5*w - 11],
[601, 601, -4/5*w^3 + 4/5*w^2 + 34/5*w + 1],
[631, 631, 2/5*w^3 + 18/5*w^2 - 27/5*w - 30],
[631, 631, 3/5*w^3 - 23/5*w^2 - 3/5*w + 15],
[661, 661, -2*w^3 + 17*w + 12],
[661, 661, -9/5*w^3 + 4/5*w^2 + 74/5*w + 1],
[691, 691, -7/5*w^3 - 18/5*w^2 + 72/5*w + 35],
[691, 691, 6/5*w^3 - 31/5*w^2 - 6/5*w + 19],
[701, 701, -1/5*w^3 + 6/5*w^2 + 16/5*w - 6],
[701, 701, 2/5*w^3 - 7/5*w^2 - 22/5*w + 8],
[709, 709, -3/5*w^3 + 13/5*w^2 + 13/5*w - 18],
[709, 709, 1/5*w^3 + 14/5*w^2 - 11/5*w - 20],
[709, 709, -2/5*w^3 - 8/5*w^2 + 17/5*w + 5],
[709, 709, 1/5*w^3 + 4/5*w^2 - 21/5*w - 10],
[719, 719, -3/5*w^3 + 8/5*w^2 + 13/5*w - 12],
[719, 719, 3/5*w^3 + 2/5*w^2 - 23/5*w + 1],
[719, 719, -7/5*w^3 + 17/5*w^2 + 32/5*w - 14],
[719, 719, 7/5*w^3 + 3/5*w^2 - 52/5*w - 8],
[739, 739, 1/5*w^3 - 1/5*w^2 - 11/5*w - 5],
[739, 739, w - 6],
[761, 761, -4/5*w^3 + 9/5*w^2 + 19/5*w - 11],
[761, 761, 4/5*w^3 + 1/5*w^2 - 29/5*w],
[769, 769, 3/5*w^3 - 13/5*w^2 - 13/5*w + 8],
[769, 769, -2/5*w^3 - 8/5*w^2 + 17/5*w + 15],
[769, 769, 7/5*w^3 - 17/5*w^2 - 37/5*w + 18],
[769, 769, -6/5*w^3 - 4/5*w^2 + 41/5*w + 5],
[809, 809, -1/5*w^3 + 6/5*w^2 - 4/5*w - 14],
[809, 809, -2/5*w^3 - 3/5*w^2 + 22/5*w - 4],
[809, 809, -6/5*w^3 + 6/5*w^2 + 41/5*w - 5],
[809, 809, w^3 - w^2 - 5*w + 4],
[821, 821, -3/5*w^3 - 17/5*w^2 + 33/5*w + 30],
[821, 821, -4/5*w^3 + 24/5*w^2 + 9/5*w - 15],
[829, 829, 7/5*w^3 - 12/5*w^2 - 52/5*w + 10],
[829, 829, -1/5*w^3 + 11/5*w^2 - 4/5*w - 19],
[829, 829, -13/5*w^3 + 3/5*w^2 + 108/5*w + 10],
[829, 829, -w^3 + 4*w^2 + 5*w - 23],
[839, 839, 4/5*w^3 - 9/5*w^2 - 24/5*w + 4],
[839, 839, 3/5*w^3 + 2/5*w^2 - 18/5*w - 8],
[841, 29, -w^3 + w^2 + 6*w - 4],
[859, 859, 6/5*w^3 - 11/5*w^2 - 36/5*w + 9],
[859, 859, w^3 - 6*w - 3],
[881, 881, 4/5*w^3 - 4/5*w^2 - 9/5*w - 1],
[881, 881, 9/5*w^3 - 24/5*w^2 - 19/5*w + 14],
[881, 881, -w^3 + 6*w + 2],
[881, 881, 1/5*w^3 - 1/5*w^2 - 1/5*w - 6],
[911, 911, w^3 - w^2 - 5*w + 1],
[911, 911, 6/5*w^3 - 6/5*w^2 - 41/5*w + 2],
[919, 919, -7/5*w^3 + 12/5*w^2 + 42/5*w - 7],
[919, 919, -6/5*w^3 + 1/5*w^2 + 36/5*w + 5],
[929, 929, -2/5*w^3 - 8/5*w^2 + 17/5*w + 16],
[929, 929, -3/5*w^3 + 13/5*w^2 + 13/5*w - 7],
[961, 31, w^3 - w^2 - 6*w + 2],
[971, 971, 8/5*w^3 - 13/5*w^2 - 48/5*w + 6],
[971, 971, -2/5*w^3 - 3/5*w^2 + 22/5*w],
[971, 971, 7/5*w^3 - 7/5*w^2 - 57/5*w + 7],
[971, 971, -1/5*w^3 + 6/5*w^2 - 4/5*w - 10],
[991, 991, 2/5*w^3 + 8/5*w^2 - 22/5*w - 7],
[991, 991, -2/5*w^3 + 12/5*w^2 + 2/5*w - 15]];
primes := [ideal<ZF | {F!x : x in I}> : I in primesArray];

heckePol := x^7 - 4*x^6 - 29*x^5 + 81*x^4 + 207*x^3 - 429*x^2 - 197*x - 12;
K<e> := NumberField(heckePol);

heckeEigenvaluesArray := [1, 1, e, e, -314/9127*e^6 + 2016/9127*e^5 + 7482/9127*e^4 - 48194/9127*e^3 - 53049/9127*e^2 + 255780/9127*e + 74861/9127, -157/9127*e^6 + 1008/9127*e^5 + 3741/9127*e^4 - 24097/9127*e^3 - 31088/9127*e^2 + 137017/9127*e + 60248/9127, 59/9127*e^6 - 30/9127*e^5 - 3208/9127*e^4 + 5335/9127*e^3 + 27437/9127*e^2 - 69977/9127*e + 380/9127, -157/9127*e^6 + 1008/9127*e^5 + 3741/9127*e^4 - 24097/9127*e^3 - 31088/9127*e^2 + 137017/9127*e + 60248/9127, 59/9127*e^6 - 30/9127*e^5 - 3208/9127*e^4 + 5335/9127*e^3 + 27437/9127*e^2 - 69977/9127*e + 380/9127, 655/9127*e^6 - 1880/9127*e^5 - 21537/9127*e^4 + 36178/9127*e^3 + 155626/9127*e^2 - 211144/9127*e - 79626/9127, 655/9127*e^6 - 1880/9127*e^5 - 21537/9127*e^4 + 36178/9127*e^3 + 155626/9127*e^2 - 211144/9127*e - 79626/9127, 445/9127*e^6 - 2392/9127*e^5 - 9964/9127*e^4 + 48128/9127*e^3 + 51317/9127*e^2 - 239596/9127*e + 5960/9127, 445/9127*e^6 - 2392/9127*e^5 - 9964/9127*e^4 + 48128/9127*e^3 + 51317/9127*e^2 - 239596/9127*e + 5960/9127, 625/9127*e^6 - 3257/9127*e^5 - 17276/9127*e^4 + 75697/9127*e^3 + 138117/9127*e^2 - 430345/9127*e - 153454/9127, 625/9127*e^6 - 3257/9127*e^5 - 17276/9127*e^4 + 75697/9127*e^3 + 138117/9127*e^2 - 430345/9127*e - 153454/9127, 1640/9127*e^6 - 6867/9127*e^5 - 44310/9127*e^4 + 130505/9127*e^3 + 284803/9127*e^2 - 630526/9127*e - 168574/9127, 1640/9127*e^6 - 6867/9127*e^5 - 44310/9127*e^4 + 130505/9127*e^3 + 284803/9127*e^2 - 630526/9127*e - 168574/9127, -1539/9127*e^6 + 5114/9127*e^5 + 46089/9127*e^4 - 91052/9127*e^3 - 313171/9127*e^2 + 473144/9127*e + 166440/9127, -1539/9127*e^6 + 5114/9127*e^5 + 46089/9127*e^4 - 91052/9127*e^3 - 313171/9127*e^2 + 473144/9127*e + 166440/9127, -285/9127*e^6 + 609/9127*e^5 + 8535/9127*e^4 - 3340/9127*e^3 - 52248/9127*e^2 - 24271/9127*e + 47048/9127, -285/9127*e^6 + 609/9127*e^5 + 8535/9127*e^4 - 3340/9127*e^3 - 52248/9127*e^2 - 24271/9127*e + 47048/9127, -730/9127*e^6 + 3001/9127*e^5 + 18499/9127*e^4 - 51468/9127*e^3 - 103565/9127*e^2 + 233579/9127*e + 132358/9127, -2115/9127*e^6 + 7882/9127*e^5 + 58535/9127*e^4 - 139114/9127*e^3 - 362756/9127*e^2 + 669175/9127*e + 161802/9127, -2115/9127*e^6 + 7882/9127*e^5 + 58535/9127*e^4 - 139114/9127*e^3 - 362756/9127*e^2 + 669175/9127*e + 161802/9127, -2157/9127*e^6 + 9605/9127*e^5 + 62675/9127*e^4 - 218867/9127*e^3 - 443856/9127*e^2 + 1267692/9127*e + 292094/9127, -586/9127*e^6 + 2309/9127*e^5 + 19951/9127*e^4 - 62270/9127*e^3 - 161903/9127*e^2 + 415028/9127*e + 110700/9127, -586/9127*e^6 + 2309/9127*e^5 + 19951/9127*e^4 - 62270/9127*e^3 - 161903/9127*e^2 + 415028/9127*e + 110700/9127, 272/9127*e^6 - 293/9127*e^5 - 12469/9127*e^4 + 14076/9127*e^3 + 99727/9127*e^2 - 168375/9127*e + 28050/9127, 272/9127*e^6 - 293/9127*e^5 - 12469/9127*e^4 + 14076/9127*e^3 + 99727/9127*e^2 - 168375/9127*e + 28050/9127, -259/9127*e^6 - 23/9127*e^5 + 7276/9127*e^4 + 20823/9127*e^3 - 19428/9127*e^2 - 177472/9127*e - 66640/9127, -259/9127*e^6 - 23/9127*e^5 + 7276/9127*e^4 + 20823/9127*e^3 - 19428/9127*e^2 - 177472/9127*e - 66640/9127, -157/9127*e^6 + 1008/9127*e^5 + 3741/9127*e^4 - 24097/9127*e^3 - 21961/9127*e^2 + 137017/9127*e - 49276/9127, -157/9127*e^6 + 1008/9127*e^5 + 3741/9127*e^4 - 24097/9127*e^3 - 21961/9127*e^2 + 137017/9127*e - 49276/9127, -612/9127*e^6 + 2941/9127*e^5 + 21210/9127*e^4 - 86433/9127*e^3 - 185596/9127*e^2 + 559102/9127*e + 114864/9127, -612/9127*e^6 + 2941/9127*e^5 + 21210/9127*e^4 - 86433/9127*e^3 - 185596/9127*e^2 + 559102/9127*e + 114864/9127, -1578/9127*e^6 + 6062/9127*e^5 + 43414/9127*e^4 - 113606/9127*e^3 - 252877/9127*e^2 + 588858/9127*e - 28108/9127, -1578/9127*e^6 + 6062/9127*e^5 + 43414/9127*e^4 - 113606/9127*e^3 - 252877/9127*e^2 + 588858/9127*e - 28108/9127, 717/9127*e^6 - 2685/9127*e^5 - 22433/9127*e^4 + 62204/9127*e^3 + 151044/9127*e^2 - 371463/9127*e + 34010/9127, 3372/9127*e^6 - 13162/9127*e^5 - 93777/9127*e^4 + 247517/9127*e^3 + 600787/9127*e^2 - 1247805/9127*e - 277462/9127, 717/9127*e^6 - 2685/9127*e^5 - 22433/9127*e^4 + 62204/9127*e^3 + 151044/9127*e^2 - 371463/9127*e + 34010/9127, 3372/9127*e^6 - 13162/9127*e^5 - 93777/9127*e^4 + 247517/9127*e^3 + 600787/9127*e^2 - 1247805/9127*e - 277462/9127, -973/9127*e^6 + 1887/9127*e^5 + 32021/9127*e^4 - 20690/9127*e^3 - 193364/9127*e^2 + 103649/9127*e - 42156/9127, -973/9127*e^6 + 1887/9127*e^5 + 32021/9127*e^4 - 20690/9127*e^3 - 193364/9127*e^2 + 103649/9127*e - 42156/9127, 2053/9127*e^6 - 7077/9127*e^5 - 57639/9127*e^4 + 113088/9127*e^3 + 367338/9127*e^2 - 517983/9127*e - 238930/9127, 2053/9127*e^6 - 7077/9127*e^5 - 57639/9127*e^4 + 113088/9127*e^3 + 367338/9127*e^2 - 517983/9127*e - 238930/9127, -1405/9127*e^6 + 3963/9127*e^5 + 45919/9127*e^4 - 79554/9127*e^3 - 310414/9127*e^2 + 526764/9127*e + 77580/9127, -1405/9127*e^6 + 3963/9127*e^5 + 45919/9127*e^4 - 79554/9127*e^3 - 310414/9127*e^2 + 526764/9127*e + 77580/9127, -2409/9127*e^6 + 10816/9127*e^5 + 69261/9127*e^4 - 241035/9127*e^3 - 492360/9127*e^2 + 1326645/9127*e + 234162/9127, -2409/9127*e^6 + 10816/9127*e^5 + 69261/9127*e^4 - 241035/9127*e^3 - 492360/9127*e^2 + 1326645/9127*e + 234162/9127, 468/9127*e^6 - 2249/9127*e^5 - 13535/9127*e^4 + 51600/9127*e^3 + 97902/9127*e^2 - 247693/9127*e - 111460/9127, 468/9127*e^6 - 2249/9127*e^5 - 13535/9127*e^4 + 51600/9127*e^3 + 97902/9127*e^2 - 247693/9127*e - 111460/9127, -419/9127*e^6 + 1760/9127*e^5 + 17832/9127*e^4 - 60473/9127*e^3 - 191910/9127*e^2 + 414967/9127*e + 227178/9127, 55/9127*e^6 - 2039/9127*e^5 - 206/9127*e^4 + 69017/9127*e^3 + 33621/9127*e^2 - 433252/9127*e - 59358/9127, 55/9127*e^6 - 2039/9127*e^5 - 206/9127*e^4 + 69017/9127*e^3 + 33621/9127*e^2 - 433252/9127*e - 59358/9127, -419/9127*e^6 + 1760/9127*e^5 + 17832/9127*e^4 - 60473/9127*e^3 - 191910/9127*e^2 + 414967/9127*e + 227178/9127, 2389/9127*e^6 - 11734/9127*e^5 - 63378/9127*e^4 + 258254/9127*e^3 + 422883/9127*e^2 - 1445398/9127*e - 259042/9127, 2389/9127*e^6 - 11734/9127*e^5 - 63378/9127*e^4 + 258254/9127*e^3 + 422883/9127*e^2 - 1445398/9127*e - 259042/9127, 812/9127*e^6 - 2888/9127*e^5 - 25278/9127*e^4 + 60275/9127*e^3 + 186714/9127*e^2 - 329907/9127*e - 121620/9127, 1185/9127*e^6 - 4934/9127*e^5 - 35968/9127*e^4 + 113804/9127*e^3 + 258073/9127*e^2 - 664791/9127*e - 59196/9127, 812/9127*e^6 - 2888/9127*e^5 - 25278/9127*e^4 + 60275/9127*e^3 + 186714/9127*e^2 - 329907/9127*e - 121620/9127, 1185/9127*e^6 - 4934/9127*e^5 - 35968/9127*e^4 + 113804/9127*e^3 + 258073/9127*e^2 - 664791/9127*e - 59196/9127, 954/9127*e^6 - 21/9127*e^5 - 31452/9127*e^4 - 55591/9127*e^3 + 167976/9127*e^2 + 422882/9127*e + 266/9127, 954/9127*e^6 - 21/9127*e^5 - 31452/9127*e^4 - 55591/9127*e^3 + 167976/9127*e^2 + 422882/9127*e + 266/9127, -2509/9127*e^6 + 6226/9127*e^5 + 80422/9127*e^4 - 91051/9127*e^3 - 547681/9127*e^2 + 486451/9127*e + 347064/9127, -223/9127*e^6 - 196/9127*e^5 + 7639/9127*e^4 + 13559/9127*e^3 - 20322/9127*e^2 - 65939/9127*e - 94872/9127, -2509/9127*e^6 + 6226/9127*e^5 + 80422/9127*e^4 - 91051/9127*e^3 - 547681/9127*e^2 + 486451/9127*e + 347064/9127, -223/9127*e^6 - 196/9127*e^5 + 7639/9127*e^4 + 13559/9127*e^3 - 20322/9127*e^2 - 65939/9127*e - 94872/9127, 356/9127*e^6 - 3739/9127*e^5 - 2495/9127*e^4 + 82312/9127*e^3 - 11883/9127*e^2 - 343185/9127*e + 41276/9127, 356/9127*e^6 - 3739/9127*e^5 - 2495/9127*e^4 + 82312/9127*e^3 - 11883/9127*e^2 - 343185/9127*e + 41276/9127, 1391/9127*e^6 - 6431/9127*e^5 - 35412/9127*e^4 + 119901/9127*e^3 + 249915/9127*e^2 - 506756/9127*e - 332298/9127, 1391/9127*e^6 - 6431/9127*e^5 - 35412/9127*e^4 + 119901/9127*e^3 + 249915/9127*e^2 - 506756/9127*e - 332298/9127, 3011/9127*e^6 - 14216/9127*e^5 - 73839/9127*e^4 + 267625/9127*e^3 + 437860/9127*e^2 - 1210400/9127*e - 178926/9127, 3011/9127*e^6 - 14216/9127*e^5 - 73839/9127*e^4 + 267625/9127*e^3 + 437860/9127*e^2 - 1210400/9127*e - 178926/9127, 1303/9127*e^6 - 4994/9127*e^5 - 42384/9127*e^4 + 115347/9127*e^3 + 349455/9127*e^2 - 695221/9127*e - 405262/9127, 1303/9127*e^6 - 4994/9127*e^5 - 42384/9127*e^4 + 115347/9127*e^3 + 349455/9127*e^2 - 695221/9127*e - 405262/9127, -2947/9127*e^6 + 9852/9127*e^5 + 89696/9127*e^4 - 182170/9127*e^3 - 618947/9127*e^2 + 925964/9127*e + 349812/9127, -2947/9127*e^6 + 9852/9127*e^5 + 89696/9127*e^4 - 182170/9127*e^3 - 618947/9127*e^2 + 925964/9127*e + 349812/9127, -2371/9127*e^6 + 7084/9127*e^5 + 77250/9127*e^4 - 143235/9127*e^3 - 541981/9127*e^2 + 912473/9127*e + 208418/9127, -2371/9127*e^6 + 7084/9127*e^5 + 77250/9127*e^4 - 143235/9127*e^3 - 541981/9127*e^2 + 912473/9127*e + 208418/9127, -4744/9127*e^6 + 17727/9127*e^5 + 137747/9127*e^4 - 345899/9127*e^3 - 907457/9127*e^2 + 1903428/9127*e + 377840/9127, -4744/9127*e^6 + 17727/9127*e^5 + 137747/9127*e^4 - 345899/9127*e^3 - 907457/9127*e^2 + 1903428/9127*e + 377840/9127, -1817/9127*e^6 + 6957/9127*e^5 + 44807/9127*e^4 - 100875/9127*e^3 - 221082/9127*e^2 + 338472/9127*e - 69868/9127, -1817/9127*e^6 + 6957/9127*e^5 + 44807/9127*e^4 - 100875/9127*e^3 - 221082/9127*e^2 + 338472/9127*e - 69868/9127, -2261/9127*e^6 + 12133/9127*e^5 + 58584/9127*e^4 - 279011/9127*e^3 - 401723/9127*e^2 + 1588432/9127*e + 235734/9127, -2261/9127*e^6 + 12133/9127*e^5 + 58584/9127*e^4 - 279011/9127*e^3 - 401723/9127*e^2 + 1588432/9127*e + 235734/9127, -442/9127*e^6 + 1617/9127*e^5 + 12276/9127*e^4 - 27437/9127*e^3 - 74209/9127*e^2 + 76238/9127*e + 16026/9127, -2378/9127*e^6 + 5850/9127*e^5 + 77940/9127*e^4 - 100244/9127*e^3 - 503778/9127*e^2 + 621286/9127*e - 128862/9127, -442/9127*e^6 + 1617/9127*e^5 + 12276/9127*e^4 - 27437/9127*e^3 - 74209/9127*e^2 + 76238/9127*e + 16026/9127, 2108/9127*e^6 - 9116/9127*e^5 - 66972/9127*e^4 + 227740/9127*e^3 + 537864/9127*e^2 - 1398458/9127*e - 371304/9127, 2108/9127*e^6 - 9116/9127*e^5 - 66972/9127*e^4 + 227740/9127*e^3 + 537864/9127*e^2 - 1398458/9127*e - 371304/9127, -4322/9127*e^6 + 15192/9127*e^5 + 122227/9127*e^4 - 255608/9127*e^3 - 774947/9127*e^2 + 1151690/9127*e + 428204/9127, -4322/9127*e^6 + 15192/9127*e^5 + 122227/9127*e^4 - 255608/9127*e^3 - 774947/9127*e^2 + 1151690/9127*e + 428204/9127, -3313/9127*e^6 + 13132/9127*e^5 + 99696/9127*e^4 - 287817/9127*e^3 - 710255/9127*e^2 + 1716321/9127*e + 387366/9127, -3313/9127*e^6 + 13132/9127*e^5 + 99696/9127*e^4 - 287817/9127*e^3 - 710255/9127*e^2 + 1716321/9127*e + 387366/9127, 173/9127*e^6 - 2099/9127*e^5 + 2505/9127*e^4 + 43179/9127*e^3 - 94045/9127*e^2 - 171618/9127*e + 342990/9127, 173/9127*e^6 - 2099/9127*e^5 + 2505/9127*e^4 + 43179/9127*e^3 - 94045/9127*e^2 - 171618/9127*e + 342990/9127, 439/9127*e^6 - 842/9127*e^5 - 14588/9127*e^4 + 6746/9127*e^3 + 115355/9127*e^2 - 40658/9127*e - 165790/9127, 439/9127*e^6 - 842/9127*e^5 - 14588/9127*e^4 + 6746/9127*e^3 + 115355/9127*e^2 - 40658/9127*e - 165790/9127, 583/9127*e^6 - 1534/9127*e^5 - 13136/9127*e^4 + 5071/9127*e^3 + 2255/9127*e^2 + 40394/9127*e + 341918/9127, -4409/9127*e^6 + 19413/9127*e^5 + 119068/9127*e^4 - 390170/9127*e^3 - 786477/9127*e^2 + 1964462/9127*e + 484262/9127, 583/9127*e^6 - 1534/9127*e^5 - 13136/9127*e^4 + 5071/9127*e^3 + 2255/9127*e^2 + 40394/9127*e + 341918/9127, -4409/9127*e^6 + 19413/9127*e^5 + 119068/9127*e^4 - 390170/9127*e^3 - 786477/9127*e^2 + 1964462/9127*e + 484262/9127, 1278/9127*e^6 - 1578/9127*e^5 - 46439/9127*e^4 + 15938/9127*e^3 + 296835/9127*e^2 - 170546/9127*e + 111258/9127, 1278/9127*e^6 - 1578/9127*e^5 - 46439/9127*e^4 + 15938/9127*e^3 + 296835/9127*e^2 - 170546/9127*e + 111258/9127, 3343/9127*e^6 - 11755/9127*e^5 - 103957/9127*e^4 + 248298/9127*e^3 + 727764/9127*e^2 - 1487993/9127*e - 331792/9127, 4520/9127*e^6 - 20707/9127*e^5 - 124794/9127*e^4 + 452958/9127*e^3 + 843046/9127*e^2 - 2505127/9127*e - 583480/9127, 4520/9127*e^6 - 20707/9127*e^5 - 124794/9127*e^4 + 452958/9127*e^3 + 843046/9127*e^2 - 2505127/9127*e - 583480/9127, 3343/9127*e^6 - 11755/9127*e^5 - 103957/9127*e^4 + 248298/9127*e^3 + 727764/9127*e^2 - 1487993/9127*e - 331792/9127, -1392/9127*e^6 + 3647/9127*e^5 + 40726/9127*e^4 - 35528/9127*e^3 - 257496/9127*e^2 + 62266/9127*e + 148514/9127, -1392/9127*e^6 + 3647/9127*e^5 + 40726/9127*e^4 - 35528/9127*e^3 - 257496/9127*e^2 + 62266/9127*e + 148514/9127, -67/9127*e^6 + 5139/9127*e^5 - 9042/9127*e^4 - 151781/9127*e^3 + 67074/9127*e^2 + 867636/9127*e - 46840/9127, -67/9127*e^6 + 5139/9127*e^5 - 9042/9127*e^4 - 151781/9127*e^3 + 67074/9127*e^2 + 867636/9127*e - 46840/9127, 483/9127*e^6 - 6124/9127*e^5 - 1975/9127*e^4 + 145928/9127*e^3 + 19950/9127*e^2 - 699403/9127*e - 257086/9127, 483/9127*e^6 - 6124/9127*e^5 - 1975/9127*e^4 + 145928/9127*e^3 + 19950/9127*e^2 - 699403/9127*e - 257086/9127, 1702/9127*e^6 - 7672/9127*e^5 - 45206/9127*e^4 + 156531/9127*e^3 + 271094/9127*e^2 - 772591/9127*e - 73192/9127, 1702/9127*e^6 - 7672/9127*e^5 - 45206/9127*e^4 + 156531/9127*e^3 + 271094/9127*e^2 - 772591/9127*e - 73192/9127, 490/9127*e^6 + 4237/9127*e^5 - 30046/9127*e^4 - 134365/9127*e^3 + 219049/9127*e^2 + 723532/9127*e - 65838/9127, 490/9127*e^6 + 4237/9127*e^5 - 30046/9127*e^4 - 134365/9127*e^3 + 219049/9127*e^2 + 723532/9127*e - 65838/9127, -5381/9127*e^6 + 24084/9127*e^5 + 145775/9127*e^4 - 495233/9127*e^3 - 972260/9127*e^2 + 2585617/9127*e + 589382/9127, 5372/9127*e^6 - 21759/9127*e^5 - 152711/9127*e^4 + 442287/9127*e^3 + 1013555/9127*e^2 - 2360226/9127*e - 454546/9127, -5381/9127*e^6 + 24084/9127*e^5 + 145775/9127*e^4 - 495233/9127*e^3 - 972260/9127*e^2 + 2585617/9127*e + 589382/9127, 5372/9127*e^6 - 21759/9127*e^5 - 152711/9127*e^4 + 442287/9127*e^3 + 1013555/9127*e^2 - 2360226/9127*e - 454546/9127, 2730/9127*e^6 - 11598/9127*e^5 - 77433/9127*e^4 + 237111/9127*e^3 + 571095/9127*e^2 - 1263857/9127*e - 546744/9127, 2730/9127*e^6 - 11598/9127*e^5 - 77433/9127*e^4 + 237111/9127*e^3 + 571095/9127*e^2 - 1263857/9127*e - 546744/9127, -66/9127*e^6 - 1204/9127*e^5 + 3898/9127*e^4 + 46783/9127*e^3 - 25742/9127*e^2 - 394623/9127*e + 136944/9127, -66/9127*e^6 - 1204/9127*e^5 + 3898/9127*e^4 + 46783/9127*e^3 - 25742/9127*e^2 - 394623/9127*e + 136944/9127, 1491/9127*e^6 - 10968/9127*e^5 - 28319/9127*e^4 + 252854/9127*e^3 + 177458/9127*e^2 - 1345930/9127*e - 153136/9127, 1491/9127*e^6 - 10968/9127*e^5 - 28319/9127*e^4 + 252854/9127*e^3 + 177458/9127*e^2 - 1345930/9127*e - 153136/9127, -690/9127*e^6 + 4837/9127*e^5 + 15860/9127*e^4 - 122414/9127*e^3 - 119770/9127*e^2 + 744895/9127*e + 182118/9127, -690/9127*e^6 + 4837/9127*e^5 + 15860/9127*e^4 - 122414/9127*e^3 - 119770/9127*e^2 + 744895/9127*e + 182118/9127, -870/9127*e^6 + 5702/9127*e^5 + 14045/9127*e^4 - 104348/9127*e^3 - 33157/9127*e^2 + 351516/9127*e - 78310/9127, -870/9127*e^6 + 5702/9127*e^5 + 14045/9127*e^4 - 104348/9127*e^3 - 33157/9127*e^2 + 351516/9127*e - 78310/9127, 79/9127*e^6 + 888/9127*e^5 + 36/9127*e^4 - 57519/9127*e^3 - 21737/9127*e^2 + 422983/9127*e - 175534/9127, 79/9127*e^6 + 888/9127*e^5 + 36/9127*e^4 - 57519/9127*e^3 - 21737/9127*e^2 + 422983/9127*e - 175534/9127, -2760/9127*e^6 + 10221/9127*e^5 + 81694/9127*e^4 - 206719/9127*e^3 - 533842/9127*e^2 + 1190688/9127*e + 126090/9127, -2760/9127*e^6 + 10221/9127*e^5 + 81694/9127*e^4 - 206719/9127*e^3 - 533842/9127*e^2 + 1190688/9127*e + 126090/9127, 5238/9127*e^6 - 20608/9127*e^5 - 152541/9127*e^4 + 412535/9127*e^3 + 1074687/9127*e^2 - 2203925/9127*e - 657750/9127, 5238/9127*e^6 - 20608/9127*e^5 - 152541/9127*e^4 + 412535/9127*e^3 + 1074687/9127*e^2 - 2203925/9127*e - 657750/9127, 907/9127*e^6 - 3091/9127*e^5 - 18996/9127*e^4 + 30965/9127*e^3 + 12463/9127*e^2 - 87557/9127*e + 343386/9127, 907/9127*e^6 - 3091/9127*e^5 - 18996/9127*e^4 + 30965/9127*e^3 + 12463/9127*e^2 - 87557/9127*e + 343386/9127, -2194/9127*e^6 + 6994/9127*e^5 + 67626/9127*e^4 - 136357/9127*e^3 - 423162/9127*e^2 + 821193/9127*e - 119014/9127, 2854/9127*e^6 - 13208/9127*e^5 - 70098/9127*e^4 + 252655/9127*e^3 + 388518/9127*e^2 - 1237669/9127*e - 9154/9127, 2854/9127*e^6 - 13208/9127*e^5 - 70098/9127*e^4 + 252655/9127*e^3 + 388518/9127*e^2 - 1237669/9127*e - 9154/9127, -2194/9127*e^6 + 6994/9127*e^5 + 67626/9127*e^4 - 136357/9127*e^3 - 423162/9127*e^2 + 821193/9127*e - 119014/9127, -5028/9127*e^6 + 21120/9127*e^5 + 140968/9127*e^4 - 433612/9127*e^3 - 933870/9127*e^2 + 2278012/9127*e + 426132/9127, -5028/9127*e^6 + 21120/9127*e^5 + 140968/9127*e^4 - 433612/9127*e^3 - 933870/9127*e^2 + 2278012/9127*e + 426132/9127, -387/9127*e^6 - 422/9127*e^5 + 21197/9127*e^4 - 4055/9127*e^3 - 186620/9127*e^2 + 208860/9127*e + 413018/9127, -291/9127*e^6 + 2159/9127*e^5 - 5216/9127*e^4 + 913/9127*e^3 + 148695/9127*e^2 - 236048/9127*e - 179464/9127, -291/9127*e^6 + 2159/9127*e^5 - 5216/9127*e^4 + 913/9127*e^3 + 148695/9127*e^2 - 236048/9127*e - 179464/9127, 2096/9127*e^6 - 6016/9127*e^5 - 67093/9127*e^4 + 108468/9127*e^3 + 446892/9127*e^2 - 580740/9127*e - 2898/9127, 2096/9127*e^6 - 6016/9127*e^5 - 67093/9127*e^4 + 108468/9127*e^3 + 446892/9127*e^2 - 580740/9127*e - 2898/9127, -3274/9127*e^6 + 12184/9127*e^5 + 84117/9127*e^4 - 183120/9127*e^3 - 451104/9127*e^2 + 697034/9127*e + 107310/9127, -3274/9127*e^6 + 12184/9127*e^5 + 84117/9127*e^4 - 183120/9127*e^3 - 451104/9127*e^2 + 697034/9127*e + 107310/9127, -486/9127*e^6 + 6899/9127*e^5 - 337/9127*e^4 - 166619/9127*e^3 + 12069/9127*e^2 + 862761/9127*e + 162084/9127, -486/9127*e^6 + 6899/9127*e^5 - 337/9127*e^4 - 166619/9127*e^3 + 12069/9127*e^2 + 862761/9127*e + 162084/9127, -562/9127*e^6 - 3891/9127*e^5 + 29320/9127*e^4 + 130639/9127*e^3 - 171626/9127*e^2 - 636280/9127*e - 334048/9127, -562/9127*e^6 - 3891/9127*e^5 + 29320/9127*e^4 + 130639/9127*e^3 - 171626/9127*e^2 - 636280/9127*e - 334048/9127, -4993/9127*e^6 + 18163/9127*e^5 + 146645/9127*e^4 - 347376/9127*e^3 - 969726/9127*e^2 + 1780769/9127*e + 378402/9127, -4993/9127*e^6 + 18163/9127*e^5 + 146645/9127*e^4 - 347376/9127*e^3 - 969726/9127*e^2 + 1780769/9127*e + 378402/9127, 2468/9127*e^6 - 10846/9127*e^5 - 63342/9127*e^4 + 200735/9127*e^3 + 382892/9127*e^2 - 931145/9127*e + 204314/9127, 3424/9127*e^6 - 14426/9127*e^5 - 96295/9127*e^4 + 304970/9127*e^3 + 620792/9127*e^2 - 1691112/9127*e - 84996/9127, -320/9127*e^6 + 3566/9127*e^5 + 2858/9127*e^4 - 89576/9127*e^3 + 10989/9127*e^2 + 491226/9127*e - 142524/9127, 3424/9127*e^6 - 14426/9127*e^5 - 96295/9127*e^4 + 304970/9127*e^3 + 620792/9127*e^2 - 1691112/9127*e - 84996/9127, -320/9127*e^6 + 3566/9127*e^5 + 2858/9127*e^4 - 89576/9127*e^3 + 10989/9127*e^2 + 491226/9127*e - 142524/9127, 5516/9127*e^6 - 22451/9127*e^5 - 151259/9127*e^4 + 431485/9127*e^3 + 973471/9127*e^2 - 2215285/9127*e - 549220/9127, 5516/9127*e^6 - 22451/9127*e^5 - 151259/9127*e^4 + 431485/9127*e^3 + 973471/9127*e^2 - 2215285/9127*e - 549220/9127];
heckeEigenvalues := AssociativeArray();
for i := 1 to #heckeEigenvaluesArray do
  heckeEigenvalues[primes[i]] := heckeEigenvaluesArray[i];
end for;

ALEigenvalues := AssociativeArray();
ALEigenvalues[ideal<ZF | {5, 5, 1/5*w^3 + 4/5*w^2 - 11/5*w - 7}>] := -1;
ALEigenvalues[ideal<ZF | {5, 5, 1/5*w^3 - 6/5*w^2 - 1/5*w + 4}>] := -1;

// EXAMPLE:
// pp := Factorization(2*ZF)[1][1];
// heckeEigenvalues[pp];

print "To reconstruct the Hilbert newform f, type
  f, iso := Explode(make_newform());";

function make_newform();
 M := HilbertCuspForms(F, NN);
 S := NewSubspace(M);
 // SetVerbose("ModFrmHil", 1);
 NFD := NewformDecomposition(S);
 newforms := [* Eigenform(U) : U in NFD *];

 if #newforms eq 0 then;
  print "No Hilbert newforms at this level";
  return 0;
 end if;

 print "Testing ", #newforms, " possible newforms";
 newforms := [* f: f in newforms | IsIsomorphic(BaseField(f), K) *];
 print #newforms, " newforms have the correct Hecke field";

 if #newforms eq 0 then;
  print "No Hilbert newform found with the correct Hecke field";
  return 0;
 end if;

 autos := Automorphisms(K);
 xnewforms := [* *];
 for f in newforms do;
  if K eq RationalField() then;
   Append(~xnewforms, [* f, autos[1] *]);
  else;
   flag, iso := IsIsomorphic(K,BaseField(f));
   for a in autos do;
    Append(~xnewforms, [* f, a*iso *]);
   end for;
  end if;
 end for;
 newforms := xnewforms;

 for P in primes do;
  xnewforms := [* *];
  for f_iso in newforms do;
   f, iso := Explode(f_iso);
   if HeckeEigenvalue(f,P) eq iso(heckeEigenvalues[P]) then;
    Append(~xnewforms, f_iso);
   end if;
  end for;
  newforms := xnewforms;
  if #newforms eq 0 then;
   print "No Hilbert newform found which matches the Hecke eigenvalues";
   return 0;
  else if #newforms eq 1 then;
   print "success: unique match";
   return newforms[1];
  end if;
  end if;
 end for;
 print #newforms, "Hilbert newforms found which match the Hecke eigenvalues";
 return newforms[1];

end function;