? ? ?    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])
% = 
[0.28334945180064027179781065475712672521 + 0.E-39*I 1 -1.2833494518006402
717978106547571267252 + 0.E-38*I]

[0.64167472590032013589890532737856336260 + 0.E-38*I -2 -0.141674725900320
13589890532737856336259 + 0.E-39*I]

[1 1 1]

? mathess(mathilbert(7))
% = 
[1 90281/58800 -1919947/4344340 4858466341/1095033030 -77651417539/8196787
326 3386888964/106615355 1/2]

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

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

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

[0 0 0 9773092581/35395807550620 -24363634138919/107305824577186620 721182
03606917/60481351061158500 55899/3088554700]

[0 0 0 0 67201501179065/8543442888354179988 -9970556426629/740828619992676
600 -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; 210, 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, 1
866]]
? mathnf(amat,2)
% = [[360360, 0, 0, 0, 0, 144144, 300300; 0, 27720, 0, 0, 0, 0, 22176; 0, 
0, 27720, 0, 0, 0, 6930; 0, 0, 0, 2520, 0, 0, 840; 0, 0, 0, 0, 2520, 0, 12
60; 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, 3080, 315, 360, 13132, 30574; 30030, 2520, 2772, 280, 315, 12
036, 27986; 27720, 2310, 2520, 252, 280, 11109, 25803], [7, 6, 5, 4, 3, 2,
 1]]
? mathnf(amat,3)
% = [[360360, 0, 0, 0, 0, 144144, 300300; 0, 27720, 0, 0, 0, 0, 22176; 0, 
0, 27720, 0, 0, 0, 6930; 0, 0, 0, 2520, 0, 0, 840; 0, 0, 0, 0, 2520, 0, 12
60; 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, 3080, 315, 360, 13132, 30574; 30030, 2520, 2772, 280, 315, 12
036, 27986; 27720, 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])
% = 
[3.1415926535897932384626433832795028841 9.4247779607693797153879301498385
086525]

[6.2831853071795864769252867665590057683 12.566370614359172953850573533118
011536]

[9.4247779607693797153879301498385086525 15.707963267948966192313216916397
514420]

? matimagecompl([1,3,5;2,4,6;3,5,7])
% = [3]
? matimagecompl(Pi*[1,3,5;2,4,6;3,5,7])
% = [3]
? 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.0806046117362794348018732148859
532074]

[-1.0806046117362794348018732148859532074 -0.16770632690571522600486354099
847562046]

[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, -23298
94066500, 3883123564320, -3803844432960, 2020775945760, -449057772020]~
? matsolvemod([2,3;5,4],[7,11],[1,4]~)
% = [[-5, -1]~, [-77, 723; 0, 1]]
? matsolvemod([2,3;5,4],[7,11],[1,4]~,1)
% = [-5, -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.2423608705752095521357284
1585070114077, 0.000012570757122625194922982397996498755027, 0.00000010827
994845655497685388772372251711485, 0.0163215213198758221243450795641915058
90, 0.00061574835418265769764919938428527140264]~, [0.74871921887909485900
280109200517845109, -0.61454482829258676899320019644273870645, 0.011144320
930724710530678340374220998541, -0.001248194084082175116939816304638783447
3, 0.24032536934252330399154228873240534568, -0.06222658815019768177515212
6611810492910; 0.44071750324351206127160083580231701801, 0.211082481678670
48675227675845247769095, -0.17973275724076003758776897803740640964, 0.0356
06642944287635266122848131812048466, -0.6976513752773701229620833504667826
5583, 0.49083920971092436297498316169060044997; 0.320696869822251901063590
24326699463106, 0.36589360730302614149086554211117169622, 0.60421220675295
973004426567844103062241, -0.24067907958842295837736719558855679285, -0.23
138937333290388042251363554209048309, -0.535476921621074865934744917509495
45456; 0.25431138634047419251788312792590944672, 0.39470677609501756783094
636145991581708, -0.44357471627623954554460416705180105301, 0.625460386549
22724457753441039459331059, 0.13286315850933553530333839628101576050, -0.4
1703769221897886840494514780771076439; 0.211530840078965246642136676739779
91959, 0.38819043387388642863111448825992418973, -0.4415366410122896622214
3649752977203423, -0.68980719929383668419801738006926829419, 0.36271492146
487147525299457604461742111, 0.047034018933115649705614518466541243873; 0.
18144297664876947372217005457727093715, 0.37069590776736280861775501084807
394603, 0.45911481681642960284551392793050866602, 0.2716054533663128693001
5536176213647001, 0.50276286675751538489260566368647786272, 0.540681563103
85293880022293448123782121]]
? 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 -83 2562]

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

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

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

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

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

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

