   realprecision = 38 significant digits
   echo = 1 (on)
? nfpol=x^5-5*x^3+5*x+25
x^5 - 5*x^3 + 5*x + 25
? qpol=y^3-y-1;un=Mod(1,qpol);w=Mod(y,qpol);p=un*(x^5-5*x+w)
Mod(1, y^3 - y - 1)*x^5 + Mod(-5, y^3 - y - 1)*x + Mod(y, y^3 - y - 1)
? p2=x^5+3021*x^4-786303*x^3-6826636057*x^2-546603588746*x+3853890514072057
x^5 + 3021*x^4 - 786303*x^3 - 6826636057*x^2 - 546603588746*x + 385389051407
2057
? fa=[11699,6;2392997,2;4987333019653,2]

[11699 6]

[2392997 2]

[4987333019653 2]

? setrand(1);a=matrix(3,5,j,k,vectorv(5,l,random\10^8));
? setrand(1);as=matrix(3,3,j,k,vectorv(5,l,random\10^8));
? nf=nfinit(nfpol)
[x^5 - 5*x^3 + 5*x + 25, [1, 2], 595125, 45, [[1, -1.08911514572050482502495
27946671612684, -2.4285174907194186068992069565359418365, 0.7194669112891317
8943997506477288225734, -2.5558200350691694950646071159426779971; 1, -0.1383
8372073406036365047976417441696636 - 0.4918163765776864349975328551474152510
7*I, 1.9647119211288133163138753392090569931 + 0.809714924188978951282940822
19556466857*I, -0.072312766896812300380582649294307897099 + 2.19808037538462
76641195195160383234877*I, -0.98796319352507039803950539735452837192 + 1.570
1452385894131769052374806001981109*I; 1, 1.682941293594312776162956161507997
6006 + 2.0500351226010726172974286983598602164*I, -0.75045317576910401286427
186094108607489 + 1.3101462685358123283560773619310445916*I, -0.787420688747
75359433940488309213323154 + 2.1336633893126618034168454610457936018*I, 1.26
58732110596551455718089553258673705 - 2.716479010374315056657802803578983483
4*I], [1, -1.0891151457205048250249527946671612684, -2.428517490719418606899
2069565359418365, 0.71946691128913178943997506477288225734, -2.5558200350691
694950646071159426779971; 1.4142135623730950488016887242096980786, -0.195704
13467375904264179382543977540672, 2.7785222450164664309920925654093065576, -
0.10226569567819614506098907018896260032, -1.3971909474085893198147151262541
540506; 0, -0.69553338995335755797766403996841143190, 1.14510982744395651299
26149974389115723, 3.1085550780550843138423672171643499922, 2.22052069130868
72788181483285734827869; 1.4142135623730950488016887242096980786, 2.38003840
20787979181834702019470475019, -1.0613010590986270398182318786558994413, -1.
1135810173202366904448352912286604471, 1.79021506332534372536778891648110361
60; 0, 2.8991874737236275652408825679737171587, 1.85282662165584876344468105
12816401036, 3.0174557027049114270734649132936867273, -3.8416814583731999185
306312841432940661], 0, [5, 2, 0, -1, -2; 2, -2, -5, -10, 20; 0, -5, 10, -10
, 5; -1, -10, -10, -17, 1; -2, 20, 5, 1, -8], [345, 0, 200, 110, 177; 0, 345
, 95, 1, 145; 0, 0, 5, 4, 3; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1], [108675, 5175, 0
, -10350, -15525; 5175, 13800, -8625, -1725, 27600; 0, -8625, 37950, -17250,
 0; -10350, -1725, -17250, -24150, -15525; -15525, 27600, 0, -15525, -3450],
 [595125, [238050, -296700, 91425, 1725, 0]~]], [-2.428517490719418606899206
9565359418365, 1.9647119211288133163138753392090569931 + 0.80971492418897895
128294082219556466857*I, -0.75045317576910401286427186094108607489 + 1.31014
62685358123283560773619310445916*I], [1, 1/15*x^4 - 2/3*x^2 + 1/3*x + 4/3, x
, 2/15*x^4 - 1/3*x^2 + 2/3*x - 1/3, -1/15*x^4 + 1/3*x^3 + 1/3*x^2 - 4/3*x - 
2/3], [1, 0, 3, 1, 10; 0, 0, -2, 1, -5; 0, 1, 0, 3, -5; 0, 0, 1, 1, 10; 0, 0
, 0, 3, 0], [1, 0, 0, 0, 0, 0, -1, -1, -2, 4, 0, -1, 3, -1, 1, 0, -2, -1, -3
, -1, 0, 4, 1, -1, -1; 0, 1, 0, 0, 0, 1, 1, -1, -1, 1, 0, -1, -2, -1, 1, 0, 
-1, -1, -1, 3, 0, 1, 1, 3, -3; 0, 0, 1, 0, 0, 0, 0, 0, 1, -1, 1, 0, 0, 0, -2
, 0, 1, 0, -1, -1, 0, -1, -2, -1, -1; 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1,
 1, 2, 1, 0, 1, 0, 0, 0, 0, 2, 0, -1; 0, 0, 0, 0, 1, 0, -1, -1, -1, 1, 0, -1
, 0, 1, 0, 0, -1, 1, 0, 0, 1, 1, 0, 0, -1]]
? nf1=nfinit(nfpol,2)
[x^5 - 2*x^4 + 3*x^3 + 8*x^2 + 3*x + 2, [1, 2], 595125, 4, [[1, -1.089115145
7205048250249527946671612684, 2.4285174907194186068992069565359418364, -0.71
946691128913178943997506477288225734, 2.555820035069169495064607115942677997
1; 1, -0.13838372073406036365047976417441696637 + 0.491816376577686434997532
85514741525107*I, -1.9647119211288133163138753392090569931 + 0.8097149241889
7895128294082219556466857*I, 0.072312766896812300380582649294307897125 + 2.1
980803753846276641195195160383234877*I, 0.9879631935250703980395053973545283
7195 + 1.5701452385894131769052374806001981109*I; 1, 1.682941293594312776162
9561615079976006 + 2.0500351226010726172974286983598602164*I, 0.750453175769
10401286427186094108607491 - 1.3101462685358123283560773619310445914*I, 0.78
742068874775359433940488309213323161 - 2.13366338931266180341684546104579360
17*I, -1.2658732110596551455718089553258673704 + 2.7164790103743150566578028
035789834835*I], [1, -1.0891151457205048250249527946671612684, 2.42851749071
94186068992069565359418364, -0.71946691128913178943997506477288225734, 2.555
8200350691694950646071159426779971; 1.4142135623730950488016887242096980786,
 -0.19570413467375904264179382543977540674, -2.77852224501646643099209256540
93065576, 0.10226569567819614506098907018896260036, 1.3971909474085893198147
151262541540506; 0, 0.69553338995335755797766403996841143191, 1.145109827443
9565129926149974389115723, 3.1085550780550843138423672171643499922, 2.220520
6913086872788181483285734827869; 1.4142135623730950488016887242096980786, 2.
3800384020787979181834702019470475019, 1.06130105909862703981823187865589944
13, 1.1135810173202366904448352912286604472, -1.7902150633253437253677889164
811036159; 0, 2.8991874737236275652408825679737171587, -1.852826621655848763
4446810512816401035, -3.0174557027049114270734649132936867271, 3.84168145837
31999185306312841432940662], 0, [5, 2, 0, 1, 2; 2, -2, 5, 10, -20; 0, 5, 10,
 -10, 5; 1, 10, -10, -17, 1; 2, -20, 5, 1, -8], [345, 0, 145, 235, 168; 0, 3
45, 250, 344, 200; 0, 0, 5, 4, 3; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1], [108675, 51
75, 0, 10350, 15525; 5175, 13800, 8625, 1725, -27600; 0, 8625, 37950, -17250
, 0; 10350, 1725, -17250, -24150, -15525; 15525, -27600, 0, -15525, -3450], 
[595125, [-238050, 296700, 91425, 1725, 0]~]], [-1.0891151457205048250249527
946671612684, -0.13838372073406036365047976417441696637 + 0.4918163765776864
3499753285514741525107*I, 1.6829412935943127761629561615079976006 + 2.050035
1226010726172974286983598602164*I], [1, x, 1/2*x^4 - 3/2*x^3 + 5/2*x^2 + 2*x
 - 1, 1/2*x^4 - x^3 + x^2 + 9/2*x + 1, 1/2*x^4 - x^3 + 2*x^2 + 7/2*x + 2], [
1, 0, -1, -7, -14; 0, 1, 1, -2, -15; 0, 0, 0, -2, -4; 0, 0, -1, -1, 2; 0, 0,
 1, 3, 4], [1, 0, 0, 0, 0, 0, -1, 1, 2, -4, 0, 1, 3, -1, 1, 0, 2, -1, -3, -1
, 0, -4, 1, -1, -1; 0, 1, 0, 0, 0, 1, 1, 1, 1, -1, 0, 1, -2, -1, 1, 0, 1, -1
, -1, 3, 0, -1, 1, 3, -3; 0, 0, 1, 0, 0, 0, 0, 0, 1, -1, 1, 0, 0, 0, 2, 0, 1
, 0, 1, 1, 0, -1, 2, 1, 1; 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, -1, -1, -2, 
1, 0, -1, 0, 0, 0, 0, -2, 0, 1; 0, 0, 0, 0, 1, 0, 1, -1, -1, 1, 0, -1, 0, -1
, 0, 0, -1, -1, 0, 0, 1, 1, 0, 0, 1]]
? nfinit(nfpol,3)
[[x^5 - 2*x^4 + 3*x^3 + 8*x^2 + 3*x + 2, [1, 2], 595125, 4, [[1, -1.08911514
57205048250249527946671612684, 2.4285174907194186068992069565359418364, -0.7
1946691128913178943997506477288225734, 2.55582003506916949506460711594267799
71; 1, -0.13838372073406036365047976417441696637 + 0.49181637657768643499753
285514741525107*I, -1.9647119211288133163138753392090569931 + 0.809714924188
97895128294082219556466857*I, 0.072312766896812300380582649294307897125 + 2.
1980803753846276641195195160383234877*I, 0.987963193525070398039505397354528
37195 + 1.5701452385894131769052374806001981109*I; 1, 1.68294129359431277616
29561615079976006 + 2.0500351226010726172974286983598602164*I, 0.75045317576
910401286427186094108607491 - 1.3101462685358123283560773619310445914*I, 0.7
8742068874775359433940488309213323161 - 2.1336633893126618034168454610457936
017*I, -1.2658732110596551455718089553258673704 + 2.716479010374315056657802
8035789834835*I], [1, -1.0891151457205048250249527946671612684, 2.4285174907
194186068992069565359418364, -0.71946691128913178943997506477288225734, 2.55
58200350691694950646071159426779971; 1.4142135623730950488016887242096980786
, -0.19570413467375904264179382543977540674, -2.7785222450164664309920925654
093065576, 0.10226569567819614506098907018896260036, 1.397190947408589319814
7151262541540506; 0, 0.69553338995335755797766403996841143191, 1.14510982744
39565129926149974389115723, 3.1085550780550843138423672171643499922, 2.22052
06913086872788181483285734827869; 1.4142135623730950488016887242096980786, 2
.3800384020787979181834702019470475019, 1.0613010590986270398182318786558994
413, 1.1135810173202366904448352912286604472, -1.790215063325343725367788916
4811036159; 0, 2.8991874737236275652408825679737171587, -1.85282662165584876
34446810512816401035, -3.0174557027049114270734649132936867271, 3.8416814583
731999185306312841432940662], 0, [5, 2, 0, 1, 2; 2, -2, 5, 10, -20; 0, 5, 10
, -10, 5; 1, 10, -10, -17, 1; 2, -20, 5, 1, -8], [345, 0, 145, 235, 168; 0, 
345, 250, 344, 200; 0, 0, 5, 4, 3; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1], [108675, 5
175, 0, 10350, 15525; 5175, 13800, 8625, 1725, -27600; 0, 8625, 37950, -1725
0, 0; 10350, 1725, -17250, -24150, -15525; 15525, -27600, 0, -15525, -3450],
 [595125, [-238050, 296700, 91425, 1725, 0]~]], [-1.089115145720504825024952
7946671612684, -0.13838372073406036365047976417441696637 + 0.491816376577686
43499753285514741525107*I, 1.6829412935943127761629561615079976006 + 2.05003
51226010726172974286983598602164*I], [1, x, 1/2*x^4 - 3/2*x^3 + 5/2*x^2 + 2*
x - 1, 1/2*x^4 - x^3 + x^2 + 9/2*x + 1, 1/2*x^4 - x^3 + 2*x^2 + 7/2*x + 2], 
[1, 0, -1, -7, -14; 0, 1, 1, -2, -15; 0, 0, 0, -2, -4; 0, 0, -1, -1, 2; 0, 0
, 1, 3, 4], [1, 0, 0, 0, 0, 0, -1, 1, 2, -4, 0, 1, 3, -1, 1, 0, 2, -1, -3, -
1, 0, -4, 1, -1, -1; 0, 1, 0, 0, 0, 1, 1, 1, 1, -1, 0, 1, -2, -1, 1, 0, 1, -
1, -1, 3, 0, -1, 1, 3, -3; 0, 0, 1, 0, 0, 0, 0, 0, 1, -1, 1, 0, 0, 0, 2, 0, 
1, 0, 1, 1, 0, -1, 2, 1, 1; 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, -1, -1, -2,
 1, 0, -1, 0, 0, 0, 0, -2, 0, 1; 0, 0, 0, 0, 1, 0, 1, -1, -1, 1, 0, -1, 0, -
1, 0, 0, -1, -1, 0, 0, 1, 1, 0, 0, 1]], Mod(-1/2*x^4 + 3/2*x^3 - 5/2*x^2 - 2
*x + 1, x^5 - 2*x^4 + 3*x^3 + 8*x^2 + 3*x + 2)]
? nfinit(nfpol,4)
[x^5 - 2*x^4 + 3*x^3 + 8*x^2 + 3*x + 2, [1, 2], 595125, 4, [[1, -1.089115145
7205048250249527946671612684, 2.4285174907194186068992069565359418364, -0.71
946691128913178943997506477288225734, 2.555820035069169495064607115942677997
1; 1, -0.13838372073406036365047976417441696637 + 0.491816376577686434997532
85514741525107*I, -1.9647119211288133163138753392090569931 + 0.8097149241889
7895128294082219556466857*I, 0.072312766896812300380582649294307897125 + 2.1
980803753846276641195195160383234877*I, 0.9879631935250703980395053973545283
7195 + 1.5701452385894131769052374806001981109*I; 1, 1.682941293594312776162
9561615079976006 + 2.0500351226010726172974286983598602164*I, 0.750453175769
10401286427186094108607491 - 1.3101462685358123283560773619310445914*I, 0.78
742068874775359433940488309213323161 - 2.13366338931266180341684546104579360
17*I, -1.2658732110596551455718089553258673704 + 2.7164790103743150566578028
035789834835*I], [1, -1.0891151457205048250249527946671612684, 2.42851749071
94186068992069565359418364, -0.71946691128913178943997506477288225734, 2.555
8200350691694950646071159426779971; 1.4142135623730950488016887242096980786,
 -0.19570413467375904264179382543977540674, -2.77852224501646643099209256540
93065576, 0.10226569567819614506098907018896260036, 1.3971909474085893198147
151262541540506; 0, 0.69553338995335755797766403996841143191, 1.145109827443
9565129926149974389115723, 3.1085550780550843138423672171643499922, 2.220520
6913086872788181483285734827869; 1.4142135623730950488016887242096980786, 2.
3800384020787979181834702019470475019, 1.06130105909862703981823187865589944
13, 1.1135810173202366904448352912286604472, -1.7902150633253437253677889164
811036159; 0, 2.8991874737236275652408825679737171587, -1.852826621655848763
4446810512816401035, -3.0174557027049114270734649132936867271, 3.84168145837
31999185306312841432940662], 0, [5, 2, 0, 1, 2; 2, -2, 5, 10, -20; 0, 5, 10,
 -10, 5; 1, 10, -10, -17, 1; 2, -20, 5, 1, -8], [345, 0, 145, 235, 168; 0, 3
45, 250, 344, 200; 0, 0, 5, 4, 3; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1], [108675, 51
75, 0, 10350, 15525; 5175, 13800, 8625, 1725, -27600; 0, 8625, 37950, -17250
, 0; 10350, 1725, -17250, -24150, -15525; 15525, -27600, 0, -15525, -3450], 
[595125, [-238050, 296700, 91425, 1725, 0]~]], [-1.0891151457205048250249527
946671612684, -0.13838372073406036365047976417441696637 + 0.4918163765776864
3499753285514741525107*I, 1.6829412935943127761629561615079976006 + 2.050035
1226010726172974286983598602164*I], [1, x, 1/2*x^4 - 3/2*x^3 + 5/2*x^2 + 2*x
 - 1, 1/2*x^4 - x^3 + x^2 + 9/2*x + 1, 1/2*x^4 - x^3 + 2*x^2 + 7/2*x + 2], [
1, 0, -1, -7, -14; 0, 1, 1, -2, -15; 0, 0, 0, -2, -4; 0, 0, -1, -1, 2; 0, 0,
 1, 3, 4], [1, 0, 0, 0, 0, 0, -1, 1, 2, -4, 0, 1, 3, -1, 1, 0, 2, -1, -3, -1
, 0, -4, 1, -1, -1; 0, 1, 0, 0, 0, 1, 1, 1, 1, -1, 0, 1, -2, -1, 1, 0, 1, -1
, -1, 3, 0, -1, 1, 3, -3; 0, 0, 1, 0, 0, 0, 0, 0, 1, -1, 1, 0, 0, 0, 2, 0, 1
, 0, 1, 1, 0, -1, 2, 1, 1; 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 0, 0, -1, -1, -2, 
1, 0, -1, 0, 0, 0, 0, -2, 0, 1; 0, 0, 0, 0, 1, 0, 1, -1, -1, 1, 0, -1, 0, -1
, 0, 0, -1, -1, 0, 0, 1, 1, 0, 0, 1]]
? nf3=nfinit(x^6+108);
? nf4=nfinit(x^3-10*x+8)
[x^3 - 10*x + 8, [3, 0], 568, 2, [[1, -0.36332823793268357037416860931988791
957, -3.1413361156553641347759399165844441383; 1, -1.76155718183189058754537
11274124874988, 2.6261980685272936133764995500786243868; 1, 3.12488541976457
41579195397367323754183, -0.48486195287192947860055963349418024846], [1, -0.
36332823793268357037416860931988791957, -3.141336115655364134775939916584444
1383; 1, -1.7615571818318905875453711274124874988, 2.62619806852729361337649
95500786243868; 1, 3.1248854197645741579195397367323754183, -0.4848619528719
2947860055963349418024846], 0, [3, 1, -1; 1, 13, -5; -1, -5, 17], [284, 76, 
46; 0, 2, 0; 0, 0, 1], [196, -12, 8; -12, 50, 14; 8, 14, 38], [568, [120, 21
0, 2]~]], [-3.5046643535880477051501085259043320579, 0.864640886695403025831
12842266613688800, 2.6400234668926446793189801032381951699], [1, 1/2*x^2 + x
 - 3, -1/2*x^2 + 3], [1, 0, 6; 0, 1, 0; 0, 1, -2], [1, 0, 0, 0, 4, -2, 0, -2
, 6; 0, 1, 0, 1, 2, 0, 0, 0, -2; 0, 0, 1, 0, 1, -1, 1, -1, -1]]
? setrand(1);bnf2=bnfinit(qpol);nf2=bnf2[7];
? setrand(1);bnf=bnfinit(x^2-x-57,,[0.2,0.2])
[Mat(3), Mat([1, 2, 1, 2, 1, 2, 1, 2, 1]), [-2.71246530518434397468087951060
61300699 + 3.1415926535897932384626433832795028842*I; 2.71246530518434397468
08795106061300699], [1.7903417566977293763292119206302198761, 1.289761953065
2735025030086072395031017 + 3.1415926535897932384626433832795028842*I, 0.701
48550268542821846861610071436900868, 1.368455532 E-48, 0.5005798036324558738
2620331339071677436 + 3.1415926535897932384626433832795028842*I, 1.088856254
0123011578605958199158508674, 1.7241634548149836441438434283070556827 + 3.14
15926535897932384626433832795028842*I, -0.3432876442770270943898878667334192
1876 + 3.1415926535897932384626433832795028842*I, 2.133629400974756470719099
7873636390948 + 3.1415926535897932384626433832795028842*I, 0.066178301882745
732185368492323164193434 + 3.1415926535897932384626433832795028842*I; -1.790
3417566977293763292119206302198761, -1.2897619530652735025030086072395031017
, -0.70148550268542821846861610071436900868, 0.E-48, -0.50057980363245587382
620331339071677436, -1.0888562540123011578605958199158508674, -1.72416345481
49836441438434283070556827, 0.34328764427702709438988786673341921876, -2.133
6294009747564707190997873636390948, -0.0661783018827457321853684923231641934
34], [[3, [0, 1]~, 1, 1, [1, 1]~], [5, [-1, 1]~, 1, 1, [2, 1]~], [11, [-1, 1
]~, 1, 1, [2, 1]~], [3, [1, 1]~, 1, 1, [0, 1]~], [5, [2, 1]~, 1, 1, [-1, 1]~
], [11, [2, 1]~, 1, 1, [-1, 1]~], [19, [1, 1]~, 1, 1, [0, 1]~], [17, [3, 1]~
, 1, 1, [-2, 1]~], [17, [-2, 1]~, 1, 1, [3, 1]~], [19, [0, 1]~, 1, 1, [1, 1]
~]], 0, [x^2 - x - 57, [2, 0], 229, 1, [[1, -8.06637297521077796359593102467
05326058; 1, 7.0663729752107779635959310246705326058], [1, -8.06637297521077
79635959310246705326058; 1, 7.0663729752107779635959310246705326058], 0, [2,
 -1; -1, 115], [229, 115; 0, 1], [115, 1; 1, 2], [229, [115, 1]~]], [-7.0663
729752107779635959310246705326058, 8.0663729752107779635959310246705326058],
 [1, x - 1], [1, 1; 0, 1], [1, 0, 0, 57; 0, 1, 1, -1]], [[3, [3], [[3, 0; 0,
 1]]], 2.7124653051843439746808795106061300699, 0.8814422512654579364, [2, -
1], [x + 7], 156], [Mat(1), [[0, 0]], [[1.7903417566977293763292119206302198
761, -1.7903417566977293763292119206302198761]]], 0]
? setrand(1);bnfinit(x^2-x-100000,1)
[Mat(5), Mat([3, 2, 1, 2, 0, 3, 0, 2, 2, 3, 1, 4, 3, 2, 2, 3, 3, 0]), [-129.
82045011403975460991182396195022419 + 6.283185307179586476925286766559005768
4*I; 129.82045011403975460991182396195022419], [-41.811264589129943393339502
258694361489 + 6.2831853071795864769252867665590057684*I, 9.2399004147902289
816376260438840931575 + 3.1415926535897932384626433832795028842*I, -11.87460
9881075406725097315997431161032 + 3.1415926535897932384626433832795028842*I,
 0.E-67, -51.051165003920172374977128302578454647 + 3.1415926535897932384626
433832795028842*I, -64.910225057019877304955911980975112095 + 3.141592653589
7932384626433832795028842*I, -29.936654708054536668242186261263200457 + 3.14
15926535897932384626433832795028842*I, -47.668319071568233997332918482707687
879 + 6.2831853071795864769252867665590057684*I, 3.8762936464778825067484824
790355076166, -6.7377511782956880607802359510546381087 + 3.14159265358979323
84626433832795028842*I, -35.073513410834255332559266307639723380 + 3.1415926
535897932384626433832795028842*I, 33.130781426597481571750300827582717074 + 
2.967364921 E-67*I, 54.878404098312329644822020875673145628 + 5.93472984 E-6
7*I, -14.980188104648613073630759189293219180 + 3.14159265358979323846264338
32795028842*I, -26.831076484481330319708743069401142309 + 3.1415926535897932
384626433832795028842*I, -19.706749066516065512488907834878146945 + 3.141592
6535897932384626433832795028842*I, -22.104515522613877880850594423816214545 
+ 3.1415926535897932384626433832795028842*I, -45.687558235607825900087984737
729869106 + 6.2831853071795864769252867665590057684*I, 47.668319071568233997
332918482707687879 + 1.186945969 E-66*I; 41.81126458912994339333950225869436
1489, -9.2399004147902289816376260438840931575, 11.8746098810754067250973159
97431161032, 0.E-67, 51.051165003920172374977128302578454647, 64.91022505701
9877304955911980975112095, 29.936654708054536668242186261263200457, 47.66831
9071568233997332918482707687879, -3.8762936464778825067484824790355076166, 6
.7377511782956880607802359510546381087, 35.073513410834255332559266307639723
380, -33.130781426597481571750300827582717074, -54.8784040983123296448220208
75673145628, 14.980188104648613073630759189293219180, 26.8310764844813303197
08743069401142309, 19.706749066516065512488907834878146945, 22.1045155226138
77880850594423816214545, 45.687558235607825900087984737729869106, -47.668319
071568233997332918482707687879], [[2, [2, 1]~, 1, 1, [1, 1]~], [5, [5, 1]~, 
1, 1, [1, 1]~], [13, [-5, 1]~, 1, 1, [6, 1]~], [2, [3, 1]~, 1, 1, [0, 1]~], 
[5, [6, 1]~, 1, 1, [0, 1]~], [7, [4, 1]~, 2, 1, [-3, 1]~], [13, [6, 1]~, 1, 
1, [-5, 1]~], [23, [7, 1]~, 1, 1, [-6, 1]~], [43, [-15, 1]~, 1, 1, [16, 1]~]
, [17, [20, 1]~, 1, 1, [-2, 1]~], [17, [15, 1]~, 1, 1, [3, 1]~], [29, [14, 1
]~, 1, 1, [-13, 1]~], [29, [-13, 1]~, 1, 1, [14, 1]~], [31, [39, 1]~, 1, 1, 
[-7, 1]~], [31, [24, 1]~, 1, 1, [8, 1]~], [41, [7, 1]~, 1, 1, [-6, 1]~], [41
, [-6, 1]~, 1, 1, [7, 1]~], [43, [16, 1]~, 1, 1, [-15, 1]~], [23, [-6, 1]~, 
1, 1, [7, 1]~]], 0, [x^2 - x - 100000, [2, 0], 400001, 1, [[1, -316.72816130
129840161392089489603747004; 1, 315.72816130129840161392089489603747004], [1
, -316.72816130129840161392089489603747004; 1, 315.7281613012984016139208948
9603747004], 0, [2, -1; -1, 200001], [400001, 200001; 0, 1], [200001, 1; 1, 
2], [400001, [200001, 1]~]], [-315.72816130129840161392089489603747004, 316.
72816130129840161392089489603747004], [1, x - 1], [1, 1; 0, 1], [1, 0, 0, 10
0000; 0, 1, 1, -1]], [[5, [5], [[2, 0; 0, 1]]], 129.820450114039754609911823
96195022419, 0.9876536979069047227, [2, -1], [379554884019013781006303254896
369154068336082609238336*x + 11983616564425078999046283595002287166517812761
1316131167], 26], [Mat(1), [[0, 0]], [[-41.811264589129943393339502258694361
489 + 6.2831853071795864769252867665590057684*I, 41.811264589129943393339502
258694361489]]], 0]
? \p19
   realprecision = 19 significant digits
