Return to linear CVS log

Up to [local] / OpenXM_contrib / pari-2.2 / src / test / 32

Imported pari-2.2.1(alpha).

   realprecision = 38 significant digits
   echo = 1 (on)
? algdep(2*cos(2*Pi/13),6)
x^6 + x^5 - 5*x^4 - 4*x^3 + 6*x^2 + 3*x - 1
? algdep(2*cos(2*Pi/13),6,15)
x^6 + x^5 - 5*x^4 - 4*x^3 + 6*x^2 + 3*x - 1
? charpoly([1,2;3,4],z)
z^2 - 5*z - 2
? charpoly(Mod(x^2+x+1,x^3+5*x+1),z)
z^3 + 7*z^2 + 16*z - 19
? charpoly([1,2;3,4],z,1)
z^2 - 5*z - 2
? charpoly(Mod(1,8191)*[1,2;3,4],z,2)
z^2 + Mod(8186, 8191)*z + Mod(8189, 8191)
? lindep(Mod(1,7)*[2,-1;1,3],-1)
[Mod(6, 7), Mod(5, 7)]~
? lindep([(1-3*sqrt(2))/(3-2*sqrt(3)),1,sqrt(2),sqrt(3),sqrt(6)])
[-3, -3, 9, -2, 6]
? lindep([(1-3*sqrt(2))/(3-2*sqrt(3)),1,sqrt(2),sqrt(3),sqrt(6)],14)
[-3, -3, 9, -2, 6]
? matadjoint([1,2;3,4])

[4 -2]

[-3 1]

? matcompanion(x^5-12*x^3+0.0005)

[0 0 0 0 -0.00049999999999999999999999999999999999999]

[1 0 0 0 0]

[0 1 0 0 0]

[0 0 1 0 12]

[0 0 0 1 0]

? matdet([1,2,3;1,5,6;9,8,7])
-30
? matdet([1,2,3;1,5,6;9,8,7],1)
-30
? matdetint([1,2,3;4,5,6])
3
? matdiagonal([2,4,6])

[2 0 0]

[0 4 0]

[0 0 6]

? mateigen([1,2,3;4,5,6;7,8,9])

[-1.2833494518006402717978106547571267252 1 0.283349451800640271797810654757
12672521]

[-0.14167472590032013589890532737856336260 -2 0.6416747259003201358989053273
7856336260]

[1 1 1]

? mathess(mathilbert(7))

[1 90281/58800 -1919947/4344340 4858466341/1095033030 -77651417539/819678732
6 3386888964/106615355 1/2]

[1/3 43/48 38789/5585580 268214641/109503303 -581330123627/126464718744 4365
450643/274153770 1/4]

[0 217/2880 442223/7447440 53953931/292008808 -32242849453/168619624992 1475
457901/1827691800 1/80]

[0 0 1604444/264539275 24208141/149362505292 847880210129/47916076768560 -45
44407141/103873817300 -29/40920]

[0 0 0 9773092581/35395807550620 -24363634138919/107305824577186620 72118203
606917/60481351061158500 55899/3088554700]

[0 0 0 0 67201501179065/8543442888354179988 -9970556426629/74082861999267660
0 -3229/13661312210]

[0 0 0 0 0 -258198800769/9279048099409000 -13183/38381527800]

? mathilbert(5)

[1 1/2 1/3 1/4 1/5]

[1/2 1/3 1/4 1/5 1/6]

[1/3 1/4 1/5 1/6 1/7]

[1/4 1/5 1/6 1/7 1/8]

[1/5 1/6 1/7 1/8 1/9]

? amat=1/mathilbert(7)

[49 -1176 8820 -29400 48510 -38808 12012]

[-1176 37632 -317520 1128960 -1940400 1596672 -504504]

[8820 -317520 2857680 -10584000 18711000 -15717240 5045040]

[-29400 1128960 -10584000 40320000 -72765000 62092800 -20180160]

[48510 -1940400 18711000 -72765000 133402500 -115259760 37837800]

[-38808 1596672 -15717240 62092800 -115259760 100590336 -33297264]