? qflll(m,7)
% = 
[-420 -420 840 630 -1092 -83 2562]

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

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

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

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

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

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

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

[0 1 3 -22 18 -24 23]

[0 1 3 -21 10 -19 13]

[0 1 3 -20 6 -14 8]

[0 1 3 -19 4 -10 5]

[0 1 3 -18 3 -7 3]

? qflllgram(mp~*mp,4)
% = [[-420, -420, 840, 630, 2982, -1092, 757; -210, -280, 630, 504, 2415, 
-876, 700; -140, -210, 504, 420, 2050, -749, 641; -105, -168, 420, 360, 17
85, -658, 589; -84, -140, 360, 315, 1582, -588, 544; -70, -120, 315, 280, 
1421, -532, 505; -60, -105, 280, 252, 1290, -486, 471; 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, -12012], [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 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,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, 757; -210, -280, 630, 504, 2415, 
-876, 700; -140, -210, 504, 420, 2050, -749, 641; -105, -168, 420, 360, 17
85, -658, 589; -84, -140, 360, 315, 1582, -588, 544; -70, -120, 315, 280, 
1421, -532, 505; -60, -105, 280, 252, 1290, -486, 471; 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, -12012], [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 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]

? 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/6
227020800, 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)
% = 
[48.999999999999999999999999999999627248 -1175.999999999999999999999999999
9851549 8819.9999999999999999999999999998567354 -29399.9999999999999999999
99999999441659 48509.999999999999999999999999998973568 -38807.999999999999
999999999999999110386 12011.999999999999999999999999999707258]

[-1175.9999999999999999999999999999839914 37631.99999999999999999999999999
9359946 -317519.99999999999999999999999999380803 1128959.99999999999999999
99999999758214 -1940399.9999999999999999999999999554891 1596671.9999999999
999999999999999613918 -504503.99999999999999999999999998728131]

[8819.9999999999999999999999999998374720 -317519.9999999999999999999999999
9348578 2857679.9999999999999999999999999368820 -10583999.9999999999999999
99999999753258 18710999.999999999999999999999999545401 -15717239.999999999
999999999999999605463 5045039.9999999999999999999999998699630]

[-29399.999999999999999999999999999342876 1128959.999999999999999999999999
9736178 -10583999.999999999999999999999999744107 40319999.9999999999999999
99999998998936 -72764999.999999999999999999999998154672 62092799.999999999
999999999999998397863 -20180159.999999999999999999999999471783]

[48509.999999999999999999999999998758142 -1940399.999999999999999999999999
9500817 18710999.999999999999999999999999515458 -72764999.9999999999999999
99999998103489 133402499.99999999999999999999999650276 -115259759.99999999
999999999999999696284 37837799.999999999999999999999998998454]

[-38807.999999999999999999999999998900612 1596671.999999999999999999999999
9557685 -15717239.999999999999999999999999570419 62092799.9999999999999999
99999998317978 -115259759.99999999999999999999999689747 100590335.99999999
999999999999999730508 -33297263.999999999999999999999999111177]

[12011.999999999999999999999999999631826 -504503.9999999999999999999999999
8517694 5045039.9999999999999999999999998559749 -20180159.9999999999999999
99999999435909 37837799.999999999999999999999998959308 -33297263.999999999
999999999999999095899 11099087.999999999999999999999999701779]

? vecindexsort([8,7,6,5])
% = [4, 3, 2, 1]
? veclexsort([[1,5],[2,4],[1,5,1],[1,4,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, 12168]
? print("Total time spent: ",gettime);
Total time spent: 1286
? \q
Good bye!