? setrand(1);sbnf=bnfinit(x^3-x^2-14*x-1,3)
[x^3 - x^2 - 14*x - 1, 3, 10889, [1, x, x^2 - x - 9], [-3.233732695981516673
, -0.07182350902743636344, 4.305556205008953036], [10889, 5698, 8994; 0, 1, 
0; 0, 0, 1], Mat(2), Mat([1, 1, 0, 1, 1, 0, 1, 1]), [9, 15, 16, 17, 39, 10, 
33, 57, 69], [2, [-1, 0, 0]~], [[0, 1, 0]~, [5, 3, 1]~], [[-4, -1, 2, 3, 10,
 3, 1, 7, 2; 1, 1, 1, 1, 5, 0, 1, 2, 1; 0, 0, 0, 0, 1, 0, 0, 0, -1], 0]]
? \p38
   realprecision = 38 significant digits
? bnrinit(bnf,[[5,4;0,1],[1,0]],1)
[[Mat(3), Mat([1, 2, 1, 2, 1, 2, 1, 2, 1]), [-2.7124653051843439746808795106
061300699 + 3.1415926535897932384626433832795028842*I; 2.7124653051843439746
808795106061300699], [1.7903417566977293763292119206302198761, 1.28976195306
52735025030086072395031017 + 3.1415926535897932384626433832795028842*I, 0.70
148550268542821846861610071436900868, 1.368455532 E-48, 0.500579803632455873
82620331339071677436 + 3.1415926535897932384626433832795028842*I, 1.08885625
40123011578605958199158508674, 1.7241634548149836441438434283070556827 + 3.1
415926535897932384626433832795028842*I, -0.343287644277027094389887866733419
21876 + 3.1415926535897932384626433832795028842*I, 2.13362940097475647071909
97873636390948 + 3.1415926535897932384626433832795028842*I, 0.06617830188274
5732185368492323164193434 + 3.1415926535897932384626433832795028842*I; -1.79
03417566977293763292119206302198761, -1.289761953065273502503008607239503101
7, -0.70148550268542821846861610071436900868, 0.E-48, -0.5005798036324558738
2620331339071677436, -1.0888562540123011578605958199158508674, -1.7241634548
149836441438434283070556827, 0.34328764427702709438988786673341921876, -2.13
36294009747564707190997873636390948, -0.066178301882745732185368492323164193
434], [[3, [0, 1]~, 1, 1, [1, 1]~], [5, [-1, 1]~, 1, 1, [2, 1]~], [11, [-1, 
1]~, 1, 1, [2, 1]~], [3, [1, 1]~, 1, 1, [0, 1]~], [5, [2, 1]~, 1, 1, [-1, 1]
~], [11, [2, 1]~, 1, 1, [-1, 1]~], [19, [1, 1]~, 1, 1, [0, 1]~], [17, [3, 1]
~, 1, 1, [-2, 1]~], [17, [-2, 1]~, 1, 1, [3, 1]~], [19, [0, 1]~, 1, 1, [1, 1
]~]], 0, [x^2 - x - 57, [2, 0], 229, 1, [[1, -8.0663729752107779635959310246
705326058; 1, 7.0663729752107779635959310246705326058], [1, -8.0663729752107
779635959310246705326058; 1, 7.0663729752107779635959310246705326058], 0, [2
, -1; -1, 115], [229, 115; 0, 1], [115, 1; 1, 2], [229, [115, 1]~]], [-7.066
3729752107779635959310246705326058, 8.0663729752107779635959310246705326058]
, [1, x - 1], [1, 1; 0, 1], [1, 0, 0, 57; 0, 1, 1, -1]], [[3, [3], [[3, 0; 0
, 1]]], 2.7124653051843439746808795106061300699, 0.8814422512654579364, [2, 
-1], [x + 7], 156], [Mat(1), [[0, 0]], [[1.790341756697729376329211920630219
8761, -1.7903417566977293763292119206302198761]]], [0, [Mat([[6, 1]~, 1])]]]
, [[[5, 4; 0, 1], [1, 0]], [8, [4, 2], [[2, 0]~, [0, 1]~]], Mat([[5, [-1, 1]
~, 1, 1, [2, 1]~], 1]), [[[[4], [[2, 0]~], [[2, 0]~], [[Mod(0, 2)]~], 1]], [
[2], [[0, 1]~], Mat(1)]], [1, 0; 0, 1]], [1], Mat([1, -3, -6]), [12, [12], [
[3, 0; 0, 1]]], [[1/2, 0; 0, 0], [1, -1; 1, 1]]]
? bnr=bnrclass(bnf,[[5,4;0,1],[1,0]],2)
[[Mat(3), Mat([1, 2, 1, 2, 1, 2, 1, 2, 1]), [-2.7124653051843439746808795106
061300699 + 3.1415926535897932384626433832795028842*I; 2.7124653051843439746
808795106061300699], [1.7903417566977293763292119206302198761, 1.28976195306
52735025030086072395031017 + 3.1415926535897932384626433832795028842*I, 0.70
148550268542821846861610071436900868, 1.368455532 E-48, 0.500579803632455873
82620331339071677436 + 3.1415926535897932384626433832795028842*I, 1.08885625
40123011578605958199158508674, 1.7241634548149836441438434283070556827 + 3.1
415926535897932384626433832795028842*I, -0.343287644277027094389887866733419
21876 + 3.1415926535897932384626433832795028842*I, 2.13362940097475647071909
97873636390948 + 3.1415926535897932384626433832795028842*I, 0.06617830188274
5732185368492323164193434 + 3.1415926535897932384626433832795028842*I; -1.79
03417566977293763292119206302198761, -1.289761953065273502503008607239503101
7, -0.70148550268542821846861610071436900868, 0.E-48, -0.5005798036324558738
2620331339071677436, -1.0888562540123011578605958199158508674, -1.7241634548
149836441438434283070556827, 0.34328764427702709438988786673341921876, -2.13
36294009747564707190997873636390948, -0.066178301882745732185368492323164193
434], [[3, [0, 1]~, 1, 1, [1, 1]~], [5, [-1, 1]~, 1, 1, [2, 1]~], [11, [-1, 
1]~, 1, 1, [2, 1]~], [3, [1, 1]~, 1, 1, [0, 1]~], [5, [2, 1]~, 1, 1, [-1, 1]
~], [11, [2, 1]~, 1, 1, [-1, 1]~], [19, [1, 1]~, 1, 1, [0, 1]~], [17, [3, 1]
~, 1, 1, [-2, 1]~], [17, [-2, 1]~, 1, 1, [3, 1]~], [19, [0, 1]~, 1, 1, [1, 1
]~]], 0, [x^2 - x - 57, [2, 0], 229, 1, [[1, -8.0663729752107779635959310246
705326058; 1, 7.0663729752107779635959310246705326058], [1, -8.0663729752107
779635959310246705326058; 1, 7.0663729752107779635959310246705326058], 0, [2
, -1; -1, 115], [229, 115; 0, 1], [115, 1; 1, 2], [229, [115, 1]~]], [-7.066
3729752107779635959310246705326058, 8.0663729752107779635959310246705326058]
, [1, x - 1], [1, 1; 0, 1], [1, 0, 0, 57; 0, 1, 1, -1]], [[3, [3], [[3, 0; 0
, 1]]], 2.7124653051843439746808795106061300699, 0.8814422512654579364, [2, 
-1], [x + 7], 156], [Mat(1), [[0, 0]], [[1.790341756697729376329211920630219
8761, -1.7903417566977293763292119206302198761]]], [0, [Mat([[6, 1]~, 1])]]]
, [[[5, 4; 0, 1], [1, 0]], [8, [4, 2], [[2, 0]~, [0, 1]~]], Mat([[5, [-1, 1]
~, 1, 1, [2, 1]~], 1]), [[[[4], [[2, 0]~], [[2, 0]~], [[Mod(0, 2)]~], 1]], [
[2], [[0, 1]~], Mat(1)]], [1, 0; 0, 1]], [1], Mat([1, -3, -6]), [12, [12], [
[3, 0; 0, 1]]], [[1/2, 0; 0, 0], [1, -1; 1, 1]]]
? rnfinit(nf2,x^5-x-2)
[x^5 - x - 2, [], [[49744, 0, 0; 0, 49744, 0; 0, 0, 49744], [3109, 0, 0]~], 
1, [], [], [[1, x, x^2, x^3, x^4], [[1, 0, 0; 0, 1, 0; 0, 0, 1], [1, 0, 0; 0
, 1, 0; 0, 0, 1], [1, 0, 0; 0, 1, 0; 0, 0, 1], [1, 0, 0; 0, 1, 0; 0, 0, 1], 
[1, 0, 0; 0, 1, 0; 0, 0, 1]]], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0;
 0, 0, 0, 1, 0; 0, 0, 0, 0, 1], [], [y^3 - y - 1, [1, 1], -23, 1, [[1, 0.754
87766624669276004950889635852869189, 1.3247179572447460259609088544780973407
; 1, -0.87743883312334638002475444817926434595 - 0.7448617666197442365931704
2860439236724*I, -0.66235897862237301298045442723904867037 + 0.5622795120623
0124389918214490937306150*I], [1, 0.75487766624669276004950889635852869189, 
1.3247179572447460259609088544780973407; 1.414213562373095048801688724209698
0786, -1.2408858979558593537192653626096055786, -0.9367170507273508470331116
4961686101697; 0, -1.0533936124468254335289038498031013276, 0.79518331180303
271004476715629658754003], 0, [3, -1, 0; -1, 1, 3; 0, 3, 2], [23, 16, 13; 0,
 1, 0; 0, 0, 1], [-7, 2, -3; 2, 6, -9; -3, -9, 2], [23, [10, 7, 1]~]], [1.32
47179572447460259609088544780973407, -0.662358978622373012980454427239048670
37 + 0.56227951206230124389918214490937306150*I], [1, y^2 - 1, y], [1, 0, 1;
 0, 0, 1; 0, 1, 0], [1, 0, 0, 0, 0, 1, 0, 1, 1; 0, 1, 0, 1, -1, 0, 0, 0, 1; 
0, 0, 1, 0, 1, 0, 1, 0, 0]], [x^15 - 5*x^13 + 5*x^12 + 7*x^11 - 26*x^10 - 5*
x^9 + 45*x^8 + 158*x^7 - 98*x^6 + 110*x^5 - 190*x^4 + 189*x^3 + 144*x^2 + 25
*x + 1, Mod(39516536165538345/83718587879473471*x^14 - 6500512476832995/8371
8587879473471*x^13 - 196215472046117185/83718587879473471*x^12 + 22990222748
0108910/83718587879473471*x^11 + 237380704030959181/83718587879473471*x^10 -
 1064931988160773805/83718587879473471*x^9 - 20657086671714300/8371858787947
3471*x^8 + 1772885205999206010/83718587879473471*x^7 + 5952033217241102348/8
3718587879473471*x^6 - 4838840187320655696/83718587879473471*x^5 + 518039072
0553188700/83718587879473471*x^4 - 8374015687535120430/83718587879473471*x^3
 + 8907744727915040221/83718587879473471*x^2 + 4155976664123434381/837185878
79473471*x + 318920215718580450/83718587879473471, x^15 - 5*x^13 + 5*x^12 + 
7*x^11 - 26*x^10 - 5*x^9 + 45*x^8 + 158*x^7 - 98*x^6 + 110*x^5 - 190*x^4 + 1
89*x^3 + 144*x^2 + 25*x + 1), -1], 0]
? bnfcertify(bnf)
1
? setrand(1);bnfclassunit(x^4-7,2,[0.2,0.2])