[12012 -504504 5045040 -20180160 37837800 -33297264 11099088]

? mathnf(amat)

[420 0 0 0 210 168 175]

[0 840 0 0 0 0 504]

[0 0 2520 0 0 0 1260]

[0 0 0 2520 0 0 840]

[0 0 0 0 13860 0 6930]

[0 0 0 0 0 5544 0]

[0 0 0 0 0 0 12012]

? mathnf(amat,1)
[[420, 0, 0, 0, 210, 168, 175; 0, 840, 0, 0, 0, 0, 504; 0, 0, 2520, 0, 0, 0,
 1260; 0, 0, 0, 2520, 0, 0, 840; 0, 0, 0, 0, 13860, 0, 6930; 0, 0, 0, 0, 0, 
5544, 0; 0, 0, 0, 0, 0, 0, 12012], [420, 420, 840, 630, 2982, 1092, 4159; 21
0, 280, 630, 504, 2415, 876, 3395; 140, 210, 504, 420, 2050, 749, 2901; 105,
 168, 420, 360, 1785, 658, 2542; 84, 140, 360, 315, 1582, 588, 2266; 70, 120
, 315, 280, 1421, 532, 2046; 60, 105, 280, 252, 1290, 486, 1866]]
? mathnf(amat,4)
[[420, 0, 0, 0, 210, 168, 175; 0, 840, 0, 0, 0, 0, 504; 0, 0, 2520, 0, 0, 0,
 1260; 0, 0, 0, 2520, 0, 0, 840; 0, 0, 0, 0, 13860, 0, 6930; 0, 0, 0, 0, 0, 
5544, 0; 0, 0, 0, 0, 0, 0, 12012], [420, 420, 840, 630, 2982, 1092, 4159; 21
0, 280, 630, 504, 2415, 876, 3395; 140, 210, 504, 420, 2050, 749, 2901; 105,
 168, 420, 360, 1785, 658, 2542; 84, 140, 360, 315, 1582, 588, 2266; 70, 120
, 315, 280, 1421, 532, 2046; 60, 105, 280, 252, 1290, 486, 1866]]
? mathnf(amat,3)
[[360360, 0, 0, 0, 0, 144144, 300300; 0, 27720, 0, 0, 0, 0, 22176; 0, 0, 277
20, 0, 0, 0, 6930; 0, 0, 0, 2520, 0, 0, 840; 0, 0, 0, 0, 2520, 0, 1260; 0, 0
, 0, 0, 0, 168, 0; 0, 0, 0, 0, 0, 0, 7], [51480, 4620, 5544, 630, 840, 20676
, 48619; 45045, 3960, 4620, 504, 630, 18074, 42347; 40040, 3465, 3960, 420, 
504, 16058, 37523; 36036, 3080, 3465, 360, 420, 14448, 33692; 32760, 2772, 3
080, 315, 360, 13132, 30574; 30030, 2520, 2772, 280, 315, 12036, 27986; 2772
0, 2310, 2520, 252, 280, 11109, 25803], [7, 6, 5, 4, 3, 2, 1]]
? mathnfmod(amat,matdetint(amat))

[420 0 0 0 210 168 175]

[0 840 0 0 0 0 504]

[0 0 2520 0 0 0 1260]

[0 0 0 2520 0 0 840]

[0 0 0 0 13860 0 6930]

[0 0 0 0 0 5544 0]

[0 0 0 0 0 0 12012]

? mathnfmodid(amat,123456789*10^100)

[60 0 0 0 30 24 35]

[0 120 0 0 0 0 24]

[0 0 360 0 0 0 180]

[0 0 0 360 0 0 240]

[0 0 0 0 180 0 90]

[0 0 0 0 0 72 0]

[0 0 0 0 0 0 12]

? 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]

? matimage([1,3,5;2,4,6;3,5,7])

[1 3]

[2 4]

[3 5]

? matimage([1,3,5;2,4,6;3,5,7],1)

[3 5]

[4 6]

[5 7]

