M7R0
I6IntegrableConnections=6"%&falseGE\[l[s%.sim_sylvesterGt+6,%"AG%"BG%"CG%"DG%"E
G%"FG%#XYG%#EFG%"cG%"hG6-%"mG%"nG%$extG%#AAG%"iG%"jG%%xyisG%$eqsG%$XY1G%$EF1G%"
sGF$F$C*@$4333-%%typeG6$9$-.%'matrixG6$%)anythingG.%'squareG-FG6$9%FJ-FG6$9&FJ-
FG6$9'FJ-%&ERRORG6#%9square~matrices~expectedG>8$-&%'linalgG6#.%'rowdimG6#FI>8%
-F[o6#FV>8,-&F\o6#FK6$Fin,$Fbo""#>8-Fgo@%5/9#""'32FcpFbp4-FG6$&9"6#""(.%&arrayG
C'>8*<#-%$seqG6$-Fdq6$&Ffo6$8(8)/Fjq;"""Fin/F[r;F^rF[p@%/Fbp"#5>8'-%&evalmG6#,&
FIF^r*&9,F^rFSF^r!"">FerFI>8+<$-Fdq6$-Fdq6$/,&-%(convertG6$7#-Fdq6$*&&Fer6$Fjq%
"kGF^r&Ffo6$F`tF[rF^r/F`tF]r%"+GF^r-Fhs6$7#-Fdq6$*&&Ffo6$Fjq,&F`tF^rFboF^rF^r&F
VFbtF^r/F`t;F^rFboFdtF\s&F^pFiqF\r/F[rF`u-Fdq6$-Fdq6$/,&-Fhs6$7#-Fdq6$*&&FSF_tF
^rFatF^rFctFdtF^r-Fhs6$7#-Fdq6$*&F[uF^r&FYFbtF^rF_uFdtF\s&F^p6$Fjq,&F[rF^rFboF^
rF\rFbu-%'assignG6#-&%+SolveToolsG6#%'LinearG6$F_sFaq@$1"")FbpC$>9*Ffo>9+F^pC$>
FfoFhw>F^pFjw@%Fbr>8./F[s9->F`x%%NULLG7$-%$mapG6$%&evalaG-%%subsG6&F`x-Fdq6$-Fd
q6$/Fau&9(FiqF\rFbu-Fdq6$-Fdq6$/Fgv&9)FiqF\rFbu-&F\o6#.%*submatrixG6%FfoF]rF`u-
Fgx6$Fix-F[y6&F`xF]yFdy-F\z6%FfoF]r;,&FboF^rF^rF^rF[pF$F$F$6#Fcr%)G_matrixGt+6'
%"KG%"xGF`t%'lambdaG%'invarsG6.F?%$nisGF6%"GG%#ccG%"zG%&zcolsGF9F:%#jjGF3%"MGF$
F$C0>Fin&Fcy6#F^r>Fbo-%'subsopG6$/F^rFdxFcy>8&-&F"6#%,mat_get_dimGF`o>Fer-Fho6%
Fc\lFc\l""!>8/-&F"6#%)mat_evalG6&FIFS,$FinF\s,&FVF^rFinF\s>F^pF^r?(F`xF^rF^r,&F
VF^rF\sF^r%%trueGC$?(FfoF^pF^r,(F^pF^r&Fbo6#F`xF^rF\sF^rFh]l?(F_sF^rF^rFc\lFh]l
>&Fer6$F_sFfo-%&coeffG6%-%'normalG6#&F]]lFa^lFS,(FinF\sF`xF^rF\sF^r>F^pFfo>FaqF
dx?(FfoF^pF^rFc\lFh]lC%>F[rFh]l?(F_sF^rF^rFc\lFh]lC$>Fjq-Fc^l6%Fe^lFS,(FinF\sFV
F^rF\sF^r@$0FjqF[]lC$>F[rF%>F`^lFjq@$F[r>Faq6$FaqFfo?&Ffo7#FaqFh]lC$?(F_sF^rF^r
Fc\lFh]l>F`^l-Fc^l6%Fe^lFSFd]l>&Fer6$FfoFfo,&Fe`lF^rFYF^r>F`x7#-Fdq6$-Fdq6$,(Fi
nF\sFfoF^rF\sF^r/F_s;F^r&Fbo6#Ffo/Ffo;F^rFg]l-%)userinfoG6&F^r.%-super_reduceG%
,pole~order:GFin-Ffal6&F^rFhal%,valuations:GF`x-Fgr6#FerF$6#%&IndexGF$6#""&%/ma
t_differenceGt+6%F*F+F^[l6(F,%#vaG%#vbGF3%#CtG%#CCGF$F$C3-%(readlibG6#.%&utilsG
@$-FG6$%0Valuation_BoundG.%%nameG>Fecl"$+">Fin-F^q6#;F^rF\q>&FinF\\l.%+differen
ceG>Fbo-&F"6#%.mat_valuationG6$FIFV>Fc\l-Fddl6$FSFV>Fer-%$minG6$-&F"6#%.mat_get
_orderGFgdl-F`elFjdl>Fjq-&F"6#%&EvalmG6$,&-F_]l6&FIFVFboFerF^r-F_]l6&FSFVFc\lFe
rF\sFV>&Fin6#F\p.%#NNG>&Fin6#""$Fjq>&Fin6#""%Fer>&Fin6#FcblFdx>&Fin6#FcpFI>&Fin
F[qFS>F[r-%+new_matrixGF$>&%0Mat_DescriptionG6#F[rFin-%'RETURNGFjglF$6$FiglFecl
F$6#FgflFg\lt+6#F][l6#Fh[lF$F$C$>Fin&FiglF`o@;/F_dl.%,transformedG-F[o6#Fifl/F_
dl.%3pseudo_transformedGFihl/F_dl.%1changed_exp_partG-Fe\l6#Fcgl/F_dl.%'copiedG
-Fe\l6#F`fl/F_dl.%,substitutedGFail/F_dl.%)ramifiedGFail/F_dl.%)splittedG-F[o6#
Fefl/F_dl.%*localizedG-F[o6#F`gl/F_dl.%.block_reducedGF]gl/F_dl.%5direct_block_
reducedGFcgl/F_dl.%(swappedG-Fe\lFgjl/F_dl.%1elementary_transGFail-Fen6$Q0Unkno
wn~symbol:F$F_dlF$6#FiglF$6#F^r%/l_moser_reduceGt+6&F+F^[lF_[l%%optsG6#%&opts1G
F$F$C$@%1F[glFbp>FinFY>Fin<"-&F"6#%/l_super_reduceG6&FIFSFV-%&unionG6$Fin<#/.%%
stepGF^rF$F$F$6#F[gl%3minimalIntegerRootGt+6%%"fG%"uG%"vG6+%#roG%#rrG%#nuG%$lcf
GF9F:%$lirG%#i0G%&foundGF$F$C*>Fer7#-%'coeffsG6$-%'expandGF`oFV>F_sF^r>FjqF\p?(
F$F^rF^rF$1Fjq-%%nopsGF_blC$@$2-%'degreeG6$&Fer6#FjqFS-Fc_m6$&Fer6#F_sFS>F_sFjq
>Fjq,&FjqF^rF^rF^r>Fbo-%&rootsGFh_m>Faq7"@$0FboFb`mC$?&FinFboFh]l@$-FG6$F_dl.%(
integerG>Faq7$F_dl-%#opG6#Faq>Faq-%%sortGF`am@%/FaqFb`m-F\hl6#%)infinityGC&>Ffo
F%>F[rF^r?(F$F^rF^rF$34Ffo1F[r-F^_mF`amC$>Fc\l&FaqFjgl@%/-%)simplifyG6#-F[y6$/F
SFc\lFIF[]l>FfoFh]l>F[r,&F[rF^rF^rF^r@%Ffo-F\hl6#Fc\lFfamF$F$F$6#Fgfl%%qtcdGt+6
-Fh[lF^[l%"pGF6F`t%%ntabG%#mkGF_[lFc[l%"TG%%invTG62%#M1G%"qG%$mk1G%#GqG%#G1G%"r
G%"YG%"bGF9F:%&ntab1G%"PG%%invPGF,F-%%invDGF$F$C0-Ffal6%Fcbl.%)ratsuperG%5enter
~qtcd-algorithmG>Fin-%%copyGF`o>FboF[]l>Fc\l-%%evalG6#Fhw>Fjq-F`em6#F[s?(F$F^rF
^rF$Fh]lC8>Fjq-F[y6$/FjwF[]l-FdemFf_m-Ffal6&FcblF[em%$q~=GFbo>Fer-&F\o6#F_z6%Fj
q;F^r,&FYF^rFboF\sFhfm>F[r-%&prankG6%-Fefm6%Fer;F^rFc\lFhfmFSFV@$2F[rFc\lC$-Ffa
l6%FcblF[em%<rows~are~linearly~dependentG[-&F"6#%-pgaussUnimodG6&-&F\o6#%*trans
poseGF_blFSFV.F]]l>Faq-&F\o6#%*subvectorG6%F]]lFhfmF^r@$/&Faq6#FifmF[]lC,-Ffal6
%FcblF[em%Epermute~matrix~and~solution~vector~YG>Ffo,(FYF^rFboF\sF\sF^r?(F$F^rF
^rF$/&FaqFbalF[]l>Ffo,&FfoF^rF\sF^r>FjhmFdim>FdimF[]l>Fin-&F\o6#%(swaprowG6%Fin
FifmFfo>Fin-&F\o6#%(swapcolGF^jm@$/Fbp"#6C$>9.-F[jm6%FijmFifmFfo>Fbx-Fajm6%FbxF
ifmFfo>Fjq-F[jm6%FjqFifmFfo>Fjq-FajmFa[n@%2F^r-Fc_m6$FVFSC$-%&gcdexG6'FjhmFVFS.
%"UG.%"VG>Faq-Fgx6&%$remG-Fgr6#*&F]\nF^rFaqF^rFVFS>Faq-Ffel6#*&FjhmF\sFaqF^r>F]
]l-&F"6#%#IdG6#FY?(FfoF^rF^rFifmFh]l>&F]]l6$FifmFfoFdim>80-&F\o6#%(inverseG6#F]
]l>Fin-Fgr6#-%#&*G6$F]]l,&-F`^n6$FinFf]nF^r-Fgx6%%%diffGFf]nFSF\s@$FejmC$>Fbx-F
gr6#-F`^n6$FbxFf]n>Fijm-Fgr6#-F`^n6$F]]lFijm@$/Fbo,&FYF^rFc\lF\sFggm-Ffal6%Fcbl
F[em%Grecomputing~pole~order,~Glambda~and~mkG@%Fejm>Fin-&F"6#%,sortColumnsG6/Fi
nFSFVFYFcyFjw.F?.F`x.%%valsG.Fc\l.FjqFbxFijm>Fin-F]`n6-FinFSFVFYFcyFjwFa`nFb`nF
c`nFe`nFf`n-Ffal6&FgflF[emFjalF?-Ffal6&F\pF[emF]bl-F_am6#-Fhs6$Fd`n.%%listG-Ffa
l6&FgflF[em%&ntab:G-Fhs6$F`xFban-Ffal6&FgflF[em%$mk:GFhw>FboFhz-Ffal6&FcblF[em%
4exit~qtcd~with~q~=~GFbo-Ffal6%FcblF[em%Aapplying~diagonal~transformationGFhwFb
o>82-&F\o6#%%diagG6%-Fdq6$FV/Ffo;F^rFhw-Fdq6$F^r/Ffo;F^r,(FYF^rFhwF\sFboF\s-Fdq
6$FV/FfoF`u>83-Ffbn6%-Fdq6$*$FVF\sF\cnF^cn-Fdq6$F\dnFecn@$FejmC$>Fbx-Fgr6#-F`^n
6$FbxFdbn>Fijm-Fgr6#-F`^n6$FgcnFijm-Fgr6#-F`^n6$Fgcn,&-F`^n6$FinFdbnF^r-Fgx6%Fg
^nFdbnFSF\sF$6#%#MMGF$6#Ffjm%)mat_swapGt+6&F*F^[lF9F:6(Fh[lF+Fh]mF3%#BtG%#BBGF$
F$C3@$Fccl>FeclFicl>Fin-F^q6#;F^rFcp>&FboF\\lF_[m>Fc\l-Fddl6$FIFS>Fer-F`elFhfn>
Fjq-F_]l6&FIFSFc\lFer>&FboFaflFc\l>&FboFfflFjq>&FboFjflFer>&FboF^gl7$FVFY>&FboF
aglFI>F[rFegl>FhglFbo>Fjq-&F\o6#.Fcjm6%FjqFVFY>Fjq-&F\o6#.F]jmF`hn>FeflFjqF[hlF
$F]hlF$6#F[gl%+randomPolyGt+6$F^[l%"dG6#%#d1GF$F$C$@%/FbpF^r>FinF[gl>FinFS-F\hl
6#-%)randpolyG6&FI/.Fe^m-%%randG6#;!#5Fcr/.%'exponsG-F[jn6#;F^rFS/.%&termsGFinF
$F$F$6#F\p%'matvalGt+6&Fh[lF^[lFfcm%$lcMG6*F5F6%#lcGF9F:%$tmpG%#v1GF?F$F$C&>Faq
Fham>Fin-&F\o6#F_oF`o>Fbo-&F\o6#%'coldimGF`o@%/FbpFgflC$?(FerF^rF^rFinFh]l?(Fjq
F^rF^rFboFh]l@$0&FI6$FerFjqF[]l>Faq-F]el6$Faq-&F"6#%%val1G6%Fc\oFSFV-F\hlF`amC&
>FY-&F\o6#FL6%FinFboF[]l?(FerF^rF^rFinFh]l?(FjqF^rF^rFboFh]l@$Fb\oC%>F_s-Fi\o6&
Fc\oFSFVFe`n>&F[rFd\o7$F_sFc\l>Faq-F]el6$FaqF_s?(FerF^rF^rFinFh]l?(FjqF^rF^rFbo
Fh]l@$/&F\^oF\\lFaq>&FYFd\o&F\^oFaflF]]oF$F$F$6#F[glFaclt+F$F$F$F$F$F$F$F$6#F[]
l%1mat_get_exp_partGt+6$F][lF^[lFahlF$F$C%@$/-&F"6#%-mat_get_typeGF`o.%*algebra
icG-F\hl6#-&F"6#%3mat_get_eigenvalueG6#Fjp>FinFdhl@6Ffhl-&F"6#F\_o6$&Fin6#""*FS
F[ilFa`oF^il-Fh^m6#,&-Fb`o6$FcglFSF^r*&F]glF^r-Fg^n6$-&F"6#%5mat_get_ramificati
onGFhfnFSF\sF^rFcil-Fb`o6$F`flFSFhil-%(collectG6%-F[y6$F]glF[aoFSFixF[jl-F[y6$/
FS*&F]glF^r)FSF`glF^rF[aoFhjl@%/-F^_mFgjlF^r-Fb`o6$&Fin6#"#9FS-Fb`o6$&Fin6#"#CF
SFcjlF[]l/F_dl.%'simpleG-Fb`o6$F`glFSF^[mF`coF$Fh[mF$6#F\p%-convert_backGt+6&F-
%"NGF^[l%#kkG6%%#iDGFh]mF*F$F$C&>Fin-&F\o6#.Fj]nF`o>Fbo-&F"6#%(mat_valG6%FinFVF
[]l>Fc\l-Fgr6#*&)FVFboF^r-Fgx6%Fhs-Fgx6&%'taylorG-&F"6#%(normalmG6#*&-F`^n6$Fin
FSF^rFjdoF\sFV,&FYF^rF^rF^r.%(polynomGF^r-Fgx6&FgaoFc\lFVFixF$F$F$6#F[gl%<mat_g
et_indicial_polynomialGt+6$F][lFg]m6#F_[oF$F$C$>Fin-&F"6#%,simple_evalG6&FIF^[l
F[]lF[]l-&F\o6#.%$detG6#-Fgr6#,&F`flF^r*&FSF^rF_dlF^rF\sF$F$F$6#F\p%*MexpsoldeG
t+6%F*F^[lF]\m6/%#t0GF_[o%$solGF_\m%#eqG%#xxG%$SolGF?F9%"RG%$IndG%%testGF`tF$F$
C,@%1FgflFbp>Fer&Fjp6#Fbp>FerFe\m@$4-FG6$Fer.%$setG>Fer<#Fer>FablF[]l>Fin-%%tim
eGF$@$4-FG6$FIFKC%@'-FG6$FI.%)equationGC$>F[r-F_am6$F^r&FjpFafl>Fbo-&F"6#%0syst
em_expsoldeG6%-%.companion_sysG6$-%$lhsG6#&FjpF\\lFajoF_joFer-FG6$-%&numerGF`oF
ieoC%>Fjq7#-%(getLODOGF^`o>F[r-F_am6$F^r&FjqFafl>Fbo--F^cl6#.%)expsoldeG6%-Fijo
6$&FjqF\\lF[\pF[rFerC$>F[rF_jo>Fbo-Fdjo6%-Fijo6$F^[pFajoF_joFer-Ffal6&F^rF`\p%2
terminated.~Time:G,&FaioF^rFinF\s-F\hl6#<#-Fdq6$*&-%$expG6#-%$IntG6$&Fc\lF\\lF[
rF^r&&Fc\lFaflF\\lF^r/Fc\lFbo>Fbo-FdjoF^`oF]]p>FaqFe\m?(FfoF^rF^r-F^_m6#FboFh]l
C(>&F_sFbal&FaalFafl>&F`xFbal-%'removeG6$f*6#%"XGF$6$%)operatorG%&arrowGF$-%$ha
sG6$FI%'RootOfGF$F$F$-%&minusG6$-F[`p6$-Fhs6$-Fgx6$F_am-Fgx6$%'indetsGFi^pF\io-
Fhs6$-Fgx6$F_am-Fgx6$Fe`pFIF\io<#FS>&F^pFbal-&F\o6#%(augmentG6#-Fdq6$-Fgx6$f*Fa
_pF$Fc_pF$-Fc^l6$FI&&T#6#T%6#F:F$F$6&F^hoF`xF9FfoFi^p/F:;F^r-F^_m6#F\_p>F]]lF^r
@$2F^rFfoC$>Ff]nF^r?(F$F^rF^rF$3/F]]lF^r1Ff]nFfimC$@$3/-Fgbm6#,&&FaalF\\lF^r&&F
bo6#Ff]nF\\lF\sF[]l-&F\o6#%'iszeroG6#,&F^apF^r&F^pFicpF\s>F]]lF[]l>Ff]n,&Ff]nF^
rF^rF^r@$F]cp>Faq-F\]m6$Faq<#7$Ffcp-Fgx6$Ff^l-Fgr6#F^apF]]oF$F`blF$6#Fgfl%)back
subsGt+6$F*F26.F^[lF`t%"eG%"SGF9F:F5F6%%varsG%)leadvarsG%)leadcolsG%&countGF$F$
C5@$4-FG6$FIFL-Fen6#%0matrix~expectedG>FaqFjn>F_s-&F\o6#.F[\oF`o>F]]lF[]l>F`x-&
F\o6#.%'vectorG6$F_sF[]l>F^pFifp?(FjqF^rF^rFaqFh]l?(F[rF^rF^rF_sFh]l@$0&FIFiqF[
]lC&>&F`xFjglFjq>&F^pFf_mF[r>F]]l,&F]]lF^rF^rF^rFggm?(FjqF^rF^r,&F_sF^rF\sF^rFh
]l@$531&F^p6#F]`mFigp0FbhpF[]l3Fdhp/FigpF[]l-Fen6#%Xmatrix~must~be~in~echelon~f
orm,~zero~rows~at~the~bottomG@$0&F`xFj_mF[]l-Fen6#%4inconsistent~systemG>Ffo-Fj
fpFj_m?(FjqF^rF^rF]hpFh]l@%0&F`xFf_mF[]l>&FfoFf_m&FinFf_m@%/FbpF\p>Fgip&FSFf_m>
Fgip-%".G6$%"tGFjq>&FfoFj_mF[]l-Ffal6&F^rF\o%'vars~=G-%&printGFbal>Fc\l-&F\o6#.
%)multiplyG6$FIFfo>Fbo-&F\o6#.%$colG6$FIF_s>Fer-Fjfp6$F]hpF[]l?(FjqF]]lF\sF^rFh
]l>&Fin6#Figp-%$rhsG6#-F_am6#-F^w6$<#/&Fc\lFf_m&FboFf_m<#F\\q?(FjqF^rF^rF]hpFh]
l@%Fdip>Fe_mFhip@%Fjip>Fe_mF\jp>Fe_mF^jp-F_amF_blF$F$F$6#F\p%(my_multGt+6$F*F+6
)F9F:%"lGF,F5F6F`tF$F$C)>F[rFjn>FjqFcfp>Faq-Fdfp6#FS@$0-F[oF\^qFjq-Fen6#%Cwrong
~dimensions~in~matrix~productG>Fer-Fa]o6$F[rFaq?(FinF^rF^rF[rFh]l?(FboF^rF^rFaq
Fh]l>&Fer6$FinFbo-Fhs6$7#-Fdq6$*&&FI6$FinFc\lF^r&FS6$Fc\lFboF^r/Fc\l;F^rFjqFdtF
^blF$F$F$6#F\p%)log_freeGt+60%%DofYGFcdmF^[lFedm%'irootsGF_[lF7Fc[l%%invGGFigoF
`tF2%$cisG%&condsG63F6F_[o%&deltaG%%cis1GFd[lFfcoFadmF[^m%#bbG%(iroots1G%%sol1G
%'conds1GF`dmF9%$lcbG%%lcb1GFicmF$F$C,-Fhs6$-Fgx6#-F_am6#Fj`pF\io>FinFjn>Ffo-Fg
x6%FgaoFSFV>F^pFcy>F`xFbx>Fer90>F]]l91>Fc\l-F]el6#-F_am6#-Fgx6&%$valG-Fhs6$FfoF
[ioFVF[]l?(F$F^rF^rF$Fh]lC'@%0Fc\lFham>Fdbn-Fgx6&Fc^lFfoFVFc\l>Fdbn-Fjfp6$FinF[
]l@$31FijmFc\l/F^pFb`mFggm>F[r-F^_m6#F^p@%32F[]lF[r2&F^pF\\lFc\lC+-Ffal6'F[gl.%
'isoldeG%*exponent:GFicq%5computing~eigenspaceG>Fbo-&F"6#%,my_linsolveG6&-F[y6$
/FjyFicq-Fdem6#FjwF\cq9/Fer-Ffal6&F[glF]dq%-dimension~isG,&-F^_m6#F_gnF^rF]_mF
\s>FjqFefn>FerF_gn>Ffo-Fgx6&&F"6#%)lc_evalaG-Fceo6#,&FfoF^r-F`^n6$-F[y6$Fidq-Fd
emF`oFjqF\sFVFhw-Ffal6&F^rF]dq%>computing~monomial,~degree~=~GFicq>F`x-Faeo6$,&
F`xF^r*&FjqF^r)FVFicqF^rF^rFV>F^p-F_\l6$Fa\lF^pC)@%/-Fix6#-%'NormalG6#-F[y6$/Fj
yFc\lFYF[]lC)>Ff]n-&F"6#%+rowechelonG6%-F[y6$FigqFjdq.84.Faq>Fgcn-&F"6#%'evalam
G6$-F`^n6$FdhqFdbnFV>F_s<#-Fdq6$&Fgcn6#81/Fdiq;,&F`bmF^rF^rF^r-F[oFicp-Ffal6(F[
glF]dqF_dqFc\l%-conditions:~GF_s>Faq-F^w6$F_s-F[`p6$-F\]m6$<#F`]q-Fe`p6$F_sFfcl
-Fhs6$-Fgx6$f*Fa_pF$Fc_pF$-F_am6#-Fe`pF`oF$F$F$FIF\io@%/FaqFdxC$-Ffal6%F[glF]dq
%,no~solutionG-F\hl6#7&%%FAILGFfoFerFiaqC)-Ffal6&F[glF]dq%(found:~GFaq>Ffo-F[y6
$F^am-FdemFbal>Fer-Fhs6$-Fe`p6#-F[y6$F^amFerFban>F`x-F[y6$F^am-FdemF]^l>F]]l-F
\]m6$F]]lFaq>Fgcn-F[y6$F^am-Fdem6#Fgcn>Fbo-Fcdq6&Ff]nFgcnF\eqFer>F^pF]gqC%>Ff]n
Fahq>FgcnFdbn>Fbo7$-Fgx6$Fix-Fgr6#-F`^n6$-F[y6$Figq-FdemFhemFgcnFer-Ffal6&F[glF
]dq%4new~indeterminates:GF_gn>FjqFefn>FerF_gn>Ffo-Fgx6'Fheq-Fgr6#,&FfoF^r-F`^n6
$-F[y6$FigqFbfqFjqF\sFVFhw,&Fc\lF^rF^rF^r-Ffal6&F^rF]dqFefqFc\l>F`x-Fgr6#,&F`xF
^r*&FjqF^r)FVFc\lF^rF^r>Fc\l-F]el6#-F_am6#-Fgx6&&F"6#FabqFbbqFVF[]l7&-FgrF]^lFf
oFerF]]lF$F$F$6#FgboF_]nt+6#F66#F9F$F$@%/FIF[]lFdx-&F\o6#.Fhbn6#-Fdq6$F^r/Fin;F
^rFIF$F$F$6#F^r%+MsylvesterGt+6)F*F+F,Fb_pF[clF2F36,F5F6F8%$CC1GF9F:%#X1G%$xisG
F<F?F$F$C*@$4FEFZ>FinFjn>FboF_^q>Faq-FhoFj^q>FerFebr@%5F]\o32FgflFbp4-FG6$&FjpF
jflF]qC(>F_s<#-Fdq6$-Fdq6$&FaqFiqF\rFbu@%/FbpF\q>Fc\l-Fgr6#,&FIF^r*&FjyF^r-F]]n
6#FinF^rF\s>Fc\lFI>Ffo<#-Fdq6$-Fdq6$/,&-Fhs6$7#-Fdq6$*&&Fc\lF_tF^r&FaqFbtF^rFct
FdtF^r-Fhs6$7#-Fdq6$*&&FaqF_tF^r&FSFbtF^rF_uFdtF\s&FerFiqF\rFbu>%#EqGFfo-F[w6#-
F^w6$FfoF_s@$1FcblFbpC$>FYFaq>FcyFerC$>FaqFY>FerFcy@%Fhcr>F^p/FjyFhw>F^pFdx-Fgx
6$Fix-F[y6%F^p-Fdq6$-Fdq6$/Fier&FVFiqF\rFbu-FdemF`amF$6#F[frF$6#F\qFgfot+6&F][l
F^[lF9F:6&F6F*FfcoFdepF$F$C&>FerFdhl>Fin-Fe\l6#&FerFagl@'52FYF[]l2FYFVC$>FinFd
\l-Fgr6#-Fho6%FinFinF[]l0&FerF\\lF^co-Fen6#%=argument~has~inproper~formatG@$2&F
erF^glFYC'-Ffal6&Fcbl.%*mat_toolsG%Acomputing~more~terms~up~to~orderG,&FYF^r&Fe
rFaflF\s>Fbo-F_]l6&FahrFS,&F[jrF\sFcirF^r,(FYF^rF[jrF\sF\sF^r>Fc\l-Faeo6#*&)FS,
&F[jrF^rF^rF^rF^r-F`^n6$Fbo&FerFfflF^r>&FerFjfl-Ffel6$,&F[[sF^rFc\lF^rFS>FcirFY
7$-Fgx6'&F"6#%)truncateGFijrFSFVFY-Fgx6'Fc[sF[[sFSFVFYF$Fh[mF$6#F[gl%*lc_evalam
Gt+6$Fh[lF^[lF$F$F$-Fgx6&Fheq-FgrF`oFS-Fdq6$&FjpF[ar/F9;FgflFbpF$F$F$6#F\pFceot
+FahlF$F$F$-Fgx6$Ff^lF^\sF$F$F$6#F^r%*mat_checkGt+F[\s6-F*F+F,F-FfcoF]hoFh]mF]d
mFicmF6F5F$F$@2/&FIF\\lFdjlC(>Fin-Fgr6#&FIFagl-Ffal6%Fcbl.%&debugG%9enter~check
~localizationG@$2,&&FIFjflF^r&FIFaflF\sF[]lC$-Ffal6%FcblFf]s%Cno~terms~computed
,~localization~OKG-F\hl6#-Fgr6#*&FinF^r)FSF]^sF^r-Ffal6&FcblFf]s%Mtesting~serie
s~expansion,~computed~terms~areG;F]^sF\^s>F[r-Fgx6$Fix-Fgr6#,&&FIFfflF^r-Fgx6&t
+6%Ff]mF^[lF`tF$F$F$-Fhs6$-F_eo6%FI/FSF[]l,&FVF^rF^rF^rFieoF$F$F$6#FgflFc]sFSF[
^sF\s@%4-&F\o6#.F]dpFjglC$>%#EEGFI-Fen6#%0in~localizationGC$-Ffal6%FcblFf]s%0lo
calization~OKGFd^s/F^]sFghlC--Ffal6%FcblFf]s%;enter~check~transformationG-Ffal6
&FcblFf]sFj^s;F]^s&FIF[q-Ffal6&FcblFf]s%(span~isGFc]s>Fin-&F"6#Fi\s6$&FIFf`oFS@
$4-FG6$F]^sFj`mC$-Ffal6%FcblFf]s%Jvaluation~not~computed,~transformation~OKG-F
\hl6#-&F"6#%*transformG6%FinFSF\^s>Faq-FddlF^bs-Ffal6&FcblFf]s%;testing~transfo
rmed~terms:G;Faq,&FfasF^rFc]sF^r>Fbo-Fgx6$Fix-Fgr6#*&)FSFaqF^r-Fgx6&t+Ff_sF$F$F
$Fg_sF$F$F$6#Fgfl-Fgx6$Ff^l-Fgr6#*&FinF^rF\dsF\sFS,(FfasF^rFc]sF^rFaqF\sF^r>Fbo
-F[cs6%FboFSF\^s>F[r-Fgx6$-%"@G6$FixFf^l-Fgr6#,&Fb_sF^rFboF\s@%F_`sC%>Ff`sFI-Fh
jp6#Ff`s-Fen6#%2in~transformationGC$-Ffal6%FcblFf]s%2transformation~OKGFjbs/F^]
sF_ilC*-Ffal6%FcblFf]s%=enter~check~exponential~partG@$5Fj]sFabsC%>Fin-Fgr6#,&&
FIF^glF\s-F\bs6$FfasFSF^r-Ffal6%FcblFf]s%Gno~terms~computed,~exponential~part~O
KG-F\hl6#-FgrF_dr>Faq-FddlF]gs>FinF\gs>Fbo-Fgx6$Ff^l-Fgr6#*&,&FinF^rF[gsF\sF^rF
g^sF\s-Ffal6&FcblFf]s%Dtesting~exponential~part,~terms~areGF[_s>F[r-Fgx6$Fix-Fg
r6#,&*&Fb_sF^rFg^sF\sF^r-Fgx6&t+Ff_sF$F$F$Fg_sF$F$F$6#FgflFboFSF[^sF\s@%F_`sC$>
Ff`sFI-Fen6#%4in~change~exp.~partGC$-Ffal6%FcblFf]s%4exponential~part~OKG-F\hl6
#-Fgr6#F]hs/F^]sFdilC$-Ffal6%FcblFf]s%1enter~check~copyG-F\hl6#-F\bs6$F]^sFS/F^
]sFiilC%-Ffal6%FcblFf]s%9enter~check~substitutionG>FinF\gs-F\hl6#-F[y6$F[gs-Fde
mF_dr/F^]sF\jlC+-Ffal6%FcblFf]s%9enter~check~ramificationG@$FefsC&-Ffal6%FcblFf
]s%Cno~terms~computed,~ramification~OKG>FinF\gs>Fbo-Ffel6$**F[gsF^rFc]sF^r)FS,&
Fc]sF^rF\sF^rF^r-F[y6$/FS*&F[gsF^r)FSFc]sF^rF][tF^rFS-F\hl6#-FgrFf^p>FaqFegs>Fi
nF\gs>FboFj[t>Fc\l-Faeo6#*&FboF^rFg^sF\s-Ffal6&FcblFf]s%@testing~ramification,~
terms~areGF[_s>F[r-Fgx6$Fix-Fgr6#,&FghsF^r-Fgx6&t+Ff_sF$F$F$Fg_sF$F$F$6#FgflFc
\lFSF[^sF\s@%F_`sC$>Ff`sFI-Fen6#%0in~ramificationGC$-Ffal6%FcblFf]s%0ramificati
on~OKGFd\t/F^]sF^coC,-Ffal6%FcblFf]s%<enter~check~simple~operatorG>Faq-Fddl6$Fc
]sFS-Ffal6&FcblFf]s%8getting~terms,~order~isG,(F]^sF\sF\sF^rF[gsF^r>Fin-F_]l6&F
c]sFSFaqF`_t>Fer-Fgr6#Fb_s-F\bsF\_t>FerFe_t>Fjq-Fgx6'Fc[s-Ffel6$*&)FS,&F]^sF^rF
^rF^rF^r,&-F`^n6$FinFerF^r-Fgx6%Fg^nFerFSF\sF^rFSFSF[]lF[gs>F[r-Fgx6$Fix-Fgx6'F
c[s-Faeo6#,&FjqF^rF\^sF\sFSF[]lF[gs@%F_`sC$>Ff`sFI-Fen6#%3in~simple~operatorGC$
-Ffal6%FcblFf]s%3simple~operator~OKGFh]l/F^]sFijlC$-Ffal6%FcblFf]s%Center~check
~block-reduced~operatorG@%/-F^_mFb]sF^rC--Ffal6%FcblFf]s%:classical~block-reduc
tionG>Fbo-F\bs6$&FIFfboFS>F_s&Fc]sF\\l-Ffal6&FcblFf]s%/getting~terms:G;,$F_sF\s
F\^s>Fin-F_]l6&FgbtFSF^ctF\^s>F^p-Fe\l6#Fgbt>F`xF[gs>Ffo-&F\o6#.%)copyintoG6&&F
I6#Ffjm-F]]n6#FcctF^r,(F^pF^rF`xF\sF^rF^r>Fin-Fgx6%Fgao-Fgx6$Fix-F]cs6%FinFSFfo
FS>F[r-Fgx6'Fc[s-Faeo6#,&-F\z6%Fin;FadtF^pFaetF^rFb_sF\sFSF^ctF\^s@'F_`sC$>Ff`s
FI-Fen6#%Qin~block--reduced~operator:~wrong~reduced~systemG42F\^s-Fbdo6%-F\z6%F
in;F^r,&F^pF^rF`xF\sFaetFSF[]lC$>Ff`sFI-Fen6#%Rin~block--reduced~operator:~orde
r~not~high~enoughGC%-Ffal6&FcblFf]s%)order~isGF\^s-Ffal6%FcblFf]s%;block--reduc
ed~operator~OKG-F\z6%-F[cs6%FboFSFfoFaetFaetC--Ffal6%FcblFf]s%<generalized~bloc
k-reductionG>Fbo-F\bs6$&FIF[coFS>F_sFibtFjbt>Fin-F_]l6&FggtFSF^ctF\^s>F^p-Fe\l6
#Fggt>F`xF[gs>Ffo-Fhct6&&FI6#"#<-F]]nFdcqF^rFadt>Fin-Fgx6%Fgao-Fgx6$Fix-F[csFhd
tFS>F[rFjdt@%F_`sC$>Ff`sFIFeetC%FfftFiftF\gtF$6$FiglFf`sF$6#F\pFa]lt+Fjgr60Fh[l
F7F6F,%#ChGFg[lFf]q%#iiGF]\nFb_pF3Fh]m%#nnG%&clistGF$F$C&@$Fa_o-F\hl6#-&F"6#%3m
at_eval_algebraicGF^`o>FinFdhl>Fc\lFd\l@:Fehr-Fgr6#Fi\lFfhlC'-Ffal6&FcblFgir%2e
valuating~terms:G;FVFY-Ffal6&FcblFgir%-lookahead~isGF`gl-Ffal6&FcblFgir%2comput
ed~order~isGFcgl@$2FcglFYC)>Fer-F_]l6&Fe`oFS,(FcglF^rF`glF^rF^rF^r,&F`glF^rFYF^
r-Ffal6%FcblFgir%<updating~transformed~matrixG>Fer-&F"6#Fb]q6$F]gl-F_\u6$FerFif
l-Ffal6%FcblFgir%;reducing~algebraic~numbersG>Fer-Fgx6&FheqFerFSFe`o>Fefl-Ffel6
$,&FeflF^rFerF^rFS>FcglFY-Fgx6'Fc[sFeflFSFVFYF[ilC'FhjtF\[uF_[u@$Fc[uC)>FerFf[u
Fj[u>FerF^\uFd\u>FerFh\u>FeflF[]u>FcglFY-Fgx6'Fc[s-&F"6#%.matrix_seriesG6%FeflF
SFYFSFVFYFcjlC%>FboFcgl@$2FiflFYC&-Ffal6%FcblFgir%Hcomputing~new~terms~of~serie
s~expansionG>Fefl-Fgx6't+6&Ff]mF^[lF`tF7F$F$F$-Fheq6%-Fhs6$-F_eo6%FIFSF\`sFieoF
SFYF$F$F$6#F[glF`glFS,&FYF^rF`flF\sFbo>FiflFY-Ffal6&FcblFgirFhftFY-Ffel6$*&)FSF
`flF^r-Fgx6'Fc[sFeflFS,&FVF^rF`flF\sFd_uF^rFSF^ilC$@$Fc^uC&-Ffal6%FcblFgir%Kupd
ating~terms~after~change~of~exponentialGFd\u>Fefl-Fgx6&Fheq-Fgr6#,&FeflF^r-F_]l
6&FcglFS,&FiflF^rF^rF^rFYF^rFSF`gl>FiflFYF_]uFcil-F_]l6&F`flFSFVFYF[jlC$@$Fc^uC
&-Ffal6&FcblFgir%?updating~ramified~terms.~Need:G;F]auFY-Ffal6&FcblFgir%&take:G
;-%%iquoG6$,(FiflF^rF\pF^rF`glF\sF`gl-F]bu6$,(FYF^rF`glF\sF^rF^rF`gl>Fefl-Ffel6
$,&FeflF^r**F]glF^rF`glF^r)FS,&F`glF^rF\sF^rF^r-F[y6$F]bo-F_]l6&FcglFSF\buF`buF
^rF^rFS>Fifl,(*&,&F^rF^rF`buF^rF^rF`glF^rF^rF`glF^r!"#F^rF_]uFhilC$@$Fc^uC%-Ffa
l6%FcblFgir%Bupdating~terms~after~substitutionG>Fefl-Fgr6#,&FeflF^r-Fgx6&Fgao-F
[y6$F]glF[auFSFixF^r>FiflFYF_]uFhjlC$-&F"6#%3eval_block_reducedG6%FinFSFYF_]uF[
[mC'-Ffal6&Fcbl.%8mat_direct_block_reduceGFjjtF[[u-Ffal6&FcblF[euFa[uFifl>Ff]n-
Fe\l6#&FinF^dt@$Fc^uC'>F]]l-Fddl6$FbeuFS-F_]l6&&Fin6#FcrFSFYFY?(F`xF]auF^rFYFh]
lC'-Ffal6&FgflF[eu%Dlifting~valuation~of~first~block~toGF`x>Ffo-F\z6%-&F"6#%*ma
t_coeffG6%FjeuFSF`x;F^rF]gl;,(Ff]nF^rFcglF\sF^rF^rFf]n>F^p-Fgr6#-F`^n6$&Fin6#Fe
wFfo?(F[rF^rF^r-FdfpFdcqFh]l?(F_sF^rF^r-F[oFdcqFh]l@$0&F^p6$F_sF[rF[]lC$-Ffal6(
FcblF[eu%;elementary~transformation:G,*Ff]nF^rF`glF\sF]glF\sF_sF^r,(Ff]nF^rFcgl
F\sF[rF^r,$*&FiguF^r)FS,&F`xF^rF]]lF\sF^rF\s>Fjeu-&F"6#%>mat_elementary_transfo
rmationG6'FjeuFSF_huF`huFahu@$0Fe`oFdxC%-Ffal6&FgflF[eu%Elifting~valuation~of~s
econd~block~toG,&F`xF^rF^rF^r>Fdiq7$-Fdq6$Faq/Faq;F^r,(Ff]nF^rF`glF\sF]glF\s-Fd
q6$Faq/Faq;,(Ff]nF^rF`glF\sF^rF^r,&Ff]nF^rFcglF\s?(F[rF^rF^rFcglFh]lC%>Ffo-F\z6
%-Fefu6%FjeuFSFaiu;,&F]glF^rF^rF^rF^ju7#,(Ff]nF^rF[rF\sF^rF^r>F^p-Fgr6#-F`^n6$F
e`oFfo?(F_sF^rF^rFfguFh]l@$0&F^p6$F_sF^rF[]lC$-Ffal6(FcblF[euF^hu&FdiqFj_mFiju,
$*&Fb[vF^rFchuF^rF\s>Fjeu-Fghu6'FjeuFSFg[vFijuFh[v>Fefl-F\z6%-F_]l6&FjeuFSFVFYF
jfuFjfu>FiflFYF_]uFb[mC&Fhjt-Ffal6&FcblFgirFa[uFifl@$Fc^uC'>FerF[au-Ffal6%FcblF
gir%Gupdating~elementary~transformed~matrixG?(FaqF^rF^r-F^_m6#F]glFh]lC&>F_s&&F
]glF`amF\\l>F[r&Fb]vFafl>Fer-&F\o6#%'addcolG6&FerF_sF[r&F`glFjgu>Fer-&F\o6#%'ad
drowG6&FerF[rF_s,$F[^vF\s>FeflF[]u>FiflFYF_]uF^[mC&FhjtFd\v@$Fc^uC(>Fer-F_]l6&F
`glFSF]auFY-Ffal6%FcblFgir%8updating~swapped~matrixG>Fer-F]hn6%Fer&F]glF\\l&F]g
lFafl>Fer-FchnF`_v>FeflF[]u>FiflFYF_]uF$Fh[mF$6#F[gl%/int_diff_splitGt+6$%#fsGF
b_p6+%%repsG%$fs1GF_[oFf]m%#f1GF6%#n1GF[inF5F$F$C'>FinFe\m>Fc\l-Fhs6$FIF[io>Fbo
Fc\l?(F$F^rF^rF$0FboFe\mC'>FerFefn>FfoF^r?&FjqFboFh]lC%>F[r-Fc_m6$F]irFS>Faq-Fc
_m6$Fe\pFS@$/F[rFaqC$>F_s-Fix6#-Fegq6#*(,&-Fc^l6%F]irFS,&F[rF^rF\sF^rF^r-Fc^l6%
Fe\pFSF`bvF\sF^r-Fc^l6%Fe\pFSF[rF\sF[rF\s@$3-FG6$F_sFj`m/-Fix6#-Fegq6#,&F]irF^r
-F[y6$/FS,&FSF^rF_sF^rFe\pF\sF[]lC$>Fc\l-F[`p6$Fc\l<#Fjq>Ffo*&FfoF^r)Fe\pF[\pF^
r>Fin-F\]m6$Fin<#7$Ffo-Fc_m6$FfoFS>FboFc\lFinF$F$F$6#F\pFbelt+F^_oFahlF$F$C%>Fi
nFdhl@$/F`flFbfl-FddlFgeo@.FfhlFcglFb[mFiflF[ilFcgl-%'memberG6$F_dl<,.%(product
G.%$sumGF_[mFijlF`dlFdjlF\jlFiilF_ilF\[mFiflFcil-F`elFeaoF]co-Fen6#%4not~yet~im
plementedGF$Fh[mF$6#F\p%'MapplyGt+6%F*F^[lFbdmF$F$F$-Fgx6$Ff^l-Fgr6#,&-F`^nFgdl
F^r-Fgx6%Fg^nFVFSF\sF$F$F$6#Fgfl%/reduceColRank4Gt+6/Fh[lF^[lFfcmF6F`tF?FgcmFd`
nFhcmFedmFfdmFicmFjcm61F\dmF]ho%#T1G%&invT1G%#TTG%&invTTGFf]qF9F:F^dmF_[o%&vals
1GFb_pF]\nF_\nF$F$C+-Ffal6%FcblF[em%5reducing~column~rankG>Fc\lF\]n>FerF\]n>Fin
F^\s@$2F^r-Fi[o6#Fbx>%#PPG-FgrFggv>F^p,&F[sF^r&Fhw6#FcyF\s?(FfoF^rF^rFfgvFh]lC0
>Ff]n-Fdhm6%Fbx;F^rF[sFfo@$/&Ff]n6#,(F[sF^rFfoF\sF^rF^rF[]lC+-Ffal6%FcblF[em%:p
ermute~solution~vector~XG>F_sFihv?(F$F^rF^rF$/&Ff]nFj_mF[]l>F_sF]hp>FghvFaiv>Fa
ivF[]l>Fbx-F[jm6%FbxF_sFihv>Fin-Fajm6%FinF_sFihv>Fin-F[jmFjiv@$/Fbp"#8C$>F\eq-F
ajm6%F\eqF_sFihv>Fgaq-F[jm6%FgaqF_sFihv@%Fe[nC$-Fj[n6'FghvFVFS.Fdiq.Fdbn>Ff]n-F
gx6&Fc\n-Fgr6#*&FdiqF^rFf]nF^rFVFS>Ff]n-Ffel6#*&FghvF\sFf]nF^r>Fbo-&F\o6#F[dt6&
-Fhs6$Ff]nFK-F]]nFhemF^rFihv>Fbo-Fh]nFf^p?(F_sF^rF^rF^pFh]l>Faiv*&FaivF^r)FV,*F
cyF^r&FjwF\\lF^r&FjwFj_mF\sF\sF^rF^r>FjqF\]n>Fjq-Fi[w6&F\\wF\]nF^rFihv>F[r-Fh]n
Ff_m@$F^jvC$>F\eq-Fgr6#-F`^n6$F\eqFjq>Fgaq-Fgr6#-F`^n6$F[rFgaq>Fc\l-Fgr6#-F`^n6
$Fc\lFjq>Fer-Fgr6#-F`^n6$F[rFer>Fin-Fgr6#-F`^n6%F[rFinFjq>Fbx-Fgr6#-F`^n6$FboFb
x>Fin-Fgr6#,&FinF^r-F`^n6$Fer-Fgx6%Fg^nFc\lFSF\sFinF$6$%$TTTGFigvF$6#F_jv%-penc
il_solveGt+6'F*F+F,F2%'invABcG6&Fd[lFb_pFbdmFicmF$F$C'@%3Fb\m4-FG6$FYFfcl>FinFY
C$>FinF[]l?(F$F^rF^rF$/-Fix6#-Fjfo6#-Fgr6#,&FIF^r*&FinF^rFSF^rF\sF[]l>Fin,&FinF
^rF^rF^r@%3/FbpFcbl-FG6$FcyFK>FerFcy>Fer-F]doF\aw@$33Feaw-FG6$FcyFfclFb`wC$>FYF
in>FcyF^bl>Fbo-Fgr6#-F`^n6$FerFV7$-FdemFf^p-&F\o6#.%*scalarmulG6$FboFinF$F$F$6#
Fcbl%+super_formGt+6(F*FcdmF^[lFdep%%invSGF_[l6(Ff]qF+F\^mF9%"LGFhgoF$F$C(>Fc\l
F[]l>Fin-&F"6#F\em6(FIFVF\dnFjyFYFcy>Fjq7$,&&F`flF\\lF^rF\sF^r-F_am6#&F`fl6#;F
\p-F^_mFgil@$0FiflF[]l>Fc\l,$-&F"6#Fijn6%FcyFVF\dnF\s@$2F[]lFc\lC$>Fcy-Fgr6#*&F
\`rF^rFcyF^r>FY-Fgr6#*&)FV,$Fc\lF\sF^rFYF^r7&-Fgr6#,&*&FVF^rF_dlF^rF^rFc\lF^r-F
gr6#*&-F`^n6$FcyFSF^rFVF^rFjq-F[y6$/Fjy,&FjyF^rFc\lF\sFeflF$F$F$6#Fcp%1mat_bloc
k_reduceGt+6'F][lF^[lF`[lFedmF_[l61F/F5F6F]dmF`tF_[o%#B0G%#D0G%#N0G%"HGFdepFicm
F+F2%#KKGF$F$C'@$Fa_o-F\hl6#-&F"6#%;mat_block_reduce_algebraicGF^`o@%/-F^_mFdoF
^rC9>Fin-F^q6#;F^r"#:>F_dlFijl>Fer&FVF\\l>Fifl,$FerF\s@%/-%(ldegreeG6$FYFcy-Fc_
mF_iw>F`flFbfl>F`flFjhw>Fc\lFd\l>FboF`iw>Ff]n-%)mat_copyGF`o>Faq-Fgx6&Fc^l-F_]l
6&Ff]nFSFjhwFjhwFSFjhw-Ffal6%F[glFgir%Ocomputing~block-reduced~form~of~leading~
matrixG>F^p-&F"6#%4pencil_block_reduceG6(Faq-F]]nFacmFYFcyFb`nFahm>Faq-Fgx6&Fc^
lF^pFcyF[]l>F]glFbo>F`gl7#Fer>Fcgl-Fgr6#*&-F\z6%Faq;F^r,&Fc\lF^rFboF\sFc[xF^r)F
SFerF\s>Fagu-F\z6%FaqFc[x;,(Fc\lF^rFboF\sF^rF^rFc\l>Fe`oFa[x>Fjeu-F\z6%FaqFi[xF
i[x>Fefl-Fgr6#*&FjeuF^rFe[xF\s>Fbeu-Fho6%Fd[xFboF[]l>Ff]n-&F"6#%.mat_transformG
6'Ff]nFSF]]l-Fgr6#*$F]]lF\sF[]l>FeboFf]n@%1FerF^r>&Fin6#FehwFdiq>Fc]xFdxCB>Fin-
F^q6#;F^r"#F>F_dlFijl>F`flFbfl>FerFhhw>Fjq,&FerF^rF_hwF\s@'/F`iwF[]lC$>F`flFjhw
>Fifl,&FerF\sF\sF^rF\iwC$>F`flFbfl>Fifl,&FjqF\sF\sF^rC$>F`flFj^x>FiflFj^x>Fc\lF
d\l>F`glFV>Ff]nFfiw>&Fin6#"#@-&F"6#%.mat_to_simpleG6%FIFSFV>F[r-Fgfo6&Fb_xFSF[]
lF[]l>F_s&F[rF\\l>Ffo&F[rFafl-Ffal6%F[glFgir%Ocomputing~block-reduced~form~of~l
eading~pencilG>F^p-Fbjw6(FfoF_sFYFcyFb`nFahm>FfoFhjw>F_s-Fgr6#,$-Fgx6&Fc^lF^pFc
yF^rF\s>FboF`iw@$/FboF[]l>Fbo,&Fc\lF^r-Fc_m6$-Fjfo6#-Fgr6#,&FfoF^r*&FcyF^rF_sF^
rF\sFcyF\s>F]glFbo>&Fin6#F_jv-F\z6%F_sFc[xFc[x>Febo-F\z6%F_sFi[xFi[x>Fc]x-F\z6%
FfoFc[xFc[x>&Fin6#"#;-F\z6%FfoFi[xFi[x>&FinFdhtFd\x>&Fin6#"#=Fd\x@$/F`flFifl>Fe
fl-Fgr6#*&FhbxF^r)FSF]`mF\s>Ff]n-Fh\x6'Ff]nFS-Fgr6#-F`^n6$&&Figl6#Fb_xFfflF]]l-
F]do6#F]dx,(FerF^rFjqF\sF\sF^r>&Fin6#"#AF`x>&Fin6#"#BF]]l>FjboFf]n@%/FjqF[]l>&F
in6#F[^xFdiq>FcexFdx>FdbnFegl>&Figl6#FdbnFinFdbnF$Fh[mF$6#Fcbl%*MpolsoldeGt+F[
\s69FicmFjcmF_[lFfcwF]hoF6Figo%#smG%#pmGFhgoFjgoFfcoF]\nF`tFfcm%#ldG%#bdGF[inF9
%#LTGFd[lFadm%$sm1GF$F$C2>F2.F2>F^pFaio>F[rFd[o@%F]\o>Ff]n&FjpFffl>Ff]n-F]gp6$F
[rF[]l-Ffal6%Fcbl.F[fx%5calling~good_form~:~G>Fer-&F"6#%*good_formG6(FIFf]nFSFi
nFboFc\l>F`x-&F\o6#F]go6#F[jr-Ffal6&FcblFagx%?indicial~equation:~at~infinityGF`
x>F_sFijr>Fjq-&F"6#%-NNIexponentsG6$F`xFc\l>F]]l-F^_mFf_m@%3/F]]lF[]l-F[dpFj_m-
F\hl6#F]gxC(>FaqF]gx>FdiqF[]l>Fgcn-&F"6#%+vectdegreeG6%F_sFSF[r>85&FgcnF\\l>Fdh
q&FgcnFafl?(F$F^rF^rF$51F[]lFjix2F[]lF]]l@%3Fajx2Fjix&FjqF[^nC*>88-&F"6#%+Mylin
solveG6(-F[y6$/Fc\lFejxF[jrF]gxF2FdiqF[rF\[w>Faq-Fgr6#,&FaqF^r*&FhjxF^r)FSFejxF
^rF^r>F_s-Fgx6$Ff^l-Fgr6#,&F_sF^r*&Ff[yF^r-F`^n6$-F[y6$F`[yF]irFhjxF^rF\s>FgcnF
dix>FjixF[jx>FdhqF]jx>F]]l,&F]]lF^rF\sF^r>FdiqFdbnC&@$3Fajx/FejxFjix>F]]lFf\y>8
7-&F"6#%,TrigonalizeG6%-F[y6$/Fc\lFjixF[jrFdhqF[r>Fhjx-Fjjx6(&F^]yF\\l&F^]yFafl
F2FdiqF[rF\[w@%0FhjxFdxC(>Faq-Fgr6#,&FaqF^r*&)FSFjixF^rFhjxF^rF^r>F_s-Fgx6$Ff^l
-Fgr6#,&F_sF^r*&Fd^yF^r-F`^n6$-F[y6$Ff]yF]irFhjxF^rF\s>FgcnFdix>FjixF[jx>FdhqF]
jx>FdiqFdbnC$>89-%&solveG6#<#-Fdq6$&F[^y6#86/F_`y;,&-&F\o6#%%rankG6#Fj]yF^rF^rF
^rF[r@%/Ff_yFdx>FjixF\sC*-F[w6#Ff_y>Fhjx-Fjjx6(Fj]y-F[y6$Ff_y-Fdem6#F[^yF2FdiqF
[rF\[w>F_s-Fgx6$Ff^l-Fgr6#,&-F[y6$Ff_y-FdemFj_mF^rF[_yF\s>Faq-Fgr6#,&-F[y6$Ff_y
FfgrF^rFc^yF^r>FgcnFdix>FjixF[jx>FdhqF]jx>FdiqFdbn>F_s-Fgx6$Fh^m-Fgx6$Fb[pF_s>%
)SM_recupG-FgrFj_m@$-&F"6#%(nonzeroGFjguC%>8:-Fhs6$-Fgx6%Fe^mF_sFSF\io>Ffo-Fh_y
6$Fgcy-%'selectG6%Fg_p-Fe`p6#FgcyF2@%/FfoFdx-F\hl6#Fe\m-F[wFbal-Fgr6#-F`^n6$Fin
FaqF$6$F2F^cyF$6#F\pFjeqt+6&Ff]mF^[lF7Fcdm6'%"gGF[in%#ddGF^[oF9F$F$C(@$32F\pFbp
/FVFe\mOFI>Fin-Fgao6%Fg^mFSFf^l@$55FjipFheyF]ar@%2FbpF[gl-F\hlF_dr-F\hl6#-Fc[s6
&FinFSFY-Fc_mFgeo@%/FbpF[gl>FboFY>Fbo-F^iwFgeo>Fc\lFgfy?(F$F^rF^rF$Fh]lC%>Fer-F
ix6#-%'ExpandG6#-Fc^l6%FinFSFbo@$550FerF[]l2Fc\lFbo3Fify1FYFbo-F\hl6#,&*&FerF^r
)FSFboF^rF^r-Fhs6$7#-Fdq6$*&-Fc^l6%FinFSFjqF^r)FSFjqF^r/Fjq;FhzFgfyFdtF^r>FboFh
zF$F$F$6#F[gl%,tensor_prodGt+Fd]q6)F9F:F,%#maG%#mbG%#naG%#nbGF$F$C)>FerFjn>FjqF
_^q>F[rFcfp>FaqF[^q>Fc\l-Fho6%*&FerF^rFjqF^r*&F[rF^rFaqF^rF[]l?(FinF[]lF^r,&Fer
F^rF\sF^rFh]l?(FboF[]lF^rF`bvFh]l>Fc\l-Fhct6&-F[cw6$FS&FI6$FbawFhzFc\l,&F^rF^r*
&FinF^rFjqF^rF^r,&F^rF^r*&FboF^rFaqF^rF^r-FgrFacmF$F$F$6#F\pFhirt+F$F$F$F$F$F$F
$F$6#F[]l%*G_matrix0Gt+6&F][lF^[lF`tF`[l6+F?Fb[lF6Fc[lF9F:Fg[lF3Fh[lF$F$C/>Fin&
FYF\\l>Fbo-F_\l6$Fa\lFY>Fc\lFd\l>FerFi\l>Ffo-F_]l6&FIFSFc]lFd_l>Ffo-Fgr6#-FgbmF
bal>FaqF^r?(F_sF^rF^rFg]lFh]lC$?(F[rFaqF^r,(FaqF^r&FboFj_mF^rF\sF^rFh]l?(FjqF^r
F^rFc\lFh]l>Fier-Fc^l6%FhqFS,(FinF\sF_sF^rF\sF^r>FaqF[r?(F[rFaqF^rFc\lFh]l?(Fjq
F^rF^rFc\lFh]l>Fier-Fc^l6%FhqFSFd_l>F_s7#-Fdq6$-Fdq6$,(FinF\sF[rF^rF\sF^r/Fjq;F
^r&FboFjgl/F[rFdalFeal-Ffal6&F^rFhalF]blF_sF^blF$F`blF$6#F[gl%'mat_lcGt+F[\s6#F
h]mF$F$C$>FinFgfn-Fefu6%FIFSFinF$F$F$6#F\pF\[yt+6(F*FcdmF2F9F6Ffcm6&FadmFh]mF:F
?F$F$C%>Fer-&F\o6#%)linsolveG6&FIFS.FinFbo>FjyFY@%30FerFdx2FinFcyC%>Fer-F[y6$-F
dq6$/&FboFacm&FV6#,&FYF^rFc\lF^r/Fc\l;F^r,&FcyF^rFinF\s-FdemF_bl>Fjy,(FYF^rFcyF
^rFinF\s-F\hl6#F\azF_azF$F$F$6#Fcp%*MratsoldeGt+6$F*F^[lF$F$F$@+Fify@%/F^cr%'cy
clicG-F\hl6#-%.cyclic_ratsolG6%FIFSF[gx-F\hl6#-%.direct_ratsolGF]bz3F]\o/F[gxFh
az-F\hl6#-F\bzFhfn3F]\o/F[gx%'directG-F\hl6#-FabzFhfnFjip-F\hl6#-&F"6#FabzFhfnF
^bzF$F$F$6#F\p%,val_ldegreeGt+6$Ff]mF^[lF$F$F$@%F]arFham-F^iwFhfnF$F$F$6#F\p%*r
andomRatGt+6%F^[l%'denomsG%&polesG6(FfcmF[inF9Fadm%$infGF5F$F$C(>FboF^r>FjqF[]l
?(Fc\lF^rF^r-F^_mF\^qFh]lC$>Fin-F[jn6#;F[]l&FVFacm@%/&FSFacmFham>FjqFjdz>Fbo*&F
boF^r)F]ez-FinF$F^r>F[r-%$maxG6$F[]l,(-Fc_m6$FboFIF^rFjqF^rFbcuF^r>Fer-Ffin6&FI
/Fiin-F[jn6#;!"&Fcbl/F`jn-F[jn6#;F[]lF[r/FfjnF\p*&FerF^rFboF\sF$F$F$6#Fgfl%/reg
ular_seriesGt+6+F-FfcoFcdmF^[lFedmFg]mF7F`tF]\m64F6Fc[l%#G0G%&invG0GFjfv%$resGF
3F9F\`qFadmF^`qF_`qF2F[in%#b1GF_[oF_\m%#mdGF$F$C9@%/FbpFg`o>FdhqF[s>FdhqFe\m>Fi
nFjn>Fbo-Fgr6#,&*&FIF^rFjyF^rF^rFSF\s>Fc\l-F[y6$/FYF[]lFibw-Ffal6%F[glF]dq%Jcom
puting~diagonal~form~of~leading~pencilG>FaqFaio>Fer-Fgr6#*$Fc\lF\s>F`xFb`m>F]]l
Fe\m>FdiqF[]l>FjixF[]l>FfoFdx?&F^p-F``m6$FcyFjyFh]l@$-FG6$FicqFj`mC%-Ffal6&F[gl
F]dq%.integer~root:GFicq>Fjix-Feez6$FjixFicq@$52FicqFjw-Fg_p6$Fdhq.%'finiteGC$>
Ffo6$FfoFicq>Fdiq,&FdiqF^r&F^pFaflF^r-Ffal6&F[glF]dq%<maximal~integer~differenc
e:GFjix>Ffo-Fcam6#7#Ffo>F[rFb`m>FdbnF\cq>F_sF[]l?(F$F^rF^rF$2-F^_mF]^lFdiqC+>Fg
cn-&F"6#Fh_q60-Fgr6#*&)FYFjyF^rFboF^r-Fgr6#,&-F`^n6$FIFdbnF\s*&FYF^rFVF^rF^rFYF
cyFfoFjyFhwFc\lFerF[]lFjwFf]nF`xF]]l@$/F[jxFg[rFggm>Fdbn-Fgr6#F[jx>F]]l&FgcnFjf
l-Ffal6&F\pF]dq%9regular~solutions~found:G-F^_m6#&FgcnFffl@$5-Fg_p6$Fdhq.%(logf
reeG4-Fa`sFdoC$>F[r-Fgx6&&F"6#%,log_collectG-Fjhq6$-FgrFiexFYFYFjw-F\hl6#-F[y6$
/Ff]n%#_CG7$-FdemFjglF]]l>F[r-F[y6$-F_amF[^n7$-Fhs6$FdbnFban-F_amFjgl@$/Fd[[lF_
][lFggm>F`xFa][l@$/F[rFb`m-F\hl6#Fb`m>F[r-Fhs6$7#-Fdq6$*(&F[rFj_mF^r-%*factoria
lG6#F]hpF\s)-%$logGF`]nF]hpF^r/F_s;F^r-F^_mFjglFdt>F[r-Fgx6&F^^[l-Fjhq6$-FgrFjg
lFYFYFjwFf^[lF$6$FablFi^[lF$6#Fg`o%+denomatrixGt+Fdaz6%F[inF9F:F$F$C%>FinF^r?(F
boF^rF^rFd[oFh]l?(Fc\lF^rF^rFh[oFh]l>Fin-%$lcmG6$Fin-%&denomG6#&FI6$FboFc\lFinF
$F$F$6#F\pFabzt+Fdaz6<F6F[inF9F:Ff]m%%lodeG%&factsG%#npGF_[lFjgo%&boundG%%minrG
Fh]mFadm%#RFGF+%#PFGF\^mFhgo%#tiGF?%#dAG%#dbGFfcm%#muGFcdmF$F$C:>F_`yFaio-Ffal6
%Fcbl.Fabz%Lcomuting~a~rational~factor~of~the~solutionsG-Ffal6%FcblFfc[l%jncomu
ting~the~denominators~of~the~matrix~system~and~of~the~r.h.sG>Ff_yF^r>FinFd[o?(F
c\lF^rF^rFinFh]l?(FerF^rF^rFinFh]l>Ff_y-F]b[l6$Ff_y-F`b[l6#&FI6$Fc\lFer>Ff_y-Ff
^l6#*&Ff_yF^r-%$gcdG6$Ff_y-Fg^n6$Ff_yFSF\s@%F]\oC%>8=F[gx>FgcyF^r?(Fc\lF^rF^rFi
nFh]l>Fgcy-F]b[l6$Fgcy-F`b[l6#&Fbe[lFacm>FgcyF^r-Ffal6&FcblFfc[l%&time=G,&FaioF
^rF_`yF\s>F^]yFaio-Ffal6%FcblFfc[l%Ofactoring~the~denominator~of~the~matrix~sys
temG>Faq7#-F_am6#&-%(factorsGF]ayFafl-Ffal6&FcblFfc[lF^f[l,&FaioF^rF^]yF\s>Fdbn
F^r?&FjqFaqFh]lC1>F^]yFaio-Ffal6&FcblFfc[l%6calling~myequaind~at:GFe\p>F^p-&F"6
#%*myequaindG6&FIFSFe\pFfo>FjixF[]l@%0F^pF[]l-Ffal6)FcblFfc[l%>indicial~equatio
n~and~height:GF^pFjixF^f[lF]g[lC(>F^]yFaio-Ffal6&FcblFfc[l%Acalling~super-irred
ucibility~at:GFe\p>F^p-F[dwFjg[l>Fjix&F^pFjfl>F^p&F^pFfflF^h[l>F^]yFaio>8<-&F"6
#Fc]m6%F^pFfoFS-Ffal6(FcblFfc[l%=the~minimal~integer~exponantGF^i[lF^f[lF]g[l@$
3/F^i[lFhamFjip-F\hl6#-F]gpF]cq>F^]yFaio-Ffal6%FcblFfc[l%Ecomputing~the~valuati
on~of~the~r.h.sG>8;-Fi\o6%FgcyFSFe\p-Ffal6(FcblFfc[l%;the~valuation~of~the~r.h.
sGFaj[lF^f[lF]g[l@$2F[]lFaj[l>Fgcy-Ff^l6#*&FgcyF^r)Fe\pFaj[lF\s>Fdbn*&FdbnF^r)F
e\p-F]el6$F^i[l,(FjixF\sFaj[lF\sF^rF^rF^r>Fdbn*&FdbnF^r-F[e[l6$Fgcy-Fg^n6$FgcyF
SF\s-Ffal6(FcblFfc[l%4the~rational~factorGFdbnF^f[lF_f[l>Fgcn-Ffel6#,&FIF^r-Ff^
l6#*&-Fg^n6$FdbnFSF^rFdbnF\sF\s>%#_BGFgcn>F_`yFaio-Ffal6%FcblFfc[l%4calling~Mpo
lsolde~:G@%F]\o>Fdhq-&F"6#F[fx6%FgcnFS-Ffel6#*&Fbe[lF^rFdbnF\s>Fdhq-F_]\l6$Fgcn
FS-Ffal6(FcblFfc[l%6the~polynomial~factorGFdhq%&Time:GF_f[l@%0FdhqFe\m-F\hl6#-F
fel6#*&FdhqF^rFdbnF^rFfdyF$6#Fg\\lF$6#F\pF_`nt+6/Fh[lF^[lFfcmF6F`tF_[lF?FgcmFd`
nFhcmFc[lFicmFjcm6-F9F:Fajp%'valtabGFddmFh]mF^[o%#s1GF5F*Ff]qF$F$C4>F^p-FfelF`o
>F_sFham>Fijm-Fa]o6$FYFY>Fer-%$tabG6#;F^rFY?(FinF^rF^rFYFh]lC&>F[r-&F"6#%'colva
lG6&-&F\o6#Fe[q6$F^pFinFSFVFehq>&FerF_drF[r?(FboF^rF^rFYFh]l>&FijmF_cw&FaqFf^p>
F_s-F]elFjgu@$1F\sF_sC$>FhwF^ct-F\hl6#Fbfq>Fjq-F^q6$;F^rF^ct7#-Fdq6$F[]l/FinFca
\l>FfoF^r>FinF_s>FboF^r?(F$F^rF^rF$Fh]lC'@$/&FerFf^pFinC,>&Fjq6#,(F^rF^rFinF^rF
_sF\s,&Fbb\lF^rF^rF^r>F^p-Fajm6%F^pFfoFbo>F^p-F[jmFhb\l>Fijm-Fajm6%FijmFfoFbo>F
ijm-F[jmF]c\l@$F^jvC$>F\eq-Fajm6%F\eqFfoFbo>Fgaq-F[jm6%FgaqFfoFbo>Fc\l&FerFbal>
Fic\lF_b\l>F_b\lFc\l>Ffo,&FfoF^rF^rF^r@$2FYFfoFggm>FboFhz@$2FYFboC$>FboFfo>FinF
baw@$/Fin,&F_sF^rFcyF^rFggm?(FinFfoF^rFYFh]l@$0Fc`\lFhd\l?(Fc\lF^rF^rFYFh]l>&Fi
jm6$Fc\lFinF[]l>FjwFjq>F[s-Fcam6#-Fhs6$FerFban>FbxFfim?(FinFfoF^rFYFh]l>&Fijm6$
FinFin,&Fie\lF^rFjyF^r>FhwF^ct-FgrFdcqF$F$F$6#F_jvFfjot+Ffgo6<FhgoF[inFf]mF_`vF
6F9F:F?FadmFg]mF]\nF+Ffcm%)expPartsGFajpF_[oF_[l%&candsG%#x0G%%solsGF7%%ext1G%#
ssGF_\mFicm%#iTGF$F$C7>FjqFd[o>FinFaio@%F]\o>Faj[lFV>Faj[lFe\m>Fhjx-F\]m6$-Fe`p
6$Fd`v.Fi_p-Fe`p6$Faj[lFcg\l>F[rF^r>Fer-Fgx6%t+6$Ff]mF76#FajpF$F$-Fdq6$F_dl/Fin
&-Fjf[lFhfnFaflF$F$F$6#F\p-Fgx6$F`b[lFd`vFhjx>Fc\lFe\m?(F[rF^rF^rF]_mFh]l@$-Fg_
p6$&FerFjglFS>Fc\l-F\]m6$Fc\l<#Fjh\l-Ffal6%F^rF`\p%Xcomputing~exponential~parts
~at~all~finite~singularitiesG>FdiqFdx>F^]yFdx>Ff_yFhjx?&Ff]nFc\lFh]lC'-Ffal6&F^
rF`\p%1singularity~is:~G-Fi_pFicp>F_`y-Fi_p6$Ff]nFS>Fgcn-&F"6#%*exp_partsG6&FI/
FSF_`yFdbn<%-F_amF]ay.%(minimalG.%%ram1G>Ff_y-F\]m6$Ff_y-Fe`p6$FgcnFcg\l>Fdiq6$
Fdiq7$Fgcn,&FSF^rF\j\lF\s>Fdiq7#Fdiq-Ffal6&F^rF`\p%'found:GFdiq@%/FdiqFb`m>Fjix
<#F[]l>Fjix-&F"6#%+candidatesG6&FdiqFSFdbnFhjx-Ffal6%F^rF`\p%Hcomputing~exponen
tial~parts~at~infinityG>F]]l-Fgr6#-&F"6#%,mat_convertG6%FIFSFham>FgcyFe\m?&Fgcn
&-&F"6#%'newtonG6%F]]lFS<$Fgj\lFij\lF\\lFh]lC)>F_s-Fb`oFg]\l>Ff_y-&F"6#%,mat_ge
t_extGFb]r>F_s-Fh^m6#,$*&FSFbcu-F[y6$/FS*$FSF\sF_sF^rF\s>F_s,&F_sF^r*&-Fc^l6%F_
sFSF\sF^rFSF\sF\s>F^i[l-Fgx6$Ff^l-F[y6$Ff^]l-Fdem6#-&F"6#%7mat_get_transformati
onGFg]\l@%0-F[hx6#F^i[lF[]l>Fbe[l-Fh]nF[`]l>Fbe[l-F[y6$Ff^]l-Fdem6#-&F"6#%;mat_
get_inv_transformationGFg]\l@$4-F\ev6$F_sFgcyC&>Fgcy-F\]m6$Fgcy<#F_s-Ffal6&F^rF
`\p%Gtesting~exponential~parts~at~infinity:GF_s-Ffal6&FcblF`\p%9checking~candid
ates,~#~=G-F^_m6#Fjix?&F^pFjixFh]lC&>F`x-Fgr6#,(FIF^rF^pF\sF_sF\s>F`x-Fgx6$Ff^l
-Fgr6#-F`^n6$Fbe[l,&-F`^n6$F`xF^i[lF^r-Fgx6%Fg^nF^i[lFSF\s-Ffal6&F^rF`\p%:check
ing~Riccati~solutionG,&F_sF^rF^pF^r?&Ff]n-&F"6#%*poly_solsG6%F`xFSFaj[lFh]l@%4-
Fa`sFicpC$-Ffal6&F^rF`\p%7Riccati~solution~foundG*&-Fh]p6#-F[^p6$F]c]lFSF^r-Fhs
6$-Fhel6$-F`^n6$F^i[lFf]nFSFbanF^r>F^]y6$F^]y7$-Fix6#-Fegq6#F]c]l-FjhqFcd]l-Ffa
l6%F^rF`\pFc[r7#F^]yF$F$F$6#FgflFdcyt+6$Fg]mF6F[arF$F$C%>FinF^r?(F$F^rF^rF$31Fi
nFS/&FIF_drF[]l>FinFbaw@%/Fin,&F^rF^rFSF^r-F\hl6#F%-F\hl6#Fh]lF$F$F$6#F\pFggxt+
6(Fh[lFfcoF^[lFicmFjcmF_[l6-F*F9%&alphaG%'xalphaGF5Fc[lFjgoFfcwF+F:FhgoF$F$C)>F
`xFaio>F_s-&F"6#Facw6(FIFSFVFYFcyFjy>Fin&F_sF\\l>Ffo&F_sFafl>F_s&F_sFffl@%2Fag]
lF[]lC*>Fc\l-Ffbn6#-Fdq6$Fag]l/F^p;F^rFcg]l>Fer-Ffbn6#-Fdq6$)FVFag]lF^h]l>FjqFc
g]l?(FboFgflF^r-F^_mFj_mFh]lC$@$0&F_sFf^pF[]lC$>Fc\l-Ffbn6$Fc\l-Fdq6$,(Fag]lF^r
FboF^rFbcuF^r/F^p;F^rF\i]l>Fer-Ffbn6$Fer-Fdq6$)FVFci]lFdi]l>Fjq,&FjqF^rF\i]lF^r
@$2Fjq-Fe[oF_drC$>Fc\l-Ffbn6$Fc\l-Fdq6$F[]l/Fbo;F^r,&F`j]lF^rFjqF\s>Fer-Ffbn6$F
er-F]]n6#Fij]l>FY-Fgr6#-F`^n6$FYFer>Fcy-Fgr6#-F`^n6$-Ffbn6#-Fdq6$)FV,$&Fc\l6$Fb
oFboF\s/Fbo;F^rF`j]lFcy>F[r-Fgr6#-F`^n6$,(FjyF^rFinF\sFc\lF^rFer>F[r-Fgr6#,&Fjy
F^rFinF\s7%F[_[l-F[y6$/FVF[]l-Fgx6$Fgbm-F[y6$/FVF\dnF[_[lF_\rF$F$F$6#FcpFgiwt+F
`hl6%Fh[lFfcoFdgwF$F$C$>FinFdhl@%FcilFIC(>Fbo-F^q6#;F^rF\p>FefnFdil>F_gnFI>Fc\l
Fegl>&FiglFacmFboFc\lF$Fh[mF$6#F^r%(mat_sumGt+FfblFgblF$F$C3F]cl@$Fccl>FeclFicl
>FinF[dl>F_dlFaev>FboFcdl>Fc\lFidl>FerF\el>Fjq-Ffel6$,&F[flF^rF]flF^rFV>F`flFbf
l>FeflFjq>FiflFer>F]glFdx>F`glFI>FcglFS>F[rFegl>FhglFinF[hlF$F]hlF$6#FgflF\gmt+
6%F*F^[lFfcm6+%#A1GF9F:F`tFf]q%&pivotGF6F5FadmF$F$C+>FaqFd[o>F_sFh[o>FinF_em@$5
/FVFS/FVFham-F\hl6#-Fd`yF_dr>FboF^r>Fc\lF^r>FfoF[]l?(F$F^rF^rF$31Fc\lF_s1FboFaq
C%>F[rF[]l?(FerFboF^rFaqFh]l@$0&Fin6$FerFc\lF[]lC&>F[rFja^l>Fin-F[jm6%FinFboFer
>FfoF]d\lFggm@%/F[rF[]l>Fc\lFd_rC(-Fj[n6'F[rFVFSF\\nF^\n>&FinFcb[lF^r?(FerFd_rF
^rF_sFh]l>&Fin6$FboFer-Fc\n6%*&F\c^lF^rF]\nF^rFVFS?(FerFhzF^rFaqFh]l@$Fia^lC$?(
FjqFd_rF^rF_sFh]l>&FinFd\o-Fc\n6%,&*&&Fin6$FboFjqF^rFja^lF^rF\sFfc^lF^rFVFS>Fja
^lF[]l>Fc\lFd_r>FboFhzFfoF$F$F$6#Fgfl%*seq_evalaGt+6'Ff]mF^[lF7F9F:6#FbeyF$F$C%
@$/-Fe`p6$FVFcg\lFe\m-F\hlF`o>Fin-Fgao6$Fg^mFS-Fhs6$7%-Fdq6$*&-Fc^l6%FinFSF3F^r
)FSF3F^r/F3;F\gy,&FYF^rF\sF^r-Fdq6$*&-Fix6#Fde^lF^rFfe^lF^r/F3;FYFcy-Fdq6$Fce^l
/F3;,&FcyF^rF^rF^rFgfyFdtF$F$F$6#Fcbl%.ramified_caseGt+Ffgo61F]dmF6F+F,F]hoFadm
Ff]mFg]mFib[lF`tF2F[inF_gzF_[oF9F$F$C+-Ffal6%F^r.%.formal_reduceG%7entering~ram
ified~caseG>Fin,$-Fddl6%FIFSF\sF\s>FboFd\l-F_]l6&FIFSFc]lF\s@%1FboFgflC$>Fjq<#F
bo-Ffal6&F[glF]g^l%7ramification~detected:GFjqC+-Ffal6%FgflF]g^l%3applying~heur
isticG>Fc\l-F_]l6&FIFSFc]l,&FinF\sF\pF^r-Ffal6%F[glF]g^l%Dcomputing~characteris
tic~polynomialG>Faq-&F\o6#.%)charpolyG6$Fc\lF_s>Ffo-&F"6#%/newton_polygonG6&Faq
FSFSF_s-Ffal6&F[glF]g^l%=algebraic~newton~polynomial:G&FfoFffl>Fjq<#-Fdq6$-F`b[
l6#F^`x/F[rFji^l>Fjq7$-F_amFf_m-F_am6#-F[`p6$<$-Fdq6$F^p/F^p;F\p,&FboF^rFbcuF^r
FboFjq-Ffal6&F[glF]g^l%=candidates~for~ramification:GFjq>F]]lFbo>Ff]nFdx?(F$F^r
F^rF$FajxC->F[rFe\p-Ffal6&F[glF]dq%Fremaining~ramified~exponential~parts:GF]]l-
Ffal6&F\pF]dq%8doing~ramification~r~=~GF[r>Fer-Faao6&FIFSF`xF[r-Ffal6&FgflF]dq%
4looking~for~slopes:G<#-Fdq6$Ff]q/Ff]q;,&*&F[rF^r,&FinF^rFbcuF^rF^rF^rF^rF^r,&*
&F[rF^r,&FinF^rF\sF^rF^rF^rF\sF^r>Fdiq-&F"6#%3ramified_reductionG6)FerFSFh\_lF[
]_lF`xF[rF_s>F]]l,&F]]lF^r&FdiqFaflF\s-Ffal6&F[glF]g^l%Bramified~exponential~pa
rts~found:GFf]_l>Ff]n6$Ff]n-Fdq6$7$-Fb`o6$&FdbnFaflFSFa^_l/Fdbn&FdiqF\\l>Fjq-F_
\l6$Fa\lFjq-Ffal6&F[glF]dq%Bremaining~possible~ramifications:GFjq7#Ff]nF$F$F$6#
Fgfl%*valuationGt+6&%"aGF^[lFfcmF^[o6$%$ordG%%tempGF$F$C'>FinF[]l>FboFI?(F$F^rF
^rF$-%'divideG6%FboFV.F]dmC$>FboF]dm>FinFbaw@$Fify@%/FinF[]l>FY-Fc\n6%FIFVFS>FY
-Fc\n6%F]dmFVFSFinF$F$F$6#F[glF[hmt+6&F*F^[lFfcmFedm6.Fa`^lFicmFgf\lF9F:F`tFf]q
Fb`^lF6F5FadmFajpF$F$C/>F^pFd[o>FfoFh[o>FinF_em@$Fh`^lC)>Fin-&F\o6#%*gausselimG
6$-F`ap6$-F^hmF_dr-Ffbn6#-Fdq6$F^r/Fer;F^rFfoFb`n>F[rF^r>FerF^r>F`xF[]l?(F$F^rF
^rF$31F[rF^p1FerFfo@%/&Fin6$FerF[rF[]l>F[rF^cmC$>Fer,&FerF^rF^rF^r>F`xFaiu@$3Fi
fy2F`xFfo>FY-F^hm6#-Fefm6%Fin;FaiuFfo;,&F^pF^rF^rF^r,&F^pF^rFfoF^r-F\hlF]^l>Fin
Fia_l>F]]lFfo>FboFja_l>FerF^r>FjqF^r>F`xF[]l?(F$F^rF^rF$31FjqF^pFfb_lC%>F_sF[]l
?(F[rFerF^rFfoFh]l@$0&Fin6$F[rFjqF[]lC'>F_sF[e_l>Fin-F[jm6%FinFerF[r>Fbo-F[jm6%
FboFerF[r>F`xFaiuFggm@%/F_sF[]l>FjqF]`mC&-Fj[n6'F_sFVFSF\\nF^\n?(F[rF^c_lF^rFfo
Fh]l@$Fjd_lC%?(FaqF]`mF^rF^pFh]l>&FinFe^q-Fc\n6%,&*(F]\nF^r&Fin6$FerFaqF^rF[e_l
F^rF\sFaf_lF^rFVFS?(FaqF^rF^rF]]lFh]l>&FboFe^q-Fh^m6#,&*(F]\nF^rF[e_lF^r&FboFgf
_lF^rF\sFjf_lF^r>F[e_lF[]l>FjqF]`m>FerF^c_l@$Fify@$Fbc_lC$>FY-F^hm6#-Fefm6%FboF
hc_lF_b_l>FY-Fgx6&Fc\n-FdemF`]nFVFSF`xF$F$F$6#F[glFgixt+6%F_\nF^[lF66(F9Fh]mFg]
mFf]qF3F:F$F$C'>FinF^r>Fc\l-F]gp6$FVF[]l?(F$F^rF^rF$31FinFVFie]l>FinFbaw@%/FinF
\`s>Fbo,$FhamF\sC'>Fbo-&F"6#%)mydegreeG6$Fje]lFS>Fer7#Fin>FinFbaw?(F$F^rF^rF$Fj
h_lC%>FjqFbi_l@&2FboFjqC$>FboFjq>FerFhi_l/FjqFbo>Fer7$F`]qFin>FinFbaw?&F[rFerFh
]l>&Fc\lFjgl-%'lcoeffG6$&FIFjglFS7$Fbo-FdemFacmF$F$F$6#Fgfl%(intDiffGt+6%F[`vFb
_p%&whichG6,F]`vF^`vF_[oFf]mF_`vF6F``vF[inF5%'idiffsGF$F$C'>FinFe\m>Fc\lFd`v>Fb
oFc\l?(F$F^rF^rF$Fh`vC(>FerFefn>FfoF[]l>F^pFdx?&FjqFboFh]lC%>F[rF_av>FaqFbav@$F
eavC$>F_sFhav@$Fgbv@$FibvC&>F^p6$F^pF_s>Fc\lFecv>Ffo,&FfoF^rF[\pF^r@'/FVF]el@$2
F_sF[]l>FerFjq/FVFeez@$2F[]lF_s>FerFjq-Fen6#%/wrong~argumentG@%F]\o>Fin-F\]m6$F
in<#7$F]irFfoC%>F^p-Fcam6#7#F^p>F^p-Fgx6%f*6$F9F:F$Fc_pF$,&FIF^rFSF\sF$F$F$F^pF
icq>Fin-F\]m6$Fin<#7$F]irF^p>FboFc\lFinF$F$F$6#FgflFhelt+FahlF$F$F$Ff\sF$F$F$6#
F^rF[]ot+6&Ff]mF^[lFfcmF^[o6%Fh]m%#l1G%#l2GF$F$C&@$F]arFfam@$Fi`^lC%>Fin,&-F^iw
6$Fa[pFSF^r-F^iw6$-F`b[lF`oFSF\s@$Fify>FY*&-%'tcoeffGF```lF^r-Fh``lFb``lF\sFbfy
@$5Fj`^l/FVFg^]lC%>Fin,&-Fc_mFb``lF^r-Fc_mF```lF\s@$Fify>FY*&-Fjj_lF```lF^r-Fjj
_lFb``lF\sFbfy@%F]\o-F\hl6#,&-&F"6#F\__l6%Fa[pFSFVF^r-F\b`l6%Fc``lFSFVF\sC&>Fin
,&-F\b`l6&Fa[pFSFV.FboF^r-F\b`l6&Fc``lFSFVFe`nF\s-Fj[n6'Fc\lFVFSF\\nF^\n>FY-Fc
\n6%*&FboF^rF]\nF^rFVFSFbfyF$F$F$6#F[glFf`]lt+F^_o6&F6F5FicmFh[lF$F$C$>FerFdhl@
6/F]irFghlC%>Fc\l-Fd`]l6$&FerFf`oFS>Fin-F[oFacm-F_\u6$-F_ar6$-F]]n6#,&FinF^r-F[
o6#FcirF\sFcirFc\l/F]irF\ilC%>Fc\lFhc`l>FinF\d`lF]d`l/F]irF_il-Fd`]l6$&FerF[qFS
/F]irFdil-Fd`]l6$F[jrFS/F]irFiil-Fgx6%Fgao-F[y6$FcirF[e`lFS/F]irF\jl-F[y6$/FS*&
FcirF^r)FSFahrF^r-Fdem6#F[e`l/F]irFijl@%/-F^_mF`hrF^rC&>Fc\l-Fd`]l6$&FerFfboFS>
Fin-Fe\l6#Fff`l>FboFcir-Fgx6%Fgao-F_\u6$-F_ar6$-F]]n6#,&F\d`lF^rFinF\s-Fhct6&-F
gr6#,$&FerF^dtF\sF^drF^r,(FinF^rFboF\sF^rF^rFc\lFSC&>Fc\l-Fd`]l6$&FerF[coFS>Fin
-Fe\l6#F_h`l>FboFcir-Fgx6%Fgao-F_\u6$-F_ar6$Fag`l-Fhct6&-Fgr6#,$&FerFdhtF\sF^dr
F^rFjg`lFc\lFS/F]irFdjl-F]]n6#Fd\l/F]irF^co-F_\u6$-Fd`]l6$FahrFSFijr/F]irF_[mFd
i`lF$Fh[mF$6#F\p%<TestIntegrabilityConditionsGt+6%'F*Fcan'F^[lFcan'FcdmFcan6'F6
F5F_hoF9F:F$F$C'>Fin-Fe[o6#-Fhs6$F^]sFL>FboFddz>Fc\lFh]l?(FerF^rF^r,&FboF^rF\sF
^rFh]l?(FjqF^c_lF^rFboFh]l@%/-F[dp6#-Fgbm6#-Fgr6#,*-Fgx6%Fg^n-Fhs6$&FIF_blFLF\j
pF^r-F`^n6$-Fhs6$&FIFf_mFLFg[alF\s-Fgx6%Fg^nF\\al&FSF_blF\s-F`^n6$Fg[alF\\alF^r
F%C$>Fc\lF%-F\hl6#QBThe~connection~is~not~integrable!F$@$3F]\o/-F[dp6#-Fgbm6#-F
gr6#,*-Fgx6%Fg^n&FVF_blF\jpF^r-F`^n6$F\\alFe]alF\s-Fgx6%Fg^n&FVFf_mFa\alF\s-F`^
n6$Fg[alFj]alF^rF%C$>Fc\lF%-F\hl6#QIThe~right-hand~sides~are~not~compatible!F$-
F\hl6#Q>The~connection~is~integrable!F$F$F$F$6#Fgfl%-l_colechelonGt+6&F*F7Ficm%
$indG6,F+%%ind0GFadmF2FcdmF9FgitF:Fg[lF3F$F$C+>FerFh[o>Fc\lFd[o>FinF_em>FV-F]]n
F_bl>FboFb`m>F_sF^r@%/FSFe\m?(F[rFc\lF\sF^r1F_sFerC)>FfoF_s?(F$F^rF^rF$31FfoFer
/&Fin6$F[rFfoF[]l>FfoF]d\l@$/FfoF^c_l\>Fbo7$F[r-F_amFf^p@$2F_sFfoC$>Fin-F]hn6%F
inFfoF_s>FV-F]hn6%FVFfoF_s?(F^pF]d\lF^rFerFh]lC%>Fjq*&&Fin6$F[rF^pF^r&Fin6$F[rF
_sF\s?(FaqF^rF^rF[rFh]l>&Fin6$FaqF^p-Ff^l6#,&FjaalF^r*&FjqF^r&FinF`^oF^rF\s?(Fa
qF^rF^rFerFh]l>&FVF[bal-Ff^l6#,&FcbalF^r*&FjqF^r&FVF`^oF^rF\s>F_s,&F_sF^rF^rF^r
?(F[rFc\lF\sF^rFg_alC)>FfoF_s?(F$F^rF^rF$F[`al>FfoF]d\l@$Fb`alFc`al>FboFe`al@$F
h`alC$>FinF[aal>FVF^aal?(F^pF]d\lF^rFerFh]lC%>Fjq-Fix6#-Fegq6#Fcaal?(FaqF^rF^rF
[rFh]l>Fjaal-Fix6#-FegqF]bal?(FaqF^rF^rFerFh]l>Fcbal-Fix6#-FegqFebal>F_sFjbal@$
Fify>FYFbo-F_amF_drF$F$F$6#F[gl%)poly_valGt+6&F___lF^[lFcf\lF^[o6%Fcdm%$lc1GFh]
mF$F$@)F]arC$@$Fify>FYF[]lFfamFj`^lC$@$Fify>FY-Fjj_l6$-FgaoFhfnFS-F\hl6#,$-Fc_m
FhfnF\sF`]^lC$@$Fify>FY-Fh``lFieal-F\hl6#FgczC(>Fin-F[y6$/FS,&FSF^rFVF^rFI>Fc\l
F[]l>Fbo-Fix6#-F[y6$F[`sFin?(F$F^rF^rF$F`axC%>Fc\lFd_r>Fin*&-Fg^nFgeoF^rFc\lF\s
>FboF]gal@$Fify>FYFboFc\lF$F$F$6#F[gl%)col_pvalGt+6&Fh]mF^[lFfcm%$lcvG6(F6F^[oF
9F_[oF`[oF?F$F$C)>FjqFham>Fin-&F\o6#%(vectdimGF`o@$F]\o-F\hl6#-F]el6#-F_am6#-Fg
x6&FabqFd`vFSFV>FY-&F\o6#F]gpF]cq?(Fc\lF^rF^rFinFh]l@$0&FIFacmF[]lC%>F[r-%%pval
G6&FgialFSFVFfb`l>&FerFacm7$F[rFbo>Fjq-F]elFiq@$0FjqFham?(Fc\lF^rF^rFinFh]l@$/&
F^jalF\\lFjq>&FYFacm&F^jalFafl-F\hlFf_mF$F$F$6#F[gl%,simple_solsGt+6%F*F^[lFfco
6)FigoF___lF9F:%#apGF?FcflF$F$C(@%F]\o>FaqFV>FaqF\p>Fin7#-Fdq6$-Fhs6$7#-Fdq6$*&
&FboFed[lF^rFe[xF^r/Fer;F[]lFaqFdt/Fc\l;F^rFjn>Fjq-Fgx6%Fe^m-Fgx6$Fb[p-&F"6#Fhe
v6%FIFSFcgsFS>F[r-Fh_y6#-Fhs6$FjqF[io@$/F[rFdx-F\hl6#F[]l-Ffel6$-F[y6$Fc_[lFinF
SF$F$F$6#FgflFcaot+F^_oFahlF$F$C$>FinFdhl@8Ffhl-FaaoFd`oF[ilF`^blF^il-FaaoF\aoF
cil-FaaoFeaoFhil-Fix6#-F[y6$F]glFa^blF[jl-F[y6$F]boFa^blFhjl@%Fabo-FaaoFdbo-Faa
oFiboFcjlFSF]co-FaaoFacoF^[mF\_bl/F_dl.%+multipliedGF\_blF$Fh[mF$6#F\pF^g^lt+Ff
go62F]dmF6F_gzF_\mFg]mF`[l%(rhopolyGF9%"wGF_[oF+Ff]mF___lFedmF7%#A0GF$F$C0>Fc\l
Fdx@%Fgey>FerFV>FerFe\m>F[rFb`m@$-Fg_p6$Fer.%+invariantsGC$>F[r-F_\q6#-F_am6#-F
`dy6%Fg_pFerF_`bl>Fer-F^_pFh`bl>FaqFb`m@$-Fg_p6$Fer.%+rhopolygonGC&>Faq-F_\q6#-
F_am6#-F`dy6%Fg_pFerF_abl>Fer-F^_pFhabl>F_s-Fe`p6$FaqFfcl@$0F_sFe\m>FjqFag]l>Fi
nFag^lFeg^l>FboFd\l>Fdbn-F\^]lFhfn>Fc\lFdx-Ffal6%F^rF]g^l%7enter~formal~reducti
onG-Ffal6(F[glF]g^l%.pole~order~isGFin%9,~dimension~of~system~isGFbo@+2FinF^rC$
-Ffal6%F[glF]g^l%3system~is~analyticG-F\hl6#7#7$F[]lFI/FinF^rC'-Ffal6%F[glF]g^l
%7system~has~simple~poleG>Fgcn-Fgx6&Fc^l-F_]l6&FIFSF\sF\sFSF\s-Ffal6%F[glF]g^l%
Vcomputing~characteristic~polynomial~of~leading~matrixG?&F]]l-&F"6#Fh_v6$&-Fjf[
l6$-F\i^l6$FgcnFjqFdbnFaflFjqFh]lC'-Ffal6&F[glF]g^l%+factor~is:G&F]]lF\\l@%2*&-
Fc_m6$FbeblFjqF^r&F]]lFaflF^rFboC$-Ffal6'F[glF]g^l%4block-reducing~withG7#F^rFb
ebl>F`x-&F"6#F\gw6'FIFSF]fblFbeblFjq>F`xFI>Ffo-Fi_p6#&&-&F"6#F`[`l6%&-Fjf[l6$Fb
eblFdbnFaflFjq.F]elF\\lF\\l-Ffal6&F[glF]g^l%3minimal~eigenvalueGFfo>Fc\l6$Fc\l7
$*&FfoF^rFSF\s-&F"6#%4mat_change_exp_partG6%F`xFSFggbl-F\hl6#7#Fc\l/Fd\lF^rC'-F
fal6%F[glF]g^l%Ereducing~first~order~scalar~equationG>FfoFeg^l>Ffo-Fgao6%&Ffo6$
F^rF^rFSFix-Ffal6&F[glF]g^l%8exponential~part~found:GFfo-F\hl6#7#7$Ffo-Figbl6%F
IFSFfoFeamC&-Ffal6%F[glF]g^l%8calling~Moser-reductionG>F^p-Fg\m6&FIFSFjqF^]m-Ff
al6&F[glF]g^l%4partial~rho-polygonGFgjz-F\hl6#-&F"6#F^g^l6%FicqFS-F\]m6$Fer<$/F
_`blF[i[l/F_ablFgjzC'-Ffal6%F[glF]g^l%8system~is~Moser-reducedG>Fdiq&&Faq6#F`bm
Fafl-Ffal6&F[glF]g^l%Ncharacteristic~polynomial~of~leading~matrix:~GFdiq?&F]]l&
-Fjf[l6$FdiqFdbnFaflFh]lC$-Ffal6(F[glF]g^l%6irreducible~factor~isGFbebl%0multip
licity~isGFhebl@%0-Fix6#-F[y6$FaexFbeblF[]l@%/FfeblF^rC$-Ffal6%F[glF]g^l%=irred
ucible~factor~is~linearG@%2F^rFheblC*-Ffal6&F[glF]g^l%4multiple~root,~m~=~GFheb
l>Ff]n-Fi_pFgebl>Ffo*&Ff]nF^r)FSFinF\s-Ffal6&F[glF]g^l%;changing~exponential~pa
rt:GFfo>F`xFbibl-Ffal6'F[glF]g^lF\fblFhi_l)FjqFhebl>F^p-Fajbl6%-F`fbl6'F`xFSFhi
_lFj]clFjqFSFer>Fc\l6$Fc\l-Fdq6$7$-Fb`o6$Fcg]lFSFcg]l/F_sF^pC%-Ffal6'F[glF]g^lF
\fblFhi_lFbebl>F`x-F`fbl6'FIFSFhi_lFbeblFjq>Fc\l6$Fc\l-F_am6#-Fajbl6%F`xFSFer@%
32F[glFbo2F[]l,&FboF^rFeeblF\sC&-Ffal6%F[glF]g^l%>isolating~algebraic~extension
G-Ffal6'F[glF]g^lF\fblFhi_l)FbeblFhebl>F`x-F`fbl6'FIFSFhi_lF_`clFjq>Fc\lF__clC+
-Ffal6%F[glF]g^l%Cconstant~field~extension~necessaryG>Ff]nF`]cl-Ffal6(F[glF]g^l
%(root~isGFf]n%1multiplicity~is:GFhebl>FfoFb]clFd]cl>F`xFbiblFh]cl>F^pF\^cl>Fc
\lFa^cl@%/FheblFboC&-Ffal6%F[glF]g^l%/nilpotent~caseG>Fdiq&&FaqF\\lFafl-Ffal6&F
[glF]g^l%3next~k-polynomial:GFdiq@%2F[]l-Fc_m6$FdiqFjqC(-Ffal6%F[glF]g^l%;separ
ating~ramified~slopesG>F`x-F`fbl6'FIFSF[rF^rFjq>Fc\l6$Fc\l-F_am6#-&F"6#Fgf^l6%F
`xFSF[r-Ffal6%F[glF]g^l%3applying~recursionG>F`x-F`fbl6'FIFSF[rFdiqFjq>Fc\lF__c
lC$-Ffal6%F[glF]g^l%7no~decoupling~possibleG>Fc\l6$Fc\l-F_am6#-F]ccl6%FIFSF[rC%
Fh]cl>F`x-F`fbl6'FIFSFhi_lFj]clFjq>Fc\lF__clF]hblF$F$F$6#Fgfl%*ReductionGt+6%F]
j`lF^j`lF]ho6)F6Fadm%+randchoiceGFedmFfdmF+F9F$F$C*>Fin-Fi[o6#F^]s>Fbo-Fi[oFdo>
Fc\l-F[jn6#;F[]lF^r>Fer-F`ap6$FV-&F\o6#%,stackmatrixG6$-Fhs6$-%'MatrixG6%Fbo,&F
inF^rFboF\s/%&shapeG%%zeroGFL-Ffbn6#-%"$G6$F^rFcfcl?(F$F^rF^rF$/-F[hxF_blF[]l>F
er-F`ap6$FV-&F\o6#%+randmatrixG6%FinFcfcl/%(entriesGFc\l>Fjq-Fh]nF_bl?(FaqF^rF^
r-F^_mF`oFh]l>&F[rF`am-Fgx6$-F^es6$%'factorGFf^l-Fgr6#-F`^n6$Fjq,&-F`^n6$&FIF`a
mFerF^r-Fgx6%Fg^nFer&FSF`amF\s-F\hl6#7#-Fdq6$-Fgr6#F^hcl/Faq;F^rF\hclF$F$F$6#Fg
fl%.l_column_rankGt+6)F*F^[lF`tF`[lF][l%'mk_newGF76.Fb[l%#K1GF6Fedm%#SSGFjfvF[g
vF9F:FhcmFajpFi^alF$F$C1>FinF\\z>Fc\lFd\l-Ffal6%F\pFhalF`gv>Fbo-Fgx6%FhsFcyFban
>Fer-&F"6#Ff^al6&-&F\o6#.F`hm6#-Fhs6$7#Ff`alFKFhwFf`nFahm>F^p,&Fc\lF^r-F^_mF[^n
F\s>F`xF^r?(F_sF^rF^rFg]lFh]lC$?(FfoF`xF^r,(F`xF^r&FinFj_mF^rF\sF^rFh]l>Fer-&F
\o6#.%'mulrowG6%FerFfo)FS,&FVF^rF_sF\s>F`xFfo>F[r-Fa]oFj\l>FfoF^r?(F_sF^rF^rFc
\lFh]l@%30F]]lFb`m/F_sFbebl>F]]l-F_\l6$Fa\lF]]lC$>&F[rFa^lF^r>FfoF]d\l?(F_sF^rF
^rFc\lFh]l?(FfoF^rF^r,&Fc\lF^rF^pF\sFh]l>&F[r6$F_sF[d_l&Fer6$F_s,*Fc\lF^rF^pF\s
F^rF^rFfoF\s>Faq-F]doFjgl>FjyF^p-Fh\x6'FIFSF[rFaqF[]lF$F$F$6#F\q%.mat_get_terms
Gt+F^_oFahlF$F$C$>FinFdhl@4Ffhl-Fgx6'Fc[sFeflFSF`flFcglF[il-Fgx6'Fc[s-F\^u6%Fef
lFSFcglFSF`flFcgl-F\ev6$F_dl<&FijlF\jlFiilF_il-Fgx6'Fc[sFeflFSF`flFiflFcil-&F"6
#Ff^dlFeaoFcjl-Fgr6#*&FeflF^rF[`uF^rF]coFdevFb[mFd_dlF^[mFd_dlF[[mFd_dlF$Fh[mF$
6#F\p%,mat_productGt+Ffbl6*F,FhblFibl%#htGF3FjblF[clF^[oF$F$C3F]cl@$Fccl>FeclFi
cl>FinF[dl>F_dlF_ev>FboFcdl>Fc\lFidl>F_s-F]el6$,(F_elF^rFboF\sF^rF^r,(FcelF^rFc
\lF\sF^rF^r>F[r-Fgx6$Ff^l-Fgr6#-F`^n6$-F_]l6&FIFVFbo,(FboF^rF_sF^rF\sF^r-F_]l6&
FSFVFc\l,(Fc\lF^rF_sF^rF\sF^r>F`flFbfl>FeflF[r>Fifl,*FboF^rFc\lF^rF_sF^rF\sF^r>
F]glF_s>F`glFI>FcglFS>FaqFegl>&FiglF`amFinF]]oF$F]hlF$6#Fgfl%*poly_pvalGt+F^__l
6&%#a1GFa__lFb__lF]dmF$F$C,>FinFjeal@$4-FG6$FIFieo-Fen6$%5polynomial~expected:G
FI@$F]arC$@$Fify>FYF[]lFfam@$5F\a`lFj`^lC$@$Fify>FY-Fjj_lFgeoF[fal@$Fi`^lC$@$Fi
fy>FY-Fh``lFgeoFcfal>FboF[]l>Fc\lFI@%/-F\]m6$-Fe`p6$FIFcg\lFhd^lFe\m?(F$F^rF^rF
$-Fh__l6%Fc\lFV.FerC$>Fc\l-Fg^n6$Fc\lFS>FboFhz?(F$F^rF^rF$-Fix6#-%'DivideGF[edl
C$>Fc\lFer>FboFhz@$Fify@%F`ax>FY-Fix6#-%$RemGFc`_l>FY-Fix6#-F`fdl6%FerFVFSFboF$
F$F$6#F[glF]cst+6%F*F^[lFicmF$F$F$-Fgx6$Ff^l-Fgr6#-F`^n6$F\dnF_fvF$F$F$6#FgflFh
_xt+6%F][lF^[lF`[l6.%#DDGF]dmF`tF6F5F9F:Fgf]lF*F-FfcoFdgwF$F$C8>FinFafn>FboFhhw
@$31FboF^r2F^rF_hw-Fen6#%;wrong~number~of~invariantsG>Fc\l-Feez6$F[]l,&FboF^rF_
hwF\s>FerFd\l-Ffal6&FcblFgir%Dconverting~to~k-simple~system,~k~=~GFc\l-Ffal6&Fc
blFgir%/invariants~areGFV@%Figdl>F_s7#-Fdq6$F[]l/F[r;F^rFerC$>Fjq,&FerF^r-Fhs6$
-F_\l6$Fa\lFVFdtF\s>F_s7$-Fdq6$-Fdq6$,&F[rF^rF_hwF\s/Faq;F^r&FV6#F^cm/F[r;F^r,&
F_hwF^rF\sF^r-Fdq6$F[]l/F[rFf_q>F^p-F_ar6#-Fdq6$)FS,$F[rF\s/F[rF_s-Ffal6&FcblFg
ir%3diagonal~matrix~isG7#F[[el-Ffal6&FcblFgir%>evaluating~system~up~to~orderG,&
Fc\lF\sF\sF^r>Ffo-F_]l6&FIFS,$FboF\sFg[el>F`x-Fceo6#*&)FSFd_rF^r-Fgr6#,&-F`^n6$
FfoF^pF^r-Fgx6%Fg^nF^pFSF\sF^r>F_dlF^co>F`flFc\l>Fefl-F_amFdcq>Fifl-F_amF]^l>F]
glF[]l>F`glFI>F]]lFegl>&FiglF[^nFinF]]lF$Fh[mF$6#FgflFei^lt+6'Ff]mF^[lFfcmFg]m%
)moreargsG64%&lode1GF6F?%#ypGF^[o%$lcsG%%dropG%&npolyGF[inFcepF9F:F`tFc[`l%&slo
peGF_gz%#dpGF_`vF$F$C5>Fbo-Fc_m6$FIFY>Fjix-Fgao6$Fg^mFY>Fin7#-Fdq6$-Fc^l6%FIFYF
`x/F`x;F[]lFbo>Fdiq.%-differentialG@&Feaw@&-FG6$&FjpF^glFj`m>Fbo,&Fe_elF^rF^rF^
r-F\ev6$Fe_el<$Ff_oF__el>FdiqFe_elFapC$@$Fc_el>FboFg_el@$-F\ev6$&FjpFaglFj_el>F
diqFb`el@%FgcdlC$>Fc\lF\s>FdhqF^rC$>Fc\lF^r>Fdhq-Fg^nFg[n>FerFdx>F[rFdx?(F`xF^r
F^rFhzFh]lC$>Fer6$Fer,&-&F"6#F[jal6&&FinF]^lFSFVFf`nF^r*&Fc\lF^r,&F`xF^rF\sF^rF
^rF\s>F[rF\e_l>FerF\[x>F[r7#F[r>Faq-F]el6#F`]q>F_sF[]l?(F]]lF^rF^rFhzFh]l@$/&Fe
rF[^nFaq@'/FdiqFf_o>F_s,&F_sF^r*&&F[rF[^nF^r)*&Fc\lF^rFYF^rFf\yF^rF^r30FVFham0F
VFg^]l>F_s,&F_sF^r*&-Fc\n6%*&)FdhqFf\yF^rF[celF^rFVFSF^r-F`ev6$,(F]celF^rFf]nF
\sF^rF^r/Ff]n;F^rFf\yF^rF^r>F_s,&F_sF^r*&F[celF^rFhcelF^rF^r>Ffo-Fc_m6$F_sFY@%2
F[]lFfo>Fgcn7$F[]l-Fgao6$-%)primpartGFbdelFY>FgcnFdx>F`xF]d\l?(F$F^rF^rF$2F`xFh
zC(>Fdbn-F]el6#-Fdq6$*&,&FebelF^r&FerF]^lF\sF^r,&F]]lF^rF`xF\sF\s/F]]l;FaiuFhz>
F_sF[]l?(F]]lF`xF^rFhzFh]l@$/,&FebelF^r*&FheelF^rFdbnF^rF\sFgeel>F_s,&F_sF^r*&F
[celF^r)F]cel*&FheelF^r-F`b[lFiexF\sF^rF^r@$3/FdiqF__elF`cel>F_s-Fc\n6%-F[y6$/F
Y*&F[aelF^rFYF^rF_sFVFS>Fgcn6$Fgcn7$FdbnFgdel>F`x,&F`xF^r*&FadelF^rFffelF^rF^r7
%FaqFY7#FgcnF$F$F$6#Fcbl%2RationalSolutionsGt+6%F]j`lF^j`l'F]\mFcan69%$optG%$pa
rGFcdm%$MatG%$varG%'newvarGF6F5F_\nFb_p%#XXGF?F]ho%#SPG%#RRG%$RRRGFcgw%%invHGFh
fv%$RecGF9F+Fc`qF$F$C$@%Fjip>Fin.%)standardG@%33F]\o/F_hwF\p/Fhhw.%&paramGC$>Fi
nFdiel>Fbo&FVFafl@%3Faiel/Fhhw.F_\qC$>FinF]jel>Fc\lFiiel@%3F]\o/F_hwF[glC$>Fin.
%*param_rhsG@%F\jelC$>Fc\lFiiel>Fbo&FVFjflC$>Fc\lF\[fl>FboFiiel-F\hl6#Q<ERROR,~
wrong~type~of~input!F$@%5553/FinF]iel/-Fji`lFhfnFd^al3/FinF]jel/-Fji`l6%FIFSFc
\lFd^al3/FinFdielFi[fl3/FinFfjelF]\flC&>FaqFcj`l>F_sF\hcl@%5F\\flFc\fl>Ffo-&F"6
#Fbaz6%Fej`l&FSF\\lF]^p>Ffo-F[]fl6$Fej`lF^]fl@%5/-F_am6#-Fhs6$FfoF\ioF[]l/FfoFe
\mOFe\mC'>F^p-Fgx6$f*Fa_pF$Fc_pF$F][rF$F$F$Ffo@%5Fh[flF\\fl>F`x-F[`p6$-F\]m6#-F
_am6#-Fhs6$F^pFcan<#-F_amF\^q>F`x-F[`p6$Fe^fl<$Ff`alF\_fl>F]]lFd[[l@%F\ix@%/F_s
F^rC$>Ff]n-Fa]o6%FaqF^r7#-Fjfcl6$F[]lFaq>FdiqF_\rC$>Ff]n-Ffbn6#Fj_fl>FdiqF_\r@%
Fh\flC$>Ff]n-&F"6#%8MatrixColumnPrimpartRatG6$-F`ap6#-Fdq6$-Fgx6$f*Fa_pF$Fc_pF$
-Fc^l6$FIF\bpF$F$6&FdhelF`xF9FhjxFfo/Fhjx;F^rF]]lF^]fl>Fdiq-F[y6$<#-Fdq6$/&F`x6
#FhjxF[]lFdaflF_\r>Ff]n-Ff`fl6$-F`ap6#-Fdq6$-Fgx6$f*Fa_pF$Fc_pF$FaaflF$F$FcaflF
foFdaflF^]fl@%Fd_fl@%Fh\flO7$-FgrFicp-Fhs6$FdiqFa]oOF]cflC*>Fdbn-F`ap6$Ff]n-Ffb
n6#-Fdq6$F^r/Fhjx;F^rFaq>Fgcn-Fca_lFiex>Fdhq-Fefm6%FgcnFeaflFeafl>Fjix-Fh]n6#Fd
hq>F_`y-Fefm6%FgcnFeafl;F[hp,&FaqF^rF]]lF^r?(FhjxF\pF^rF_sFh]lC$>&Ff_yF^bfl-Fgx
6$Fgbm-Fgr6#-F`^n6$-F`^n6$FjixF_`y,&-F`^n6$-Fhs6$&FIF^bflFLFf]nF^r-Fgr6#-Fgx6%F
g^nF]cfl&FSF^bflF\s@$Fh\fl>&FgcyF^bfl-Fgx6$Fgbm-Fgr6#-F`^n6$Fbefl-Fgr6#,(&Fc\lF
^bflF^r-Fgx6%Fg^n-FgrFciqF^fflF\s-F`^n6$FgeflFdiqF^r@%Fc\fl>F^]y-&F"6#Fjgel6%7#
-Fdq6$F[efl/Fhjx;F\pF_s7#-Fdq6$F^fflFjgfl7&FdielFboF]jel7#-Fdq6$FafflFjgfl@%Fa
\fl>F^]y-Fdgfl6%FggflF\hfl7$FdielFbo@%F\\fl>F^]y-Fdgfl6%FggflF\hfl7$F]jelF`hfl>
F^]y-Fdgfl6$FggflF\hfl@%/F^]yFe\mOFe\m@%Fh\flO7$-Fgx6$Ff^l-Fgr6#-F`^n6$Ff]n-F_a
m6$F^rF^]y-Fgx6$Ff^l-Fgr6#,&-Fgx6$Ff^l-Fgr6#-F`^n6$Ff]n-F_am6$F\pF^]yF^rF^cflF^
rO-Fgx6$Ff^l-Fgr6#-F`^n6$Ff]nF^]y-F\hl6#QIERROR,~the~connection~is~not~integrab
le!F$F$F$F$6#Fgfl%(getargsGt+6#%)someargsG6%F9F]\mF_helF$F$C&@$2FbpF\p-Fen6#%:w
rong~number~of~argumentsG>FboFdx?(FinFgflF^rFbpFh]lC&>Fc\l&FjpF_dr@$4-FG6$Fc\lF
[jo-Fen6#%Eoptions~must~be~an~equation~sequenceG@$4-F\ev6$-F\[pFacm<&.%&constG.
%'methodG.%'outputG.%'placesG-Fen6$%0unknown~option:GFa]gl>FboFcb[l@'-FG6$F^[pF
K7&-FgrF][pFajoF[]l7#Fbo-FG6$F^[pFban7%-%(getLodoGF\]pF[]lFd^gl7$-%'getOdeGF\]p
Fd^glF$F$F$6#F^r%-randomMatrixGt+6'F5F^[lF\dzF]dz%#raG6(Fh[lF9F:FadmF[^m%#nrGF$
F$C$?(F$F^rF^rF$Fh]lC'>Fin-Fa]o6%FIFIF[]l@%Feaw>Fer-F[jn6#;F[]lFcy>Fer-F[jn6#;F
[]lF\p@%/FVFb`mC$>F[r-F[jn6#;!"$Fgfl?(FboF^rF^rFIFh]l?(Fc\lF^rF^rFIFh]lC$>Fjq-F
erF$@$Faex>Fib^l-F[rF$?(FboF^rF^rFIFh]l?(Fc\lF^rF^rFIFh]lC$>FjqF_agl@$Faex>Fib^
l-&F"6#Ficz6%FSFVFY@$2F[glFIFggm@$0-Fgbm6#-F[hxF_drF[]lFggmFcgsF$F$F$6#FcblF\]]
lt+6&F*F^[lFcf\lFc[`l6&F6Fh[lFh]mFfcoF$F$C4-Ffal6%FcblFgir%2enter~mat_convertG>
Fbo-F^q6#;F^rFew>FinFjn>FefnFdjl>FagnFjhr>FegnFV>&FboFbguF__el@$3Fify/FYFf_o>Fe
cglFf_o-Ffal6&FcblFgir%.localizing~atGFV@%F_cel@%-Fg_pFid^l>Fhgn-Fgx6%t+FeczF$F
$F$-Ff^l6#*&Fa[pF^r-Fheq6%Fc``lFS<#T$F\sF$F$6$Fcf\lFV6#F\p-F[y6$FifalFbfqFS>Fhg
n-Fgx6%Fheq-Faeo6#F]eglFS>Fhgn-Fgx6$Fix-Faeo6#,$*&FSFbcu-F[y6$Ff^]lFbfqF^rF\s>F
c\l-F]el6#-F_am6#-Fgx6&Fabq-Fhs6$FhgnF[ioFSF[]l>Fhgn-Faeo6#*&FhgnF^r)FSFc\lF\s>
F_gnFc\l>Fcgn,&Fc\lF^rF\sF^r>&FboF[q-F\]m6$Fhd^lFag\l>FerFegl>&FiglF_blFboFerF$
Fh[mF$6#F[glFa\]lt+6&%)exppartsGF^[lFajpF76)%#p1G%#p2GFfcm%&partsGF_gzF9Fef\lF$
F$C'@$/F\hclF^r-F\hl6#<#-Fdq6$-Fix6#-%,evala/TraceG6%-F[y6$/FV&F^]sFafl&&&F^]sF
\\lFjglFaflFe\mFY/F[r;F^r-F^_m6#F`igl>Fer-F_\]l6&-F_\l6$Fa\lFIFSFVFY>FjqFe\m?&F
inF`iglFh]l?&FboFerFh]lC$>Faq-Fe`p6$FinFcg\l>Fjq-F\]m6$Fjq<#-Ff^l6#-Fix6#,&FboF
^r-%&TraceG6%-F[y6$F\iglF`flFe\mFYF^rFjqF$F$F$6#F[glFgfut+6%Fh[lF^[lF9F$F$F$-Fg
x6&Fc^l-F_]l6&FIFSFVFVFSFVF$F$F$6#FgflFb]yt+6%F*FcdmF6FahlF$F$C$>Fin-Fca_l6#-F`
apFhfn7$-Fefm6%Fin;F^rFVFa\hl-F_`\l6$FinF\`sF$F$F$6#Fgfl%.gen_sylvesterGt+6*F*F
+F,F-F.F][lFfcwF26-F5F6F7F8FicmF^jcl%#L1GFf`sF9F:Fb_pF$F$C)@$FabrFZ>FinFjn>FboF
_^q@%5Feaw32FcblFbp4-FG6$Fb`elF]qC&@%/FbpFew>Fer-Fgr6#,&FIF^r*&FjwF^rFSF^rF\s>F
erFI>Fjq-Fgr6#,&-Faiy6$-F^[dlF\^qFIF^r-Faiy6$-F^[dlF`]nFVF\s>F[r-F]hq6$FjqFehq@
$1F\qFbpC$>FjyFaa[l>Fhw-FgrF`amC$>F[rFjy>FaqFhw>F_s-Fjfp6#7#-Fdq6$-Fdq6$&FcyFe
\el/FfoF]r/F^pF`u@%Fg]hl>F`x-&F"6#F_ep6$-&F\o6#.Fbap6$F[r-Fgr6#-F`^nF`^o&FjpFbg
u>F`x-Fb`hl6#Fe`hl-Fho6#7#-Fdq6$7#-Fdq6$&F`x6#,&FfoF^r*&F^pF^rFinF^rF^r/F^p;F[]
lFjj`lF]`hlF$F$F$6#FewFfdut+6%Fh[lF^[lFgco68F6F5F3F`tFf]qF7F]dmFicmFdepFgf\lF*%
%A_21G%#XhG%#YhGF-%%D_21GFfco%%N_21G%#QhG%#RhGF_[oF2F$F$C'@$1FVF\^sFjd^l>Fc\lF
\^s>FaqFibt>FboF[gs@%F^btC3-Ffal6&F[glFgir%Bblock~reducing~system~up~to~orderGF
V>FinFcct>F[r-F\^]lFdct@$Fabs>F]^s,$FaqF\s>Ff_y&FIFd]x>F_s-Fi[w6&F]dtF^drF^rFjg
`l>F^p-Fi[w6&-Fgr6#,$F]dtF\sF^drF^rFjg`l-Ffal6&F[glFgir%Aupdating~terms~before~
splitting:G;Fd_rFV>F`x-Ffel6$-F_\u6$-F_\u6$F^p-F_]l6&FgbtFSFd_rFVF_sFS>Ffas-Fgr
6#,&FfasF^r-F\z6%F`x;F^rFcfclFiehlF^r>F]]l-F\z6%-F_]l6&FgbtFSFjchlFV;Fjg`lFinFi
ehl>&FIFbgu-Fgr6#,&FafhlF^r-Fgx6&Fc[s-F\z6%F`xFiehlF_fhlFSFd_rF^r>Fb_s-Fgr6#,&F
b_sF^r-F\z6%F`xF_fhlF_fhlF^r?(F$F^rF^rF$Fh]lC/>Fc\l-F]el6#-F_am6#-Fhs6$-Fgx6%&F
"6#FcczFafhlFSF[io@$2FVFc\lFggm-Ffal6&F[glFgir%9block~reducing,~order~isGFc\l>F
_`y-Fgx6&Fc^lFafhlFSFc\l-Ffal6%F[glFgir%;solving~sylvester~equationG>Ff]n-Fgr6#
*&)FS,&Fc\lF^rFaqF^rF^r-&F"6#%*sylvesterG6)F_bs&FIF[fu-Fgr6#,$F_`yF\s&FI6#"#7&F
IF^bxFf_yF\ihlF^r-Ffal6%F[glFgir%4transforming~blocksG>F]dt-Fgr6#,&F]dtF^rFf]nF
^r>Fb_s-Fgx6&Fheq-Ffel6$,&Fb_sF^r-F_\u6$F]]lFf]nF^rFSFSF[r-Ffal6%F[glFgirFf\u>F
afhl-Fgx6'Fheq-Ffel6$,*FafhlF^r-F_\u6$FfasFf]nF^r-F_\u6$Ff]nFb_sF\sFe^nF\sFSFSF
[rFd_rFjihl>Ffas-Ffel6$,&FfasF^r-F_\u6$Ff]nF]]lF\sFS>F\^sFV>Fb_s-Fgx6(&F"6#Fad^
lFb_sFSF[rF]^sF\^sFIC=>FinF]ht>F[r-F\^]lF^ht@$FabsC$>F]^sF\^s@$2F]^sFjchl>F]^sF
jchl>Fer,&FaqF^rF_btF\s>Ff_y&FIFdex@%/F]^sFjchl>FjqFaq>FjqFer-Ffal6&F[glFgir%Kb
lock~reducing~simple~operator~up~to~orderG,(FjqF^rFVF^rF^rF^r>Fhjx-Fefo6&&FIFc_
xFSF[]lFf]il>Fdbn-F_\u6$&FIFidx-F_\u6$&FhjxF\\l&FIF]ex>Fdhq-F_\u6$F^^il-F_\u6$&
FhjxFaflFb^il>Ffo-Fi[w6&FchtF^drF^rFjg`l>F^p-Fi[w6&-Fgr6#,$&FIFacxF\sF^drF^rFjg
`l-Ffal6%F[glFgir%Lupdating~terms~before~generalized~splittingG>Fdbn-Fgx6&Fheq-
F_\u6$F^p-F_\u6$FdbnFfoFSF[r>Fdhq-Fgx6&Fheq-F_\u6$F^p-Fgr6#,&-F_\u6$FdhqFfoF^r-
F_\u6$Fdbn-Fgr6#*&)FSF^c_lF^r-Fgx6%Fg^nFfoFSF^rF\sFSF[r>Ffas-F\z6%FdbnFiehlFieh
l>Fgcn-F\z6%FdbnF_fhlFiehl>Fafhl-F\z6%FdbnFiehlF_fhl>F_bs-F\z6%FdbnF_fhlF_fhl>F
bihl-F\z6%FdhqFiehlFiehl>Fjix-F\z6%FdhqF_fhlFiehl>F]dt-F\z6%FdhqFiehlF_fhl>Ffih
l-F\z6%FdhqF_fhlF_fhl?(F$F^rF^rF$Fh]lC6>Fc\l-F]el6$-Fbdo6%FafhlFSF[]l-Fddo6%F]d
tFSF[]l@$2Ff]ilFc\lFggm>F_`yFbhhl>F^]y-Fgx6&Fc^lF]dtFSFc\lF^hhlFdhhl>Fhjx-&F"6#
F'6,F\dhlFiihl&FIFibxFgbt-Fgr6#,$F^]yF\sFcihl&FI6#"#>&FI6#"#?Ff_yFc\lFjihl>Ff]n
-Fgr6#*&FjfglF^rFa^ilF^r>Fdiq-Fgr6#*&FjfglF^rFh^ilF^r>Fcht-Fgr6#,&FchtF^rFf]nF^
r>Fb_il-Fgr6#,&Fb_ilF^rFdiqF^r>F_bs-Ffel6$,&F_bsF^r-F_\u6$FgcnFf]nF^rFS-Ffal6%F
[glFgir%@reducing~algebraic~numbers~in~DG>Fafhl-Fgx6&Fheq-Fgr6#,(FafhlF^rFa[ilF
^r-F_\u6$FdiqF_bsF\sFSF[r>Ffihl-Ffel6$,(FfihlF^r-F_\u6$FjixFf]nF^r*&F\ailF^r-F_
\u6$FgcnFe^nF^rF\sFS-Ffal6%F[glFgir%@reducing~algebraic~numbers~in~NG>F]dt-Fgx6
&Fheq-Fgr6#,*F]dtF^r-F_\u6$FbihlFf]nF^r-F_\u6$FdiqFfihlF\s*&F\ailF^r-F_\u6$Ffas
Fe^nF^rF\sFSF[r>Ffas-Fgr6#,&FfasF^r-F_\u6$FdiqFgcnF\s>Fbihl-Fgr6#,&FbihlF^r-F_
\u6$FdiqFjixF\s>Fb_s-Fgx6$Fh^m-Fgr6#*&-&F"6#Fcco6&F_bsFfihlFSFf]ilF^rF\ailF\s>F
\^sFVFIF$Fh[mF$6#FgflFe[st+6&Ff]mF^[lF9F:6#F`tF$F$@%F]arF[]l@%Fify-Fhs6$7#-Fdq6
$*&-Fc^lF\_zF^rFc]clF^r/FinF[[uFdt-Fhs6$7#-Fdq6$Fajil/Fin;FVF^falFdtF$F$F$6#F[g
lFg_]lt+F^_oFac`lF$F$C$>FerFdhl@6Fec`lC%>Fc\l-Fe_]lFic`l>FinF\d`l-F_\u6$Fc\l-F_
ar6$-F]]n6#,&FinF^r-F[o6#F[[sF\sF[[sFfd`lC%>Fc\lFb[jl>FinF\d`lFd[jlFjd`l-Fe_]lF
\e`lF^e`l-Fe_]lF`e`lFae`l-Fgx6%Fgao-F[y6$FcirF`\jlFSFfe`l-F[y6$Fie`l-Fdem6#F`\j
lF^f`l@%F`f`lC&>Fc\l-Fe_]lFef`l>FinFhf`l>FboFcir-Fgx6%Fgao-F_\u6$Fc\l-F_ar6$Fag
`l-Fhct6&Fig`lF^drF^rFjg`lFSC&>Fc\l-Fe_]lF^h`l>FinFah`l>FboFcir-Fgx6%Fgao-F_\u6
$Fc\l-F_ar6$Fag`l-Fhct6&F_i`lF^drF^rFjg`lFSF`i`lFai`lFci`l-F_\u6$-Fe_]lFgi`lFij
rFhi`lFe^jlF$Fh[mF$6#F\p%1mat_colvaluationGt+6%F*F^[lF9F$F$F$-Fen6#%=should~not
~be~called~anymoreGF$F$F$6#FgflF[hblt+6%F][lF^[lFd_bl6%FfcoFdgwF7F$F$C0@$Fa_o-F
\hl6#-&F"6#%6mat_change_eigenvalueGFi_x@$/-FixFdoF[]lFjd^l>FinF[dl>F_dlF_il>F`f
lFbfl>Fc\l-F\^]lF`o>Fefl-Fgx6&Fheq-Fgr6#,&-Fg_dlFhfnF^rFVF\sFSFc\l>FiflFjfn>F]g
lFV>F`gl-F\]m6$Fc\lFhd^l>FcglFI>FboFegl>&FiglFf^pFinFboF$Fh[mF$6#Fgfl%)reg_sols
Gt+6$F*Ff]el6AF6FcdmF`gzFc`qFg]mF_[o%#k1GF+F,FdepFedmFf]mF3F9F-FfcoFicmFgf\lF_
\mF7Fef\lF^[lFcf\lFajpFhgoF_gzFadmF[^mF_[lF_`q%#coGF$F$C?>F^i[lFaio>FablF[]l>F_
`yFe\m@$2FbpFgfl-Fen6#%9wrong~numer~of~argumentsG@%5-FG6$FajoFban-FG6$FajoF\gpC
'>Fbo-Fgr6#Fajo@%2F[glFbp>FaqFe_el>Faq%&OrderG@$Fa]hl>F_`y-F\]m6$F_`yFb`el@'-FG
6$F[gxFfclC$>Ff_yF[gx>FgcyF[]l-FG6$F[gxF[joC$>Ff_y-F\[p6#F[gx>Fgcy-F_\qFadjl-Fe
n6#%Msecond~argument~must~be~a~equation~or~a~nameG>Faj[lF^crC'>Fbo-Fjfp6$FjnF[]
l@%Fjbr>FaqF^cr>FaqFacjl@$F^cjl>F_`y-F\]m6$F_`yFe_el@'-FG6$FajoFfclC$>Ff_yFajo>
FgcyF[]l-FG6$FajoF[joC$>Ff_y-F\[pF\cjl>Fgcy-F_\qF\cjlFddjl>Faj[lF[gx>F_s-Fj\]l6
%FIFf_yFgcy>F^]y-F\^]lFj_m@%0FgcyFham>Fc\l-F[y6$/Ff_y,&Ff_yF^rFgcyF^rFibw>Fc\l-
Fgr6#,$*&Ff_yFbcu-F[y6$/Ff_y*$Ff_yF\sFibwF^rF\s>Fin-Fe\lFj_m>F[r-Fg\m6%F_sFf_yF
jq>FfoF^`x>Fdhq-Fe_]l6$FfoFf_y>Fjix-Fd`]lF^hjl>F`x&&F``xF\\lFafl@$/-Fc_m6$F`xFj
qF[]lFi_[l-Ffal6&F[glF]dq%4indicial~polynomialGF`x>F^p-Ff_x6%F^`xFf_y&F[rFffl>F
f]n,(F^rF^rFaqF^r-F]el6$F[]l-Fbdo6%FboFf_yF[]lF\s-Ffal6&F[glF]dq%Gevaluating~si
mple~operator~up~to~orderGFf]n>F[r-Fefo6&F^pFf_yF[]lFf]n@$/%'_debugGFh]l-F\bs6$
F^pFf_y>FdbnF^`x-Ffal6&F[glF]dq%1diagonal~matrix:G7#-Fdq6$&Fdbn6$FdiqFdiq/FdiqF
]r>FgcnF``x>Fbe[lFdx>F[r-Fjfbl6%&-Fjf[l6$F`xF^]yFaflFjqF`gbl?&F]]lF[rFh]lC'>8@F
`]cl>Fhjx-F\]m6$F^]y-Fe`p6$Fe[[mFcg\l-Ffal6&F[glF]dq%2minimal~exponent:GFe[[m@%
4-FG6$Fe[[mFj`mC'-Ffal6%F[glF]dq%<right~hand~side~must~vanishG>Fer-Fhs6$Fc\lF[i
o>8B-F^w6$Fer-F[`p6$-Fe`p6$FerFfcl<#Ff_y-Ffal6&F[glF]dq%1conditions~foundGFj\[m
@%0Fj\[mFdxC%-Ffal6%F\pF]dq%Etrying~to~solve~homogeneous~equationG>8>-&F"6#Fifz
6+Fdbn-Fgr6#,&FgcnF^r*&Fe[[mF^rFdbnF^rF\sF\cqFf_y-F[y6$/Fjq,&FjqF^rFe[[mF^rF`xF
jqFhjxFaqF_`y>8AFj\[m>F\^[mFb`mC$>F\^[m-F^^[m6+FdbnFa^[m-Fgr6#*&-F`^n6$FjixFc\l
F^r)Ff_yFe[[mF\sFf_yFe^[mFjqFhjxFaqF_`y@$0F\^[mFb`mC$>Fj^[m-F`dy6%t+6$F2%#ibGF$
F$F$5-%&evalbG6#0Fe\m-%*intersectG6$-Fe`p6$-F\[pF`oFfclFS-Fa`[m6#0Fe\m-Fe`[m6$-
Fe`p6$-F_\qF`oFfclFSF$F$F$6#F\p&F\^[mFafl-Fe`p6#-Fhs6$Fc\lFban>Fj^[m-F`dy6$t+6#
F2F$F$F$-Fa`[m6#0Fi`[mFaa[mF$F$F$6#F^rFj^[m@$Fg_[mC$>8?-Fgx6&F^^[l-Fjhq6$-F`^n6
$F\iq&F\^[mF\\lFf_yFf_yFaq>Fbe[l6$Fbe[l7$-F[y6$/Ff_yFaj[l*&Fe_[mF^r-Fdem6#Fdb[m
F^rFj^[m-Ffal6&F^rF]dq%Ftime~for~computing~regular~solutions:G,&FaioF^rF^i[lF\s
7#Fbe[lF$6$FablF]jjlF$6#F\pF\`\lt+F]halF_halF$F$C%>FjqFham>FinFchal@%F]\oC$?(Fc
\lF^rF^rFinFh]l@$Ffial>Fjq-F]el6$Fjq-Fi\o6%FgialFSFVF[[blC&>FYFaial?(Fc\lF^rF^r
FinFh]l@$FfialC%>F[r-Fi\oF\jal>F^jalF_jal>FjqFajal@$Fcjal?(Fc\lF^rF^rFinFh]l@$F
fjal>FijalFjjalF[[blF$F$F$6#F[glFh`flt+6$F]hoF^[l6)F9%#R1GF,F2F[inF:FfhelF$F$C$
?(FinF^rF^rFh[oFh]lC(>Fbo-Fefm6%FI;F^rFd[o;FinFin>Fc\l-Fgbm6#-Fgx6%%(contentG-F
gx6$Fb[pFboFS>Fjq-Fhf[m6$-F`b[l6#&FboFjhblFS>Fer-Ff^l6#&Fc\lFjhbl?(F[rF\pF^r-Fe
[oFacmFh]lC$@%0-F`b[l6#&Fbo6$F[rF^rF^r>Fjq-F[e[l6$Fjq-Fhf[m6$Fjg[mFSF$>Fer-F[e[
l6$Fer-Ff^l6#&Fc\lF]h[m>&FaqF_dr-Fgx6$Ff^l-Fgr6#*(FjqF^rFerF\sFboF^r-F\hl6#-F`a
p6#-Fdq6$Fjh[m/Fin;F^rFh[oF$F$F$6#F\pFig[lt+6&Fh[lF^[lFfcmF_[l6*F9Fgf]lF^^elF*F
6F^[oFh]mFe_blF$F$C(>FjqFd[o>Fc\l*&FVF^rF[aelF\s>Fer-Fgr6#*&Fc\lF^rFIF^r>F_sFdx
?(FinF^rF^rFjqFh]lC&>Faq-Fj_\l6&-&F\o6#%$rowG6$FerFinFSFV.F[r@$2F[]lFaq>F[r-Fbi
al6$FjqF[]l>F_s-FjeclFjgu@$1F[]lFaq>&F_sFje\l,&Fj[\mF^rFYF\s-Fc\n6%-F[hxFj_mFVF
SF$F$F$6#F[glFfhxt+6$FjgoF_[l6%F^^mF[^mFj]mF$F$C&>Fbo-F``mFhfn>FinFb`m@$Fd`mC$?
&Fc\lFboFh]l@$3-FG6$F]^pFj`m2F\sF]^p>Fin7$F]^pFjdal>Fin-FcamF_drFinF$F$F$6#F\p%
*EigenringGt+6%F]j`lF^j`l'FeielFcan6'F6F5F9F+%$eigGF$F$@%Fi[flC'>FinFcj`l>FboF
\hcl?(Fc\lF^rF^rFboFh]l>F^jal-Fgr6#,&-Fhs6$-&%.LinearAlgebraG6#%1KroneckerProdu
ctG6$-Fhs6$FgialFafcl-Fhs6$-Ffbn6#-Fjfcl6$F^rFinFafclFLF^r-Ffbn6#-Fjfcl6$-F^hm6
#FgialFinF\s@%Fjip>Fjq-Fdgfl6$7#-Fdq6$F^jal/Fc\lF`uFS>Fjq-Fdgfl6%F]`\mFSFVO7#-F
dq6$-Fa]o6%FinFin-F_`\l6$FjqFc\l/Fc\l;F^r-Fi[oFf_mFb[glF$F$F$6#FgflFabqt+6&Ff]m
F^[lFcf\lF^[oFf_`lF$F$C&@$F]\o-F\hl6#,&-&F"6#F\ealF^b`lF^r-Fha\mF`b`lF\s>Fin,&-
Fha\mFeb`lF^r-Fha\mFhb`lF\s@$Fify>FY-Fix6#*&FboF^rFc\lF\sFinF$F$F$6#F[gl%5Matri
xColumnPrimpartGt+Fhe[m6)F9Fje[mF,F[inF2F:FfhelF$F$C$?(FinF^rF^rFh[oFh]lC(>FboF
_f[m>Fc\l-Fgbm6#-Fgx6%Fhf[mFboFS>Fer-F`b[lFcg[m>Fjq-Ff^l6#*&Fdg[mF^rFbc\mF^r?(F
[rF\pF^rFfg[mFh]lC$@%0-F`b[lFgh[mF^r>Fer-F[e[l6$FerF[d\mF$>Fjq-F[e[l6$Fjq-Ff^l6
#*&Fhh[mF^rF[d\mF^r>Fjh[m-Fgx6$Ff^l-Fgr6#*(FerF^rFjqF\sFboF^rF`i[mF$F$F$6#F\pFb
c]lt+Ffgo6<Fd_blF`tF_[oF+F,F-FfcoFfcmFdepFedmFicmFgf\l%#chGF?F3F9F:F6F7F[^mFc[l
FadmF_gzFhajl%$CisGFg]mF$F$CC>FjixFjn@$-Fa`sF`o-F\hl6#7#-Fjfp6$Fjix7#-Fdq6$&Fi^
[lFb]r/Fgcn;F^rFjix>F_`yFag\l>FerFi\]l>Fc\l-Fg\m6%FerFSFbe[l>Ffo-Ff_x6%F]^pFS&F
c\lFffl>F^p&F^^pFafl>Ff_y-F`dy6%FG-Fgx6$t+6#Ff]mF$F$F$F^]sF$F$F$6#F^r-F``m6$F^p
Fbe[lFj`m@$/Ff_yFb`m-F\hl6#7#Fjdjl>F_`y-F\^]l6#F]^p>F`x-Fe_]l6$F]^pFS>Faj[l,$-F
]el6$F[]lFfj\lF\s-Ffal6&F[glF]dqFgijl,&Faj[lF^rF^rF^r>Fc\l-Fefo6&FfoFSF[]lFeh\m
@$F\jjl-Fi\sFadv>F[rF]^p-Ffal6&F[glF]dqFcjjl7#-Fdq6$&F[r6$FgcnFgcnF]f\m>FaqF_^p
>Fhjx-Fhel6$,&*&F[rF^rFbe[lF^rF^rFaqF\sFS-Ffal6%F[glF]dq%Dcomputing~polynomial~
solution~spaceG-Ffal6&F[glF]dq%3bound~on~degree~isGFaj[l>Fc\l-F^^[m6+F[r-Fgr6#,
&FaqF^r*&Faj[lF^rF[rF^rF^r-Fjfp6$FjixF[]lFS-F[y6$/Fbe[l,&Fbe[lF^rFaj[lF\sF^pFbe
[lFef\lFeh\m<#Ff][l@$/Fc\lFb`mFfg\m>Fc\l-Fgr6#-F`^n6$-F`^n6$F`xF[rF]^p>F_s-Fhel
6$-Fgr6#*&)FSFaj[lF^r-F[y6$Ff^]lF^[`lF^rFS>F^i[l-F`dy6%Fg_p-Fe`p6#-Fhs6$F_sFban
Fi^[l>Ff]n-Fgx6%Fgao-Fgx6$Fb[p-Fgr6#,&-Fgx6%Fg^nF_sFSF^r-F`^nFf[qF\sFS>Fdiq-F^w
6$-Fhs6$-Fgx6%Fe^mFf]nFSF[ioF^i[l>F_s-F[y6$-F_amFciqF]by>Fdbn-F]el6#-F_am6#-Fgx
6&Fabq-Fhs6$F_sF[ioFSF[]l@$2FdbnF[]lC$>Fdiq-Fh_y6#-Fhs6$-Fgx6%Fe^m-Fgx6'Fc[sF_s
FSFdbnF\sFSF[io>F_sFi]]m>F^i[l7#-F_amF[`]l@%3Fgey-Fg_p6$FV.%&basisG7#-Fdq6$-F[y
6&-Fdq6$/&F^i[lFbdflF[]l/Fdhq;F^r,&FgcnF^rF\sF^r/&F^i[lFb]rF^r-Fdq6$/&F^i[lFf^p
F[]l/Fbo;,&FgcnF^rF^rF^r-F^_mF[`]lF]by/Fgcn;F^rF`a]m7#F_cyF$6%FablFi^[lF]jjlF$6
#FgflFfdlt+6%F][lF^[lFjb[l6*Fh[lFh]mF,%&sigmaGF7FcdmFhblFiblF$F$C'@$Fa_o-F\hl6#
-&F"6#%8mat_valuation_algebraicGF^`o@$4-FG6$FeclFj`m>Fecl"$+#@%F]\o>F[rFV>F[rFe
cl>FinFdhl@%5-FG6$F`flFj`m/F`flFhamF`flC$>FjqFdbbl@=FfhlC(-Ffal6%FcblFgir%6reco
mputing~valuationG>FerF`gl>Fbo-F]el6#-F_am6#-Fgx6&Fd`r-Fhs6$FeflF[ioFSF[]l?(F$F
^rF^rF$32FcglFbo2FboF[rC&>Fc\l-F_]l6&Fe`oFS,(FcglF^rFerF^rF^rF^r,&FboF^rFerF^r>
Fefl-Fgx6&Fheq-Fgr6#,&FeflF^r-F_\u6$F]gl-F_\u6$Fc\lFiflF^rFSFagu>FcglFbo>Fbo-F]
el6#-F_am6#-Fgx6%FjghlF`d]mFS@%3F]\o1F[rFbo>F`flFbfl>F`flFboFbo/F_dlFaevC'Fec]m
@%4-Fa`sFbjl>FboFhe]mC$>FboFifl?(F$F^rF^rF$31FboFeclFgf]mC%>FboFhz>Fefl-Ffel6$,
(FeflF^r-F_]l6&F`glFSFboFboF^r-F_]l6&FcglFSFboFboF^rFS-Ffal6&FcblFgir%3checking
~valuationGFbo>F`flFbo>FiflFboFbo/F_dlF`dlC'Fec]m@%Fff]m>FboFhe]mC$>FboFifl?(F$
F^rF^rF$F\g]mC%>FboFhz>Fefl-Ffel6$,(FeflF^rFdg]mF^rFfg]mF\sFSFhg]m>F`flFbo>Fifl
FboFbo/F_dlF_evC(Fec]m>Faq-FddlFaco>F_s-FddlF\ao@%4-Fa`s6#-Fgx6'Fc[sFeflFS,&Faq
F^rF_sF^rFifl>FboFhe]mC%>FboFifl?(F$F^rF^rF$3F]g]m-Fa`s6#-Fgx6'Fc[sFeflFSFhi]mF
boC'>FboFhz>F]glFgju>Fefl-Ffel6$,*FeflF^r-F_\u6$-F_]l6&F`glFS,&FboF^rF_sF\sF][^
m-F_]l6&FcglFSF_s,(FboF^rF\sF^rFaqF\sF^r-F_\u6$-F_]l6&F`glFSFaq,(FboF^rF\sF^rF_
sF\s-F_]l6&FcglFS,&FboF^rFaqF\sFh[^mF^r-F_\u6$F[[^mFf[^mF^rFS-Fhjp6$Q&M[3]=F$Fe
flFhg]m>FiflFbo>F`flFboFbo/F_dl.%&applyGC'Fec]m@%Fff]m>FboFhe]mC$>FboFifl?(F$F^
rF^rF$F\g]mC%>FboFhz>Fefl-Ffel6$,(FeflF^r-F_]l6&F]glFSFboFboF^rFdg]mF^rFSFhg]m>
F`flFbo>FiflFboFboF[ilC(Fec]m>FerF`gl>FboFjc]m?(F$F^rF^rF$Fcd]mC&>Fc\lFhd]m>Fef
lF]e]m>FcglFbo>FboFjc]m@%F_f]m>F`flFbfl>F`flFboFboFcjl-Fen6#%Kundetermined~valu
ation~of~localized~matrixGF^ilC'Fec]m>FboFhe]m?(F$F^rF^rF$2FiflFboC'@$/FboFham>
FboF]au>Fefl-Fgx6&Fheq-Fgr6#,&FeflF^r-F_]l6&FcglFSF]auFboF^rFSFjq>FiflFbo>FboFh
e]m@$2F[rFiflFggm@%F_f]m>F`flFbfl>F`flFbo-F]el6$FboF[rFcil-Fddl6%F`flFSF[rF[jl-
Fen6#%Hundefined~valuation~for~ramified~matrixGFhil-Fen6#%Kundefined~valuation~
for~substituted~matrixGFhjlC(Fec]m>FboFifl?(F$F^rF^rF$Fh]lC%@%2FboF[]l>Fbo-F]bu
6$FboF\p>FboFhz-FdduFggy@$5Fff]mF`f]mFggm>Fbo-F]el6$F[rFje]m@%F_f]m>F`flFbfl>F`
flFboFboF[[mC(Fec]m>FboFifl?(F$F^rF^rF$Fh]lC%-F_]l6&FIFSFboFbo@$F_a^mFggm>FboFh
z>FboFaa^m@%F_f]m>F`flFbfl>F`flFboFbo-Fen6#QCvaluation~of~unknow~data~structure
F$F$F]hlF$6#FgflF\emt+6(F*F^[lFfcmF_[lFicmFjcm68FhgoFcf\lF9F:F?FjgglF^^elF`tF6F
hcm%&mktabGFgcmFd`nFh[l%%subGGFc[l%&thetaG%#rkGFedm%#LmGF\^elF\^mF$F$C:>Ff_yF[]
l>F8F^\s@$5555FafyFeaw4-FG6$FSFfclFa`w4-FG6$FI-FK6$-%(ratpolyG6$FNFSFP-Fen6#%Bw
rong~number~or~type~of~argumentsG-Ffal6'F^rF[em%2enter~ratsuper~atGFV%3reductio
n~step~=~1G>FinFaio@'F\a`lC%>Fdiq-Fgx6$Fgbm-FgrFhegl>F[rFS>FboFg^]l3/Ff[nF^r0FV
FSC%>Fbo-Fi_pFg[n>Fdiq-Fgx6$Fgbm-F[y6$/FS,&FSF^rFboF^rFbfq>F[rFSC%>FboF[]l>F[rF
V>FdiqF_em>FfoFd[o@$FapC$>Fcy-F]]nFbal>Fjy-FgrF^hv>F_sF^r@%/F[rFham>FaqF^r>Faq-
Fg^n6$F[rFS-Ffal6%FcblF[em%Ecomputing~pole~order,~Glambda~and~mkG@%Fify>Fdiq-F]
`n6-FdiqFSF[rFfoF_sFYFf`nFahm.Ff]n.F^p.Fgcn>Fdiq-F]`n6/FdiqFSF[rFfoF_sFYFf`nFah
mFhg^mFig^mFjg^mFcyFjy-Ffal6&FgflF[emFjalFjq-Ffal6&F\pF[emF]bl-F_am6#-Fhs6$Ff]n
Fban-Ffal6&FgflF[emFfan-Fhs6$F]]lFban-Ffal6&FgflF[emF[bnF^p>&F`xF\\lF^p>FhjxFdx
?(F$F^rF^rF$32F^rFjq2F_sFjqC.>%#GGG-FgrFb]r>FeenF^gfl>Fdhq-Fh^m6#-Fc\n6%-F[hxFb
]rF[rFS-Ffal6(FcblF[em%0reduction~step:GF_s%'theta:GFdhq@%/FdhqF[]lC$>Fjix-Figm
6&-Fefm6%FgcnF_b_l;F^rF^pFSF[r.F_`y@%2FjixF^pC%-Ffal6%F[glF[em%?columns~are~lin
early~dependentG-Ffal6&FcblF\em%.kernel~dim~=~G,&F^pF^rFjixF\s@%Fify>Fdiq-&F"6#
Fdfv6-FdiqFSF[rFfoF_s,$FjqF\sF]]lFf]nF^pF_`y%#iPG>Fdiq-Fi[_m6/FdiqFSF[rFfoF_sF
\\_mF]]lFf]nF^pF_`yF]\_mFcyFjyC%-Ffal6%F[glF[em%Acolumns~are~linearly~independe
ntG>Ff_y,&Ff_yF^rF^rF^r@%Fify>Fdiq-&F"6#Fccm6+FdiqFSF[rFfoF_sF]]lF^pFYFgcn>Fdiq
-Fj\_m6-FdiqFSF[rFfoF_sF]]lF^pFYFgcnFcyFjyC)@$2F_s,&FjqF^rF\sF^r>Fdhq-Fh^m6#*&F
dhqF^r)FY-F^iw6$FdhqFYF\s@$3Fd_fl0FjqF^rC%>Fc\l-&F\o6#F_i^l6$-Fefm6%Fgcn;F^rFbe
blFf^_mFY>Fhjx7$Fc]_m-Fjdel6$-Fc\n6%-F[y6$/FY*&FaqF^rFYF^rFc\lF[rFSFY-Ffal6(F\p
F[em%<slope~(Katz~approximation):GFc]_m%+invariant:G&-F_am6$F^r7#FhjxFafl>Fhjx6
$7$,(FjqF^rF_sF\sF\sF^r-Fjdel6$-Fc\n6%-F[y6$/FY,$F`__mF\sFdhqF[rFSFYFhjx-Ffal6(
F\pF[em%'slope:GF\`_mFd__mFe__m@$/F_sFc]_mFggm>F_sFjbal-Ffal6&F^rF[em%=increase
~reduction~step~to:~GF_sFg_n@%Fify>FdiqFfg^m>FdiqF\h^mF^h^mF`h^mFfh^mFjh^m>F\ip
F^p@'FaexC%-Ffal6&F\pF[em%6matrix~is~analytic~atGF[r@%/FboFg^]l>Fjq,&FjqF\sF\pF
^r>FjqF\\_m-F\hl6#7&F^\s7$FjqFfo)FYFfoF[]l/FjqF^rC(-Ffal6&F\pF[em%>matrix~is~re
gular~singular~atGF[r-F]ew6&FdiqFSF[rFjg^m>Fdhq-Fa^_m6$FgcnFY>Fhjx7$F[]l-Fjdel6
$-Fc\n6%-F[y6$F___mFdhqF[rFSFY-Ffal6(F\pF[emFg`_mF[]l%3indicial~equation:GFe__m
>FbeblFfo>&F]]lFf_m,&FfoF^rF^pF\s@&Fha_mC%@$F^cjlC$>Fcy-Fgx6$Fgbm-F[y6$Ff^]l-Fd
emF^hv>Fjy-Fgx6$Fgbm-F[y6$Ff^]l-Fdem6#Fjy>Fdiq-Fgx6$Fgbm-Fgr6#,$*&FSFbcu-F[y6$F
f^]l-FdemFciqF^rF\s>Fjq,&FjqF^rFbcuF^r0FboF[]lC$@$F^cjlC$>Fcy-Fgx6$Fgbm-F[y6$/F
S,&FSF^rFboF\sFcd_m>Fjy-Fgx6$Fgbm-F[y6$F`f_mFid_m>Fdiq-Fgx6$Fgbm-F[y6$F`f_mFde_
m-Ffal6&F^rF[em%'time:~GF`]p>FhjxFh__m@%/FhjxFb`m7&F^gfl7$F\\_m-F_am6#Fhh^mFb`m
Ff_y7&F^gflFcg_m&Fa^ilFaflFf_yF$6'%#PTG%&subGTGFei^mFeenF8F$6#FcpFi\mt+F\\m65Fh
goF9F?FicmFjcmF`t%#k2GF6FhcmF][lF_\mF`[lF*Fc[lF[c^m%'theta0GFc_blF7Fi^alF$F$C<@
$Fcb]m>FeclFgb]m>FinFaio@$555F_bjlF^cjlFfc^m4-FG6$FVFfclFad^m>FaqFham>Ff]nFfiw>
F_s-Fe\lFicp>Fer-F]]nFj_m>FjqF`i_m>F[rF^r>Fc\lFham>F]]lFb`m@%FifyC%@$-Fg_p6$FYF
_`blC$>F]]l-F_\q6#-F_am6#-F`dy6%Fg_pFYF_`bl@$F^]dlC$>Fc\l,$-FddlF]j\lF\s>F[rFg[
dl@$-Fg_p6$FYF`]m>Faq-F_\q6#-F_am6#-F`dy6%Fg_pFYF`]m>F`xFY>F`xFe\m>FinFaio>Fc\l
Fej_m>Fjix-F\^]lFicp-Ffal6&F^rFhal%2reduction~step~=~GF[r>FdhqFdx>FfoF[]l>FdbnF
[]l@$/F]]lFb`m>Fc\lFham?(F$F^rF^rF$31F[r-F]el6$F]gglFaq2FfoF_sC%?(F$F^rF^rF$3/F
dbnF[]l2F^rFc\l@%/F[rF^rC0@$F_\`mC$>Fc\lFej_m-Ffal6&F[glFhalF[cblFc\l@$1Fc\lF^r
Fggm>Fdiq-&F"6#Fd[z6&Ff]nFSF[rF_hbl-Fjjcl6&FdiqFjixF\edlF[[_m>Ffo-F^_mF^`y>Fjq-
Fgx6$Fix-Fgr6#*$FerF\s-Ffal6&F^rFhal%6Moser~reduction:~r~=~GFfo>Ff]n-Fh\x6'Ff]n
FSFerFjqF[]l>Fc\lFej_m>Fdiq-Fgx6$Ff^l-&F"6#Fjz6'Ff]nFSF[rFVF_hbl>Fdbn-Fix6#-%)P
rimpartG6$-FjfoFciqFV-Ffal6&FcblFhalFaj^mFdbn@$0FdbnF[]lFggm>Ff]n-&F"6#%'l_qtcd
G6+Ff]nFSF[rFfoF_hblFVFdiqFjixF`xC*>Fdiq-Fh]`m6&Ff]nFSF[rF]]l>F^p-&F\o6#.%'kern
elGFciq@%0F^pFe\mC%-Ffal6%F[glFhalFa[_m-Ffal6&FcblFhalFd[_mFccq>Ff]n-&F"6#Fiicl
6*Ff]nFSF[rF]]lF^p.FfoFjixF`xC$>FfoF_s-Ffal6%F[glFhalFd\_m>Fdiq-F`_`m6'Ff]nFSF[
rFVF]]l>FdbnFd_`mFj_`m@$F]``mFggm>Ff]n-F```m6+Ff]nFSF[rFfoF]]lFVFdiqFjixF`x@%F]
]`m@%Fe]`mC(>Fgcn-Fix6#-Fg_`m6$-F\i^l6$-Fefu6%Ff]nFSF\sFVFV-Ffal6%F^rFhal%;syst
em~is~regular~singularG>Fdhq7$F[]lFgcn-Ffal6(F^rFhalFg`_mF[]lFec_mFgcn>F]]lF]fb
lFggmC'>Fgcn-Fix6#-Fg_`m6$-F\i^l6$-Fefu6%Ff]nFSF\fwFVFV>Fdhq6$7$,&Fc\lF^rFbcuF^
r-Fix6#-Fg_`m6$-F[y6$/FV,$FVF\sFdbnFV7$F]ggl-Fix6#-Fg_`m6$FgcnFV-Ffal6(F^rFhalF
g`_mF]gglFd__m&-F_am6$F\p7#FdhqFafl-Ffal6(F^rFhal%'slope~GFhd`mFd__m&-F_am6$F^r
F[f`mFafl>F]]l7$Fc\lFfoC%>Fdhq6$7$,(Fc\lF^rF[rF\sF\sF^rFid`mFdhq-Ffal6(F^rFhalF
^f`mFhf`mFd__mF_f`m>F]]l7$F__[l,&FfoF^r-Fhs6$7#-Fdq6$&F]]lFf^p/Fbo;F\pFg[dlFdtF
\s@$F[]`mC%>F[rF^cm>FdbnF[]l-Ffal6&F^rFhalF]a_mF[r@$/FfoF_s?(FboF[rF^rF]gglFh]l
C%-Ffal6&F^rFhalF]a_mFbo>Fdhq6$7$,(Fc\lF^rFboF\sF\sF^rF^rFdhq-Ffal6(F^rFhalFg`_
mFeh`mFd__mF_f`m-Ffal6&F^rFhalF^g_mF`]p@$1Fc\lF[]lC$-Ffal6%F^rFhal%7system~is~n
on-singularG-F\hl6#7%Ff]n7#7$F[]l)FVF_s7#F[]l-Ffal6&FcblFhal%-rho-polygon:GF[f`
m7%Ff]nF[f`mF]]lF$6#FeclF$6#F[glFbj\lt+6&F*%&xeqx0GFajpF]\m6+F^[lFcf\lF_\mF_[o%
%tmp1G%#epGF+F7FhgoF$F$C->FfoFaio@$Fdio-Fen6#%Efirst~argument~must~be~square~ma
trixG@'Fgc^mC$>FinFS>FboF[]l-FG6$FSF[joC$>Fin-F\[pF\^q>Fbo-F_\qF\^qFddjl@$Fhh_m
-Fen6#%>third~argument~must~be~a~nameG>Fc\lFe\m@$Fb\m@%-FG6$FYF[io>Fc\lFY-Fen6#
%;4th~argument~must~be~a~setG>F_s-F\]m6$-F\]m6$Fag\l-Fe`p6$FboFcg\l-Fe`p6$Fc\lF
cg\l@%4-Fg_p6$Fc\l.%&MoserGC)>Faq-Fj\]l6%FIFinFbo>Fjq-Fb]]l6%FaqFinFc\l@$30F[\p
Fjn4-Fg_p6$Fc\lFij\l-Fen6#%Pnot~all~solutions~computed.~Please~report~bug~!G>Fj
qFe\p@$F\jjl?&F[rFjqFh]l-F\bs6$F[rFinFjq>Fer7#-Fdq6$7$/Fin-F[y6$/FinFV-FaaoFd^a
m-F[y6$F]_am-Fb`oFd^am/F[rFjqC(>Faq-&F"6#%)localizeGFc]am-F^cl6#.%0RationalExpP
artG>F_s-Fe`p6$-Fhs6$FaqF[ioFcg\l>Fer-&F"6#F\`am6&FaqFinFVF_s>Fjq7#-Fdq6$7$F^`x
-Fh^m6#**-Fg^n6$F``xFVF^r-Fjj_l6$-F_\qF`j^lFVF\s-Fc_mFbaamF\s)FV,&FdaamF^rF\sF^
rF\s/F[rFer>FerFjq@%Fh^^m>Fer-F[y6$/Fin*$FinF\sFer>Fer-F[y6$/FinFcfclFer-Ffal6&
F^rF]dq%Ftime~for~computing~exponential~parts:G,&FaioF^rFfoF\sFerF$6#F]jjlF$6#F
[gl%4PolynomialSolutionsGt+F\helF^helF$F$C$@%Fjip>FinF]iel@%F`ielC$>FinFdiel>Fb
oFiiel@%F[jelC$>FinF]jel>Fc\lFiiel@%FbjelC$>FinFfjel@%F\jelC$>Fc\lFiiel>FboF\[f
lC$>Fc\lF\[fl>FboFiielF`[fl@%Fd[flC&>FaqFcj`l>F_sF\hcl@%Fh\fl>Ffo-Fgbm6#-F_]\lF
]]fl>Ffo-Fgbm6#-F_]\lFa]fl@%Fc]flOFe\mC'>F^p-Fgx6$f*Fa_pF$Fc_pF$F][rF$F$F$Ffo@%
Fa^fl>F`xFc^fl>F`xF^_fl>F]]lFd[[l@%F\ix@%Fd_flC$>Ff]nFg_fl>FdiqF_\rC$>Ff]nF_`fl
>FdiqF_\r@%Fh\flC$>Ff]n-Ff`fl6$-F`ap6#-Fdq6$-Fgx6$f*Fa_pF$Fc_pF$FaaflF$F$FcaflF
foFdaflF^]fl>FdiqFgafl>Ff]n-Ff`fl6$-F`ap6#-Fdq6$-Fgx6$f*Fa_pF$Fc_pF$FaaflF$F$Fc
aflFfoFdaflF^]fl@%Fd_fl@%Fh\flOF\cflOF]cflC*>FdbnFccfl>FgcnF\dfl>FdhqF^dfl>Fjix
Fadfl>F_`yFddfl?(FhjxF\pF^rF_sFh]lC$>F[eflF\efl@$Fh\fl>FafflFbffl@%Fc\fl>F^]y-&
F"6#FibamFfgfl@%Fa\fl>F^]y-FihamFfhfl@%F\\fl>F^]y-FihamF[ifl>F^]y-FihamF_ifl@%F
aiflOFe\m@%Fh\flOFeiflOF\[glFb[glF$F$F$6#FgflFjhqt+F[\sF$F$F$-Fgx6$t+F`g\mF$F$F
$*&-Fgao6%Fa[pFjdglFixF^r-Fgao6%Fc``lFjdglFixF\sF$F$6$F^[lFS6#F^rF^\sF$F$F$6#F
\pFdjwt+6(F*F+FedmF_[lFdepFicm60F6F2%&dummyGF,F-FigvF][lF/%#iFGF`tF5F_jclFjfvFi
^alF$F$C'>FinFjn@$/-Fc_m6$FVFYF[]lC&>F`x,&FinF^r-Fc_m6$-Fjfo6#-Fgr6#,&FIF^r*&FY
F^rFSF^rF\sFYF\s-Fbjw6(FSFI)FYF`xFYFahmFhg^m@$FapC$-F[w6$FcyF]]l-F[w6$FjyFf]n-F
\hl6#-Fgx6$Fix-F_\u6$-F_\u6$F]]lFg[bmFf]n@%0-Fix6#-FjfoF\^qF[]lC)>Fer-F]doF\^q>
Faq-Fgr6#-F[y6$/FY-Fgr6#-F`^n6$FerFIFV@$4-FG6$FaqFK>Faq-Fgr6#*&F^drF^rFaqF^r>F_
s-F^[dl6#-Fhs6$7#-Fdq6$-Fhs6$F^pFban/F^p-Fj``mF`amFK>F_s-&F"6#%+colechelonG6%F_
sFc\lFdiq>Ff]n-Ff`hl6$-%+complementG6$F_sFdiqF_s>F]]l-Fgx6$Fix-F_\u6$-F]doFicpF
erC1>FboF[]l?(F$F^rF^rF$/-Fix6#-Fjfo6#-Fgr6#,&FIF^r*&FboF^rFSF^rF\sF[]l>FboFhz>
Fer-F]doF[abm>Fjq-Fgr6#-F`^n6$FSFer>F[r-Fb[p6#-F[y6$/FY,&FboF^r*$FYF\sF^rFV>Faq
-Fgr6#-F[y6$/FYFjqF[r@$F_^bm>FaqFc^bm>F_sFg^bm>F_sFc_bm>F_s-Ff`hl6$-&F"6#F\`bmF
]`bmF_s>Ffo-Fgr6#*$F_sF\s>Fjq-Fgr6#-F`^n6$-F`^n6$FfoFjqF_s>F`x,&FinF^rF^[bmF\s>
F]]l-Fgx6$Fix-Fgr6#-F`^n6$-F_ar6$-F]]nF]^l-F]do6#-F\z6%Fjq;FaiuFinF[ebmFfo>Ff]n
-Fgx6$Fix-Fgr6#-F`^n6$FerF_s@$FapC$F`\bmFb\bm-Fgx6$Fix-F_\u6$F]]l-F_\u6$Fg[bmFf
]nF$F$F$6#FcpFb]_lt+6*F*F^[lFhajlF^h_mF2FadmFg]mF]\m69F]dmF6F_gzF_\mF`[lFc_blF9
Fd_blF_[oF+Ff]mF___lFedmF]hoF7Fe_blF]\nF_\nFbeyFcdmFd[lF8F`^mF$F$C2>Fc\lFdx@%2F
\qFbp>FerFjw>FerFe\m>FjqFb`m@$F]`blC$>FjqFc`bl>FerFj`bl>F[rFb`m@$F]ablC$>F[rFca
bl>FerFjabl>Fin,$FgfnF\s@$2F]]_lFV-F\hl6#7$Fb`mF[]lFeg^l>FboFd\l>FdbnFdbbl>Fc\l
Fdx>FgcyF[]l-Ffal6%F^rF]g^l%Jenter~formal~reduction~after~ramificationGFibbl@%F
h_[lC&Feibl>Ffo-Fg\m6&FIFSFhwF^]m-Ffal6&F[glF]g^lF]jbl&FfoFafl-F\hl6#-F`]_l6*&F
foF\\lFSFVFYFcyFjyFhw-F\]m6$Fer<$/F_`blFji^l/F_ablFdhbmC(Fjjbl>Ff]n&&F[r6#F[a[l
Fafl-Ffal6&F[glF]g^lFc[clFf]n@$/-%%igcdG6$F]]_lFjyF^r?&F`x&-Fjf[l6$-Ff^l6#*&Ff]
nF^r)Fhw-F^iw6$Ff]nFhwF\sFdbnFaflFh]lC0-Ffal6(F[glF]g^lF[\clF]i^mF\\cl&F`xFafl>
Fdiq-Ff^l6#-F[y6$/Fhw)Fhw*$FjyF\s-F[y6$/Fhw*&FjyF^rFhwF^rF]i^m-Ffal6&F[glF]g^l%
3rational~slope~k~=G*&F]]_lF^rFjyF\s-%'igcdexG6&F]]_lFjyFchq.Fjix>F^]y-Fi_p6$-F
[y6$/FcyF^rFdiqFhw>Fhjx)F^]yFdhq-Ffal6&F[glF]g^l%8smart~ramification,~c~=GFhjx>
Fdiq-Fix6#-F[y6$/FcyFhjxFdiq-Ffal6&F[glF]g^l%0reduced~factor:GFdiq>&F_`yF\\l,&F
hwF^r*&)F^]yFjixF^rFjyF^rF\s>&F_`yFaflFhjbm>Ff_y-&F"6#%)mat_subsG6%F]]cmFIFS>Fg
cy,&FgcyF^r*(-Fc_m6$FdiqFhwF^rFg]cmF^rFjyF^rF^r@%/-Fc_m6$Fb]cmFhwF^rC$Ff\cl@%2F
^rFg]cmC*-Ffal6&F[glF]g^lF^]clFg]cm>F]]l-Fi_pFf^cm>F_s*&F]]lF^rFc]clF\s-Ffal6&F
[glF]g^lFf]clF_s>F^p-Figbl6%Ff_yFSF_s-Ffal6'F[glF]g^lF\fblFhi_l)FhwFg]cm>Ffo-Fa
jbl6%-F`fbl6'F^pFSFhi_lFh_cmFhwFSFer>Fc\l6$Fc\l-F_amFbalC%-Ffal6'F[glF]g^lF\fbl
Fhi_lFb]cm>F^p-F`fbl6'Ff_yFSFhi_lFb]cmFhw>Fc\l6$Fc\l-F_am6#-Fajbl6%F^pFSFer@%3F
f_cl2F[]l,&FboF^r*&Fg]cmF^rFe^cmF^rF\sC&Fj_cl-Ffal6'F[glF]g^lF\fblFhi_l)Fb]cmFg
]cm>F^p-F`fbl6'Ff_yFSFhi_lFeacmFhw>Fc\lFh`cmC+Fe`cl>F]]lF^_cm-Ffal6(F[glF]g^lF[
aclF]]lF\aclFg]cm>F_sF`_cmFa_cm>F^pFd_cmFf_cm>FfoFj_cm>Fc\l6$Fc\l-Fdq6$7$-Fb`o6
$&FaqFaflFSFhbcm/FaqFfo@$2F[]l-Fc_mFdjbm@&/FcjbmF\ccmC&Fdacl>Ff]n&F^`xFafl-Ffal
6&F[glF]g^lF\bclFf]n@%F[ccmC'F`ccl>F^p-F`fbl6'FIFSFjqFf]nFhw>Ffo-F`]_l6*F^pFSFV
FYFcyFjyFhwFer>Fgcy,&FgcyF^rFdhbmF^r>Fc\l6$Fc\l-F_am6#FihbmC$FhcclFcgbm2F[]lFcj
bmC'-Ffal6'F[glF]g^lF\fblFhi_lFbjbm>F^p-F`fbl6'FIFSFhi_lFbjbmFhw>FfoFjccm>FgcyF
]dcm>Fc\lF_dcm-F\hl6#7$F_hblFgcyF$F$F$6#Few%:HyperexponentialSolutionsGt+6$F]j`
lF^j`l65F6F5F_\nF`tF9FcepF:Ff]mF+F]hoF?FfhelFghelFcgwFhhelFicmF8FihelFigoF$F$@%
Fi[flC&>FinFcj`l>FboF\hcl>Fc\l-&F"6#FdgoFa]fl@%/Fc\lFe\mOFe\mC$@%/FboF^r?(FerF^
rF^r-F^_mFacmFh]l>&F_`yF_bl-Fgbm6#-Fh^m6#*&-Fh]p6#-%$intG6$&&Fc\lF_blF\\lF^]flF
^r&FbgcmFaflF^r?(FjqF^rF^rFdfcmFh]lC.>&F[rFf_m-Fh]p6#-F_gcm6$&Fg\qF\\lF^]fl>Fgi
p7#-Fdq6$-Fgr6#,&-Fhs6$F[iclFLF^r*&-Ff^l6#*&-Fg^n6$FggcmF^iclF^rFggcmF\sF^rF__
\mF^rF\s/Faq;F\pFbo>Figp-Fgr6#&Fg\qFafl>Feip-Fi[oF]\q>Fhc_m-F`ap6$FigpF__\m>&Ff
]nFf_m-Fca_l6#Fhc_m>&FdiqFf_m-Fefm6%Fhicm;F^rFeipF_jcm>&FdbnFf_m-Fh]n6#F\jcm>&F
gcnFf_m-Fefm6%FhicmF_jcm;,&FeipF^rF^rF^r,&FinF^rFeipF^r>&FdhqFf_m7#-Fdq6$-Fgx6$
Fgbm-Fgr6#-F`^n6$-F`^n6$FajcmFejcm,&-F`^n6$&Fgip6#,&FaqF^rF\sF^rFigpF^r-Fgr6#-F
gx6%Fg^n-FgrF]\qF^iclF\sF\icm>Fjix-&F"6#Faecm6$F\[dm7#-Fdq6$F^iclF\icm@%/FjixFe
\mOFe\m>&F_`yFf_m-Fgx6$Ff^l-Fgr6#-F`^n6$*&FggcmF^rFigpF^rFjixO-F`ap6#-Fdq6$-Fgb
m6#-Fgx6$Fh^m-Fgr6#-Fgbm6#F_]dm/Fjq;F^rFdfcmFb[glF$F$F$6#F\pFe_ot+F`hlFahlF$F$C
$>FinFdhl@4Ffhl-Fc_o6#Fe`oF[ilF[_dmF^il-Fc_oFbilFcil-Fc_oFgilFhilF]_dmF[jlF]_dm
Fhjl@%Fabo-Fc_o6#Febo-Fc_o6#FjboFcjlFaguF]co-Fc_oFgjlF$Fh[mF$6#F^r%)mat_pvalGt+
F[[o6-Fh]mF2Fd[lF[inF9F:F6F5F]\nF_\nF___lF$F$C+@$F]\o-F\hl6#-F]el6#-F_am6#-Fgx6
&FdaelFd`vFSFV>F_sFjn>FaqFcfp>FinFham>Fc\lFdx?(FjqF^rF^rF_sFh]l?(F[rF^rF^rFaqFh
]lC$>Fer-Fdael6%FdgpFSFV@&2FerFinC$>FinFer>Fc\l7$FjqF[r/FerFin>Fc\l6$Fc\lFbadm>
FY-Fa]o6%F_sFaqF[]l?&FboF_hblFh]lC%>F`x-Ff^l6#*&)FVFc]lF^r&FI6$FefnF_gnF^r-Fj[n
6'-F`b[lF]^lFVFSFja`mFig^m>&FYFabdm-Fc\n6%*&-Fb[pF]^lF^rFfoF^rFVFSFinF$F$F$6#F[
glF[jalt+Fe_`l6'Fh]mFg_`lFh_`lF]\nF_\nF$F$C&@%Fify>Fin,&-&F"6#Febdl6&-Fh^m6#Fa[
pFSFVFfb`lF^r-Fccdm6&-Fh^m6#Fc``lFSFVFe`nF\s>Fin,&-Fccdm6%FfcdmFSFVF^r-Fccdm6%F
jcdmFSFVF\s@$55Fi`^lF\a`lFj`^lC$@$Fify>FYFcb\mFbfy@$FifyC$-Fj[n6'Fc\lFVFSF\edlF
f`n>FY-Fc\n6%*&FboF^rFerF^rFVFSFinF$F$F$6#F[gl%(l_superGt+F$F$F$F$F$F$F$F$6#F[]
lFb``mt+6*F+F^[lF`tFhcmF`[lF_[lFc[lF763F6Fei^mF]dmF^dmF_dmFadmF2F9F:F][lFedmF-F
gdmF]\n%%invUGFi^alF*F$F$C0>FdhqFI-Ffal6%F\pFhalF]em>Fin-Fe\lFbdfl>Fc\lF[]l>Fer
-F`emF`]n?(F$F^rF^rF$Fh]lC6>Fbo-F[y6$/FjyF[]l-Fhel6$FhwFS-Ffal6&FcblFhalFbfmFc
\l>Fjq-F\z6%Fbo;F^r,&FinF^rFc\lF\sF^gdm@%/FjwFe\m>F[r-Fd`y6#-F\z6%FjqF_idlF^gdm
C$-Fjjcl6&-F^hmFdgdmFjw.FhjamFjg^m>F[r-F^_mFb]r@$2F[rFerC$-Ffal6%FcblFhalFfgmFg
gm>F^p-Fgx6%Fhs-&F\o6#F]a`m6#-F^[dlFf_mFban@$/F^pFe\mFggm>F`x-Fjjcl6&-F^hm6#-Fh
s6$7#F[]elFKFjwF[hdmFjg^m>F`x-F^hmF]^l>Faq-F[oF]^l-Ffal6&FcblFhalFd[_mFaq?(F_sF
aqF\sF^rFh]l@$/&F`x6$F_s,*FinF^rFc\lF\sFaqF\sF_sF^rF[]lC'>FfoF`jdm?(F$F^rF^rF$/
&F`xFa^lF[]l>FfoFfim-Ffal6'FcblFhal%-permute~rowsGF`jdmFfo>F`x-F]hn6%F`xF`jdmFf
o>Fdhq-&F"6#Fgen6&FdhqFSF`jdmFfo>FdbnF^dr>Fdbn-Fhct6&F`xFdbn,*FinF^rFc\lF\sFaqF
\sF^rF^rF^r>Fdiq-F]doFiex>Fdhq-Fh\x6'FdhqFSFdiqFdbnF[]l@$/Fc\l,&FinF^rFerF\sFgg
m-Ffal6%FgflFhal%5recomputing~G-matrixG>Fhw-Fgx6$Ff^l-F`_`m6'FdhqFSFVFjyFcy>Fc
\lF\ihl>Fbo-Fgx6$Ff^lFdfdm-Ffal6&FcblFhalF_bnFc\l-Ffal6(FgflFhal%Bapplying~diag
onal~transformation:GFY,(FinF^rFYF\sFc\lF\sFc\lFYFc\l>F]]l-F_ar6%-Fdq6$FS/F_sFe
_\l-Fdq6$F^r/F_s;F^rF`]em-Fdq6$FS/F_sF`gm>Ff]n-F_ar6%-Fdq6$Fg^]lFf]emFg]em-Fdq6
$Fg^]lF]^em-Fh\x6'FdhqFSF]]lFf]nF^rF$F$F$6#FewF^^]lt+F`hlFahlF$F$C$>FinFdhl@8Ff
hl-F\^]lF\_dmF[ilF\_emF^ilF`glFcil-F\^]lFgilFhilF`glF[jlFaguFhjl@%Fabo-F\^]lFa_
dm-F\^]lFc_dmFcjlFcglF]co-F\^]lFgjlF^[mFa_emFb[mFa_emF$Fh[mF$6#F^rFd]]lt+Ffgo6F
F_[oFaj`mF+F,Fg]mF_gzFc_blF\^elF`[lFedmF`tF9Ff]mFd_blF7%%rootGF_\m%&opts2G%&opt
s3G%(invars1G%'step_kG%$mulG%%rmulGF[inFf]qFadm%&rlistG%)r_testedG%*slope_setG%
+slope_set1GF]\nF_\nFhajlF^h_mF`^mFceyF$F$C9>F[rFdx>F^i[lFham@%Fgey>FdhqFV>Fdhq
Fe\m>FhjxF[]l@$-Fg_p6$FdhqF`]mC$>Fhjx-F_\q6#-F_am6#-F`dy6%Fg_pFVF`]m>Fdhq-F^_p6
%Fg_pFdhqF`]m@%-Fg_p6$Fdhq.%'slopesGC$>Fe[[m-F_\q6#-F_am6#-F`dy6%Fg_pFVFhaem>Fd
hq-F^_p6%Fg_pFdhqFhaem>Fe[[m<#-Fdq6$F]]l/F]]l;F[]l-Feez6$F[]l,&--F^cl6#.FfdlFhf
nF\sF\sF^r>Fj^[mFe[[m>FfoFb`m@$-Fg_p6$FdhqF_`blC$>Ffo-F_\q6#-F_am6#-F`dy6%Fg_pF
VF_`bl>Fdhq-F^_p6%Fg_pFdhqF_`bl>Ff_yFd\l@$-Fg_p6$Fdhq.%)max_solsGC$>Ff_y-F_\q6#
-F_am6#-F`dy6%Fg_pFVFfdem>Fdhq-F^_p6%Fg_pFdhqFfdem>Fbe[lF^r@$-Fg_p6$Fdhq.%)test
_ramG>Fbe[l-F_\q6#-F_am6#-F`dy6%Fg_pFVFgeem-Ffal6'F^rF]dq%5enter~system~Newton:
GFf_y%=exponential~parts~to~computeG>Fj^[mFe[[m>FaqFb`m>F_sFdx>FdbnFb`jl>8FF[]l
>Fc\l-&F"6#FgiwF`o?(F$F^rF^rF$Fh]lC+@$FeamC)>Fhjx,&FhjxF^rF^rF^r-Ffal6&F\pF]dq%
Jcalling~super~reduction,~reduction~step~=GFhjx>Fin-Fg\m6&Fc\lFSFjq-F\]m6$-F\]m
6$Fdhq<#/F`]mFhjx<#/F_`blFfo>Fc\lF_dl>FaqF`fl>FfoFefl>Fe[[m<#-Fdq6$F]]l/F]]l;-F
]el6#-F_am6#Fe[[m-F]el6$-FeezF\iem-Feez6$F[]l,&-FddlF`edlF\sF\sF^r-Ffal6&F[glF]
dq%5partial~rho-polygon:GFaq-Ffal6&F^rF]dq%:searching~for~slopes~k~inGFe[[m>F`x
&F_[clF\\l>F^pF^[cl>F_s6$7$F`xF^pF_s@$-F\ev6$F`xFe[[mC%>Fe[[m-F[`p6$Fe[[m<#F`x-
Ffal6&FcblF]dq%+slope~is:~GF`x@%/F`xF[]l@$2F[]l-Fc_m6$F^pFjqC'-Ffal6(F[glF]dq%)
slope~=~GF`xFjhjlF^p>Fin&-Fjf[l6$F^pFdbnFafl@$-Fg_p6$FdhqFgj\l>Fin-Fjfbl6%FinFj
qF`gbl>8G-Fc_m6$-FaaoF`edlFS?&Ff]nFinFh]lC*>Fgcn-Fi_p6$&Ff]nF\\lFjq@$2F^r-Fc_mF
j\fm>Fdbn-F\]m6$Fdbn<#Fgcn-Ffal6&F[glF]dqF]\[mFgcn-Ffal6&F[glF]dq%9exponential~
parts~found:G*(&Ff]nFaflF^rF^]fmF^rFb\fmF^r>Fer-Figbl6%Fc\lFS*&FgcnF^rFSF\s-Ffa
l6%F[glF]dq%<converting~to~simple~systemG>F[r6$F[r-Ff_x6%FerFSFfo>Fifem,&FifemF
^rFh]fmF^rC$>F^p-Fix6#-Fegq6#*&F^pF^r)Fjq-F^iwFb[fmF\s@$F`[fmC%@%3-Fg_p6$F^pF2/
-Fiibm6$Fbe[lF`xF^rC&>F^p-Ff^l6#-F[y6$/Fjq*(FjqF^rFbe[lF^r)F2*&F`xF^rFbe[lF\sF
\sF^p>F^p-Ff^l6#-F[y6$/Fjq)Fjq*$Fbe[lF\sF^p-Ffal6(F[glF]dq%2rational~slope~=~GF
b`fm%<reduced~Newton~polynomial~=GF^p>Fjix-F^_p6%Fg_pFdhqFgeemC$-Ffal6(F[glF]dq
%1integer~slope~=~GF`x%4Newton~polynomial~=GF^p>FjixFdhq@$4-Fg_p6$FjixFgeemC(-F
j[cm6&F`xFbe[l.Fj\[m.8C>F^]y-F_\l6$/-F^_mFbalFdxFfo@$/-F^_m6#F^]yF^r>F^]yFb`m>F
inFdx>Fb\fmFc\fm?&Ff]nFh[fmFh]lC.>F_`y-F\]m6$Fjix<%F_hem/F_`blF^]y/FfdemFh]fm@%
-Fg_p6$Fdhq.%(regularG>F_`y-F\]m6$F_`y<#/FhaemF\\]l>F_`y-F\]m6$F_`y<#/Fhaem<#-F
dq6$F]]l/F]]l;F[]lFiael>FgcnFi\fm@%2F^rFbe[lC$-Ffal6&F[glF]dq%7smart~ramificati
on~c~=G)FgcnFj\[m>Fer-Fj]cm6%/F2FaefmFc\lFS>Fer-F\gemFacm>Fdiq-Fgx6$Fix-Ff^l6#*
(Fbe[lF^r)FgcnFabfmF^r)FSFaiuF\s>Fer-Figbl6%FerFSFdiq@$F]]fm>FdbnF`]fm-Ffal6&F[
glF]dq%>completing~exponential~part:~GFdiq>Fbo-Fb]]l6%FerFSF_`y>Fin6$Fin-F_am6#
Fefn>Fdbn-F\]m6$Fdbn<#-Fdq6$-F_am6#-F\^]lF[`]l/F^i[lFefn>FifemFf^fm>F[rFd^am@$5
5F^[fm1Ff_yFifem/Fe[[mFe\mFggm>Faq-F_\l6$/F`bmFdxFaq@$332FifemFf_y4Feeem4-Fg_p6
$FdhqFij\lC$>F_s7#F_s?(F]]lF\pF^rFhh]lFh]lC%>Faj[l,&-F^iw6$&&F_sF[^nFaflFjqF^r-
Fc_m6$&&F_s6#Ff\yFaflFjqF\s@$330Faj[lF[]l0Faj[lF^r5-F\ev6$Ff\yFj^[m-F\ev6$,&F]]
lF^rFbcuF^rFj^[mC+>8D&FjifmF\\l>8E&FfifmF\\l-Ffal6&F[glF]dq%@upper~bound~for~ra
mification~isGFaj[l-Ffal6(F[glF]dq%3bounds~for~slope:~GFijfm%&<~k~<GF\[gm>F\^[m
7#-F_am6#-Fcam6#<$-Fdq6$F^i[l/F^i[l;F\p,&Faj[lF^rFbcuF^rFaj[l>Fdb[mFe\m>FinFdx-
Ffal6&F[glF]dq%8possible~ramifications:GF\^[m?(F$F^rF^rF$Fhj[lC2>Fbe[lF[c[m-Ffa
l6&F[glF]dqFi[_l,&Ff_yF^rFifemF\s-Ffal6&F\pF]dqF\\_lFbe[l>Fer-Faao6&Fc\lFSF2Fbe
[l-Ffal6&FgflF]dqFb\_l<#-Fdq6$F^i[l/F^i[l;,&*&Fbe[lF^rFijfmF^rF^rF^rF^r,&*&Fbe[
lF^rF\[gmF^rF^rF\sF^r>Fin-Fb]]l6%FerFS-F\]m6$Fdhq<%/FhaemFc]gm/FfdemF[]gm/Fgeem
Fbe[l>F[r6$F[r-F_am6#&Fin6%F^rF^rF^r>FgcyF`fl>Fdbn-F\]m6$Fdbn<#-Fdq6$Fegfm/F^i[
lF_dl-Ffal6&F[glF]dqFi]_lFgcy>Fifem,&FifemF^rFgcyF^r>Faj[l,&Faj[lF^rFgcyF\s>Fdb
[m-F\]m6$Fdb[m<#Fbe[l>F\^[m7#-F_am6#-Fcam6#-F[`p6$F[\gmFdb[m-Ffal6&F[glF]dqFi^_
lF\^[m@$F^hfmFggm@$F^hfmFggm7$F]belFifemF$F\b[mF$6#FgflFei_lt+FeczF$F$F$@'F]ar-
F\hl6#F_i_l-FG6$FIFjeo-F\hl6#F^fal-F\hl6#,&Faa`lF^rF`a`lF\sF$F$F$6#F\pF^^ut+6%F
h[lF^[lF:6#FfcoF$F$@%Fce\mF^\sC$>Fin-Fgx6&t+6%F___lF^[lFg[lFh^zF$F$C$>Fin-Fabq6
%FIFSF[]l@%1Fd]lF\sF[]l-Fh^m6#*&Fc]clF^r-Fhs6$-%'seriesG6%-Ff^l6#*&FIF^r)FSFc]l
F^rFS,(FVF^rFinF\sF^rF^rFieoF^rF$F$F$6#FgflFIFSFVFcgsF$F$F$6#FgflFh_amt+6%Fh[lF
^[lFcf\lF$F$F$@%F_cel@%F_dgl-Fgx6&t+6%Ff]mF^[lFcf\lF$F$F$-Ff^l6#*&Fa[pF^r-Fheq6
%FjcdmFS<#FVF\sF$F$F$6#FgflF]eglFSFV-FgrFceglF^e^mF$F$F$6#Fgfl%-mat_lead_monGt+
F[\sFh^zF$F$C$>FinFgfn-F_]l6&FIFSFinFinF$F$F$6#F\pFj\xt+6)F][lF^[lFicmFjcmFia]m
%$phiGFa`q6&Fh[lFfcoF,F7F$F$C&@$Fa_o-F\hl6#-&F"6#%8mat_transform_algebraicGF^`o
>FinFdhl>Fer-F\]m6$-F\]m6$-F\^]lF_dr-Fe`p6$-Fhs6$FVF[ioFcg\l-Fe`p6$-FhsF\\amFcg
\l@%0F_dlFghlC4>Fin-F^q6#;F^rFfjm>F_dlFghl>F`flFbfl>FeflFi`jl-Ffal6%FcblFgir%4t
ransforming~systemG>Fc\l-F_\u6$FY-Fgr6#,&-F_\u6$FeflFVF^rFafvF\sFd\u>Fefl-Fgx6&
FheqFc\lFSFer>Fifl-FgrFdo>F]gl-Fgr6#F\dn>F`glFcy>Fcgl,&FjfnF^rFcyF\s>FaguFer>Fe
`oFI@$FhcrC$>FjeuFjy>FbeuFhw>FboFegl>FbajlFinFboC/>F`flFbfl>F`gl,&F`glF^rFcyF^r
@$2F[]lFcy-Ffal6&FcblFgir%6updating~span,~is~nowGF`gl>FaguFer>Fcgl,&FcglF^rFcyF
\s-Ffal6%FcblFgir%?update~transformation~matricesG>Fifl-Fgx6&Fheq-F_\u6$FiflFVF
SFer>F]gl-Fgx6&Fheq-F_\u6$FYF]glFSFerFbggm>Fc\lFfggmFd\u>FeflF^hgmFIF$Fh[mF$6#F
\qFf_bmt+6%F*FicmFi^al6.F7F+F[_alFadmF2FcdmF9FgitF:Fg[lF3FjfvF$F$C->FjqFcfp>Fer
Fjn>FinFag\l>FboF_em>F]]l-F]]nFf_m>Fc\lFb`m>FfoF^r@%/FinFe\m?(FaqFerF\sF^r1FfoF
jqC)>F^pFfo?(F$F^rF^rF$31F^pFjq/&FboF[balF[]l>F^pFjc_l@$/F^pF]`mFc`al>Fc\l7$Faq
-F_amFacm@$2FfoF^pC$>Fbo-F]hn6%FboF^pFfo>F]]l-F]hn6%F]]lF^pFfo?(F`xFjc_lF^rFjqF
h]lC%>F[r*&&Fbo6$FaqF`xF^r&FboF\`lF\s?(F_sF^rF^rFaqFh]l>&Fbo6$F_sF`x-Ff^l6#,&F[
^hmF^r*&F[rF^r&FboFa^lF^rF\s?(F_sF^rF^rFjqFh]l>&F]]lF\^hm-Ff^l6#,&Fd^hmF^r*&F[r
F^rFh^lF^rF\s>FfoF]d\l?(FaqFerF\sF^rF[\hmC)>F^pFfo?(F$F^rF^rF$F_\hm>F^pFjc_l@$F
e\hmFc`al>Fc\lFg\hm@$Fj\hmC$>FboF]]hm>F]]lF`]hm?(F`xFjc_lF^rFjqFh]lC%>F[r-Fix6#
-Fegq6#Fe]hm?(F_sF^rF^rFaqFh]l>F[^hm-Fix6#-FegqF^^hm?(F_sF^rF^rFjqFh]l>Fd^hm-Fi
x6#-FegqFf^hm>FfoF]d\l@$F]\o>FVFc\l@$1F\pFbp-F[w6$FSF]]lFf`alF$F$F$6#FgflFddot+
6&Fh[lF^[lFcf\lF\[o6+Fh]mF2Fd[lF6F5F[inF9F:F___lF$F$C+@$F]\oFhhal>FjqFjn>FerFcf
p>FinFham>Fc\lFdx?(FaqF^rF^rFjqFh]l?(F_sF^rF^rFerFh]lC$>F[r-Fabq6%&FIF`^oFSFV@&
2F[rFinC$>FinF[r>Fc\l7$FaqF_s/F[rFin>Fc\l6$Fc\lFcbhm>FY-Fho6%FjqFerF[]l?&FboF_h
blFh]lC$>Ffo-Ff^l6#*&),&FSF^rFVF\sFc]lF^rF`bdmF^r>Ffbdm-F[y6$F[`sFfoFinF$F$F$6#
F[glF\^cmt+6%F?F][lF^[l6$FfcoFdgwF$F$C->FinF[dl>F_dlFiil>F`fl-FddlF\^q>Fefl-Fgx
6&Fgao-F[y6$FI-Fg_dlFjdlFVFix>FiflFcel>F]glFI>F`gl-F\]m6$-F\^]lF\^qFgddl>FcglFS
>FboFegl>FbajlFinFboF$Fh[mF$6#FgflF`^[lt+Ff_s6(FbeyF9F:Fd[l%#ttGF_gzF$F$C(>Fin-
Fgao6$Fg^m-Fh`[lF\^q>Fer7#-Fe^m6%FinFeehmFjq>Fjq7#Fjq>F[rF[]l?(FboF^rF^rFihxFh]
lC$>Fin-Fgao6$F_b\lFS>F[r,&F[rF^r*&-Fhs6$7#-Fdq6$*&-Fc^l6%FinFSFc\lF^rFjfglF^r/
Fc\l;F\gyFg]lFdtF^r&FjqFf^pF^rF^rF[rF$F$F$6#Fgfl%,denomvectorGt+6$FcdmF^[l6$F[i
nF9F$F$C%>FinF^r?(FboF^rF^r-FfhalF`oFh]l>Fin-F]b[l6$Fin-F`b[l6#&FIFf^pFinF$F$F$
6#F\pFedqt+6&F*FcdmF2Fiit6'F_gzFhjamF,%$dimGF9F$F$C)>Fin-&F\o6#.Ff_z6&FIFSFfb`l
Fc\l@$/FinFdx-F\hl6#Fdx>Fin-Fhs6$FinFban>Fer-F^_m6#-F`dy6%Fg_p-Fe`pF_drFc\l>Fin
-F[y6$-Fdq6$/Fg\q&FV6#,&FablF^rFjqF^r/FjqF_idlFin>Fabl,&FablF^rFerF^r7$Fcgs7$-F
_amF`]n-Fdq6$Fj]al/Fjq;,(FablF^rFerF\sF^rF^rFablF$F`blF$6#F[glFfglt+F$FiagmF$F$
C&@$-FG6$%)Mat_KeysGFfcl>Fb[imFe\m>Fin(F*,&-F^_m6#Fb[imF^rF^rF^r>Fb[im-F\]m6$Fb
[im<#FinFinF$Fh[imF$6#F[]lF\`bmt+6$F*Fi^al6*F+Ficm%%ind1GF,F6F3F:F9F$F$C)@%F`in
>Fin-Fd_bm6%FIFfb`lFe`nC$>FinF^\s>Fc\lFS>FjqFjn>Fer-Fho6%Fjq,&FjqF^rFdfcmF\sF[]
l>F[rF^r>FaqF^r?(F_sF^rF^rFjqFh]l@%30Fc\lFb`m/F_sF]^p>Fc\l-F_\l6$Fa\lFc\lC$>&Fe
rFbwF^r>Faq,&FaqF^rF^rF^rF^blF$F$F$6#F\pF_hqt+F][hm62F7F5F6F]hoF+FjfvF9Fa\imF:F
[_alFadmF2FcdmFgitFg[lF3F$F$C&>FboFjn>Fc\lFcfp>FinFag\l@%Fi[hmC(>Fer-&F\o6#.Fea
_l6#-Ff`hl6$FI-F]]nFf^p>Fjq-F\z6%FerF`uF`gm>F[r-F\z6%FerF`u;Fd_r,&Fc\lF^rFboF^r
@$Fj`hm-F[w6$FSF[r@$F]\oC&>FaqF^r>F_sFdx?(FfoF^rF^rFc\lFh]lC$@$0&FjqF\`lF[]lC$>
F_sFa^l>FaqF]^im@$2FboFaqFggm>FVF^ifmFdj^lC,>F]]lFcfp>F`xFjn>FjqF_em>F[rFfdbm>F
^pFb`m>FaqF^r?(FfoF^rF^rF]]l1FaqF`xC)>FdiqFaq?(F$F^rF^rF$31FdiqF`x/-Fix6#&Fjq6$
FdiqFfoF[]l>Fdiq,&FdiqF^rF^rF^r@$/FdiqFaiuFc`al>F^p7$FfoF[]el@$2FaqFdiqC$>Fjq-F
chn6%FjqFdiqFaq>F[r-Fchn6%F[rFdiqFaq?(FgcnF^bimF^rF`xFh]lC&>Ff]n-Fix6#-Fegq6#*&
&Fjq6$FgcnFfoF^rFc`imF\s>Fdcim-Fix6#-Fegq6#,&FdcimF^r*&Ff]nF^rFc`imF^rF\s?(Fdbn
F]d\lF^rF]]lFh]l>&Fjq6$FgcnFdbn,&F_dimF^r*&Ff]nF^r&Fjq6$FaqFdbnF^rF\s?(FdbnF^rF
^rF`xFh]l>&F[rF`dim-Fix6#-Fegq6#,&FgdimF^r*&Ff]nF^r&F[rFddimF^rF\s>FaqF]^im@$F]
\o>FVF^p@$Fj`hmFi_imFdj^lF$F$F$6#FgflF$