[x^4 - 7]

[[2, 1]]

[[-87808, 1]]

[[1, x, x^2, x^3]]

[[2, [2], [[3, 1, 2, 1; 0, 1, 0, 0; 0, 0, 1, 0; 0, 0, 0, 1]]]]

[14.229975145405511722395637833443108790]

[1.1211171071527562293]

? setrand(1);bnfclassunit(x^2-x-100000)
  ***   Warning: insufficient precision for fundamental units, not given.

[x^2 - x - 100000]

[[2, 0]]

[[400001, 1]]

[[1, x - 1]]

[[5, [5], [[2, 0; 0, 1]]]]

[129.82045011403975460991182396195022419]

[0.9876536979069047227]

[[2, -1]]

[[;]]

[-27]

? setrand(1);bnfclassunit(x^2-x-100000,1)

[x^2 - x - 100000]

[[2, 0]]

[[400001, 1]]

[[1, x - 1]]

[[5, [5], [[2, 0; 0, 1]]]]

[129.82045011403975460991182396195022419]

[0.9876536979069047227]

[[2, -1]]

[[379554884019013781006303254896369154068336082609238336*x + 119836165644250
789990462835950022871665178127611316131167]]

[26]

? setrand(1);bnfclassunit(x^4+24*x^2+585*x+1791,,[0.1,0.1])