? matimage(Pi*[1,3,5;2,4,6;3,5,7])

[9.4247779607693797153879301498385086525 15.70796326794896619231321691639751
4420]

[12.566370614359172953850573533118011536 18.84955592153875943077586029967701
7305]

[15.707963267948966192313216916397514420 21.99114857512855266923850368295652
0189]

? matimagecompl([1,3,5;2,4,6;3,5,7])
[3]
? matimagecompl(Pi*[1,3,5;2,4,6;3,5,7])
[1]
? matindexrank([1,1,1;1,1,1;1,1,2])
[[1, 3], [1, 3]]
? matintersect([1,2;3,4;5,6],[2,3;7,8;8,9])

[-1]

[-1]

[-1]

? matinverseimage([1,1;2,3;5,7],[2,2,6]~)
[4, -2]~
? matisdiagonal([1,0,0;0,5,0;0,0,0])
1
? matker(matrix(4,4,x,y,x/y))

[-1/2 -1/3 -1/4]

[1 0 0]

[0 1 0]

[0 0 1]

? matker(matrix(4,4,x,y,sin(x+y)))

[1.0000000000000000000000000000000000000 1.080604611736279434801873214885953
2074]

[-1.0806046117362794348018732148859532074 -0.1677063269057152260048635409984
7562046]

[1 0]

[0 1]

? matker(matrix(4,4,x,y,x+y),1)

[1 2]

[-2 -3]

[1 0]

[0 1]

? matkerint(matrix(4,4,x,y,x*y))

[-1 -1 -1]

[-1 0 1]

[1 -1 1]

[0 1 -1]

? matkerint(matrix(4,4,x,y,x*y),1)

[-1 -1 -1]

[-1 0 1]

[1 -1 1]

[0 1 -1]

? matkerint(matrix(4,6,x,y,2520/(x+y)),2)

[3 1]

[-30 -15]

[70 70]

[0 -140]

[-126 126]

[84 -42]

? matmuldiagonal(amat,[1,2,3,4,5,6,7])

[49 -2352 26460 -117600 242550 -232848 84084]

[-1176 75264 -952560 4515840 -9702000 9580032 -3531528]

[8820 -635040 8573040 -42336000 93555000 -94303440 35315280]

[-29400 2257920 -31752000 161280000 -363825000 372556800 -141261120]

[48510 -3880800 56133000 -291060000 667012500 -691558560 264864600]

[-38808 3193344 -47151720 248371200 -576298800 603542016 -233080848]

[12012 -1009008 15135120 -80720640 189189000 -199783584 77693616]

? matmultodiagonal(amat^-1,%)

[1 0 0 0 0 0 0]

[0 2 0 0 0 0 0]

[0 0 3 0 0 0 0]

[0 0 0 4 0 0 0]

[0 0 0 0 5 0 0]

[0 0 0 0 0 6 0]

[0 0 0 0 0 0 7]

? matpascal(8)

[1 0 0 0 0 0 0 0 0]

[1 1 0 0 0 0 0 0 0]

[1 2 1 0 0 0 0 0 0]

[1 3 3 1 0 0 0 0 0]

[1 4 6 4 1 0 0 0 0]

[1 5 10 10 5 1 0 0 0]

[1 6 15 20 15 6 1 0 0]

[1 7 21 35 35 21 7 1 0]

[1 8 28 56 70 56 28 8 1]

? matrank(matrix(5,5,x,y,x+y))
2
? matrix(5,5,x,y,gcd(x,y))

[1 1 1 1 1]

[1 2 1 2 1]

[1 1 3 1 1]

[1 2 1 4 1]

[1 1 1 1 5]

? matrixqz([1,3;3,5;5,7],0)

[1 1]

[3 2]

[5 3]

? matrixqz([1/3,1/4,1/6;1/2,1/4,-1/4;1/3,1,0],-1)

[19 12 2]

[0 1 0]

[0 0 1]

? matrixqz([1,3;3,5;5,7],-2)

[2 -1]

[1 0]

[0 1]