[x^4 + 24*x^2 + 585*x + 1791]

[[0, 2]]

[[18981, 3087]]

[[1, -10/1029*x^3 + 13/343*x^2 - 165/343*x - 1135/343, 17/1029*x^3 - 32/1029
*x^2 + 109/343*x + 2444/343, -11/343*x^3 + 163/1029*x^2 - 373/343*x - 4260/3
43]]

[[4, [4], [[7, 2, 4, 0; 0, 1, 0, 0; 0, 0, 1, 0; 0, 0, 0, 1]]]]

[3.7941269688216589341408274220859400302]

[0.8826018286655581300]

[[6, 10/1029*x^3 - 13/343*x^2 + 165/343*x + 1478/343]]

[[4/1029*x^3 + 53/1029*x^2 + 66/343*x + 111/343]]

[151]

? setrand(1);bnfclgp(17)
[1, [], []]
? setrand(1);bnfclgp(-31)
[3, [3], [Qfb(2, 1, 4)]]
? setrand(1);bnfclgp(x^4+24*x^2+585*x+1791)
[4, [4], [[7, 4, 5, 6; 0, 1, 0, 0; 0, 0, 1, 0; 0, 0, 0, 1]]]
? bnrconductor(bnf,[[25,14;0,1],[1,1]])
[[5, 4; 0, 1], [1, 0]]
? bnrconductorofchar(bnr,[2])
[[5, 4; 0, 1], [0, 0]]
? bnfisprincipal(bnf,[5,2;0,1],0)
[1]~
? bnfisprincipal(bnf,[5,2;0,1])
[[1]~, [7/3, 1/3]~, 157]
? bnfisunit(bnf,Mod(3405*x-27466,x^2-x-57))
[-4, Mod(1, 2)]~
? \p19
   realprecision = 19 significant digits