? matsize([1,2;3,4;5,6])
[3, 2]
? matsnf(matrix(5,5,j,k,random))
[741799239614624774584532992, 2147483648, 2147483648, 1, 1]
? matsnf(1/mathilbert(6))
[27720, 2520, 2520, 840, 210, 6]
? matsnf(x*matid(5)-matrix(5,5,j,k,1),2)
[x^2 - 5*x, x, x, x, 1]
? matsolve(mathilbert(10),[1,2,3,4,5,6,7,8,9,0]~)
[9236800, -831303990, 18288515520, -170691240720, 832112321040, -23298940665
00, 3883123564320, -3803844432960, 2020775945760, -449057772020]~
? matsolvemod([2,3;5,4],[7,11],[1,4]~)
[-5, -1]~
? matsolvemod([2,3;5,4],[7,11],[1,4]~,1)
[[-5, -1]~, [-77, 723; 0, 1]]
? matsupplement([1,3;2,4;3,6])

[1 3 0]

[2 4 0]

[3 6 1]

? mattranspose(vector(2,x,x))
[1, 2]~
? %*%~

[1 2]

[2 4]

? norml2(vector(10,x,x))
385
? qfgaussred(mathilbert(5))

[1 1/2 1/3 1/4 1/5]

[0 1/12 1 9/10 4/5]

[0 0 1/180 3/2 12/7]

[0 0 0 1/2800 2]

[0 0 0 0 1/44100]

? qfjacobi(mathilbert(6))
[[1.6188998589243390969705881471257800712, 0.2423608705752095521357284158507
0114077, 0.000012570757122625194922982397996498755027, 0.0000001082799484565
5497685388772372251711485, 0.016321521319875822124345079564191505890, 0.0006
1574835418265769764919938428527140264]~, [0.74871921887909485900280109200517
845109, -0.61454482829258676899320019644273870645, 0.01114432093072471053067
8340374220998541, -0.0012481940840821751169398163046387834473, 0.24032536934
252330399154228873240534568, -0.062226588150197681775152126611810492910; 0.4
4071750324351206127160083580231701801, 0.21108248167867048675227675845247769
095, -0.17973275724076003758776897803740640964, 0.03560664294428763526612284
8131812048466, -0.69765137527737012296208335046678265583, 0.4908392097109243
6297498316169060044997; 0.32069686982225190106359024326699463106, 0.36589360
730302614149086554211117169622, 0.60421220675295973004426567844103062241, -0
.24067907958842295837736719558855679285, -0.23138937333290388042251363554209
048309, -0.53547692162107486593474491750949545456; 0.25431138634047419251788
312792590944672, 0.39470677609501756783094636145991581708, -0.44357471627623
954554460416705180105301, 0.62546038654922724457753441039459331059, 0.132863
15850933553530333839628101576050, -0.41703769221897886840494514780771076439;
 0.21153084007896524664213667673977991959, 0.3881904338738864286311144882599
2418973, -0.44153664101228966222143649752977203423, -0.689807199293836684198
01738006926829419, 0.36271492146487147525299457604461742111, 0.0470340189331
15649705614518466541243873; 0.18144297664876947372217005457727093715, 0.3706
9590776736280861775501084807394603, 0.45911481681642960284551392793050866602
, 0.27160545336631286930015536176213647001, 0.502762866757515384892605663686
47786272, 0.54068156310385293880022293448123782121]]
? m=1/mathilbert(7)

[49 -1176 8820 -29400 48510 -38808 12012]

[-1176 37632 -317520 1128960 -1940400 1596672 -504504]

[8820 -317520 2857680 -10584000 18711000 -15717240 5045040]

[-29400 1128960 -10584000 40320000 -72765000 62092800 -20180160]

[48510 -1940400 18711000 -72765000 133402500 -115259760 37837800]

[-38808 1596672 -15717240 62092800 -115259760 100590336 -33297264]

[12012 -504504 5045040 -20180160 37837800 -33297264 11099088]

? mp=concat(m,matid(7))

[49 -1176 8820 -29400 48510 -38808 12012 1 0 0 0 0 0 0]