? bnfmake(sbnf)
[Mat(2), Mat([1, 1, 0, 1, 1, 0, 1, 1]), [1.173637103435061715 + 3.1415926535
89793239*I, -4.562279014988837911 + 3.141592653589793239*I; -2.6335434327389
76050 + 3.141592653589793239*I, 1.420330600779487358 + 3.141592653589793239*
I; 1.459906329303914335, 3.141948414209350544], [1.246346989334819161 + 3.14
1592653589793238*I, 0.5404006376129469728 + 3.141592653589793239*I, -0.69263
91142471042844 + 3.141592653589793239*I, -1.990056445584799713 + 3.141592653
589793239*I, -0.8305625946607188641 + 3.141592653589793239*I, 0.E-57, 0.0043
75616572659815433 + 3.141592653589793239*I, -1.977791147836553954, 0.3677262
014027817707 + 3.141592653589793239*I; 0.6716827432867392934 + 3.14159265358
9793238*I, -0.8333219883742404172 + 3.141592653589793238*I, -0.2461086674077
943078, 0.5379005671092853266, -1.552661549868775854, 0.E-57, -0.87383180430
71131265, 0.5774919091398324093, 0.9729063188316092380; -1.91802973262155845
4, 0.2929213507612934446, 0.9387477816548985923, 1.452155878475514387, 2.383
224144529494719, 0.E-57, 0.8694561877344533111, 1.400299238696721544, -1.340
632520234391008], [[3, [-1, 1, 0]~, 1, 1, [1, 1, 1]~], [5, [-1, 1, 0]~, 1, 1
, [0, 1, 1]~], [5, [2, 1, 0]~, 1, 1, [1, -2, 1]~], [5, [3, 1, 0]~, 1, 1, [2,
 2, 1]~], [13, [19, 1, 0]~, 1, 1, [-2, -6, 1]~], [3, [10, 1, 1]~, 1, 2, [-1,
 1, 0]~], [11, [1, 1, 0]~, 1, 1, [-3, -1, 1]~], [19, [-6, 1, 0]~, 1, 1, [6, 
6, 1]~], [23, [-10, 1, 0]~, 1, 1, [-7, 10, 1]~]]~, 0, [x^3 - x^2 - 14*x - 1,
 [3, 0], 10889, 1, [[1, -3.233732695981516673, 4.690759845041404812; 1, -0.0
7182350902743636345, -8.923017874523549404; 1, 4.305556205008953037, 5.23225
8029482144592], [1, -3.233732695981516673, 4.690759845041404812; 1, -0.07182
350902743636345, -8.923017874523549404; 1, 4.305556205008953037, 5.232258029
482144592], 0, [3, 1, 1; 1, 29, 8; 1, 8, 129], [10889, 5698, 8994; 0, 1, 0; 
0, 0, 1], [3677, -121, -21; -121, 386, -23; -21, -23, 86], [10889, [1899, 51
91, 1]~]], [-3.233732695981516673, -0.07182350902743636345, 4.30555620500895
3037], [1, x, x^2 - x - 9], [1, 0, 9; 0, 1, 1; 0, 0, 1], [1, 0, 0, 0, 9, 1, 
0, 1, 44; 0, 1, 0, 1, 1, 5, 0, 5, 1; 0, 0, 1, 0, 1, 0, 1, 0, -4]], [[2, [2],
 [[3, 2, 0; 0, 1, 0; 0, 0, 1]]], 10.34800724602768001, 1.000000000000000000,
 [2, -1], [x, x^2 + 2*x - 4], 1000], [Mat(1), [[0.E-57, 0.E-57, 0.E-57]], [[
1.246346989334819161 + 3.141592653589793238*I, 0.6716827432867392934 + 3.141
592653589793238*I, -1.918029732621558454]]], [[-4, -1, 2, 3, 10, 3, 1, 7, 2;
 1, 1, 1, 1, 5, 0, 1, 2, 1; 0, 0, 0, 0, 1, 0, 0, 0, -1], 0]]
? \p38
   realprecision = 38 significant digits
? bnfnarrow(bnf)
[3, [3], [[3, 0; 0, 1]]]
? bnfreg(x^2-x-57)
2.7124653051843439746808795106061300699
? bnfsignunit(bnf)

[-1]

[1]

? bnfunit(bnf)
[[x + 7], 156]
? bnrclass(bnf,[[5,4;0,1],[1,0]])
[12, [12], [[3, 0; 0, 1]]]
? bnr2=bnrclass(bnf,[[25,14;0,1],[1,1]],2)
[[Mat(3), Mat([1, 2, 1, 2, 1, 2, 1, 2, 1]), [-2.7124653051843439746808795106
061300699 + 3.1415926535897932384626433832795028842*I; 2.7124653051843439746
808795106061300699], [1.7903417566977293763292119206302198761, 1.28976195306
52735025030086072395031017 + 3.1415926535897932384626433832795028842*I, 0.70
148550268542821846861610071436900868, 1.368455532 E-48, 0.500579803632455873
82620331339071677436 + 3.1415926535897932384626433832795028842*I, 1.08885625
40123011578605958199158508674, 1.7241634548149836441438434283070556827 + 3.1
415926535897932384626433832795028842*I, -0.343287644277027094389887866733419
21876 + 3.1415926535897932384626433832795028842*I, 2.13362940097475647071909
97873636390948 + 3.1415926535897932384626433832795028842*I, 0.06617830188274
5732185368492323164193434 + 3.1415926535897932384626433832795028842*I; -1.79
03417566977293763292119206302198761, -1.289761953065273502503008607239503101
7, -0.70148550268542821846861610071436900868, 0.E-48, -0.5005798036324558738
2620331339071677436, -1.0888562540123011578605958199158508674, -1.7241634548
149836441438434283070556827, 0.34328764427702709438988786673341921876, -2.13
36294009747564707190997873636390948, -0.066178301882745732185368492323164193
434], [[3, [0, 1]~, 1, 1, [1, 1]~], [5, [-1, 1]~, 1, 1, [2, 1]~], [11, [-1, 
1]~, 1, 1, [2, 1]~], [3, [1, 1]~, 1, 1, [0, 1]~], [5, [2, 1]~, 1, 1, [-1, 1]
~], [11, [2, 1]~, 1, 1, [-1, 1]~], [19, [1, 1]~, 1, 1, [0, 1]~], [17, [3, 1]
~, 1, 1, [-2, 1]~], [17, [-2, 1]~, 1, 1, [3, 1]~], [19, [0, 1]~, 1, 1, [1, 1
]~]], 0, [x^2 - x - 57, [2, 0], 229, 1, [[1, -8.0663729752107779635959310246
705326058; 1, 7.0663729752107779635959310246705326058], [1, -8.0663729752107
779635959310246705326058; 1, 7.0663729752107779635959310246705326058], 0, [2
, -1; -1, 115], [229, 115; 0, 1], [115, 1; 1, 2], [229, [115, 1]~]], [-7.066
3729752107779635959310246705326058, 8.0663729752107779635959310246705326058]
, [1, x - 1], [1, 1; 0, 1], [1, 0, 0, 57; 0, 1, 1, -1]], [[3, [3], [[3, 0; 0
, 1]]], 2.7124653051843439746808795106061300699, 0.8814422512654579364, [2, 
-1], [x + 7], 156], [Mat(1), [[0, 0]], [[1.790341756697729376329211920630219
8761, -1.7903417566977293763292119206302198761]]], [0, [Mat([[6, 1]~, 1])]]]
, [[[25, 14; 0, 1], [1, 1]], [80, [20, 2, 2], [[2, 0]~, [4, 2]~, [-2, -2]~]]
, Mat([[5, [-1, 1]~, 1, 1, [2, 1]~], 2]), [[[[4], [[2, 0]~], [[2, 0]~], [[Mo
d(0, 2), Mod(0, 2)]~], 1], [[5], [[6, 0]~], [[6, 0]~], [[Mod(0, 2), Mod(0, 2
)]~], Mat([1/5, -14/5])]], [[2, 2], [[4, 2]~, [-2, -2]~], [1, 0; 0, 1]]], [1
, -12, 0, 0; 0, 0, 1, 0; 0, 0, 0, 1]], [1], Mat([1, -3, -6, 0]), [12, [12], 
[[3, 0; 0, 1]]], [[1, -18, 9; -1/2, 10, -5], [-2, 0; 0, -10]]]
? bnrclassno(bnf,[[5,4;0,1],[1,0]])
12
? lu=ideallist(bnf,55,3);
? bnrclassnolist(bnf,lu)
[[3], [], [3, 3], [3], [6, 6], [], [], [], [3, 3, 3], [], [3, 3], [3, 3], []
, [], [12, 6, 6, 12], [3], [3, 3], [], [9, 9], [6, 6], [], [], [], [], [6, 1
2, 6], [], [3, 3, 3, 3], [], [], [], [], [], [3, 6, 6, 3], [], [], [9, 3, 9]
, [6, 6], [], [], [], [], [], [3, 3], [3, 3], [12, 12, 6, 6, 12, 12], [], []
, [6, 6], [9], [], [3, 3, 3, 3], [], [3, 3], [], [6, 12, 12, 6]]
? bnrdisc(bnr,Mat(6))
[12, 12, 18026977100265125]
? bnrdisc(bnr)
[24, 12, 40621487921685401825918161408203125]
? bnrdisc(bnr2,,,2)
0
? bnrdisc(bnr,Mat(6),,1)
[6, 2, [125, 14; 0, 1]]
? bnrdisc(bnr,,,1)
[12, 1, [1953125, 1160889; 0, 1]]
? bnrdisc(bnr2,,,3)
0
? bnrdisclist(bnf,lu)
[[[6, 6, Mat([229, 3])]], [], [[], []], [[]], [[12, 12, [5, 3; 229, 6]], [12
, 12, [5, 3; 229, 6]]], [], [], [], [[], [], []], [], [[], []], [[], []], []
, [], [[24, 24, [3, 6; 5, 9; 229, 12]], [], [], [24, 24, [3, 6; 5, 9; 229, 1
2]]], [[]], [[], []], [], [[18, 18, [19, 6; 229, 9]], [18, 18, [19, 6; 229, 
9]]], [[], []], [], [], [], [], [[], [24, 24, [5, 12; 229, 12]], []], [], [[
], [], [], []], [], [], [], [], [], [[], [12, 12, [3, 3; 11, 3; 229, 6]], [1
2, 12, [3, 3; 11, 3; 229, 6]], []], [], [], [[18, 18, [2, 12; 3, 12; 229, 9]
], [], [18, 18, [2, 12; 3, 12; 229, 9]]], [[12, 12, [37, 3; 229, 6]], [12, 1
2, [37, 3; 229, 6]]], [], [], [], [], [], [[], []], [[], []], [[], [], [], [
], [], []], [], [], [[12, 12, [2, 12; 3, 3; 229, 6]], [12, 12, [2, 12; 3, 3;
 229, 6]]], [[18, 18, [7, 12; 229, 9]]], [], [[], [], [], []], [], [[], []],
 [], [[], [24, 24, [5, 9; 11, 6; 229, 12]], [24, 24, [5, 9; 11, 6; 229, 12]]
, []]]
? bnrdisclist(bnf,20,,1)
[[[[matrix(0,2), [[6, 6, Mat([229, 3])], [0, 0, 0], [0, 0, 0], [0, 0, 0]]]],
 [], [[Mat([12, 1]), [[0, 0, 0], [0, 0, 0], [0, 0, 0], [12, 0, [3, 3; 229, 6
]]]], [Mat([13, 1]), [[0, 0, 0], [0, 0, 0], [12, 6, [-1, 1; 3, 3; 229, 6]], 
[0, 0, 0]]]], [[Mat([10, 1]), [[0, 0, 0], [0, 0, 0], [0, 0, 0], [0, 0, 0]]]]
, [[Mat([20, 1]), [[12, 12, [5, 3; 229, 6]], [0, 0, 0], [0, 0, 0], [24, 0, [
5, 9; 229, 12]]]], [Mat([21, 1]), [[12, 12, [5, 3; 229, 6]], [0, 0, 0], [24,
 12, [5, 9; 229, 12]], [0, 0, 0]]]], [], [], [], [[Mat([12, 2]), [[0, 0, 0],
 [0, 0, 0], [0, 0, 0], [0, 0, 0]]], [[12, 1; 13, 1], [[0, 0, 0], [12, 6, [-1
, 1; 3, 6; 229, 6]], [0, 0, 0], [24, 0, [3, 12; 229, 12]]]], [Mat([13, 2]), 
[[0, 0, 0], [0, 0, 0], [0, 0, 0], [0, 0, 0]]]], [], [[Mat([44, 1]), [[0, 0, 
0], [0, 0, 0], [12, 6, [-1, 1; 11, 3; 229, 6]], [0, 0, 0]]], [Mat([45, 1]), 
[[0, 0, 0], [0, 0, 0], [0, 0, 0], [12, 0, [11, 3; 229, 6]]]]], [[[10, 1; 12,
 1], [[0, 0, 0], [0, 0, 0], [0, 0, 0], [0, 0, 0]]], [[10, 1; 13, 1], [[0, 0,
 0], [0, 0, 0], [0, 0, 0], [0, 0, 0]]]], [], [], [[[12, 1; 20, 1], [[24, 24,
 [3, 6; 5, 9; 229, 12]], [0, 0, 0], [0, 0, 0], [48, 0, [3, 12; 5, 18; 229, 2
4]]]], [[13, 1; 20, 1], [[0, 0, 0], [24, 12, [3, 6; 5, 9; 229, 12]], [24, 12
, [3, 6; 5, 6; 229, 12]], [48, 0, [3, 12; 5, 18; 229, 24]]]], [[12, 1; 21, 1
], [[0, 0, 0], [24, 12, [3, 6; 5, 9; 229, 12]], [0, 0, 0], [48, 0, [3, 12; 5
, 18; 229, 24]]]], [[13, 1; 21, 1], [[24, 24, [3, 6; 5, 9; 229, 12]], [0, 0,
 0], [48, 24, [3, 12; 5, 18; 229, 24]], [0, 0, 0]]]], [[Mat([10, 2]), [[0, 0
, 0], [12, 6, [-1, 1; 2, 12; 229, 6]], [12, 6, [-1, 1; 2, 12; 229, 6]], [24,
 0, [2, 36; 229, 12]]]]], [[Mat([68, 1]), [[0, 0, 0], [12, 6, [-1, 1; 17, 3;
 229, 6]], [0, 0, 0], [0, 0, 0]]], [Mat([69, 1]), [[0, 0, 0], [12, 6, [-1, 1
; 17, 3; 229, 6]], [0, 0, 0], [0, 0, 0]]]], [], [[Mat([76, 1]), [[18, 18, [1
9, 6; 229, 9]], [0, 0, 0], [0, 0, 0], [36, 0, [19, 15; 229, 18]]]], [Mat([77
, 1]), [[18, 18, [19, 6; 229, 9]], [0, 0, 0], [36, 18, [-1, 1; 19, 15; 229, 
18]], [0, 0, 0]]]], [[[10, 1; 20, 1], [[0, 0, 0], [0, 0, 0], [0, 0, 0], [0, 
0, 0]]], [[10, 1; 21, 1], [[0, 0, 0], [0, 0, 0], [0, 0, 0], [0, 0, 0]]]]]]
? bnrisprincipal(bnr,idealprimedec(bnf,7)[1])
[[9]~, [112595/19683, 13958/19683]~, 192]
? dirzetak(nf4,30)
[1, 2, 0, 3, 1, 0, 0, 4, 0, 2, 1, 0, 0, 0, 0, 5, 1, 0, 0, 3, 0, 2, 0, 0, 2, 
0, 1, 0, 1, 0]
? factornf(x^3+x^2-2*x-1,t^3+t^2-2*t-1)