[-1176 37632 -317520 1128960 -1940400 1596672 -504504 0 1 0 0 0 0 0]

[8820 -317520 2857680 -10584000 18711000 -15717240 5045040 0 0 1 0 0 0 0]

[-29400 1128960 -10584000 40320000 -72765000 62092800 -20180160 0 0 0 1 0 0 
0]

[48510 -1940400 18711000 -72765000 133402500 -115259760 37837800 0 0 0 0 1 0
 0]

[-38808 1596672 -15717240 62092800 -115259760 100590336 -33297264 0 0 0 0 0 
1 0]

[12012 -504504 5045040 -20180160 37837800 -33297264 11099088 0 0 0 0 0 0 1]

? qflll(m)

[-420 -420 840 630 -1092 757 2982]

[-210 -280 630 504 -876 700 2415]

[-140 -210 504 420 -749 641 2050]

[-105 -168 420 360 -658 589 1785]

[-84 -140 360 315 -588 544 1582]

[-70 -120 315 280 -532 505 1421]

[-60 -105 280 252 -486 471 1290]

? qflll(m,7)

[-420 -420 840 630 -1092 757 2982]

[-210 -280 630 504 -876 700 2415]

[-140 -210 504 420 -749 641 2050]

[-105 -168 420 360 -658 589 1785]

[-84 -140 360 315 -588 544 1582]

[-70 -120 315 280 -532 505 1421]

[-60 -105 280 252 -486 471 1290]

? qflllgram(m)

[1 1 27 -27 69 0 141]

[0 1 4 -22 34 -24 49]

[0 1 3 -21 18 -24 23]

[0 1 3 -20 10 -19 13]

[0 1 3 -19 6 -14 8]

[0 1 3 -18 4 -10 5]

[0 1 3 -17 3 -7 3]

? qflllgram(m,7)

[1 1 27 -27 69 0 141]

[0 1 4 -22 34 -24 49]

[0 1 3 -21 18 -24 23]

[0 1 3 -20 10 -19 13]

[0 1 3 -19 6 -14 8]

[0 1 3 -18 4 -10 5]

[0 1 3 -17 3 -7 3]

? qflllgram(m,1)

[1 1 27 -27 69 0 141]

[0 1 4 -23 34 -24 91]

[0 1 3 -22 18 -24 65]

[0 1 3 -21 10 -19 49]

[0 1 3 -20 6 -14 38]

[0 1 3 -19 4 -10 30]

[0 1 3 -18 3 -7 24]

? qflllgram(mp~*mp,4)
[[-420, -420, 840, 630, 2982, -1092, -83; -210, -280, 630, 504, 2415, -876, 
70; -140, -210, 504, 420, 2050, -749, 137; -105, -168, 420, 360, 1785, -658,
 169; -84, -140, 360, 315, 1582, -588, 184; -70, -120, 315, 280, 1421, -532,
 190; -60, -105, 280, 252, 1290, -486, 191; 420, 0, 0, 0, -210, 168, 35; 0, 
840, 0, 0, 0, 0, 336; 0, 0, -2520, 0, 0, 0, 1260; 0, 0, 0, -2520, 0, 0, -840
; 0, 0, 0, 0, -13860, 0, 6930; 0, 0, 0, 0, 0, 5544, 0; 0, 0, 0, 0, 0, 0, -12
012], [0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 
0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0
; 1, 0, 0, 0, 0, 0, 0; 0, 1, 0, 0, 0, 0, 0; 0, 0, 1, 0, 0, 0, 0; 0, 0, 0, 1,
 0, 0, 0; 0, 0, 0, 0, 1, 0, 0; 0, 0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 0, 1]]
? qflll(m,1)

[-420 -420 840 630 -1092 -83 2982]

[-210 -280 630 504 -876 70 2415]

[-140 -210 504 420 -749 137 2050]

[-105 -168 420 360 -658 169 1785]

[-84 -140 360 315 -588 184 1582]

[-70 -120 315 280 -532 190 1421]

[-60 -105 280 252 -486 191 1290]

? qflll(m,2)

[-420 -420 -630 840 1092 2982 -83]

[-210 -280 -504 630 876 2415 70]

[-140 -210 -420 504 749 2050 137]

[-105 -168 -360 420 658 1785 169]

[-84 -140 -315 360 588 1582 184]

[-70 -120 -280 315 532 1421 190]

[-60 -105 -252 280 486 1290 191]

? qflll(mp,4)
[[-420, -420, 840, 630, 2982, -1092, -83; -210, -280, 630, 504, 2415, -876, 
70; -140, -210, 504, 420, 2050, -749, 137; -105, -168, 420, 360, 1785, -658,
 169; -84, -140, 360, 315, 1582, -588, 184; -70, -120, 315, 280, 1421, -532,
 190; -60, -105, 280, 252, 1290, -486, 191; 420, 0, 0, 0, -210, 168, 35; 0, 
840, 0, 0, 0, 0, 336; 0, 0, -2520, 0, 0, 0, 1260; 0, 0, 0, -2520, 0, 0, -840
; 0, 0, 0, 0, -13860, 0, 6930; 0, 0, 0, 0, 0, 5544, 0; 0, 0, 0, 0, 0, 0, -12
012], [0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 
0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0; 0, 0, 0, 0, 0, 0, 0
; 1, 0, 0, 0, 0, 0, 0; 0, 1, 0, 0, 0, 0, 0; 0, 0, 1, 0, 0, 0, 0; 0, 0, 0, 1,
 0, 0, 0; 0, 0, 0, 0, 1, 0, 0; 0, 0, 0, 0, 0, 1, 0; 0, 0, 0, 0, 0, 0, 1]]
? qflll(m,3)

[-420 -420 840 630 -1092 -83 2982]

[-210 -280 630 504 -876 70 2415]

[-140 -210 504 420 -749 137 2050]

[-105 -168 420 360 -658 169 1785]

[-84 -140 360 315 -588 184 1582]

[-70 -120 315 280 -532 190 1421]

[-60 -105 280 252 -486 191 1290]

? qfminim([2,1;1,2],4,6)
[6, 2, [0, -1, 1; 1, 1, 0]]
? qfperfection([2,0,1;0,2,1;1,1,2])
6
? qfsign(mathilbert(5)-0.11*matid(5))
[2, 3]
? aset=Set([5,-2,7,3,5,1])
["-2", "1", "3", "5", "7"]
? bset=Set([7,5,-5,7,2])
["-5", "2", "5", "7"]
? setintersect(aset,bset)
["5", "7"]
? setisset([-3,5,7,7])
0
? setminus(aset,bset)
["-2", "1", "3"]
? setsearch(aset,3)
3
? setsearch(bset,3)
0
? setunion(aset,bset)
["-2", "-5", "1", "2", "3", "5", "7"]
? trace(1+I)
2
? trace(Mod(x+5,x^3+x+1))
15
? Vec(sin(x))
[1, 0, -1/6, 0, 1/120, 0, -1/5040, 0, 1/362880, 0, -1/39916800, 0, 1/6227020
800, 0, -1/1307674368000]
? vecmax([-3,7,-2,11])
11
? vecmin([-3,7,-2,11])
-3
? concat([1,2],[3,4])
[1, 2, 3, 4]
? concat(Mat(vector(4,x,x)~),vector(4,x,10+x)~)

[1 11]

[2 12]

[3 13]

[4 14]

? vecextract([1,2,3,4,5,6,7,8,9,10],1000)
[4, 6, 7, 8, 9, 10]
? vecextract(matrix(15,15,x,y,x+y),vector(5,x,3*x),vector(3,y,3*y))

[6 9 12]

[9 12 15]

[12 15 18]

[15 18 21]

[18 21 24]

? (1.*mathilbert(7))^(-1)

[49.000000000000000000000000000000045975 -1176.00000000000000000000000000000
20892 8820.0000000000000000000000000000216289 -29400.00000000000000000000000
0000087526 48510.000000000000000000000000000164477 -38808.000000000000000000
000000000145051 12012.000000000000000000000000000048237]