[Mod(1, t^3 + t^2 - 2*t - 1)*x + Mod(-t, t^3 + t^2 - 2*t - 1) 1]

[Mod(1, t^3 + t^2 - 2*t - 1)*x + Mod(-t^2 + 2, t^3 + t^2 - 2*t - 1) 1]

[Mod(1, t^3 + t^2 - 2*t - 1)*x + Mod(t^2 + t - 1, t^3 + t^2 - 2*t - 1) 1]

? vp=idealprimedec(nf,3)[1]
[3, [1, 0, 1, 0, 0]~, 1, 1, [1, -1, -1, -1, 0]~]
? idx=idealmul(nf,matid(5),vp)

[3 2 1 0 1]

[0 1 0 0 0]

[0 0 1 0 0]

[0 0 0 1 0]

[0 0 0 0 1]

? idealinv(nf,idx)

[1 0 0 2/3 0]

[0 1 0 1/3 0]

[0 0 1 1/3 0]

[0 0 0 1/3 0]

[0 0 0 0 1]

? idy=idealred(nf,idx,[1,5,6])

[5 0 0 0 2]

[0 5 0 0 2]

[0 0 5 0 1]

[0 0 0 5 2]

[0 0 0 0 1]

? idx2=idealmul(nf,idx,idx)

[9 5 7 0 4]

[0 1 0 0 0]

[0 0 1 0 0]

[0 0 0 1 0]

[0 0 0 0 1]

? idt=idealmul(nf,idx,idx,1)

[2 0 0 0 0]

[0 2 0 0 0]

[0 0 2 0 0]

[0 0 0 2 1]

[0 0 0 0 1]

? idz=idealintersect(nf,idx,idy)

[15 10 5 0 12]

[0 5 0 0 2]

[0 0 5 0 1]

[0 0 0 5 2]

[0 0 0 0 1]

? aid=[idx,idy,idz,matid(5),idx]
[[3, 2, 1, 0, 1; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1]
, [5, 0, 0, 0, 2; 0, 5, 0, 0, 2; 0, 0, 5, 0, 1; 0, 0, 0, 5, 2; 0, 0, 0, 0, 1
], [15, 10, 5, 0, 12; 0, 5, 0, 0, 2; 0, 0, 5, 0, 1; 0, 0, 0, 5, 2; 0, 0, 0, 
0, 1], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0,
 0, 1], [3, 2, 1, 0, 1; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0
, 0, 1]]
? bid=idealstar(nf2,54,1)
[[[54, 0, 0; 0, 54, 0; 0, 0, 54], [0]], [132678, [1638, 9, 9]], [[2, [2, 0, 
0]~, 1, 3, [1, 0, 0]~], 1; [3, [3, 0, 0]~, 1, 3, [1, 0, 0]~], 3], [[[[7], [[
1, 0, 1]~], [[1, 0, -27]~], [[]~], 1]], [[[26], [[2, 2, 2]~], [[-25, 2, 2]~]
, [[]~], 1], [[3, 3, 3], [[1, 3, 0]~, [1, 0, 3]~, [4, 0, 0]~], [[1, -24, 0]~
, [1, 0, -24]~, [-23, 0, 0]~], [[]~, []~, []~], [0, 1/3, 0; 0, 0, 1/3; 1/3, 
0, 0]], [[3, 3, 3], [[1, 9, 0]~, [1, 0, 9]~, [10, 0, 0]~], [[1, -18, 0]~, [1
, 0, -18]~, [-17, 0, 0]~], [[]~, []~, []~], [0, 1/9, 0; 0, 0, 1/9; 1/9, 0, 0
]]], [[], [], [;]]], [468, 469, 0, 0, -48776, 0, 0, -36582; 0, 0, 1, 0, -3, 
-6, 0, 0; 0, 0, 0, 1, -1, 0, -6, -3]]
? vaid=[idx,idy,matid(5)]
[[3, 2, 1, 0, 1; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1]
, [5, 0, 0, 0, 2; 0, 5, 0, 0, 2; 0, 0, 5, 0, 1; 0, 0, 0, 5, 2; 0, 0, 0, 0, 1
], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 
1]]
? haid=[matid(5),matid(5),matid(5)]
[[1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1]
, [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1
], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 
1]]
? idealadd(nf,idx,idy)

[1 0 0 0 0]

[0 1 0 0 0]

[0 0 1 0 0]

[0 0 0 1 0]

[0 0 0 0 1]

? idealaddtoone(nf,idx,idy)
[[3, 2, 1, -3, 1]~, [-2, -2, -1, 3, -1]~]
? idealaddtoone(nf,[idy,idx])
[[-5, 0, 0, 0, 0]~, [6, 0, 0, 0, 0]~]
? idealappr(nf,idy)
[-1, -1, 2, 4, -3]~
? idealappr(nf,idealfactor(nf,idy),1)
[-1, -1, 2, 4, -3]~
? idealcoprime(nf,idx,idx)
[-1/3, 1/3, 1/3, 1/3, 0]~
? idealdiv(nf,idy,idt)

[5 0 5/2 0 1]

[0 5/2 0 0 1]

[0 0 5/2 0 1/2]

[0 0 0 5/2 1]

[0 0 0 0 1/2]

? idealdiv(nf,idx2,idx,1)

[3 2 1 0 1]

[0 1 0 0 0]

[0 0 1 0 0]

[0 0 0 1 0]

[0 0 0 0 1]

? idf=idealfactor(nf,idz)

[[3, [1, 0, 1, 0, 0]~, 1, 1, [1, -1, -1, -1, 0]~] 1]

[[5, [-1, 0, 0, 0, 2]~, 4, 1, [2, 2, 1, 2, 1]~] 3]

[[5, [2, 0, 0, 0, -2]~, 1, 1, [2, 0, 3, 0, 1]~] 1]

? idealhnf(nf,vp)

[3 2 1 0 1]

[0 1 0 0 0]

[0 0 1 0 0]

[0 0 0 1 0]

[0 0 0 0 1]

? idealhnf(nf,vp[2],3)

[3 2 1 0 1]

[0 1 0 0 0]

[0 0 1 0 0]

[0 0 0 1 0]

[0 0 0 0 1]