[-1176.0000000000000000000000000000007015 37632.0000000000000000000000000000
36155 -317520.00000000000000000000000000039285 1128960.000000000000000000000
0000016298 -1940400.0000000000000000000000000031060 1596672.0000000000000000
000000000027521 -504504.00000000000000000000000000091794]

[8819.9999999999999999999999999999987063 -317520.000000000000000000000000000
01369 2857680.0000000000000000000000000004729 -10584000.00000000000000000000
0000002587 18711000.000000000000000000000000005552 -15717240.000000000000000
000000000005216 5045040.0000000000000000000000000017929]

[-29399.999999999999999999999999999970929 1128959.99999999999999999999999999
90570 -10583999.999999999999999999999999992003 40319999.99999999999999999999
9999971163 -72764999.999999999999999999999999949359 62092799.999999999999999
999999999957242 -20180159.999999999999999999999999986112]

[48509.999999999999999999999999999911823 -1940399.99999999999999999999999999
68289 18710999.999999999999999999999999971121 -72764999.99999999999999999999
9999890954 133402499.99999999999999999999999980291 -115259759.99999999999999
999999999983068 37837799.999999999999999999999999944464]

[-38807.999999999999999999999999999899366 1596671.99999999999999999999999999
62508 -15717239.999999999999999999999999965108 62092799.99999999999999999999
9999866538 -115259759.99999999999999999999999975693 100590335.99999999999999
999999999979026 -33297263.999999999999999999999999931034]

[12011.999999999999999999999999999960320 -504503.999999999999999999999999998
49528 5045039.9999999999999999999999999858501 -20180159.99999999999999999999
9999945550 37837799.999999999999999999999999900488 -33297263.999999999999999
999999999913962 11099087.999999999999999999999999971679]

? vecsort([8,7,6,5],,1)
[4, 3, 2, 1]
? vecsort([[1,5],[2,4],[1,5,1],[1,4,2]],,2)
[[1, 4, 2], [1, 5], [1, 5, 1], [2, 4]]
? vecsort(vector(17,x,5*x%17))
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]
? vecsort([[1,8,5],[2,5,8],[3,6,-6],[4,8,6]],2)
[[2, 5, 8], [3, 6, -6], [4, 8, 6], [1, 8, 5]]
? vecsort([[1,8,5],[2,5,8],[3,6,-6],[4,8,6]],[2,1])
[[2, 5, 8], [3, 6, -6], [1, 8, 5], [4, 8, 6]]
? vector(10,x,1/x)
[1, 1/2, 1/3, 1/4, 1/5, 1/6, 1/7, 1/8, 1/9, 1/10]
? setrand(1);a=matrix(3,5,j,k,vectorv(5,l,random\10^8))

[[10, 7, 8, 7, 18]~ [17, 0, 9, 20, 10]~ [5, 4, 7, 18, 20]~ [0, 16, 4, 2, 0]~
 [17, 19, 17, 1, 14]~]

[[17, 16, 6, 3, 6]~ [17, 13, 9, 19, 6]~ [1, 14, 12, 20, 8]~ [6, 1, 8, 17, 21
]~ [18, 17, 9, 10, 13]~]

[[4, 13, 3, 17, 14]~ [14, 16, 11, 5, 4]~ [9, 11, 13, 7, 15]~ [19, 21, 2, 4, 
5]~ [14, 16, 6, 20, 14]~]

? setrand(1);as=matrix(3,3,j,k,vectorv(5,l,random\10^8))

[[10, 7, 8, 7, 18]~ [17, 0, 9, 20, 10]~ [5, 4, 7, 18, 20]~]

[[17, 16, 6, 3, 6]~ [17, 13, 9, 19, 6]~ [1, 14, 12, 20, 8]~]

[[4, 13, 3, 17, 14]~ [14, 16, 11, 5, 4]~ [9, 11, 13, 7, 15]~]

? getheap
[111, 12331]
? print("Total time spent: ",gettime);
Total time spent: 180
? \q