? ideallist(bnf,20)
[[[1, 0; 0, 1]], [], [[3, 0; 0, 1], [3, 1; 0, 1]], [[2, 0; 0, 2]], [[5, 4; 0
, 1], [5, 2; 0, 1]], [], [], [], [[9, 6; 0, 1], [3, 0; 0, 3], [9, 4; 0, 1]],
 [], [[11, 10; 0, 1], [11, 2; 0, 1]], [[6, 0; 0, 2], [6, 2; 0, 2]], [], [], 
[[15, 9; 0, 1], [15, 4; 0, 1], [15, 12; 0, 1], [15, 7; 0, 1]], [[4, 0; 0, 4]
], [[17, 15; 0, 1], [17, 3; 0, 1]], [], [[19, 0; 0, 1], [19, 1; 0, 1]], [[10
, 8; 0, 2], [10, 4; 0, 2]]]
? ideallog(nf2,w,bid)
[1388, 6, 8]~
? idealmin(nf,idx,[1,2,3])
[[-2; 1; 1; 0; 1], [2.0885812311199768913287869744681966009 + 3.141592653589
7932384626433832795028842*I, 1.5921096812520196555597562531657929785 + 4.244
7196639216499665715751642189271112*I, -0.79031915447583185468082063233076160
203 + 2.5437460822678889883600220330800078854*I]]
? idealnorm(nf,idt)
16
? idp=idealpow(nf,idx,7)

[2187 1436 1807 630 1822]

[0 1 0 0 0]

[0 0 1 0 0]

[0 0 0 1 0]

[0 0 0 0 1]

? idealpow(nf,idx,7,1)

[2 0 0 0 0]

[0 2 0 0 0]

[0 0 2 0 0]

[0 0 0 2 1]

[0 0 0 0 1]

? idealprimedec(nf,2)
[[2, [3, 0, 1, 0, 0]~, 1, 1, [0, 0, 0, 1, 1]~], [2, [12, -4, -2, 11, 3]~, 1,
 4, [1, 0, 1, 0, 0]~]]
? idealprimedec(nf,3)
[[3, [1, 0, 1, 0, 0]~, 1, 1, [1, -1, -1, -1, 0]~], [3, [1, 1, 1, 0, 0]~, 2, 
2, [0, 2, 2, 1, 0]~]]
? idealprimedec(nf,11)
[[11, [11, 0, 0, 0, 0]~, 1, 5, [1, 0, 0, 0, 0]~]]
? idealprincipal(nf,Mod(x^3+5,nfpol))

[6]

[1]

[3]

[1]

[3]

? idealtwoelt(nf,idy)
[5, [2, 2, 1, 2, 1]~]
? idealtwoelt(nf,idy,10)
[-1, -1, 2, -1, 2]~
? idealstar(nf2,54)
[[[54, 0, 0; 0, 54, 0; 0, 0, 54], [0]], [132678, [1638, 9, 9]], [[2, [2, 0, 
0]~, 1, 3, [1, 0, 0]~], 1; [3, [3, 0, 0]~, 1, 3, [1, 0, 0]~], 3], [[[[7], [[
1, 1, 1]~], [[1, -27, -27]~], [[]~], 1]], [[[26], [[3, 2, 0]~], [[3, 2, 0]~]
, [[]~], 1], [[3, 3, 3], [[1, 3, 0]~, [1, 0, 3]~, [4, 0, 0]~], [[1, -24, 0]~
, [1, 0, -24]~, [-23, 0, 0]~], [[]~, []~, []~], [0, 1/3, 0; 0, 0, 1/3; 1/3, 
0, 0]], [[3, 3, 3], [[1, 9, 0]~, [1, 0, 9]~, [10, 0, 0]~], [[1, -18, 0]~, [1
, 0, -18]~, [-17, 0, 0]~], [[]~, []~, []~], [0, 1/9, 0; 0, 0, 1/9; 1/9, 0, 0
]]], [[], [], [;]]], [2106, -77, 0, -11102, 2184, 0, 6006, -13104; 0, 0, 1, 
-3, 0, -6, 0, 0; -27, 1, 0, 142, -28, 0, -78, 168]]
? idealval(nf,idp,vp)
7
? ideleprincipal(nf,Mod(x^3+5,nfpol))
[[6; 1; 3; 1; 3], [2.2324480827796254080981385584384939684 + 3.1415926535897
932384626433832795028842*I, 5.0387659675158716386435353106610489968 + 1.5851
760343512250049897278861965702424*I, 4.2664040272651028743625910797589683173
 - 0.0083630478144368246110910258645462996115*I]]
? ba=nfalgtobasis(nf,Mod(x^3+5,nfpol))
[6, 1, 3, 1, 3]~
? bb=nfalgtobasis(nf,Mod(x^3+x,nfpol))
[1, 1, 4, 1, 3]~
? bc=matalgtobasis(nf,[Mod(x^2+x,nfpol);Mod(x^2+1,nfpol)])

[[3, -2, 1, 1, 0]~]

[[4, -2, 0, 1, 0]~]

? matbasistoalg(nf,bc)

[Mod(x^2 + x, x^5 - 5*x^3 + 5*x + 25)]

[Mod(x^2 + 1, x^5 - 5*x^3 + 5*x + 25)]

? nfbasis(x^3+4*x+5)
[1, x, 1/7*x^2 - 1/7*x - 2/7]
? nfbasis(x^3+4*x+5,2)
[1, x, 1/7*x^2 - 1/7*x - 2/7]
? nfbasis(x^3+4*x+12,1)
[1, x, 1/2*x^2]
? nfbasistoalg(nf,ba)
Mod(x^3 + 5, x^5 - 5*x^3 + 5*x + 25)
? nfbasis(p2,0,fa)
[1, x, x^2, 1/11699*x^3 + 1847/11699*x^2 - 132/11699*x - 2641/11699, 1/13962
3738889203638909659*x^4 - 1552451622081122020/139623738889203638909659*x^3 +
 418509858130821123141/139623738889203638909659*x^2 - 6810913798507599407313
4/139623738889203638909659*x - 13185339461968406/58346808996920447]
? da=nfdetint(nf,[a,aid])

[90 70 35 0 65]

[0 5 0 0 0]

[0 0 5 0 0]

[0 0 0 5 0]

[0 0 0 0 5]

? nfdisc(x^3+4*x+12)
-1036
? nfdisc(x^3+4*x+12,1)
-1036
? nfdisc(p2,0,fa)
136866601
? nfeltdiv(nf,ba,bb)
[584/373, 66/373, -32/373, -105/373, 120/373]~
? nfeltdiveuc(nf,ba,bb)
[2, 0, 0, 0, 0]~
? nfeltdivrem(nf,ba,bb)
[[2, 0, 0, 0, 0]~, [4, -1, -5, -1, -3]~]
? nfeltmod(nf,ba,bb)
[4, -1, -5, -1, -3]~
? nfeltmul(nf,ba,bb)
[50, -15, -35, 60, 15]~
? nfeltpow(nf,bb,5)
[-291920, 136855, 230560, -178520, 74190]~
? nfeltreduce(nf,ba,idx)
[1, 0, 0, 0, 0]~
? nfeltval(nf,ba,vp)
0
? nffactor(nf2,x^3+x)

[x 1]

[x^2 + 1 1]

? aut=nfgaloisconj(nf3)
[-x, x, -1/12*x^4 - 1/2*x, -1/12*x^4 + 1/2*x, 1/12*x^4 - 1/2*x, 1/12*x^4 + 1
/2*x]~
? nfgaloisapply(nf3,aut[5],Mod(x^5,x^6+108))
Mod(-1/2*x^5 + 9*x^2, x^6 + 108)
? nfhilbert(nf,3,5)
-1
? nfhilbert(nf,3,5,idf[1,1])
-1
? nfhnf(nf,[a,aid])
[[[1, 0, 0, 0, 0]~, [-2, 0, 0, 0, 0]~, [0, 0, 0, 0, 0]~; [0, 0, 0, 0, 0]~, [
1, 0, 0, 0, 0]~, [0, 0, 0, 0, 0]~; [0, 0, 0, 0, 0]~, [0, 0, 0, 0, 0]~, [1, 0
, 0, 0, 0]~], [[10, 4, 5, 6, 7; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0;
 0, 0, 0, 0, 1], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0
; 0, 0, 0, 0, 1], [3, 2, 1, 0, 1; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 
0; 0, 0, 0, 0, 1]]]
? nfhnfmod(nf,[a,aid],da)
[[[1, 0, 0, 0, 0]~, [-2, 0, 0, 0, 0]~, [0, 0, 0, 0, 0]~; [0, 0, 0, 0, 0]~, [
1, 0, 0, 0, 0]~, [0, 0, 0, 0, 0]~; [0, 0, 0, 0, 0]~, [0, 0, 0, 0, 0]~, [1, 0
, 0, 0, 0]~], [[10, 4, 5, 6, 7; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0;
 0, 0, 0, 0, 1], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 0
; 0, 0, 0, 0, 1], [3, 2, 1, 0, 1; 0, 1, 0, 0, 0; 0, 0, 1, 0, 0; 0, 0, 0, 1, 
0; 0, 0, 0, 0, 1]]]
? nfisideal(bnf[7],[5,2;0,1])
1
? nfisincl(x^2+1,x^4+1)
[-x^2, x^2]
? nfisincl(x^2+1,nfinit(x^4+1))
[-x^2, x^2]
? nfisisom(x^3+x^2-2*x-1,x^3+x^2-2*x-1)
[x, -x^2 - x + 1, x^2 - 2]
? nfisisom(x^3-2,nfinit(x^3-6*x^2-6*x-30))
[-1/25*x^2 + 13/25*x - 2/5]
? nfroots(nf2,x+2)
[Mod(-2, y^3 - y - 1)]
? nfrootsof1(nf)
[2, [-1, 0, 0, 0, 0]~]
? nfsnf(nf,[as,haid,vaid])
[[2562748315629757085585610, 436545976069778274371140, 123799938628701108220
1405, 2356446991473627724963350, 801407102592194537169612; 0, 5, 0, 0, 2; 0,
 0, 5, 0, 1; 0, 0, 0, 5, 2; 0, 0, 0, 0, 1], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 0
, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1], [1, 0, 0, 0, 0; 0, 1, 0, 0, 0; 
0, 0, 1, 0, 0; 0, 0, 0, 1, 0; 0, 0, 0, 0, 1]]
? nfsubfields(nf)
[[x^5 - 5*x^3 + 5*x + 25, x], [x, x^5 - 5*x^3 + 5*x + 25]]
? polcompositum(x^4-4*x+2,x^3-x-1)
[x^12 - 4*x^10 + 8*x^9 + 12*x^8 + 12*x^7 + 138*x^6 + 132*x^5 - 43*x^4 + 58*x
^2 - 128*x - 5]
? polcompositum(x^4-4*x+2,x^3-x-1,1)
[[x^12 - 4*x^10 + 8*x^9 + 12*x^8 + 12*x^7 + 138*x^6 + 132*x^5 - 43*x^4 + 58*
x^2 - 128*x - 5, Mod(-279140305176/29063006931199*x^11 + 129916611552/290630
06931199*x^10 + 1272919322296/29063006931199*x^9 - 2813750209005/29063006931
199*x^8 - 2859411937992/29063006931199*x^7 - 414533880536/29063006931199*x^6
 - 35713977492936/29063006931199*x^5 - 17432607267590/29063006931199*x^4 + 4
9785595543672/29063006931199*x^3 + 9423768373204/29063006931199*x^2 - 427797
76146743/29063006931199*x + 37962587857138/29063006931199, x^12 - 4*x^10 + 8
*x^9 + 12*x^8 + 12*x^7 + 138*x^6 + 132*x^5 - 43*x^4 + 58*x^2 - 128*x - 5), M
od(-279140305176/29063006931199*x^11 + 129916611552/29063006931199*x^10 + 12
72919322296/29063006931199*x^9 - 2813750209005/29063006931199*x^8 - 28594119
37992/29063006931199*x^7 - 414533880536/29063006931199*x^6 - 35713977492936/
29063006931199*x^5 - 17432607267590/29063006931199*x^4 + 49785595543672/2906
3006931199*x^3 + 9423768373204/29063006931199*x^2 - 13716769215544/290630069
31199*x + 37962587857138/29063006931199, x^12 - 4*x^10 + 8*x^9 + 12*x^8 + 12
*x^7 + 138*x^6 + 132*x^5 - 43*x^4 + 58*x^2 - 128*x - 5), -1]]
? polgalois(x^6-3*x^2-1)
[12, 1, 1]
? polred(x^5-2*x^4-4*x^3-96*x^2-352*x-568)
[x - 1, x^5 - x^4 - 6*x^3 + 6*x^2 + 13*x - 5, x^5 - x^4 + 2*x^3 - 4*x^2 + x 
- 1, x^5 - x^4 + 4*x^3 - 2*x^2 + x - 1, x^5 + 4*x^3 - 4*x^2 + 8*x - 8]
? polred(x^4-28*x^3-458*x^2+9156*x-25321,3)

[1 x - 1]

[1/115*x^2 - 14/115*x - 327/115 x^2 - 10]

[3/1495*x^3 - 63/1495*x^2 - 1607/1495*x + 13307/1495 x^4 - 32*x^2 + 216]

[1/4485*x^3 - 7/1495*x^2 - 1034/4485*x + 7924/4485 x^4 - 8*x^2 + 6]

? polred(x^4+576,1)
[x - 1, x^2 - x + 1, x^2 + 1, x^4 - x^2 + 1]
? polred(x^4+576,3)

[1 x - 1]

[1/192*x^3 + 1/8*x + 1/2 x^2 - x + 1]

[-1/24*x^2 x^2 + 1]

[-1/192*x^3 + 1/48*x^2 + 1/8*x x^4 - x^2 + 1]

? polred(p2,0,fa)
[x - 1, x^5 - 2*x^4 - 62*x^3 + 85*x^2 + 818*x + 1, x^5 - 2*x^4 - 53*x^3 - 46
*x^2 + 508*x + 913, x^5 - 2*x^4 - 13*x^3 + 37*x^2 - 21*x - 1, x^5 - x^4 - 52
*x^3 - 197*x^2 - 273*x - 127]
? polred(p2,1,fa)
[x - 1, x^5 - 2*x^4 - 62*x^3 + 85*x^2 + 818*x + 1, x^5 - 2*x^4 - 53*x^3 - 46
*x^2 + 508*x + 913, x^5 - 2*x^4 - 13*x^3 + 37*x^2 - 21*x - 1, x^5 - x^4 - 52
*x^3 - 197*x^2 - 273*x - 127]
? polredabs(x^5-2*x^4-4*x^3-96*x^2-352*x-568)
x^5 - x^4 + 2*x^3 - 4*x^2 + x - 1
? polredabs(x^5-2*x^4-4*x^3-96*x^2-352*x-568,1)
[x^5 - x^4 + 2*x^3 - 4*x^2 + x - 1, Mod(2*x^4 - x^3 + 3*x^2 - 3*x - 1, x^5 -
 x^4 + 2*x^3 - 4*x^2 + x - 1)]
? polredord(x^3-12*x+45*x-1)
[x - 1, x^3 - 363*x - 2663, x^3 + 33*x - 1]
? polsubcyclo(31,5)
x^5 + x^4 - 12*x^3 - 21*x^2 + x + 5
? setrand(1);poltschirnhaus(x^5-x-1)
x^5 - 15*x^4 + 88*x^3 - 278*x^2 + 452*x - 289
? aa=rnfpseudobasis(nf2,p)
[[[1, 0, 0]~, [0, 0, 0]~, [0, 0, 0]~, [-2, 0, 0]~, [3, 1, 0]~; [0, 0, 0]~, [
1, 0, 0]~, [0, 0, 0]~, [2, 0, 0]~, [0, -1, 0]~; [0, 0, 0]~, [0, 0, 0]~, [1, 
0, 0]~, [1, 0, 0]~, [-5, -2, 0]~; [0, 0, 0]~, [0, 0, 0]~, [0, 0, 0]~, [1, 0,
 0]~, [-2, 0, 0]~; [0, 0, 0]~, [0, 0, 0]~, [0, 0, 0]~, [0, 0, 0]~, [1, 0, 0]
~], [[1, 0, 0; 0, 1, 0; 0, 0, 1], [1, 0, 0; 0, 1, 0; 0, 0, 1], [1, 0, 0; 0, 
1, 0; 0, 0, 1], [1, 0, 2/5; 0, 1, 3/5; 0, 0, 1/5], [1, 0, 22/25; 0, 1, 8/25;
 0, 0, 1/25]], [416134375, 202396875, 60056800; 0, 3125, 2700; 0, 0, 25], [-
1275, 5, 5]~]
? rnfbasis(bnf2,aa)

[[1, 0, 0]~ [0, 0, 0]~ [0, 0, 0]~ [-26/25, 11/25, -8/25]~ [0, -4, 7]~]

[[0, 0, 0]~ [1, 0, 0]~ [0, 0, 0]~ [53/25, -8/25, -1/25]~ [-6/5, 41/5, -53/5]
~]

[[0, 0, 0]~ [0, 0, 0]~ [1, 0, 0]~ [-14/25, -21/25, 13/25]~ [16/5, -1/5, -7/5
]~]

[[0, 0, 0]~ [0, 0, 0]~ [0, 0, 0]~ [7/25, -2/25, 6/25]~ [-2/5, 2/5, -11/5]~]

[[0, 0, 0]~ [0, 0, 0]~ [0, 0, 0]~ [9/25, 1/25, -3/25]~ [-2/5, 7/5, -6/5]~]

? rnfdisc(nf2,p)
[[416134375, 202396875, 60056800; 0, 3125, 2700; 0, 0, 25], [-1275, 5, 5]~]
? rnfequation(nf2,p)
x^15 - 15*x^11 + 75*x^7 - x^5 - 125*x^3 + 5*x + 1
? rnfequation(nf2,p,1)
[x^15 - 15*x^11 + 75*x^7 - x^5 - 125*x^3 + 5*x + 1, Mod(-x^5 + 5*x, x^15 - 1
5*x^11 + 75*x^7 - x^5 - 125*x^3 + 5*x + 1), 0]
? rnfhnfbasis(bnf2,aa)

[[1, 0, 0]~ [0, 0, 0]~ [0, 0, 0]~ [6/5, 4/5, -2/5]~ [3/25, -8/25, 24/25]~]

[[0, 0, 0]~ [1, 0, 0]~ [0, 0, 0]~ [-6/5, -4/5, 2/5]~ [-9/25, -1/25, 3/25]~]

[[0, 0, 0]~ [0, 0, 0]~ [1, 0, 0]~ [-3/5, -2/5, 1/5]~ [-8/25, 13/25, -39/25]~
]

[[0, 0, 0]~ [0, 0, 0]~ [0, 0, 0]~ [-3/5, -2/5, 1/5]~ [4/25, 6/25, -18/25]~]

[[0, 0, 0]~ [0, 0, 0]~ [0, 0, 0]~ [0, 0, 0]~ [-2/25, -3/25, 9/25]~]

? rnfisfree(bnf2,aa)
1
? rnfsteinitz(nf2,aa)
[[[1, 0, 0]~, [0, 0, 0]~, [0, 0, 0]~, [-26/25, 11/25, -8/25]~, [29/125, -2/2
5, 8/125]~; [0, 0, 0]~, [1, 0, 0]~, [0, 0, 0]~, [53/25, -8/25, -1/25]~, [-53
/125, 7/125, 1/125]~; [0, 0, 0]~, [0, 0, 0]~, [1, 0, 0]~, [-14/25, -21/25, 1
3/25]~, [9/125, 19/125, -13/125]~; [0, 0, 0]~, [0, 0, 0]~, [0, 0, 0]~, [7/25
, -2/25, 6/25]~, [-9/125, 2/125, -6/125]~; [0, 0, 0]~, [0, 0, 0]~, [0, 0, 0]
~, [9/25, 1/25, -3/25]~, [-8/125, -1/125, 3/125]~], [[1, 0, 0; 0, 1, 0; 0, 0
, 1], [1, 0, 0; 0, 1, 0; 0, 0, 1], [1, 0, 0; 0, 1, 0; 0, 0, 1], [1, 0, 0; 0,
 1, 0; 0, 0, 1], [125, 0, 22; 0, 125, 108; 0, 0, 1]], [416134375, 202396875,
 60056800; 0, 3125, 2700; 0, 0, 25], [-1275, 5, 5]~]
? nfz=zetakinit(x^2-2);
? zetak(nfz,-3)
0.091666666666666666666666666666666666667
? zetak(nfz,1.5+3*I)
0.88324345992059326405525724366416928890 - 0.2067536250233895222724230899142
7938845*I
? setrand(1);quadclassunit(1-10^7,,[1,1])
  ***   Warning: not a fundamental discriminant in quadclassunit.
[2416, [1208, 2], [Qfb(277, 55, 9028), Qfb(1700, 1249, 1700)], 1, 1.00257481
62993077507]
? setrand(1);quadclassunit(10^9-3,,[0.5,0.5])
[4, [4], [Qfb(3, 1, -83333333, 0.E-48)], 2800.625251907016076486370621737074
5514, 0.9849577285369119736]
? sizebyte(%)
172
? getheap
[200, 122407]
? print("Total time spent: ",gettime);
Total time spent: 910
? \q
