MVR5
I'desolv=6"%&falseGE\[l]t%)diff/DerGRF$6(%$funG%&ldvarG%$varG%'lindexG%"pG%#iiG
F$F$C'>8$-%$DerG6#&9"6#;"""""$>8%7#-%$seqG6$&F76#8)/FD;""%9#>8'-%#opG6$""#F2?&8
&F=%%trueG@%-%'memberG6%FP-FL6$F;F28(>FJ-%'subsopG6$/FX,&&FJ6#FXF:F:F:FJ-%'RETU
RNG6#""!-Fen6$/FNFJF2F$F$F$%-desolv/splitGR6$%$eqnG%$lfnG6%%'resultG%%teqnG%$fa
cGF$F$C'-%1desolv/storecaseG6#-%&denomG6#9$>F=*&-%'factorG6#-%&numerGFapF:-%)ic
ontentGFap!""@'-%%typeG6$F=%"*GC&>F2%%NULLG?&FPF=FQ@'4-%$hasG6$FPF4-F]p6#FP-F_q
6$FP%"^GC$-F]p6#-FL6$FNFP>F26$F2FP>F2Ffr@$/F2Fdq-%&ERRORG6$%.contradictionG/Fbp
F^o>F27#F2-F_q6$F=F_r@%-Fiq6$-FL6$F:F=F4C$-F]p6#-FL6$FNF=>F27#FfsFjr>F27#F=@$2F
:-%%nopsG6#F2@$3-%)assignedG6#&%*infolevelG6#%(pdesolvG2F^oF[uC$-%&printG6#-%$c
atG6%%6Braching~occurs~with~GFct%+~factor(s)G-Fbu6#-%2desolv/convertD2mGFetF2F$
F$F$%2desolv/intlintestGR6&Feo%%funcG%$ldfG%%xvarG6#%$posGF$F$@+4-Fiq6$Fbp-FL6$
F^o9%%&nofctG4-Fiq6$F[w9'%'noxdepG0-%'degreeG6$Fbp9&F:%'nonlinG2F^o-%(convertG6
$-FL6$FNFfw%"+GC$-FT6%F`w-FL6$F;FfwF2@%2F^o-FL6$F2F\xF%%'noxderG%&noderGF$F$F$%
/desolv/sortHDTGR6$%%leqnGFfo6)%'numeqnG%%tempGF/%#jjG%'lfnameG%&nleqnG%)lastsw
apGF$F$C(>F2-FdtFap@$2F2FN-F\oFap>FX-%$mapG6$R6#%"fGF$6$%)operatorG%&arrowGF$@%
-F_q6$Fbp%)functionG-FL6$F^oFbpFbpF$F$F$F[w?(FPF:F:F2FQ>&FDF\r&FbpF\r?(FPF2F\qF
NFQC%>8*F:?(FJF:F:,&FPF:F\qF:FQ@$2F^o-%/desolv/orderDTG6%-FL6#&&FD6#FJ6#FN-FL6#
&&FD6#,&FJF:F:F:F_\lFXC&>F=Fc\l>Fc\lF]\l>F]\lF=>Fb[lFJ@$2Fb[lFd[l>FPFb[l-Fjw6$F
D%%listGF$F$F$%3desolv/rank_factorGR6%%%lfacGFfo%%lvarG6'%(_moduleG%&tlfacG%"iG
FeoFhoF$F$C(@$2FgyFNFjy>F27'%(onetermG%'gensepG%(solveDEG%-integratepdeG%)odeso
lveG>FXFdq>F=Fbp?(FPF:F:""&FQ?&FJF=FQC%@$/-&F2F\r6%FJF[wFfw%%FAILG%%nextG>FX6$F
XFJ>F=-%%subsG6$/FJFdqF=7$FX-FL6#F=F$F$F$%/desolv/isderivGR6$F`v%%dvarG6#%$derG
F$F$@)4FfzF%-Fiq6$Fiz%"DG-FT6$-FL6#FizF[w-Fiq6$Fiz%%diffGC%>F2Fbp?(F$F:F:F$-Fiq
6$F2F\al>F2-FL6$F:F2-FT6$-FL6$F^oF2F[wF%F$F$F$%.desolv/lderivGR6%FeoF`vFfo6&%&f
nameG%%LderGF``l%%lderGF$F$C)@$43-F_q6$F[wFhzFgvF[o>F2Fiv>FJ-%'selectG6%Fiq-%'i
ndetsGFgzF2>F=-FL6$F:FJ>FJ-Fen6$/F:FdqFJ?&FPFJFQ@$2F^o-%.desolv/comDerG6%FPF=Ff
w>F=FPF=F$F$F$%0desolv/sortDer2GR6$F_blFfo6*%$resG%%backGF/FayF`yFdyFby%'numder
GF$F$C)>Fb[l-F]z6$RF`zF$FbzF$@%FfzFizFbpF$F$F$F[w@%2FNFH>F=&F76#F;>F=F:>8+Fgy@$
2F_elFNFjy?(FPF:F:F_elFQ>F]_lF^[l?(FPF_elF\qFNFQC%>FDF:?(FJF:F:Fd[lFQ@$2*&-Fh[l
6'&F2F^\lF:&F2Fd\lF:Fb[lF:F=F:F^oC&>FXF]fl>F]flF^fl>F^flFX>FDFJ@$2FDFd[l>FPFD-F
jw6$F2F`]lF$F$F$%.desolv/gcdDerGR6$%'Hterm1G%'Hterm2G6+%&acof1G%&acof2G%$ex1G%$
ex2G%&sfac1G%&sfac2GFjo%%tfacG%&scindGF$F$C0>FP-Ffp6#-%)simplifyGFap>FJ-Ffp6#-F
\hl6#F[w@%-F_q6$FPFaq>FX<#-FLF\r>FX<#FP@%-F_q6$FJFaq>FD<#-FLF^\l>FD<#FJ>FX-F]z6
$R6#%"aGF$FbzF$@%-F_q6$FbpF_r7$-FL6$F:Fbp-FL6$FNFbp7$FbpF:F$F$F$FX>FD-F]z6$RFfi
lF$FbzF$@%FiilF[jlF`jlF$F$F$FD>8,-%*intersectG6$-F]z6$RFfilF$FbzF$&Fbp6#F:F$F$F
$FX-F]z6$RFfilF$FbzF$F^[mF$F$F$FD>F2Fdq>F=Fdq?&Fb[lFXFQ@%-FT6$&Fb[lF_[mFgjlC&>F
_el-Fibl6%R6$Fgil%"bGF$FbzF$@%/F^[mF[wFQF%F$F$F$FDFi[m>F_el-FL6$F:F_el@%33-F_q6
$&Fb[lF_\l%(numericG-F_q6$&F_elF_\lF\]m2F_]mF[]m>F=6$F=)Fi[m,&F[]mF:F_]mF\q>F26
$F2)Fi[m,&F_]mF:F[]mF\q>FD-%&minusG6$FD<#F_el>F=6$F=)Fi[mF[]m>FD-F]z6$RFfilF$Fb
zF$)F^[m&FbpF_\lF$F$F$FD>F2*&-Fjw6$<#F2FaqF:-Fjw6$FDFaqF:>F=-Fjw6$<#F=Faq7$F2F=
F$F$F$%-desolv/separGR6%FeoFfoFe]l6'%$eqlG%'eqlistG%&DIfacG%'tDIfacG%'mulfacGF$
F$@%43-F_q6$Fbp%)fractionG2F^o-Fdt6#-%/desolv/fmemberG6$-Ff`m6$FbpF[wFfwC&>F2-%
-desolv/eqsepG6%7$F:FbpF[wFfw>F=Fdq?(F$F:F:F$0F27"C&>FP&F2F_[m>F2-Fen6$FcclF2?&
FJF2FQC&>FX*&&FJF_\lF:&FPF_\lF\q@$2F^o-Fdt6#-Ff`m6$FX7$-FL6#Ffw-FLFahlF`_l>F2-F
e_l6$Fg_lF2>FP7$,&&FPF_[mF:*&FXF:&FJF_[mF:F:F_bm>F=6$F=FPF`t7#F_amF$F$F$%2desol
v/checkICminGR6#%)ICvertexG6%F]bl%$pICG%%flagGF$F$C'>F2F\jl@$4-Fit6#&%'_ICminGF
et>FcdmFcam?&F=FcdmFQC$>FP-%0desolv/isderivDG6$F^jlF[t@&5/FPFQ/FP.%&equalG-F\o6
#F%/FP.%'negderG>Fcdm-Fe_l6$/F=FdqFcdm>Fcdm7$-FLFbdmFbpFQF$6#FddmF$%.desolv/new
arbGR6#Fe]l6#%(arbnameGF$F$C&@$4-Fit6#%1desolv_arbitraryG>FifmF^o>Fifm,&FifmF:F
:F:>F2-Feu6$%$AAAGFifm@%/FgyF^oF2-F26#-FLFapF$FhfmF$F\vR6#Feo6&%%DeqnGF_bl%&lar
bfGFhoF$F$C&>F=F[cl>FP-Fibl6%FiqF=%%DerAG>F=-Fibl6%FiqF=F4@%/-FdtFj_lF^oFbpC&>F
=-F]z6$RF_`lF$F$F$/Fbp-%/desolv/Der2mapGFapF$F$F$F=>FJ-Fe_l6$F=Fbp@$2F^o-FdtF\r
C$>FP-F]z6$RF_`lF$F$F$/Fbp-%0desolv/Der2mapAGFapF$F$F$FP>FJ-Fe_l6$FPFJFJF$F$F$%
(con2mapGR6$%$vecGFe]l6#F,F$F$C%@$F^sF[o>F2-F]z6%FjwF[w%'symbolG-Fjw6$-F]z6$R6#
%#zzGF$FbzF$*&F^[mF:-F\al6$-%$_PPG6#-FL6#T#-FL6#Ff^mF:F$F$6$F,F2-Fe_l6$-F]z6$R6
#%"xGF$FbzF$/Fbp-Fjw6$FbpFjjmF$F$F$F[wFbpF^xF$F$F$%/directseparateGRFe_m6(%$Eqn
GFho%%lofnG%&lovarG%&lsvarG%&lDfacGF$F$C'>F2-%'expandGFgp>FP-Ff`m6$F2F[w>FJ-Ff`
m6$F2Ffw>FX-%/desolv/stardepG6%F2FPFJ@%5/FXF__l0&FXF_\lF^oF__lC(>F=-Fc_mF]^n>FD
-Fjw6$-F]z6$R6#%#exGF$FbzF$F^[mF$F$F$F=%$setG>F=-F]z6$RF\_nF$FbzF$Ff^mF$F$F$F=@
$3-Fit6#&F\u6#%+separationG2F^oFg_nC)-FbuF$-Fbu6#%JDirect~separation~applied~to
~the~equationG-Fbu6#-F\vFap-Fbu6#%Zwith~assumption~that~expressions~are~linearl
y~independentG-FbuFC-Fbu6#%&givesG-Fbu6#-F\vFj_l@%-Fit6#%0desolv_liassumeG>F`an
<$FD-FLF_an>F`an<#FD7%F=FcamF:F$F_anF$Fc^lRFe_m6#FioF$F$C%>F2-%2desolv/convertm
2DGFi`m>F2-%*odesolveDG6%F2F[wFfw@%/F2F__lF__l-F\v6#FfamF$F$F$%.desolv/mindifGR
6%%&derivGFeo%'maxderG6$%%lordGF+F$F$C%>F2-%$zipG6%R6$Fd\n%"yGF$FbzF$,&FbpF:F[w
F\qF$F$F$FfwF^jl>F=-F_cn6%R6$%%varxG%%ordxGF$F$F$-%"$GFi`mF$F$F$-FL6$F;FbpF2@%F
ehmF[w-F\al6$F[wFi_lF$F$F$%0desolv/subsfunDGR6$Feo%'subeqnG6'F+FhoF``lF_bl%%sfu
nGF$F$C(>FJF[cl>FJ-Fibl6%FiqFJ-FL6$F:-%$lhsGFahl>FXF`en?&FPFJFQC%>F2-F_cn6%RFbc
nF$FbzF$FdcnF$F$F$Fcr-FL6$FNFX>F2-F_cn6%F]dn-FL6$F;FXF2@%/FctF^o>FJ-Fe_l6$/FPF[
wFJ>FJ-Fe_l6$/FP-F\al6$F[w-FLFetFJ>F=-F\hl6#-Fc]n6#-Fe_l6$FJFbp@%3-Fit6#%0desol
v_simplifyGFhgn-%.desolv/mysimpGFj_lF=F$FggnF$%.desolv/xderivGR6%%%der1G%%der2G
Fbv6#%%dordGF$F$C$>F2,&-Fjw6$F^jlF^xF:-Fjw6$-FL6$FNF[wF^xF\q@%32F^oF2/Fbp-F\al6
$F[w-F]dn6$FfwF2FQF%F$F$F$%/desolv/sortDerGRF]dlF^dlF$F$C)@%Fidl>F=F[el>F=F:>Fb
[l-F]z6$RF`zF$FbzF$@%FfzFizFbpF$F$F$F[w>F_elFgy@$FaelFjy?(FPF:F:F_elFQ>F]_lF^[l
?(FPF_elF\qFNFQC%>FDF:?(FJF:F:Fd[lFQ@$2*&-Fh[l6%-FL6#F]fl-FL6#F^flFb[lF:F=F:F^o
C&>FXF]fl>F]flF^fl>F^flFX>FDFJ@$Fefl>FPFDFgflF$F$F$%-desolv/rmfacGRFdo6'Fio%%te
rmGFd]lFjoFhoF$F$C%F\p>F2Fdp@'-F_q6$F2FaqC'>FPFdq?&FJF2FQC%@$/-Fdt6#-Ff`m6$FJF[
wF^oF`_l@$-F_q6$FJF_r>FJF^cl>FPF_jm>FX-Fjw6$7#FPFaq-F]p6#*&F2F:FXF\q@%/FXF:FjrF
X-F_q6$F2F_r@%2F^o-Fdt6#-Ff`m6$FcalF[wC$-F]p6#FKFcal-F[s6$%Icontradiction:~Deso
lv~can't~handle~'eqn'GF^sF2F$F$F$%/desolv/der2strGR6#Fhbn6'F]blFe]lF+F-F/F$F$@)
Fb`lFbpFc`lC'>F=7#Fegm>F2Fh`l>FJ7#-FL6#-FL6$F^oFiz>FP7#-F@6$&F=6#Fin/FX;F:-FdtF
^\l-Feu6%F2%"_GFghlFj`lC'>F2Fbp>FPFdq?(F$F:F:F$F`alC$>FP6$FPFK>F2Fcal>F2Fgal-Fe
u6%F2F_`oFPFizF$F$F$%.desolv/SEPAR2GR6$Feo%%svarG6'%&lsfacGFh]lFho%%dfacG%%ifac
GF$F$@%-F_q6$FbpF^xC&>F2-F]z6%%1desolv/sepfactorG-Fjw6$FbpF`]lF[w>FPFdq?(F$F:F:
F$FbamC&>FJFcal>F=-Fibl6%RF_\mF$F$F$@%/Ff^mF[wFQF%F$F$F$F2F^bm>F2-Fibl6%RF_\mF$
F$F$@%FfboF%FQF$F$F$F2F^bm@%/FfhmF:>FPF_jmC$>FX-Fjw6$-F]z6$RFfilF$F$F$F^[mF$F$F
$F=F^x>FP6$FP7$FXF^bmFa]o7#-FjaoFi`mF$F$F$%0desolv/normformGR6%%&oleqnGFfoFe]l6
.F/Fay%#kkGFeoF^bl%&nLderGFbyF]y%&eqnumG%%lHDTG%%narbGF_yF$F$C*@$3-Fit6#&F\u6#%
)decoupleG2F^oFjdo-Fbu6#%EReduce~the~system~into~a~normal~formG>Fb[l-F]z6$RF`zF
$FbzF$@%FfzFizFbpF$F$F$F[w>F_el-Fjx6$FbpFb[l>F2F:>8/-Fdt6#F_el?(F$F:F:F$2F2F\fo
C0>F2,&F2F:F:F:@$2F\foF[fo@$FgdoC$>F[foF\fo-Fbu6#-Feu6$%/#~equations~:~GF[fo>Fg
jl&&F_elFetF_[m>FX&F_goF_\l>FJ&F_goF\el>FJ-%.desolv/reduceG6%FJ7#-F@6$&F_elF\r/
FP;F:,&F2F:F\qF:F[w@$/FJF^oC%>F_el-Fen6$FirF_el>F2F^hoF`_l>FJ-Ff[oFh\o>FD-%1des
olv/diffparseGFh\o@$2%0desolv_sublimitG-%'lengthG6#&FDF\elC'-Fbu6$%3before~repl
acementG-F_ioF^\l>8.-F_fm6#-Ff`m6$FaioFfw>%+desolv_subG<$-FL6#F^jo/FhioFaio>FJ-
Fc]n6#,&FJF:*&,&FhioF:FaioF\qF:)&FDF_[m&FDF_\lF:F:@$3Fhdo2FNFjdo-Fbu6#-Feu6$%UA
rbitrary~replacement~;~The~coefficient's~length~is~GF^io>FJ-FjgnFh\o>FD-Fen6%/F
;Fdq/FGFdqFD>FJ7%FgjlFDFJ@%/FDFX>F_el-Fen6$/F2FJF_elC%?(F=F:F:F^hoFQ@$2F^o-Fh[l
6%-FL6#&&F_elFj_lF_\l-FLFCFb[l%&breakG@%/F=F\fo>F_elF`\pC$>F_elFcho>F_el-Fen6$/
F=6$FJF\]pF_el>F2F=@$3Fhdo2F:Fjdo-Fbu6#-Feu6&F\goF\fo%3~;~total~length~:~G-F_io
F]foF_elF$6$F^joF]ioF$F]amR6%Fg_mFfoFe]l6)%$EqlG%'sepgrpGFhoF^aoF]y%%olfnG%&olv
arGF$F$C&>F27#Fbp>FPFdq?(F$F:F:F$FbamC(>F=Ffam>F2Fham>FD-Ff`m6$&F=F_\lF[w>Fb[l-
Ff`m6$Ff_pFfw>FJ-Fibl6%R6$F,%#lfGF$F$F$4-Fiq6$F[wFbpF$F$F$Fb[lFD@%/F\`oF^o>FP6$
FPF=C%>FX-%.desolv/sumsepG6$Ff_pFJ>FX-F]z6%RF_\mF$FbzF$7$*&F[wF:F^[mF:Ff^mF$F$F
$FX&F=F_[m>F27$-FLFjnF\gnFa]oF$F$F$%2desolv/lineartestGR6%%%expnGFhbnFbv6#%(var
_numGF$F$@%33-Fiq6$F[wFfw/-FdwFi`mF:2F^oFghnC$-FT6%Ffw-FL6$F;F[wF2@%2F^o-FL6$F2
FihnFQF%F%F$F$F$F]pR6#Fiap6%FhoFjoFioF$F$@$4-F_q6$FbpF\]mC%>FP-FfpFap@%FchlC$>F
2Fdq?&F=FPFQ@$4-F_q6$F=F\]m>F26$F2F=>F2FP@%-Fit6#%0desolv_nzassumeG>Fedp-%&unio
nG6$FedpF[_m>FedpF[_mF$FddpF$%)separateGRFe_mF$F$F$-F]z6$RF\_nF$F$F$Ff^mF$F$F$-
Fc_m6%FbpF[wFfwF$F$F$%(findunkGRFdo6'F_bl%&lnderG%&ltermGFi[o%&lnameGF$F$C(>FX-
F]z6$RF`zF$F$F$@%FfzFizFbpF$F$F$F[w>F2Fdq>F=Fdq>FP<#Fbp?(F$F:F:F$0FP<"C%>FJ-FL6
$F:FP>FP-Fen6$FcclFP@*-FTFh\o>F=6$F=FJ-F_q6$FJFhz@$3-Fiq6$-FL6$F^oFJ<$F\alFe`l-
Fiq6$FJFX>F26$F2FJ5-F_q6$FJF^xF[il>FP<$F_ilFghl3Fj\o-F_q6$-FL6$FNFJF\]m>FP<$Fgh
lF^cl7$-Fjw6$F[_mF`]l-Fjw6$Fa_mF`]lF$F$F$%0desolv/stardep2GR6$FfoFe]l6)%%freqG%
%min_GFhoF`v%%ldepGFjcnF/F$F$C+>FXFdq?&FJFbpFQ@$F`gp>FX6$FXF_il>FX7#FX>F27#-F@6
$F^o/Fb[l;F:-FdtFahl?&FDFXFQC$-FT6%FDF[w.Fb[l>F2-Fen6$/Fb[l,&&F26#Fb[lF:F:F:F2>
F=-%$minG6#F\gn>FPFdq?(Fb[lF:F:FcjpFQ@$/F^[qF=>FP6$FP&F[wF_[qFa]oF$F$F$FjdmR6$F
^hnF_hn6$F-%'linnegGF$F$C%>F2-F_cn6%RFbcnF$FbzF$FdcnF$F$F$FbpF[w>F2-Fe_l6$/F^oF
dqF2@%FafnF`emC$>F=-Fibl6%F_qF2%)negativeG@'FehmFQ/FfhmFctFeemF%F$F$F$%*diff/De
rAGRF$F)F$F$C'>F2-F`hmF5>F=F>>FJFK?&FPF=FQ@%FS>FJFZF[oF_oF$F$F$%3desolv/collect
termGR6#Fd]l6&Fho%&lxfunGFjoFh]lF$F$C&>F2Fdq>F=Fdq?&FPFbpFQ@$4-FT6$F`cmF`tC&>F=
6$F=F`cm>FJ-Fibl6%R6$F]_nFazF$F$F$@%Fb\mFQF%F$F$F$FbpF`cm>FJ-Fjw6$-F]z6$RF\_nF$
F$F$Ff^mF$F$F$FJF^x>F26$F27$F`cmFJF`sF$F$F$%3desolv/clearderdepGR6%FeoF]yFfoF$F
$F$-%0desolv/reducecnG6%Fbp-F]z6$RF\_nF$FbzF$F\_oF$F$F$F[wFfwF$F$F$%.desolv/sim
pscGRF^cpF_cpF$F$C%>F2Fdq@$FacpC$>FP-Ffp6#-Fc]n6#*&FbpF:FjpF\q@'Fchl?&F=FPFQ@$F
\dp>F2F`dp3F]r-F_q6$FcrF\]m>F2Fgfp>F2FPF2F$F$F$%/desolv/sortfunGR6#Ffo6'%'num_f
nGF[cnF/%%num_G%&lndepGF$F$C*>F2Fgy@$FiyFjy>FX-F]z6$RF`zF$FbzF$@%FfzFgyF^oF$F$F
$Fbp>FJ)"#5-F_ioFet>F=7#-F@6$,&*&&FXF\rF:FJF:F:FPF:/FP;F:F2>F=-%%sortGFj_l>F=-F
]z6%RF_\mF$F$F$-%%iremGFi`mF$F$F$F=FJ@%32F:FH/&F7F_\lF\q7#-F@6$&Fbp6#&F=F\rFacq
7#-F@6$&Fbp6#&F=6#,&F2F:FPF\q/FP;F^oF^hoF$F$F$%.desolv/newunkGR6$Fe]lFfoF$F$F$@
%FbgmC&@$4-Fit6#%0desolv_constantG>FjeqF^o>Fjeq,&FjeqF:F:F:?(F$F:F:F$-Fiq6$F[w-
Feu6$%#C_GFjeq>FjeqF]fq-F46%FafqFcamFcamC&@$4-Fit6#%0desolv_functionG>F\gqF^o>F
\gq,&F\gqF:F:F:?(F$F:F:F$-Fiq6$F[w-Feu6$%#F_GF\gq>F\gqF_gq-F46%Fcgq7#-F]dn6$F^o
FgyFbpF$6$F\gqFjeqF$%0desolv/SEPARATEGRFe_m6(F]y%)new_eqnsGF/%'occfnsG%(occvars
G%%seqnGF$F$C&>F27#-Fc]nFap>F=Fdq?(F$F:F:F$2F^oFctC$>FPF^o?&FDF2FQC&>FP,&FPF:F:
F:>FJ-Ff`m6$FDF[w>FX-Fibl6&R6%FgilF`\m%"cGF$F$F$3Fa`p4-Fiq6$FfwFbpF$F$F$FfwFDFJ
@%/-FdtFjnF^oC$>F=6$F=FD@$1FctFP>F2-Fen6$/FPFdqF2C$>F27$-FL6#Fgjq-FL6#-FL6$F:-%
.desolv/sepsumG6$FbpFXF^]pFa_mF$F$F$%0desolv/sortDer3GR6#F_bl6)F_dlF`dlF/FayF`y
FdyFadlF$F$C(@%F^dq>F=F`dq>F=F:>Fb[lFgy@$2Fb[lFNFjy?(FPF:F:Fb[lFQ>F]_lF^[l?(FPF
b[lF\qFNFQC%>FDF:?(FJF:F:Fd[lFQ@$2*&,&-Fjw6$-FL6$FNF]flF^xF:-Fjw6$-FL6$FNF^flF^
xF\qF:F=F:F^oC&>FXF]fl>F]flF^fl>F^flFX>FDFJ@$Fefl>FPFDFgflF$F$F$%/desolv/backin
tGR6&Feo%&lvlcoGFfoFe]l65F\dmFio%%nlfnG%$lcoGF,%#coGF\]n%%fnewG%%ieqnG%&fnew1GF
*%#tfG%%fvarGFho%%ilcoG%%var1G%$tcoG%(reverseG%%tlfnGF$F$C*>85R6#%"lG6#Fi]lF$F$
7#-F@6$&Fbp6#,(FgyF:F2F\qF:F:/F2;F:FgyF$F$F$>F2FQ>F=Fbp>FPFdq>86Ffw?&FJF[wF2C%>
FXFbcm>FJF^bm?&FDFJFQC%>Fb[l-Ff`m6$F=F_`r>Fgjl-%,integrateDEG6%F=Fb[lF\jp@%0Fgj
lF__lC)>F_`r&FgjlF_\l>Fgjl&&FgjlF_[mF_[m>F_el-Fibl6%R6$FazF_`pF$F$F$4-FTFi`mF$F
$F$F_`rFb[l>F2FQ>81-Fc]n6#*&FDF:FgjlF:?&FhioF_elFQC&>F[foFhio>80-F[^m6$<#-FL6#F
`w<#-FL6#Fhio?&82-F^_rFahlFQC$>83&F`crF_[m@$-Fiq6$FhioFdcrC$>F`cr&F`crF_\l?&84F
`crFQ@$2F^o-Fdt6#-Ff`m6$F]drFfbr>F[fo-F\al6$*&F[foF:F]drF\qFdcr@$0FhioF[fo>F^br
-Fe_l6$/FhioF[foF^br>FP6$FP-FLF]fo>F2F%>FPFa]o@%F27$F^brFP7$FbpFPF$F$F$%'comtab
GR6$%%lvecGFe]l62F/Fay%)gamcountG%%commGFeoFio%&leqnsG%'lambdaG%$psiG%%solnG%'n
umvecG%%ceqnG%(llambdaG%&gammaG%'lgammaG%*commatrixGF$F$C/>FhioFgy@$3-Fit6#&F\u
6#Ffer2F^oF]grC$-Fbu6#%2The~vectors~are~:G?(F2F:F:FgyFQ-Fbu6#/&FgjlFet-FL6$F2Fb
p>Ffbr7#-F@6$F_go/F2;F:Fhio>F[fo-Fjw6$7#-F@6$*&F_goF:FhgrF:F_hrF^x>Fdcr7#-F@6$7
#-F]dn6$F^oFhioF_hr>Fdcr-Fjw6$Fdcr%'matrixG>F`crFdq>FPF^o>FD-Fjw6$-F_cn6%RFbcnF
$FbzF$*&FbpF:F[wF:F$F$F$Ffbr-F]z6%F`jmFbpF[wF^x?(F2F:F:,&FhioF:F\qF:FQ?(F=FbfoF
:FhioFQC$>FJ-%+liebracketG6%&FbpFet&FbpFj_lF[w@$0FJF^oC'>FX-F\hl6#-Fc]n6#-Fip6#
,&FDF:FJF\q>Fb[l-F]hq6%FXFcamF[w@$/Fb[lF%F`_l>8--%&solveG6$Fb[l-Fjw6$FfbrF^_n@%
/Fi[sFdqC'>FPF`iq>FJ-%1desolv/conDif2opGF^\l>F`cr6$F`cr/&F^brF\rFJ>&FdcrF`dpFi\
s>&Fdcr6$F=F2,$Fi\sF\qC$>F[]s-Fe_l6$Fi[sF[fo>F]]s,$F[]sF\q>F`cr7#F`cr@$FjfrC%-F
bu6#%9The~commutator~table~is:G-Fbu6#Fdcr@$2F^oFPC$-Fbu6#%&whereG-Fbu6#F`cr@%F`
^s7$-%%copyGF^^sF`crFi^sF$F$F$%'genvecGR6%Ffo%%lparGFe]l6&%$parG%$GenG%(Lievecs
G%'symvecGF$F$C'>F=-F_cn6%R6$F]_nF,F$F$F$@%-F_q6$Fbp%"=G*&-%$rhsGFapF:-%%DiffG6
$Ff[nF[wF:*&FbpF:Fa`sF:F$F$F$FbpFfw>F=-Fc]n6#-Fjw6$F=F^x>FPFdq?&F2F[wFQC*@%4-F_
q6$F2Fhz>FJ-Ffp6#-%3desolv/coefficientGF^]s>FJ-Ffp6#-%/desolv/deptermGF^]s@$F`h
oF`_l@$F[il>FJ-Fibl6%FiqFJFfw>FJ-Fc]nF^\l>F=,&F=F:*&FJF:F2F:F\q>FJFd\s>FPF_jm@$
3-Fit6#&F\u6#F[_s2F^oF[csC$-Fbu6$%*parameterGF2-FbuF^\lFa]oF$F$F$%/desolv/eqsor
tDGR6#F]y6'%'sum_ntG%'ot_ordGF/%(num_eqnGFeoF$F$C.>FJFgy@$2FJFNFbp>F2Fdq>F=Fdq>
FPF^o?&FXFbpFQC$>FPF`iq@%-F_q6$FXF^x>F26$F2,&*&F`jqF:)FhbqFfioF:F:FPF:>F=6$F=,&
*&-Fjw6$FienF^xF:F[esF:F:FPF:>F2-Fecq6#F`s>F2-F]z6%RF_\mF$F$F$-F[dq6$Fbp)Fhbq-F
_ioFahlF$F$F$F2FJ>F=-Fecq6#F`t>F=-F]z6%RF_\mF$F$F$FiesF$F$F$F=FJ>F27$Fi_lF\gn7#
-F@6$&Fbp6#F]_l/FP;F:FctF$F$F$%1desolv/splithnihGRFe_m6)F]]nF\]n%'maxvarG%%lmfn
G%(hompartG%*inhompartG%"sGF$F$C+>F2-Ff`mFew>F=Fh`m>FPFct>FJ-F]z6%R6$Faz%"nGF$F
bzF$@%3Ffz/FgyF[wFizFdqF$F$F$F=FP@$Fd`p-F\o6#7$F^oFbp>FXFdq>FDFdq@'Fdao?&Fb[lFb
pFQ@%-Fiq6$Fb[lFJ>FX6$FXFb[l>FD6$FDFb[l-Fiq6$FbpFJ>FXFb[l>FDFb[l7$-Fjw6$F\jpF^x
-Fjw6$7#FDF^xF$F$F$Fd[rR6$%(sep_eqnG%(sep_varG6&%(indlistG%(deplistG%(temlistGF
i[oF$F$@%FdaoC&>F2Fcam>F=Fcam?&FJFbpFQC%>FP-%0desolv/colterm2G6%F2F=-FjaoFh\o>F
2F`cm>F=F_bmFb_mC$>FPFjco7$7#F`cm7#F_bmF$F$F$%3desolv/indsepTypeIGR6&FeoFfoFe]l
F^ao61F\]n%%lnfnG%%nunkG%%xfunGF]]n%%fdepG%'fnodepG%(nunkdepG%'divfacG%,implici
tsumG%,explicitsumG%&exsumGF^jsF`hqFioF$F$C:>F2Fh`m>FXFggs>FJ-Fibl6%FiqF2F`w@$0
F\`oF:-F\o6#F__l>FJF^cl>FD7#F_il>Fb[l-Fibl6%R6$Fd\n%#lxGF$FbzF$FjarF$F$F$FXFD>F
`crFghq>Fi[s-Fibl6%FiqF`crFfgp>Fhio-%%evalG6#,&F`crF:Fi[sF\q>Fi[s-Ffp6#Fi[s>Fgj
lF:@$-F_q6$Fi[sFaq>Fgjl-Fibl6%FiqFi[sFb[l@$0FgjlF:C&@$2F^o-Fdt6#-Ff`m6$FgjlF2F\
]t-F]p6#Fgjl>Fi[s*&Fi[sF:FgjlF\q>Fhio-F\hl6#-Fc]n6#*&FhioF:FgjlF\q@$-Fiq6$Fi[sF
b[lF\]t>Fhio-Fj`p6$FhioF`w>FfbrFdq>F^brFdq>F=Fdq?&F[foFhioFQ@'2F^o-Fdt6#-Ff`m6$
&F[foF_[mFb[l>F^br6$F^br&F[foF_\l-Fiq6$F`atFb[lC(>F_el-Ff`m6$F`atFD>FP-Faeq6$F_
elF[w>F^br6$F^br,&FPF:F`atF\q>Ffbr6$Ffbr*&FPF:F]atF:>FP-F^imF\r>F=Fdcm>Ffbr6$Ff
br*&F]atF:F`atF:>Ffbr-Fc]n6#,&Fi[sF:-Fjw6$7#FfbrF^xF:>F^br-F\hl6#7$FfbrF^br@$Fd
_nC*F\`n-Fbu6#%FIndirect~separation~type~I~applied~toG-Fbu6#-Feu6$%5the~equatio
n~w.r.t.~GF`wF``nFg`n-Fbu6#-F\v6#F^br-Fbu6#%=with~new~unknown~function(s)G-FbuF
_fs7$F^brF`tF$F$F$%.desolv/subsOTGRFedn6'F]blF``lF_bl%%lindG%%lsubGF$F$C)>F2-FL
6$F:F[w>FP-Fibl6%FiqF[clF4>FP-Fibl6%FiqFPF2>FXFdq?&F=FPFQC%>FJ-F_cn6%RFbcnF$Fbz
F$FdcnF$F$F$F[tFihn>FJ-F]z6%F_qFJF]]q@$4-Fiq6$FJFQ>FX6$FX/F=F^o>FXF\jp@%F_jqFbp
-Fc]n6#-Fe_l6$FXFbpF$F$F$%)classifyGRF$6%%%lunkGFgepF/F$F$C(@%/FHF^o>F2Fcam>F27
#-F@6$&F7F\r/FP;F:FH@$4F^an>F`anFdfp@$4Fcdp>FedpFdfp>F=-F]z6$RF\_nF$FbzF$@%FfzF
izFbpF$F$F$F2>Fedp-F]z6$Fb`qFedp@%2F^oFfhmC&-Fbu6#%7non-zero~assumptions~:G?(FP
F:F:-FdtFddpFQ@$-Fiq6$-FL6$FPFedpF=-Fbu6#Fbit-Fbu6#%Clinearly~independent~assum
ptions~:G?(FPF:F:-FdtF_anFQ@$-Fiq6$-FL6$FPF`anF=-Fbu6#F^jtC&Fjht?(FPF:F:F^itFQF
ditFfit?(FPF:F:FjitFQF`jtF$6$F`anFedpF$%'gendefGR6%F]yF^`l%%ivarG6.%'linvarG%(l
dfeqnsGFeoF/%%exprG%%lfndG%(inf_fnsGFe]lF`v%&lzetaG%&sleqnG%'symconGF$F$C7-%2de
solv/checkinputGFaep>F2-%-desolv/gendlG6%FfwF[wF^o>F_elFfbm>Fb[l-F@6$-&%#xiG6#&
FfwF^\l6#F`er/FJ;F:-FdtFhbm>Fb[l7$Fb[l-F@6$-&%$etaG6#&F[wF^\lFc\u/FJFbjp>F=Fdq>
FD-F]z6%R6$FazFa_rF$FbzF$-Fbp6#FibmF$F$F$F[wFfw>Fhio-F]z6$RF\_nF$F$F$-Fc]n6#-Fi
p6#,&-FaenFapF:F_`sF\qF$F$F$Fbp>Fhio-F]z6%%/desolv/LDsolveGFhioFD@$3-Fit6#&F\u6
#Ffjt2F^oF[_u-Fbu6#%=Calculate~defining~equationsG?&FPFhioFQC+@$Fh^u-Fbu6#-Feu6
$%,equation~:~GFP>FX-%(prolongG6%-Fc]n6#-Fip6#,&-FaenF\rF:-F``sF\rF\qFfwF[w>FX-
Fe_l6$/F\gnF`erFX>FX-Fi_q6%FXFhioFD>FX-%.desolv/conderG6%FXF[wFfw>F[fo-Fj`p6#Fe
ap>F[fo-F]z6$RFfilF$FbzF$Ff^mF$F$F$F[fo@$Fh^u-Fbu6#-Feu6$%8#~defining~equations
~:~G-Fdt6#F[fo>F=6$F=-FLF]bu>F=Fa_m>F=-F]z6%Ff[oF=Fb[l@$Fh^u-Fbu6#-Feu6$%A#~def
ining~equations~in~total~:~GFfhm>F=-F]z6&F\bnF=Fb[lFcam>F=-%1desolv/reductionG6
$F=Fb[l>F=-F]z6$F\vF=>F=Fcbu>F=-%.desolv/eqsortG6#-Fjw6$F=F`]l@$Fh^u-Fbu6#-Feu6
$%H#~defining~equations~after~reduction~:~GFfhm7%F=Fb[lF_elF$F$F$%1desolv/depal
lvarGR6%FiapFbvFe]l6$%(set_facG%&livarGF$F$C%@%-F_q6$FbpFaq>F2<#Fegm>F2Fafp>F=-
Ff`m6$-Fibl6%FiqF2F[wFfw@%/F=FfwFQF%F$F$F$%/desolv/stepoutGR6(%%snumG%(eqn_numG
F`frF]y%$ofnGFfo6&%%llenG%'maxlenG%'minlenG%)totallenGF$F$@$3Fht2F:F[uC&-Fbu6#-
Feu6&%(module~GFbp%6~applies~to~equation~GF[w@$2FNF[uC%-Fbu6#-F\v6#&FfwF_[m-Fbu
6#-Feu6$%'gives~G-F\v6#&FfwF_\l@$2F;F[uC(-Fbu6#%9The~current~system~is~:~G-Fbu6
#-F\vF[cr-Fbu6#%<The~current~unknowns~are~:~G-Fbu6#9)-Fbu6#%=The~current~soluti
ons~are~:~G-Fbu6#-F\v6#9(@%2F^o-FdtF[crC&>F2-F]z6$F_ioF`w>FPFa[q>F=-%$maxGFc[q>
FJ-Fjw6$F2F^xC%>FPF^o>F=F^o>FJF^o-Fbu6#-Feu6-F\goFeiu%0,~#~unknowns~:~G-Fdt6#-F
f`m6$F`wFjhu%2,~system~size~:~[GFP%",GF=F_[vFJ%"]GF$F$F$%2desolv/transform1GR6&
FeoFfo%%lstrGF^`l6+FhoF/F_blF``l%%tderG%'lfdvarGF*%'lsubfnG%$lfdGF$F$C->F_el-F_
cn6%RFbcnF$FbzF$/Fbp-Fjw6$F[wFjjmF$F$F$F[wFfw>F2-Fe_l6$F_elFbp>Fgjl-F\clF`as>FD
-Fibl6%R6$Faz%#ldGF$FbzF$-FT6$FizF[wF$F$F$FgjlF`w>Fgjl-F[^m6$FgjlFD>F_el-F]z6$R
F`zF$FbzF$/FbpFizF$F$F$FD>F2-Fe_l6$F_elF2>Fgjl-Fe_l6$F_elFgjl>Fgjl-Fibl6$RF`zF$
FbzF$33Fc`l-%&evalbG6#/FgyF:-F_q6$Fegm%%nameGF$F$F$Fgjl>F_el-F]z6$RF_`l6$%'ordd
erGF`vF$F$C%@%-F_q6$Fb_oFhz>F2-FL6$FNFb_o>F2F:>F=-Fh`lFdgm/Fbp-F\al6$F=-F]dn6$F
i_lF2F$F$F$FgjlFh]vF$F$F$%(coltermGR6#%%lgrpG6(F/%$lDIG%%tgrpG%&grpDIG%$temG%%t
sumGF$F$C(@$FbgmFjy>FP-Fjw6$FbpF^_n>F2F^o>F=-%&tableGF$?(F$F:F:F$FcfpC'>F2Fbfo>
FJFgfp>FPFjfp>FX-Fibl6%R6$%%gdi1G%%gdi2GF$F$F$@%/F^[m&F[wF_[mFQF%F$F$F$FPFJ@%/F
XFdfp>&F=FetFJC%>FD,&-Fjw6$-F]z6$R6#%$gdiGF$F$F$Ff^mF$F$F$FXF^xF:F^bmF:>Fdbv7$F
bcmFD>FP-F[^m6$FPFXFjcuF$F$F$%/desolv/diffvarGRFg^o6'Fe]lF+F-%%ind_GF`vF$F$@'-F
iq6$FbpFe`lC'>F2F\_o>FPF__o>F=Fdq?&FJFPFQ>F=6$F=F]flF=-Fiq6$FbpF\alC&>FXFbp>F=F
dq?(F$F:F:F$-Fiq6$FXF\alC$>F=6$F=Fien>FX-FL6$F:FXF=F:F$F$F$%.desolv/findicGRF]d
o68F/Fay%'neweqnGFibn%'maxindGFicsF\dm%&newicGFby%'flagICG%'eqnumiG%%HDTiG%%eqn
iG%'coeffiG%'eqnumjG%%HDTjG%%eqnjG%'coeffjG%'acofijGF]y%(eqnHeadGFddoF$F$C-@$Fg
do-Fbu6#%ICalculate~integrability~conditions~(ICs)G>Fgjl-F]z6$RF`zF$FbzF$@%FfzF
izFbpF$F$F$F[w>FddmF`av>87Fbp>FD-Fdt6#F]gv>F_elFdq>Fb[lFQ?(F$F:F:F$Fb[lC$>Fb[lF
%?(F2FNF:FD4Fb[lC&>Fhio&&F]gvFetF_[m>F[fo&F[hvF_\l>Ffbr&F[hvF\el?(F=F:F:F^hoFgg
vCC>F`cr&&F]gvFj_lF_[m>Fdcr&FdhvF_\l>F]dr&FdhvF\el@$0-FL6$F:F]at-FL6$F:&FdcrF_[
mF`_l>FX-F_cn6%F]ju-FL6$FNF]at-FL6$FNF_iv>FJ-Fen6$/FNFXF]at>Fi[s-FfcmF^\l@$4Fi[
sF`_l@%-FT6$<$FhioF`cr%)_hasdoneGF`_l>Fcjv<$Fbjv-FL6#Fcjv>Ffbr-Febn6%F]atFfbrFX
>F]dr-Febn6%F_ivF]drFX>F^br-Feas6$FfbrFJ>F^_r-Feas6$F]drFJ>F_`r-Fifl6$F^brF^_r-
F]p6#&F_`rF_[m-F]p6#&F_`rF_\l>FP-Fc]n6#,&*&Fi[wF:FfbrF:F:*&F\\wF:F]drF:F\q@$/FP
F^oF`_l@$%.desolv_reduceGC$@$F][p-Fbu6#-Feu6$%+Reduce~IC~GFJ>FP-Ffgo6%FPF]gvF[w
@$Fd\wF`_l@$F][p-Fbu6$%8New~IC~whose~length~is~G-F_ioF\r>FP-Ff[o6$FPF[w>FD,&FDF
:F:F:>%(_ecountG,&F]^wF:F:F:>88-FjhoFi]w@$2F]io-F_io6#&F`^wF\elC&>89-F_fm6#-Ff`
m6$FPFfw>F^jo<$F`jo/Fi^wFf^w>FP-Fc]n6#,&FPF:*&,&Fi^wF:Ff^wF\qF:)&F`^wF_[m&F`^wF
_\lF:F:@$F][p-Fbu6#-Feu6$Fc[pFd^w>FP-FjgnFi]w>F`^w-Fen6%Fi[pFj[pF`^w>FP7%F]^wF`
^wFP>F_el6$F_elFP>F]gv7$FP-FLF`gv>F]gv-Fjx6$F]gvFgjl>Fb[lFQ>F_el7#F_el@$F[^pC$-
Fbu6#-Feu6&%4~#~new~equations~:~GF\foFa^pFb^p@$F^[p-Fbu6#-Feu6$%+~length~:~G-Fj
w6$-F]z6$F_ioF]gvF`]lF_elF$6(FcjvFddmF]^wFf\wF^joF]ioF$FjgnRF^cp6'F_blFhoF``l%&
stermGFi[oF$F$C)>F2F^et@%F^dq>F2-F[dl6$F2F`dq>F2-Ff[rFet>FJFghq@%F^hp>FJ-Fjw6$F
JF^_n>FJFail>F=Fdq?&FPF2FQC%>FX-Fibl6%FiqFJFP>FJ-F[^mFjgp>F=6$F=-Ffp6#-FjwFfds-
Fjw6$Fa_mF^xF$F$F$%2desolv/convertodeGR6&FeoF`vFhbnFbv6$F_bl%'lowordGF$F$C'>F2F
^et>F2-Fibl6%FiqF2Fiv>F=Fiw>F2-F]z6&R6%F``lF,F[cn6#F__vF$F$C$>F2,&FehnF:FfwF\q@
%/F2F^o/Fbp-%&ODE_YGFahl/Fbp-F\al6$Fjew-F]dn6$F[wF2F$F$F$F2F`wF=-Fe_lFjgrF$F$F$
Fj`pR6$FeoF,6$%%numeG%%denoGF$F$C%>F=-%.desolv/facsepG6$F_pF[w>F2Fb]n@%-F_qF`ju
C%>F2-F]z6%FjfwFgflF[w>F2-F_`vFet-F]z6%R6$%$numG%$denGF$FbzF$7$*&F^[mF:F`bvF\q*
&Ff^mF:&F[wF_\lF\qF$F$F$F2F=C$>F2-FjfwFf]n7#7$*&FfamF:FbapF\q*&&F2F_\lF:Ff_pF\q
F$F$F$%(int/DerGRF$6(F-F.F`v%(DerpartG%(nonDerpGFhoF$F$C%>FPF7-FT6%.%#_XG-FL6$F
;FP.F=@%2F^o-FL6$F=FcrC%>F2Fcr>F2-Fen6$/F=,&&F2Fj_lF:F\qF:F2-Fen6$/FNF2FPF__lF$
F$F$%0desolv/islinearGR6$FeoF`v6$Fi[o%&nopsfGF$F$C$>F=R6$Fi[o%#fnGF$F$F$-Fdt6#-
Fibl6%FiqF[clF[wF$F$F$@'FdaoC$?&F2FbpFQ@$2F:-F=6$F2FivFbemFQ2F:-F=FhvF%FQF$F$F$
%1indirectseparateGRFe_m64F\]nF]]nF]yF^]n%&tlvarGF_^rF[_r%%ldcoGF_]n%$vl1G%#L1G
%%elfnGF[]nFioF,FffwF^aoFhoF$F$C(>FfbrFb]n>F2-Ff`m6$FfbrF[w>F=-Ff`m6$FfbrFfw>FJ
-F\^n6%FfbrF2F=>F_elFdq@)/FJF__lF__l/F^bmF^o-F\o6#-Fh\nF]]x/F^bmF:C$?&F]drFbcmF
QC$>F^_r-Ff[t6&FfbrF2F=F]dr@$0F^_rF__l-F\o6#7%&F^_rF_[m7$Fibm-FL6#&F^_rF_\lFNF\
]tC->FPFcam>FhioF__l>FJFbcm>F^brFfbr?(F$F:F:F$2F^oF\`oC.>F`crFbcm>FJFacl>F[fo-F
ibl6%FiqF2F`cr>F[fo-Ff`m6$F^brF[fo>FhioF__l>FXF__l>Fhio-%-desolv/felimG6&F^brF2
F=F[fo@$/FhioF__lF\]t>F_el6$F_el-FL6#&FhioF_\l>F^br&FhioF_[m>FX-F\^n6%F^br-Ff`m
6$F^brF2F=@$5F`^n/Fb^nF^o>FJFcam@%/-F`pF^dtF:>Fb[lFcamC%>FdcrF`ax>F^br-Fc]n6#-F
ipF^dt>Fb[lFh`x@%/Fb[lFcam>Fi[sF=>Fi[s-Fibl6%R6$Fd\nFe^rF$F$F$F``pF$F$F$F=Fb[l@
%330FXF__lF\ax0F^brF^o>Fhio-F[ep6%F^brF2Fi[s>Fhio7#F^br@$-FitF^^s>Fhio-F]z6%R6$
F]_n%%DenoGF$F$F$*&FbpF:F[wF\qF$F$F$FhioFdcr@$Fd_nC'-Fbu6#%<Indirect~separation
~type~IIG-Fbu6#%3Separated~equationG-FbuFap-Fbu6#%4Resultant~equationsG-FbuF^cr
7%FhioFcamF;F$F_anF$%1prepare_decoupleGRF\y6$FicsF/F$F$C$>F2Fgy7#-F@6$7%F=-Fjho
6%FgjrF[wF:Fgjr/F=FbcqF$F$F$%4desolv/explicitpartGR6%FeoFfoFbvF$F$F$@'Fdao-Fibl
6&R6%F]_nF_`p%#xvGF$FbzF$4-FiqFg`mF$F$F$FbpF[wFfwFgexFbpF^oF$F$F$%0desolv/condo
perGRF_`l6&Ff^rF]bl%#l1G%(ord_derGF$F$@%Fc`lC$>FPFb_o@%-Fiq6$FPFe`lC'>F2Fegm@&F
_as>F2Fgal4-F_q6$F2Fi^vFjy>F=Fh`l@%2F:Fdim>FJFcr>FJF:-F\al6$-F=Fet-F]dnF\hpFbpF
bpF$F$F$%/desolv/lderivxGR6%FeoF`vFbv6(F_blF``lF^bl%)ord_derxGF]fx%(var_posGF$F
$C+@$FbblFbem>F2-Fibl6%FiqF[clFiv>FPFcal>F2Fham-FT6%FfwFbiw.FD>FX-Fjw6$FcrF^x>F
J-FL6$FDFcr?&F=F2FQ@$52FJ-FL6$FDF[t3/F^ixFJ2FJ-Fjw6$F[tF^xC$>FXFcix>FJF^ixFPF$F
$F$%2desolv/transform2GR6%FeoFfoFd[v6-FhoF_blF``l%$memGF/%&func1G%&func2GFh[v%&
fvarsGFf^rF\dmF$F$C*>F_el-F_cn6%RFbcnF$FbzF$/Ff\nF[wF$F$F$FfwF[w>F2Fc\v>F=Ff\v>
F=-Fibl6$RF`zF$FbzF$Fj`lF$F$F$F=>F=-%/desolv/sortderGFj_l>F_el-F_cn6%RFbcnF$Fbz
F$/FbpF[wF$F$F$FfwF[w?&FPF=FQC)>FDFP?(F$F:F:F$-Fiq6$-FL6$F^oFDF\al>FD-FL6$F:FD>
Fgjl7#F]]p>Fb[lFP>FhioF%?&Fi[sFgjlFQ@$-FT6%Fi[sF[w.FJC$>Fb[l-Fe_l6$/Fi[sFb\uFb[
l>FhioFQ@$FhioC%>Fb[l-Fjw6$Fb[lFe`l>Fb[l-Fe_l6$F_elFb[l>F2-Fe_l6$/FPFb[lF2F2F$F
$F$%0desolv/isderiv2GRF\\q6#F-F$F$@%F_bv@%2Ff^mF^hwF%C$>F2-F_cn6%R6$%#a1G%#a2G6
#%"dGF$F$C$>F2Fdcn@%2F2F^oF%F2F$F$F$&FbpF\el&F[wF\el@%-Fiq6$F2F%F%FQF%F$F$F$%(f
rechetGRFhjt6/%%var_GFeoFbfr%&tleqnG%(fmatrixGF/FayFfoF``l%(ope_eleG%%cof_G%(su
bexpnGF+F$F$C)@%-F_q6$F^[mF]`s>FJ-F]z6$RF\_nF$F$F$Fa^uF$F$F$Fbp>FJFbp?&F2F[wFQC
$>F[fo/-F26#Fgbm,&Fi`yF:*&%%eps_GF:-(%"QGF2Fj`yF:F:>FJ-%/desolv/subsfunG6$FJF[f
o>FJ-F\al6$FJF]ay>FJ-F\hl6#-Fe_l6$/F]ayF^oFJ>FX-%&arrayG6$Fj_rFbjp?(FDF:F:F\`oF
QC$>F=&FJFC?(Fb[lF:F:FcjpFQC)>F_el-F\cl6$F=Fhz>F_el-Fibl6%FiqF_el(F`ayFj[q>F_el
-F^[y6$-Fjw6$F_elF`]lF\q>FP-%(collectG6$F=F_el>Fi[sFdq?&FgjlF_el0FPF^oC'>Fhio-%
&coeffG6%FPFgjlF:@$/FhioF^oF`_l>FP,&FPF:*&FhioF:FgjlF:F\q>Ffbr-FdcvFd_t@%/FfbrF
:>Fi[s6$Fi[s7$FhioFcam>Fi[s6$Fi[s7$FhioF\ct>&FXFais7#Fi[s-Fj^sFjnF$F$F$F\[tR6%F
ajsFbjs%(add_facG6#%(ele_numGF$F$@'Fbgm7$7#Fcgu7#F[hu-FT6%F[huF[w.F27$-Fen6$/F2
,&FfjrF:FcguF:FbpF[w7$7$FegmFcgu7$FibmF[huF$F$F$FhcuRFecsFfcsF$F$C.>FJFgy@$F]ds
Fbp>F2Fdq>F=Fdq>FPF^o?&FXFbpFQC$>FPF`iq@%Feds>F2Fhds>F=6$F=,&*&-%.desolv/orderd
GFjnF:F[esF:F:FPF:>F2Fces>F2-F]z6%RF_\mF$F$F$FiesF$F$F$F2FJ>F=F^fs>F=-F]z6%RF_\
mF$F$F$FiesF$F$F$F=FJ>F2FefsFffsF$F$F$%/desolv/isfunofGR6$F`vFgep6#F][xF$F$@'Fj
`lC%>F2F\jl?(F$F:F:F$-Fiq6$FgalF\al>F2FcalFealFc`lFf`lF]]vF$F$F$F^[yRFh[r6&%(nu
m_derGF[cnF/F\bqF$F$C)>F2Fgy@$FiyFjy>FJFgbq>F=7#-F@6$,&*&-Fhgy6#F^[lF:FJF:F:FPF
:Facq>F=Fdcq>F=-F]z6%RF_\mF$F$F$FjcqF$F$F$F=FJ@%F]dqFgdqFadqF$F$F$%1desolv/high
coeffGRFfjw6*F]blF``lF_bl%(highderG%*highpowerG%)maxorderG%*der_orderGFhoF$F$C+
@%Fdbl>F2Fiv>F2F[w>FP-Fibl6%RF^_qF$FbzF$3-Fiq6$FbpF4-FiqFi`mF$F$F$-F\clFapF2@$/
FdimF^oF[o>FDF\q?&F=FPFQC$>Fb[lFcix@$2FDFb[lC$>FDFb[l>FJF=>FX-FdwFcis>F_el-Fibl
6%FiqFghq)FJFX>F_el-Fc]n6#*&F_elF:Fh\zF\qF_elF$F$F$%*genfactorGRFhjt6/FfoF^[uFg
_y%&deropG%'lhsderGFbjs%)deplist2G%'lopeqnGF/%%operGFeoFhoF]blF$F$C4>FP-F]z6$RF
\_nF$FbzF$F]^uF$F$F$Fbp>FJ-%)adjointDG6$-Fc_y6%FPF[wFfwFfw>FJ-Fjw6$FJ%)listlist
G>FDFg[u>F27#-F@6$-&F^frFd_t6#F]]p/FgjlFbjp>F_el-F]z6%RF_\mF$FbzF$-Fjw6$-F_cn6%
%/desolv/appoperGFbpF[wF^xF$F$F$FJF2>Fb[l-F]z6$Fe^oFD>F_el-Fe_l6$/F]]p-FLF_[qF_
el>F=-F]z6%RFf]uF$F$F$Fg]uF$F$F$F[wFfw>F_el-F]z6&Fi_qF_elFbpF=>F_el-F]z6&F]auF_
elF[wFfw>F[fo-F]z6$RFfilF$FbzF$-FL6#-%3desolv/separation2G6%F^[mFcamFf^mF$F$F$F
_el>F2-Fe_l6$F\`zF2>F[fo-F]z6%Ff[oF[foF2>F[fo-F]z6&F\bnF[foF2Fcam>F[fo-F`cu6$F[
foF2>F[fo-F]z6$F\vF[fo7%-FhcuF]buF2Fb[lF$F$F$Fe[uRFhjt6#Ff_yF$F$C&@$4-F_qF\bo-F
[s6#%BThe~first~argument~must~be~a~listG@$4-F_q6$F[wF`]l-F[s6#%CThe~second~argu
ment~must~be~a~listG?&F2F[wFQ@$F_as-F[s6#%YEach~element~of~the~second~argument~
cannot~be~a~functionG@$4-F_q6$FfwF`]l-F[s6#%BThe~third~argument~must~be~a~listG
F$F$F$%&polyDGR6&F,%%degpGFe]lFfo6&F]gtFhoF/%(unknownGF$F$C&>F2Fdq>F=Fdq?(FPF:F
:F[wFQC%>FJ-Faeq6$-Fe_l6$/FbpFdqFfwF`w>F=6$F=*&)FbpFd[lF:FJF:@%/-Fdt6#-FL6$F;FJ
F^o>F26$F2F^cl>F26$F2-F^cl6#-FLFfez7$-Fjw6$F`tF^xF`sF$F$F$Fh[uR6%FijtF^`lF__v6(
FbvFbjs%&lder1G%&lder2G%&tlderGF/F$F$C&>F=Fdq>FP-F]z6%R6$F]blFe]lF$F$F$Fg]uF$F$
F$F[wFbp?(FDF:F:FfwFQC&>FJFdq?&F2FbpFQC%>FX-F]z6%F\alFPF2>FX-Fibl6%R6$FgilFa_rF
$F$F$@%F[brF%FQF$F$F$FX7#FJ>FJ6$FJFeap>F=6$F=Fghl>FPF^hz7%FegmF=FghlF$F$F$%/des
olv/map2DerGR6#F*6*F``lFgdt%&lind2GF+Fe]l%&orderGF/%&D2DerGF$F$C$>F_elRFghz6&%&
lind1GFihz%%fordGF/F$F$@'-F_q6$Fb_o%(indexedGC&>F2F__o>FPFct>F=7#-F@6$,&FPF:-Fd
t6#F[cmF\q/FJFj_r-F46%Fh`lF=F\_oFb_v-F46%Fh`l7#Fe_vF\_o-F46%Fh`l7#F:F\_oF$F$F$@
)Fb`l-F46%FbpFcamFcamFc`l-F_elFapFj`lC&>FJFdq>F2Fbp?(F$F:F:F$F\iyC$>FJ6$FJFK>F2
Fcal@%-Fiq6$FgalFe`l-F\al6$-F_elFetFJC'>FJF^hz>FDF\`o>FX7#F\gn>F=7#-F@6$,&FDF:-
Fdt6#-Fe_l6$/&FXF_[qFdqFJF\q/Fb[lF\gs-F46%FgalF=FX-F46%FizFigqF\_oF$F$F$%(doche
ckGR6&%)org_eqnsG%&lsolnGF]gtFe]l6'%$solGF]y%%zeleGF_dl%%nsolGF$F$C*>F=F[bo?&F2
F[wFQ@$-F_q6$F2F]`s>F=-F]z6%R6$FeoFf][lF`zF$F$C$@%-F_q6$F`enFhz>F2-FL6$F^oF`en>
F2F`en@%-Fiq6$FbpF2-FcayFi`mFbpF$F$F$F=F2?(F$F:F:F$-FT6%F^oF=.FP>F=-Fen6$FijqF=
>F=-F]z6%Ff[oF=Ffw>FJ-F^u6%F=FfwF`w>FXF[w?&F2F^bmFQ@$F]^[l>FX-F]z6%RFc^[lF`zF$F
$C$@%Ff^[l>F2Fi^[l>F2F`en@%F]_[l/Fb^u-Fcay6$F_`sF[wFbpF$F$F$FXF27%FbcmFX&FJF\el
F$F$F$%/desolv/conm2opGR6$FiapFfo6'FhoFioF``lFj_yF+F$F$C'>F2Fdq>F=Fbp?&FPF[w0F=
F^oC$>FJ-F^dy6%F=FPF:@$FijrC%>F=,&F=F:*&FJF:FPF:F\q>FX-FdcvF\r@%Ff]o>F26$F27$FJ
Fcam>F26$F27$FJ7#Fab[l@$Ffa[l>F26$F27$F=FcamF`sF$F$F$%1desolv/gintorderGR6%FeoF
^aoFe]l6(%'deneqnGF_yF`aoF^^q%&lpfacG%&lffacGF$F$C(>F2F_p>F=Fb]n>FJ-Fibl6%FiqFi
byF[w@%/FJFdfp>FP7$F:F=C%>FX-%'coeffsG6%F=FJFchx@$-Fiq6$-F]z6&FcduF[jsF[wFfwFQ-
F\o6#Fcam>FP-F_cn6%R6$%%fac1G%%fac2GF$F$F$7$FbpF[wF$F$F$F[jsF\jp@$0F2F:@%-FiqFf
]n>FP-F]z6%R6$F]_nFjgwF$F$F$7$F^[m*&Ff^mF:F[wF\qF$F$F$FPF2>FP-F]z6%RF[f[lF$F$F$
7$*&F^[mF:F[wF\qFf^mF$F$F$FPF2-F]z6%R6$FjoF,6#%&lopowGF$F$@%-F_q6$Ff^mF^xC$>F2-
%(ldegreeG6$Ff^mF[w@%F]in7$*&F^[mF:)F[wF2F:*&Ff^mF:Feg[lF:FbpFbpF$F$F$FPF[wF$F$
F$%2desolv/initializeGRF$6#F/F$F$C+>F\gqF^o>FjeqF^o@%FcdpC$>%1desolv_pnzassumeG
Fedp>FedpFdfp>FedpFdfp@%-Fit6#%.desolv_moduleG@$2""(Ffh[l>Ffh[lFih[l>Ffh[lFih[l
>%/desolv_modstatG7#-F]dn6$F^oFih[l>F`anFdfp@$4-Fit6#%-desolv_NtermG>Ffi[lF:@$4
-Fit6#%/desolv_efactorG>F\j[l"%+5@$4Ffgn>FhgnF%F$6,F\gqFjeqF`h[lFedpF\j[lFfh[lF
`anF]i[lFfi[lFhgnF$%.desolv/findscGRFggm6&FfoF*Fa^rFioF$F$C$>F2F^et?&F=F2FQC%>F
P-F^dy6$FbpF=@$5/FPF:FijqF`_lF[rF$F$F$%-desolv/RmfacGRFdo6$Fd]lFjoF$F$@%F[euC%>
F2Fdq?&F=FbpFQC&@$/-Fdt6#-Ff`m6$F=F[wF^oF`_l@$Fas>F=Ffs@$3-F_qFjby-FTF]\\lF`_l>
F2F`dp-Fjw6$F`sFaqFbpF$F$F$%-desolv/infinGR6%F`vFijtF^`l6(F/Fho%)zeta_km1GF]blF
e]lF+F$F$C$>FX-FL6#-Fh[u6%F[wFfwF^o@'4-F_q6$FbpFeiz@%F[br-&F`\uFapFjn-&F]]uFapF
jnFf^vC&>FJFiz>FDF\_o>F=7#-F@6$*&-F\al6$-FJFh]u&F[wFetF:-F\al6$-&F`\u6#Ff^\lFjn
FegmF:/F2Fbjp,&-F\al6$-&F]]uF^\lFjnFegmF:Ff`sF\qC'>FJFiz>FDF\_o>FP-Fg\\l6%-Fen6
$FcclFbpF[wFfw>F=7#-F@6$*&-F\al6$Fe^\l-FL6#-Fen6$/F:Ff^\lFDF:-F\al6$Fi^\lF]\yF:
F\_\l,&-F\al6$FPF]\yF:Fh`sF\qF$F$F$F^uR6%F]yF\]nFe]l6,FcyFfo%&nlofnGFho%$sysG%%
lsysG%%lsolG%&tlinqG%%linqG%*starttimeGF$F$C4@$4-F_q6$Fbp<$F^_nF`]l-F[s6#%TINPU
T~ERROR~!!~first~argument~must~be~a~set~or~listG@$4-F_q6$F[wFja\l-F[s6#%UINPUT~
ERROR~!!~second~argument~must~be~a~set~or~listG@$4-F_q6$FfwFja\l-F[s6#%TINPUT~E
RROR~!!~third~argument~must~be~a~set~or~listG@%1FGFH>Fgjl&F76#FG>FgjlFcam>Fi[s-
%%timeGF$-Fgg[lF$>FP-F]z6&F\bnF[wF[wFcam>FP-F]z6$RFghzF$F$F$/FbpFbpF$F$F$FP>F=F
[w>F2-F]z6&F\bnFbpF=Fcam>F2Fgfl>F_el-F]z6&F\bnFgjlF=Fcam>FD7#7&F2F_elFPF=>Fb[lF
dq?(F$F:F:F$2F^o-FdtFCC'@$Fgt-Fbu6#-Feu6$%,#~systems~:GF\e\l>FXFjjo>FD-Fen6$Fcc
lFD>FJ-%0desolv/pdesolveG6$FeapFfw@'F`]xF`_l/F\`oF:>Fb[l6$Fb[lFbcm>FD7$F_ilF]]p
>Fb[l-F\v6#7#Fb[l@$FgtC$-Fbu6$%6Modules'statistics~:~GF]i[l-Fbu6$%6Computationa
l~time~:~G,&Fcc\lF:Fi[sF\q@%/-FdtF_[qF:F]`zFb[lF$F$F$F_^lRFe_mFhanF$F$C%>F2F[bn
>F2-%)onetermDGF`bn@%FbbnF__lFcbnF$F$F$FjfwR6$Fi[oF,6#FaaoF$F$@'F[euC$>F2-Fibl6
%FiqFbpF[w7$F2*&FbpF:F2F\qFf[zF`jlF_amF$F$F$%0desolv/eqsortD2GR6&F]y%%lhilG%&Nl
eqnG%&NlhilG6,%%ssysG%'unssysGF/Fay%&nlhilG%(luninfoG%'lsinfoG%&info1G%&info2G%
%infoGF$F$C.>Fi[sRFggmF$F$F$@%Fdao7%Fgy-F_ioFapF^o7%F:Fgi\lFehnF$F$F$>F2Fbp>F=F
dq?(FPFcjpF\qF:FQ@$/&F[wF\rF^oC$@$0F]_lF^o>F=6$F=F]_l>F2Fgjq>F=F`t>FD-F]z6$Fi[s
F=>Fb[l-F]z6$Fi[sF2>FX-Fe_l6$Fg\qF[w>FJF^o?&F_elFDFQC'>FJFe\l@$5-FT6$&F=F^\lF2-
FT6$,$Fg[]lF\qF2F`_l>FPF^o?&FgjlFb[lFQC$@&2&F_elF_[mFdarC&>Fb[l-Fen6$/F`iq-FL6#
7$F_elFgjlFb[l>F2-Fen6$/F`iq-FL6#7$Fg[]l&F26#F`iqF2>FX-Fen6$/F`iq-FL6#7$F^o&FXF
a]]lFXF^]p/F`\]lFdar@&/F`\]lF:@$2&F_elF\el&FgjlF\elC&>Fb[lFc\]l>F2Fj\]l>FXFc]]l
F^]p2F_]mFaarC&>Fb[lFc\]l>F2Fj\]l>FXFc]]lF^]p>FPF`iq@$/FPFctC%>Fb[l7$F]`zF_el>F
27$F\gnFg[]l>FX7$FeapF^o>FfwF2>F`wFXF$F$F$F]`qRF[`q6.F_blFhoFfdnF``lF[_r%&lLder
GF/%&ltrigG%(indetfnG%+reducediffGFgepFhdtF$F$C,>Fi[sR6%F``lF]yFgep6$FchqF^blF$
F$C$?&F2F[wFQC$>F=-FehzFdbn@&0F\jlFfsF`_l-FT6$Fidm<$FQFaem-F\o6#/F]im-Febn6%F=F
fhwF^jlFdqF$F$F$>FgjlR6$%%expxG%$lneGFiaqF$F$C$@%F[`s>F2-Fibl6%F[`l-F\cl6$F_`sF
hzF[w>F2-Fibl6%F[`lF[clF[w-FL6#-F]z6$FehzF2F$F$F$@$/FcjpF^oFjy>Fhio-F]z6$RF`zF$
FbzF$@%FfzFizFbpF$F$F$Ffw>F=Fghq>FD-F]z6$RF\_nF$FbzF$-Fehz6#F^[mF$F$F$F[w>F2<#-
Fgjl6$F=Fhio?(F$F:F:F$FjhqC%>F[fo-F]z6&Fi[sF2F[wFhio@%/F\buF:>F2<#-Fgjl6$-FL6$F
:F[foFhio>F2-F]z6%FgjlF[foFhio>F=-Fe_l6$F[foF=>F_el7:%$sinG%$cosG%$tanG%$secG%$
cscG%$cotG%%sinhG%%coshG%%tanhG%%sechG%%cschG%%cothG%'arcsinG%'arccosG%'arctanG
%'arcsecG%'arccscG%'arccotG%(arcsinhG%(arccoshG%(arctanhG%(arcsechG%(arccschG%(
arccothG@%-FiqFgcy-F\hl6$-Fc]n6#-FipFj_l%%trigGFff]lF$F$F$Fe\sRF^cp6#%%lsepGF$F
$C%@$F^sF[o>F2-Fd[r6$FbpFf[n-F_cn6%RF_\mF$FbzF$7$Fbp7#FihnF$F$F$FfamFfhwF$F$F$F
jhoRFdo6*FioF_bl%%HderGF``l%'HpowerG%'HcoeffGFbyF\dmF$F$C*@%Fidl>F_elF[el>F_elF
^o>Fb[l-F]z6$RF`zF$FbzF$@%FfzFizFbpF$F$F$F[w>F2Fghq>F=F^et>FPFfs?&FJF=FQ@$2F^o-
Fh[l6'FJF:FPF:Fb[l>FPFJ>FX-FdwFfr@%/F_elF^oC'>FD-Feas6$F2)FPFX@%F^gw>F=-Fjw6$F2
F^_n>F=F[_m>F=-F]z6$RF\_nF$FbzF$@%F[euFegmFbpF$F$F$F=>F=-Fjw6$FbhmF`]l7&FPFXFDF
=7$FPFXF$F$F$%.desolv/approlGR6$%&prolnGFeo6#%%eqn2GF$F$C%@%-F_q6$F[wF]`s>F2,&F
`enF:-F``sFahlF\q>F2F[w-F]z6%R6$Fd]zFiapF$F$F$*&F^[mF:-F\al6$F[wFf^mF:F$F$F$Fbp
F2-Fc]n6#-Fjw6$%"%GF^xF$F$F$F_azRFe_m6%F`hq%)fn_occurGF^]nF$F$C&>F2F%>F=Fh`m>FP
-Fibl6%R6$Fd\nF_`pF$F$F$F``pF$F$F$FfwF=@%3Fcim-Fiq6$FbpFP-FL6$F:-Fd[r6%FghqFPFf
wF]_pF$F$F$Fb^lRFe_mFhanF$F$C%>F2F[bn>F2-%.integratepdeDG6&F2F[wFfwFQ@%FbbnF__l
-FL6#FcbnF$F$F$Fh]^lR6%FiapFfoFe]l63%(lfn_occGFhduFijt%%lintG%(num_intGF/FcyF\d
m%)int_expnG%'varnumG%&expn1GF]y%)int_consGF_^r%(genflagGF]]nFagsF$F$C1>F2Fh`m>
FdcrFggs>F]dr-F]z6%RF^hsF$FbzF$@%FahsFbpFdqF$F$F$F2-FdtF^^s@$/-Fdt6#F]drF^oF\]t
@%2F;FH>F`crF_c\l>F`crF%@%F`cr>F=Ffw>F=-%3desolv/intminderivG6%FbpF2Ffw>FgjlFgh
q>F_elFQ>FJFdq>Fb[lFdq>F=-Fibl6%R6$%"vGF]_nF$F$F$Fa`pF$F$F$F=Fbp?&FPF=FQ@%-FiqF
frC&>FXF^o>FhioFgjl?(F$F:F:F$0FhioF^oC$>Fhio-%.desolv/intpdeG6&FgjlF2FdcrFP@$F`
a^l@'/Fb`xF^oC$>FgjlFd`x>FX,&FXF:F:F:4F`cr>FhioF^oC%>F[fo-%/desolv/partintG6&Fb
`xF2FdcrFP@$0F[foF^oC'>Fgjl,&Fd`xF:&F]atF_[mF:>FXF\b^l>F2F`at>Fb[l6$Fb[l-FL6#-F
en6$FcclF]at>F_elF%>FhioF^o>FJ6$FJFej]lC$>Fhio-Fc]n6#-%$intG6$FgjlFP@$4-Fiq6$Fh
ioF[d^lC$>FgjlFhio>FJ6$FJ7$FPF:>FJF^hz>FX-Fjw6$-F]z6$RF\_nF$FbzF$Ff^mF$F$F$FJF^
x@'/FXF^oF__l-Fiq6$FgjlF[d^lC$@$3-Fit6#&F\u6#%*integrateG2F^oFfe^l-Fbu6$%2canno
t~integrate:GFgjlF__lC&>F^brFdq?&Fi[sFJFQ@$2F^o&Fi[sF_\lC%>Ffbr-F^dz6%-FLFb^tFd
crF[w>F^br6$F^br-FL6#&FfbrF_\l>Fgjl,&FgjlF:&FfbrF_[mF:>F2-Fjw6$<$F\gnFibmF`]l7$
7$FgjlFb[l7$F\gnF^brF$F$F$%'ordervGR6$FiapFijt6'%*deriv_varG%'fn_varGF+F/F`vF$F
$C(>F=F\_o@$0FfhmFcjp-F[s6#%TINPUT~ERROR~!!!~Check~input~functions~and~variable
sG>FPF__o>FX-(%#D_GFh`lFh]u>F2-F@6$&F=6#&FPF^\l/FJ;F:Fdim-F\al6$FXF2F$F$F$%1des
olv/gdefiningGR6&F]y%%liscGF^`lFijt6.F[[uF\[uFeoF/F][uF^[uFe]lF`vFa[uF_[u%'tsle
qnGFb[uF$F$C+>F2-Fh[u6%F`wFfwF^o>Fb[l7$FjbrFgbm>F=Fdq>FD-F]z6%RFf]uF$FbzF$Fg]uF
$F$F$FfwF`w>Fgjl-F]z6$RF\_nF$F$F$F]^uF$F$F$Fbp>Fgjl-F]z6&Fi_qFgjlF[wFD>Fgjl-F]z
6%Ff^uFgjlFD?&FPFgjlFQC'>FX-F[`u6%F]`uF`wFfw>FX-Fe_l6$/F\gnF]`zFX>FX-Fi_q6%FX7$
Fibm-FLFd_tFD>FX-F]au6%FXFfwF`w@%/-Fdt6#Fb^nF^o>F=6$F=&FXF_[mC%>F[foF`au>F[fo-F
]z6$RFfilF$FbzF$Ff^mF$F$F$F[fo>F=F_buF`tF$F$F$F[arRFe_m6$Fc^r%%oeqnGF$F$C$>F2-F
b^l6&FbpF[wFfwF%@%FbbnC$>F=-Fc^lFaep@%/F=F__lF__l7$7#FbapFf_pC$>F=-Fc^l6%&FfamF
_[mFfhwFfw@%Fi]_lF2Fj]_lF$F$F$Fda^lR6&FiapFfoFe]lFijt6-F\dm%&iflagGF[_rFeoF`vF_
blF`^r%)int_lderGFho%*temp_lderG%(ex_partGF$F$C)>F2F%>FJFghq>FgjlFdq>FPFh`m>FP-
Fibl6%FiqF[wF`w?(F$F:F:F$3Fcim4F2C'>FXFgfp>FD-Fdgx6%FJFXF`w>F=-F]v6&FJFXFDF`w?(
F$F:F:F$3Fc__l/F=F%C%>Fb[l-%'lcoeffG6$FJFD@%-%4desolv/intcoefftestG6&-Fial6%Fb[
lFXF[wFDF`wF[wC)>F_el-F[d^l6$FDF`w>Fb[l-F[d^l6$-Fe_l6$/F_el%%DUM_GFb[lFda_l>Fb[
l-Fe_l6$/Fda_lF_elFb[l>Fgjl6$FgjlFb[l>FJ-F\hl6#,&FJF:-F\al6$Fb[lF`wF\q>Fi[s-Fia
l6%FJFXF[w@%/Fi[sFD>F2FQ>FDFi[s>F2FQ>F=Fj__l@)4-Fiq6$FJ-FL6$F^oFX>FP-Ff`m6$FJFP
2F`jqFf\uC$>FPFjfp>F2F%50F=F%F2>F2FQ>F2F%@%3Fc__l0FgjlFdqC&>Fhio-F[d^l6$-F^ex6%
FJF[wF`wF`w>Fgjl6$FgjlFhio>FJ,&FJF:-F\alFb`tF\q7$-Fjw6$7#FgjlF^x-F\hlF^\lF^oF$F
$F$%/desolv/subleqnGR6)F]yF\fuFfh\lFfdnFcy%%nofnGF\i\l6(F/F]blFicsFeo%&tlhilGFg
_yF$F$C->F=-FL6$F:-FaenF[cr>FPFf\u>FXFfw>FDFbp>FX-F]z6%R6$%#eqGF]blF$FbzF$@%Ff[
zF^oF:F$F$F$FbpF=>FX-F_cn6%RF_\mF$FbzF$@%/F[wF^oF^oFbpF$F$F$FfwFX>FD-Fcdn6$FbpF
`w?(F2FPF\qF:FQ@$/&FDFetF^oC$>FD-Fen6$FirFD>FX-Fen6$FirFX>FbiuFD>Fjhu-F]z6&R6%F
eoFchqFazF$F$F$@%-FiqFew/Fb^u-FcdnF[a[lFbpF$F$F$F[wF`wF=>9*FXF$F$F$%(gensepDGRF
e_m65F[]nF\]nF]]nF]yF,FioF]^rF[_rF^]nFa\xF^\xFffwF_]nFb\xF_\xF`\xF_^rF^aoFhoF$F
$C(>F2Fb]n>F=Fe]n>FPFh]n>Fgjl-F\^n6%F2F=FP>F`crFdq@'/FgjlF__lF__l/FaarF^oC(>F_`
r-Fc_mF\i_l>Ffbr-F]z6$RF\_nF$FbzF$F^[mF$F$F$F_`r>F_`r-F]z6$RF\_nF$FbzF$Ff^mF$F$
F$F_`r@$Fd_nC)F\`nF]`nF``nFc`n-Fbu6#FfbrFg`n-Fbu6#-F\v6#F_`r>F`an<$FcanF]\s7%F_
`rFcamF:@%/FaarF:C$?&F^_rFdarFQC$>F_`r-Ff[t6&F2F=FPF^_r@$0F_`rF__l-F\o6#7%Fi[w7
$Fibm-FLF[\wFNF\]tC->FJFcam>Fi[sF__l>FgjlFdar>FDF2?(F$F:F:F$2F^o-FdtFd_tC.>FXFd
ar>Fgjl-Fen6$FcclFgjl>F^br-Fibl6%FiqF=FX>F^br-Ff`m6$FDF^br>Fi[sF__l>FhioF__l>Fi
[s-Fj_x6&FDF=FPF^br@$/Fi[sF__lF\]t>F`cr6$F`cr-FL6#Fbf^l>FD&Fi[sF_[m>Fhio-F\^n6%
FD-Ff`m6$FDF=FP@$5F]`xFha^l>FgjlFcam@%/-F`pFCF:>F_elFcamC%>F[foF`^`l>FD-Fc]n6#-
FipFC>F_elFi]`l@%/F_elFcam>FdcrFP>Fdcr-Fibl6%RF`bxF$F$F$F``pF$F$F$FPF_el@%330Fh
ioF__lFha^l0FDF^o>Fi[s-F[ep6%FDF=Fdcr>Fi[sF[js@$-FitF]bu>Fi[s-F]z6%RFacxF$F$F$F
ccxF$F$F$Fi[sF[fo@$Fd_nC'FfcxFicxF\dxF]dx-FbuFb^t7%Fi[sFcamF;F$F_anF$%0desolv/o
decheckGRFe_m6.%&funcxGF_blF[cn%'lowderG%&hiderGFgdw%&hiordGFdv%*order_ODEGFbvF
f^r%'findexGF$F$C$?&Fi[sFfwFQC8>F2-Fibl6%FiqF[wFi[s@$0FctF:F`_l>F2Fcal>F=-Fibl6
%FiqF[clFgal@$F]coF`_l>F=Fjcu>FP-F]z6$RF_`lF$F$F$FehnF$F$F$F=>FD-Fb[q6#Fghl>Fb[
l-F]juF`b`l-FT6%FDFP.F_el>FJ-FL6$F_elF=>F=-Fen6$/F_elFdqF=>FP-Fen6$F\c`lFP>Fhio
Fgez>F[foFehp-FT6%Fi[sFhioFeb`l>F[fo-Fen6$Fai]lF[fo>FJ-Fen6$/FNF[foFJ>FD-Fjw6$F
[foF^x-FT6%Fb[lFPFeb`l>FXFgb`l@$4-Fiq6$-F]z6&F[hnF=FJFi[sF%C$>Fgjl,&F`esF:-Fjw6
$FehpF^xF\q-F\o6#7&F2FJFi[sFgjlF%F$F$F$F`^lRFe_mFhanF$F$C%>F2F[bn>F2-Fch_lF`bn@
%FbbnF__lF[^^lF$F$F$F\^zR6$%*op_matrixGFijt6%%)admatrixGF/FayF$F$C%>F2-&%'linal
gG6#%*transposeGFap?(F=F:F:-&F^f`l6#%'rowdimGFetFQ?(FPF:F:-&F^f`l6#%'coldimGFet
FQ>&F2Fdcm-%1desolv/adjointD2G6$F\g`lF[w-Fj^sFetF$F$F$%2desolv/gintorder1GR6$Fi
apF^ao6&%*psvar_facG%*fsvar_facGFjo%(fac_eqnGF$F$C'@%4F[eu>FJFafp>FJFbp>F2Fdq>F
=Fdq?&FPFJFQ@%554-FiqFi]w/F[wFP33F]r/FgfpF[wFcaq>F2Ffr>F=Fdcm7$-Fjw6$F`tFaqFe\\
lF$F$F$Fa^lRFe_m6.F_bl%%lfunGF]bq%'maxdepG%'maxfnsGF`v%%fsolGF/F[h]lFio%%ldegGF
dvF$F$C7>F2-Fjw6$F^etF`]l>FP-F]z6$RF`zF$F$F$-Fdt6#F^dnF$F$F$F2>FJFbb`l@$2FJ-Fdt
6#FggsF\]t>FXFdq?(F_elF:F:FdimFQ@$/-FLFg`wFJ>FX6$FX-FLFi]v>FXF\jp>F2Fdq>F=Fdq?&
FDFXFQ@%2F^o-Fjw6$-FL6$FNFDF^x>F26$F2FD>F=Fcjq>F2F`s>F=F`t@$FbtF\]t>F=-F]z6%RF_
\mF$FbzF$@%-Fiq6$F[wF\jlFdqFbpF$F$F$F=F2>Fhio-F]z6%RF_\mF$FbzF$-FdwFb`pF$F$F$F=
Fbp-FT6%-Fb[q6#F]crFhioF[fo>FD-FLFgd]l>Fb[l-F[\s6$FbpFD@$5/Fb[lFdq2F:-Fdt6#<#Fb
[lF\]t>Fgjl-Fb`_lF`]al@%4-Fiq6$FgjlF4C$@$F[_tFc_t7#/FDFb[lC$>Fi[s-Fc]n6#,&FbpF:
*&FgjlF:)FD-FdwF`]alF:F\q7$Fa^al7$FgjlFi[sF$F$F$%/desolv/map2indGRF_`l6'F+F`vF-
Fe]lF/F$F$C$@%/FizF\alC&>F2Fdq>F=Fbp?(F$F:F:F$/-FL6$F^oF=F\alC$>F26$F2F[t>F=Ffs
>F=Fg_alC&>FPF__o>FJF\_o>F2-F@6$&FJ6#&FPFjn/FXF`i^l>F=Fh`lFdbvF$F$F$%/desolv/re
ducelGRF[`q6,F``lF_blFgfzFhoF/FayFhdtFchqF_dl%*reduceDerGF$F$C(>Fi[sR6$F``lF]y6
$FchqF_dlF$F$C$?&F2F[wFQC$>F=-%/desolv/subsDerG6%FbpF2F:@$Ffa[l-F\oFj_lFdqF$F$F
$@$Fdb]lFjy>FJFghq>F=-Fibl6%Fiq-F\clFagpF4?(F$F:F:F$FhhtC%>Fb[l-F]z6%Fi[sF=F[w@
%F`g\l>F=-Fibl6%Fiq-F\cl6$-F``s6#-FL6$F:Fb[lFhzF4>F=-F]z6$RF\_nF$FbzF$-FL6#-Fib
l6%FiqFja]lF4F$F$F$Fb[l>FJ-Fe_lF]is-Fc]n6#-FipF^\lF$F$F$FhclR6%F^hnF_hnFfo6(Fho
F/%'fnord1G%'fnord2G%&countGFgepF$F$C&>FD-F]z6$RF`zF$FbzF$@%FfzFizFbpF$F$F$Ffw>
F2Fdhn@$FhewC%>FX-Fb[q6$-Fdt6#F^jl-Fdt6#Fihn?(F=F:F:FXFhew>F2,&-FL6$F=F^jlF:-FL
6$F=FihnF\q@$FhewC%-FT6%F\jlFDFc_[l-FT6%F[etFDFg\y>F2,&FJF:FPF\qF2F$F$F$Fje\lR6
'F]y%&lineqGFaa\lFfoFe]l6?%(lmoduleG%*flag_mod6GFg]l%,startmoduleGFcyF_^r%&nlso
lG%'nlineqG%'tlineqGFfh\lFhoFeo%%eqnsGFioFfdnFd]lF`a\l%'nleqn2GF]bl%%sys1G%%sys
2GF/Fay%'eqn_noGF^fuF_fuF`fuFafu%+flag_delayGF$F$C0>F27)Fgg\lFh\nFa^lF[\xFh]^lF
_bn%*decoupleDG>FDF`w>Fb[lFfw>F_elF[w>FXFbp?(Fi^wF:F:F`jqFQC$>Fdcr-Fbo6$&FX6#Fi
^wFD@%/F`_^lF:>FX-Fen6$/Fi^w-FLF^^sFXC(>FX-Fen6$/Fi^wFdqFX>F]drFdq?(8:F:F:F`_^l
FQC$>Fgjl7$-FL6#-Fen6$/F_ialFdqFdcrF`er>F]dr6$F]dr7&7$&Fdcr6#F_ialFeapFgjlFb[lF
D>F]dr7#F]dr@$Fgt-Fbu6#%8The~splitting~systems~:G-F\oFd_^l>Fi[s7#-F]dn6$F\qF`jq
@$Fhgn>FX-F]z6$FjgnFX>FX-FccsFjn>F=F^o>FJF:?(F$F:F:F$2F^o-FdtFb^tC(?(Fi^wF:F:Fc
[blFQ@$2&Fi[sF`halF:C$>FdcrF]hal@%Fbhal>FXFdhalC(>FXFjhal>F]drFdq?(F_ialF:F:F`_
^lFQC$>FgjlFbial>F]drFiial>F]drF_jal@$FgtFajalFdjal-Fch\l6&FXFi[s.FX.Fi[s>8;F:>
FhioF__l?(Fi^wF:F:Fc[blFQC$@%54-F_q6$F_halF^x2-Fdt6#F_halFfi[lC%>F[]blFi^w>F[fo
F_hal?(FPF:F:Ffh[l2FPFih[l@$2Fh[blFPC%@$/FPFg^l@$/-Fdt6#-Ff`m6$F[foFDF:C$>Fhio-
Fc^l6%F[foFDFbiu@$Fc_`lC$>FP""'F^]p>Fhio-F]_lFj^bl@%F]`x>Fi[s-Fen6$/Fi^wFPFi[sF
^]p?(FPF:F:Ffh[lF[^blC$?(F_ialF:F:Fc[blFQ@$2&Fi[sF]jalFPC'>F[]blF_ial>F[fo&FXF]
jal@$F`^bl@$Fb^blC$>FhioFi^bl@$Fc_`lC$>FPF^_blF^]p>FhioF`_bl@%F]`x>Fi[s-Fen6$/F
_ialFPFi[sF^]p@$Fc_`lF^]p@$Fc_`lF^]p@'Fc_`lC%@.F_[\lC->F]i[l-Fen6$/F:,&&F]i[lF_
[mF:F:F:F]i[l>8@F%>F_`r-FL6$F:-Faen6#Fd`x>FX-Fen6$/F[]blFdqFX>Fi[s-Fen6$FcbblFi
[s>F_el-F]z6%FcdnF_elFd`x@$-FT6$F^oF_elF\]t>F_el-Fibl6$RF\_nF$FbzF$FacpF$F$F$F_
el>Fb[l-F]z6%R6$F]_n%#seGF$FbzF$F_h_lF$F$F$Fb[lFd`x?(F_ialFc[blF\qF:FQ@$-Fiq6$F
_`blF_`rC$>F^br-Fcdn6$F_`blFd`x@'/F^brF^oC$>FX-Fen6$FgialFX>Fi[s-Fen6$FgialFi[s
2*&-F_io6#F_`blF:F\j[lF:-F_ioF^dt>FjablFQC$>FX-Fen6$/F_ialF^brFX>Fi[s-Fen6$/F_i
alF^oFi[s@%FjablC%>FX7$Feap,&F^bblF:-F``sF_bblF\q>Fi[s7$Fgf^lFhbq>FDFb`x>FD-Fe_
l6$/-F^im6#F^bblFdqFb`x/FPFNC&>F]i[l-Fen6$/FN,&&F]i[lF_\lF:F:F:F]i[l@%Fegn>Fhio
-F]z6$FjgnFd`x>FhioFd`x>FX7$F]cr-FL6#Fabbl>Fi[s7$-F]dn6$F^o-FdtF^cr-FL6#Febbl/F
PF;C%>F]i[l-Fen6$/F;,&&F]i[lF\elF:F:F:F]i[l>FjablF%@%2F:F[hblC->F^_r-Fen6$/F[]b
lF``xFX>F`^w7&F^_rF_elFb[lFD>F_elFhbbl@$F[cbl-F\o6#7#F`^w>F_el7$F`er&Fb`xF_[m>F
XFabbl>FX-F]z6%FcdnFXFd`x>Fb[l-F]z6%RFecblF$FbzF$F_h_lF$F$F$Fb[lFd`x>FD-Fe_l6$F
cfblFD>F]gv7&FXF_elFb[lFD-F\o6#7$F]gvF`^wC+>F_`rF\bbl>FXFabbl>Fi[sFebbl>F_elFhb
bl@$F[cblF\]t>F_el-Fibl6$RF\_nF$FbzF$FacpF$F$F$F_el>Fb[l-F]z6%RFecblF$FbzF$F_h_
lF$F$F$Fb[lFd`x?(F_ialFc[blF\qF:FQ@$FicblC$>F^brF]dbl@'F`dblC$>FXFcdbl>Fi[sFfdb
lFhdbl>FjablFQC$>FXF`ebl>Fi[sFdebl@%FjablC$>FX7$Feap-Fc]n6#-Fip6#F[fbl>Fi[sF^fb
l>FDFajbl/FPFGC'>F]i[l-Fen6$/FG,&&F]i[lF`c\lF:F:F:F]i[l@%Fegn>FfbrF`gbl>FfbrFd`
x>FX7$-FLF_j_lFegbl>Fi[s7$-F]dn6$F^o-FdtF_j_lF\hbl@&/&FhioF\elFN>FDFb`x/Fa^clF;
C$>FX7$FeapF[fo>Fi[s7$Fgf^l"#6F`^blC'>F]i[l-Fen6$/Fg^l,&&F]i[l6#Fg^lF:F:F:F]i[l
@%Fegn>FfbrF`gbl>FfbrFd`x>FXFh]cl>Fi[sF[^cl>FDFb`x/FPF^_blC$>F]i[l-Fen6$/F^_bl,
&&F]i[l6#F^_blF:F:F:F]i[l@%2F;F[hblC)>F^_r-Fen6$/F[]bl-FL6#Fa^clFX>F`^wF_ibl>FX
-Fen6$/F[]blFd`xFX>F_el7$F`er&Fa^clF_[m>FDFb`x>F]gvFdjblFejblC&@%Fegn>F^br-Fjgn
F_bbl>F^brFd`x>FX7$F^brFegbl>Fi[s7$F^oF\hbl>FDFb`x-Fgeu6)FPF[]bl7$F[foFd`xFXFb[
lFDFi[s@$F[cblF\]t-Fiq6$Fi[sFhbqC&-FT6%FhbqFi[s.F[]bl>F[fo&FX6#F[]bl>Fhio-Fa^lF
j^bl@$Fc_`lC'-Fbu6#%6Do~delay~substitutionG>FXFabbl>Fi[sFebbl>F_`rF\bbl?(F_ialF
c[blF\qF:FQ@$FicblC$>F^brF]dbl@%F`dblC$>FXFcdbl>Fi[sFfdblC$>FXF`ebl>Fi[sFdeblC,
@$52Fc[blFN2Ffh[lFih[lF^]p@%2F=FN>F=,&F=F:F:F:F^]p>F]i[l-Fen6$/Fih[l,&&F]i[l6#F
ih[lF:F:F:F]i[l@$Fcfu-Fbu6#%?module~7~applies~to~the~systemG@%Ffdcl>%.desolv_to
rderGF:>FfeclFN>Fhio-F]_l6%FXFDFbiu>FXFhio>Fi[sFfjal>FJF:@$FcfuC'>8<-F]z6$F_ioF
X>8>-Fb[q6#-FL6#Fafcl>8=-F]juFgfcl>8?-Fjw6$FafclF^x-Fbu6#-Feu6-F\goFc[blFiju-Fd
t6#-Ff`m6$FXFDF^[vFefclF_[vF[gclF_[vF^gclF`[v@$FgtC$-Fbu6#%6Complete~one~soluti
onG@$F]guC&-Fbu6#%7Unsolved~equations~areG-Fbu6#-F\vFjn-Fbu6#%.Solutions~areG-F
bu6#-F\vF_[q7#7&FXF_elFb[l-Ff`m6$7%FXF_elFb[lFDF$6(Ffh[lF\j[lF]i[lFfi[lFfeclFhg
nF$%0desolv/diffcofsGR6$%&Lder1G%&Lder2G6%FdevF`izFihzF$F$C&>F2-F_cn6%F]juF^jlF
ihn>F=-F_cn6%RFbcnF$FbzF$FdcnF$F$F$F2F^jl>FP-F_cn6%RFbcnF$FbzF$FdcnF$F$F$F2Fihn
7$F=FPF$F$F$Fe_zR6$Fd]zF`vF$F$F$@'Fdao-F]zFd_z/Ff^mF:Faap*&F^[mF:-F\al6$F[wF[\n
F:F$F$F$Fgg\lRFe_m6(FhoF`vF_^r%(int_conGFhdu%$lvnGF$F$C$@$-F_q6$FghqF^xF\]t@%/F
ehnF^o7$F^sF[wC(>FX-F_cn6%R6$F_hsFd\nF$F$F$7$F[wFbpF$F$F$F^jlF^dn>FPFdq>F2Fdq?&
FDFXFQ@$2F^oF[[pC%>FJ-F^dz6%F]]pF^dnF[w>FP6$FP-FL6#F^bm>F26$F2Fbcm>F=-Fen6$/FN7
#-F]dn6$F^oF]j`lFbp7$/F=-Fjw6$F`sF^x7$FibmFPF$F$F$%(solvedeGRFe_mFhanF$F$C%>F2F
[bn>F2-Fa^lF`bn@%FbbnF__lFcbnF$F$F$FeasR6$FeoFi[oFhanF$F$C$>F2Fghq@'F^gwC$>F2Fc
eu-Fc]n6#*&F2F:F[wF\qFfe[lFh^dlF^oF$F$F$Fbb^lR6&FeoFfoFe]lFijt6,Fd]lFjo%)new_co
ndGF/FayF_^rFho%&newfnG%&hipowG%#PxGF$F$C(>F2-F^c[l6%FghqF`wFfw>FDFdq>FPFdq>Fb[
lFdq?&F=F2FQC&>F_el-Faeq6$-Ff`m6$FbapFfwF[w>FD6$FD-F^imF]fo>Fgjl-Fdw6$Ff_pF`w@'
F_i_lF[o/FgjlF^oC$>FP6$FP,&FbapF:-F\al6$F_elF`wF\q>Fb[l6$Fb[l*&Ff_pF:F_elF:C&>F
P6$FP,&FbapF:-F\al6$F_el-F]dn6$F`w,&FgjlF:F:F:F\q>Fb[l6$Fb[l*&Ff_pF:-F\al6$F_el
-F]dn6$F`wFgjlF:>Fb[l6$Fb[l-F@6$*()F\qFXF:-F\al6$Ff_p-F]dn6$F`wFXF:-F\al6$F_el-
F]dn6$F`w,&FgjlF:FXF\qF:/FX;F:,&FgjlF:F\qF:>Fb[l6$Fb[l*()F\qFgjlF:F_elF:-F\al6$
Ff_pF_bdlF:@%Fc]alF^o7$7$-Fc]n6#-Fjw6$Fef\lF^xFP7$FibmFDF$F$F$FjaoR6$%(sep_facG
F_js6#%(dep_facGF$F$@'F[euC$>F2F_h\l7$Fbh\lF2Ff[zF_amF`jlF$F$F$F\eoR6%F]yFfoFe]
l6#Fg_yF$F$C%>F2-F]z6%F\bnFbpF[w>F2-FegalF`bnF[vF$F$F$F]auR6%FeoF^`lFijt6)Ffo%(
lfn_strGFhoF/Fh[v%(selDderGF_blF$F$C->FDRF`zF$F$F$33FfzFc`lFf^vF$F$F$>F2F[cl>Fb
[l-Fjw6$-Fibl6%F[`lF2F[wF`]l>F=-F]z6$Fe^oFb[l>FX-F_cn6%R6$F][x%$fnsGF$F$F$Fc[yF
$F$F$Fb[lF=>FX7$Feap-F@6$/-F_]uFj`yF_]uF`]u>FPFhft>F2-Fe_lFbi^l>F2-Fibl6$FDF2>F
X-F]z6$RF`zF$F$F$/Fbp-FiexFapF$F$F$Fgfl7$-Fe_l6$FXFPF=F$F$F$Fj_xR6&FeoFfoFe]lFb
\x6-F^`l%&lfnv1G%&lfnv2GF\g]lFhoF_\xF`^rFio%&oelfnGF]]nFa^rF$F$C(>Fi[sFggs>Fgjl
F`w>FXFbp>FDFdq?(F$F:F:F$F\\`lC/>F2-FL6$F:-F^ip6$FgjlFi[s>F=-Fibl6%R6$FazFh`^lF
$FbzF$4Ff[zF$F$F$F[wF2>FP-Fibl6%RFgidlF$FbzF$FhidlF$F$F$FgjlF2>FP-F]z6$RF`zF$Fb
zF$FizF$F$F$FP@$2F^o-Fdt6#-Ff`m6$-F`pFjnFPF\]t>FJ-Fc_m6%FXF=F`s>Fb[l-F]z6%RF^_q
F$FbzF$@%-FiqFag[lF^[mFdqF$F$F$FJFP>Fb[l-%/desolv/sortlenGF_[q?(F$F:F:F$2F^oFag
\lC(>FhioF^cal>Fb[l-Fen6$FcclFb[l@$FadyF`_l>FD6$FhioFD>Fb[l-F]z6&R6%F]_n%$co1G%
#v1GF$FbzF$-F\hl6#-F\al6$FccxFfwF$F$F$Fb[lFhioF2>FJ-F]z6&RFb\el6#%$texGF$F$C$>F
2-F\hl6#-F\al6$Fcf[lFfw@%FhewFdq-Fen6$/F:F2FbpF$F$F$FJFhioF2>FJ-F]z6%FjwFJFaq>F
X-FjwF_hp>Fi[s-Ff`m6$FXFi[s>Fgjl-Ff`m6$FXFgjl7$FXF[jsF$F$F$%)printDerGRF^cp6#%%
lDerGF$F$C&>F2F^et>F2-F]z6$RF\_n6#F+F$F$C$>F2-F_cn6%R6$Fd\nF_hsF$F$F$F\dnF$F$F$
F^dnF^jl@%/F2Fcam/FbpF\jl/Fbp&F\jl6#-FeuFc[qF$F$F$F2-Fbu6#FafwFdqF$F$F$F\bnRFdo
6)F_bl%&lderAGFgep%'lnameAG%&lfindG%'lfindAG%%larbGF$F$C*@%Fidl>Fb[lF[el>Fb[lFc
am>FP-F]z6$RF`zF$FbzF$@%FfzFizFbpF$F$F$F[w>FX-Fibl6%FiqFg[zFP>FX-F]z6$RFfilF$Fb
zF$/Fbp-FehzFapF$F$F$FX>FJ-F]z6$RF`zF$FbzF$@%FfzFizFbpF$F$F$Fb[l>FD-Fibl6%FiqFg
[zFJ>FD-F]z6$RFfilF$FbzF$FfaelF$F$F$FD-F]^t6#-Fe_l6$7$FeapF]]pFbpF$F$F$%&genncG
RFhjt61F[[uF\[uFeoF/F_[uFe]lF`vFfi^l%&gliscG%'geniscG%&lleqnG%%lconG%&tlconG%&l
caseG%%caseGF$F$C.>Fi[sR6&%%UvarGF^`lFijtFe]l6$F/%'isceqnGF$F$C$>F=-F@6$*&-&F`\
u6#&FfwFet6#FjbrF:-F\al6$-FbpFj`yFbdelF:/F2Fe\u,&FafzF:-Fj]\lFcdelF\qF$F$F$@%Ff
_^lC$>F^br7#-FL6#F_c\l>F^br-F]z6%RFe_elF$FbzF$,(F[wF:FbpF\qF:F:F$F$F$F^brFf\u>F
^br7#-F@6$,&Ff\uF:FJF\q/FJ;F^o,&Ff\uF:F\qF:Fd[u>%)spe_caseGFdfp>F2Fg[u>FDFfbm>F
X-F@6$-F_\uF[_zFd\u>FX7$FX-F@6$-F\]uF[_zF`]u>Fgjl-F]z6'Fi[sF[wF[wFfwF2>FhioFdq?
&F`crF^brFQC2>F[fo/-&F`\u6#&FfwFf^sFc[qF:?(FJ,&F`crF:F:F:F:Ff\uFQ>F[fo6$F[fo/-F
_\uFc[qF^o>F[fo7#F[fo>F_el-Fe_l6$F[foFgjl?(FJF:F:FcjpFQC$>Fb[l-F\al6$FagdlFegel
>F_el-Fen6$/FJ/Fb[l-F[\s6$&F_elF^\lFb[lF_el>F[fo-Fe_l6$/F\gnF]]pF[fo>F=-Fci^l6&
FbpF_elF[wFfw>F=-F]z6&F\bnF=FXFcam>Ffbr-F]z6%R6$F]_nFfoF$FbzF$7$-F\bn6%Fb^uF[wF
camF_`sF$F$F$F[foFX>F[fo-F]z6$RF\_nF$F$F$Fa^uF$F$F$F[fo>F=-F]z6&Fi`alF=FfbrFX>F
=-F`cu6$F=FX>F=Fccu>F=-F]z6%Ff[oF=FX>F=-Fhcu6#7$F`buFi_l>Fhio6$FhioF=7%7#FhioFX
FDF$6#F^felF$FgaalR6%F]jxFfdn%'statusG6#F\dmF$F$C%@$0F\jl-FL6$F:F`bvF[o>F2-Fjdm
6$F^jl-FL6$FNF`bv@)/F2FQ/Fbp-Febn6$FibmF^jl/F2Faem/F`bvF^hw/FfwF:F^oFbpF$F$F$Ff
`_lR6&F^hnF_hnFbvFfoF$F$F$@)F^sFQFhf_lF%2-Fhcl6%-F[d^lF`bpFbpF`wF^oF%FQF$F$F$F^
imRF_`l6$F+Fe]lF$F$@%/F_ealF^oF\jlC%>F=F^dn>F2-F_cn6%RFicnF$F$F$F\dnF$F$F$F=F^j
l@%Fafn-F\jl6#Fi_l-F\al6$F_^flF\gnF$F$F$Ff`mRFdoF$F$F$-Fibl6%R6$FazF]_nF$F$F$@%
Ffz-Fiq6$F[wFizFa`pF$F$F$F[wFbpF$F$F$%.desolv/rmfac2GRFdo6%FioFi[oFhoF$F$C%F\p>
F2Fdp@'F]\oC'>FPFdq?&F=F2FQ@$2F^oFj[\l>FPFf`p>FPF_]o-F]p6#*&F2F:FPF\q@%F_[\lFjr
FPFg]o@%Fj]oC$F`^oFcalFjrF2F$F$F$Ff^uRFdo6)%&resfnGF_bl%&LDderGF/%(new_eqnGF``l
FgepF$F$C'>Fb[l-F]z6$RF`zF$F$F$@%FfzFizFbpF$F$F$F[w>FXFghq>F=-F\cl6$FXFhz>F=-Fi
bl6%FdhyF=Fb[l@%Fehm-F[s6$F]sFbpC)>F=-F]z6$FehzF=>FPFfs>F=-Fen6$FcclF=?&FDF=FQ@
$2-Fhcl6%FPFDFb[lF^o>FPFD>FPFabt>F2-F[\sF`hdl@%Fir-F[s6$%=This~is~not~a~polynom
ial~PDEGFX/FP-Fc]nFetF$F$F$FfgoRF[`q63F_blFho%(result2GF/F`do%$DTkG%%pDTkGFhev%
&pHDTiGFbdal%#EiGFjev%'coeffrG%'acofriG%)exponentGFcdoF``lF$F$C(@$F^sF[o@$Fdb]l
Fjy@%Ff_^l>Fi[sF_c\l>Fi[sFcjp>F=Fghq?(FJFi[sF\qF:FQC)>F_el&F_]uF_\l>FgjlF_]m>F_
elF`\]l>F2-Fjho6$F=Ffw>F2-F]z6$RF\_nF$FbzF$@%FiilF[jlF`jlF$F$F$&F2F`c\l>F2-Fcin
Fi]n?&F]drF2FQC.@$4-Fiq6$F=)&F]drF_[m&F]drF_\lF`_l>FDFcefl>Fb[lFdefl@$0F]\yFd\m
F`_l@'/FDF_elC%@$2Fb[lFgjlF`_l>Fhio&F_]uF\el>F[fo-Feas6$Fhio)F_elFgjl/-Fjdm6$Fd
[al-FL6$FNF_elFQC$>Fhio-Febn6%F_elF_fflFd[al>F[fo-FeasF]\elF`_l>Ffbr-Feas6$F=)F
DFb[l>F^br-Fifl6$FfbrF[fo@%Fjefl>F`cr,&Fb[lF:FgjlF\q>F`cr,&Fb[lF:F\qF:>F=-F\hl6
#-Fc]n6#,&*&&F^brF_[mF:F=F:F:*(&F^brF_\lF:)FDF`crF:FhioF:F\q-F]p6#Fbhfl-F]p6#Fd
hfl@$FcftF[o-Fc]n6#*&F=F:-F[qFj_lF\qF$F$F$FegalR6%F]yFi^pFe]l60F/FcyFeev%%loopG
F^fuF\dmFeoF`vFfoFdfal%'hcoeffGF`fvFddo%'lanameGF$F$C5@$4-Fit6#Ffecl>FfeclF:@$4
-Fit6#Ff\w>Ff\wFQ>F]^wFgy@$2F]^wFNFjy>F^joFdfp@$4-Fit6#F]io>F]io"'++5>FcjvFdfp>
FddmF`av>Fgjl-FgaqFahl?&F_elFgjlFQ@%-F_q6$F_elFhz>&Fddm6#-FLF\cblFcam>&FddmF]fo
Fcam>F=-F]z6%Ff[oFbpFgjl?(F2F:F:F]^wFQC$>F[fo-Fjho6$FdbvFgjl@%2F]io-F_io6#&F[fo
F\elC(@$F][p-Fbu6#-Feu6$Fc[pFb\gl>Ffbr-F_fm6#-Ff`m6$Fd\glFfw>F^jo<$F`jo/FfbrFd\
gl>Fb[l-Fe_l6$/Fd\glFfbrFdbv>F[fo-Fen6%Fi[pFj[pF[fo>F=-Fen6$/F27%F2F[foFb[lF=C$
>F[foFh]gl>F=-Fen6$/F27%F2F[foFdbvF=?&Fb[lF=FQC$>Fhio-Feas6$&Fb[lF\el)&F[]mF_[m
&F[]mF_\l@%-Fiq6$FhioF4>Fi[s<$FhioFgf^l-F]pF^cr>FJF^o>FDFQ?(F$F:F:F$FDC*>FDF%>F
JFe\l@$FgdoC%-Fbu6#-Feu6$%&loop~GFJ>FX7#-F@6$-F_io6#Fdbv/F2;F:Ffhm@%/FXFcam-Fbu
6#-Feu6%%0~#~equations~:~GF`jq%5~;~lengths~:~[0,0,0]G-Fbu6#-Feu6+F`aglF`jq%/~;~
lengths~:~[G-Fb[qFaauF_[v-F]juFaauF_[vF`dwF`[v>F=-F[do6%F=FgjlFfw@$2F^o-FdtFajo
C%>F^br-F]z6$FaenF^jo>F^br-F]z6$RF`zF$FbzF$@%FfzFizFbpF$F$F$F^br?(F2F:F:FfhmFQ@
$-Fiq6$FdbvF^brC%>F]^wF^^w>Fb[l-Fen6$/F:F]^wFdbv>F=-Fen6$/F2Fb[lF=>FP-F`evF[bgl
@$FcimC$>F=7$FghlFi_l>FDFQ@$F]bglC'>F^brFabgl>F^br-F]z6$RF`zF$FbzF$@%FfzFizFbpF
$F$F$F^br?(F2F:F:FfhmFQ@$-Fiq6$F=F^brC%>F]^wF^^w>Fb[lF_cgl>F=Fccgl>F=-Fe_l6$F^j
oF=>F^joFdfp@$Fgdo-Fbu6#-Feu6&%0|*#~equations~:~GFfhmFa^p-F_ioFj_l@$F]bgl>F=F^e
gl-F]z6$RF\_nF$FbzF$F^_yF$F$F$F=F$6)F]^wFcjvFddmFfeclFf\wF^joF]ioF$FhgyR6#F`v6$
F]fxF``lF$F$@'FbdvC&>F2F^o>F=Fbp?(F$F:F:F$-Fiq6$F=F\alC$>F2Fbfo>F=FfsF2Fc`l@'Fc
iz-FdtFa_oFb_vFe_vF:F^oF$F$F$F[`uR6%FeoFijtF^`l6)FfoF/Fc`el%'eqnsmbG%%prolGF][u
FhoF$F$C->F2-Fibl6%F[`lF[clFfw>F2-F^[y6$FgflF\q>FP-F]z6$F\_alF2>FJ-Fa[v6&FbpF2F
PFfw>FX-F@6$7$-%&ZETA_G6%&F[wFj_lF[wFfwFfhgl/F=Fbjp>FX6$FX-F@6$7$-Fdhgl6%&FfwFj
_lF[wFfwF_igl/F=Fe\u>FX7$FX-F@6$7$-Fdhgl6%&FPFj_lF[wFfw-Fjw6$FhiglFjjm/F=F`i^l>
Fb[l-F]z6%R6$Fd]zF]_nF$F$F$Fh[^lF$F$F$FXFJ>Fb[l-Fjw6$Fb[lF^x>Fb[l-Fc]n6#-F]^t6#
-Fe_l6$/FdhglFg\\lFb[l>Fb[l-Fhix6%Fb[lF2FPF$F$F$%*diffSpolyGR6&%%eqn1GF[[^lFfoF
e]l6,%'lHead1G%'lHead2G%&Head1G%&Head2G%%acofG%&SpolyGF+%%dcofG%&deqn1G%&deqn2G
F$F$C)>F2-FjhoFew>F=-FjhoF`bp>FP*&)FfamFfhwF:&F2F\elF:>FJ*&)FbapFf_pF:&F=F\elF:
>FX-FiflF_jm@%/FfamFbap@%3/Fe\_lFP/Fb^nFJ>FDF^o>FD-Fc]n6#,&*&Fe\_lF:F[wF:F:*&Fb
^nF:FbpF:F\qC->F_el-Fbicl6$FfamFbap>Fb[l-F_cn6%F]dn-FL6$F;FfamF`\]l@%/Fag\lF^o>
FgjlFbp>Fgjl-F\al6$FbpF]`z>Fb[l-F_cn6%F]dn-FL6$F;FbapF_]m@%Fc^hl>Fi[sF[w>Fi[s-F
\al6$F[wF]`z>F2-Fjho6$FgjlFfw>F=-Fjho6$Fi[sFfw>FPFc\hl>FJFg\hl>FXF[]hl>FD-Fc]n6
#,&*&Fe\_lF:Fi[sF:F:*&Fb^nF:FgjlF:F\q-FjgnFCF$F$F$FjasR6$FiapF`vF$F$F$@'Fdao-Fi
bl6%R6$%"tGF]blF$F$F$Ff[zF$F$F$FbpFivFgvFbpF^oF$F$F$FcayR6$FiapFfdn6(Fe]lF+FhoF
/F``lF_blF$F$@%4Ff^[l-F\hl6#-Fc]n6#-Fe_lFb`pC)>FDF[cl>FD-Fibl6%FiqFDFi^[l>FD-F[
^m6$FD-Fibl6%FiqFDF`en>FD-Fibl6$RF`zF$F$F$Fc`lF$F$F$FD>FPFbp?&FXFDFQC&>F27#Feap
>F=7#-FL6#-FL6$F^oF^c_l>F=-F@6$&F26#Fg[]l/FJFi`gl>FP-Fe_l6$/FX-F\al6$Fb[^lF=FP-
F\hl6#-Fc]n6#-Fe_l6$F[wFPF$F$F$F^g`lR6$%&doperGFijt6%Fiap%&AdjopGFfoF$F$C(>F=R6
$Fd]zFe]lF$F$F$@%/Ff^mFcam*&F^[mF:-%$QQ_GFh]uF:*&)F\q-FdtF\\nF:-F\al6$FjdhlF[\n
F:F$F$F$@%Fgbz>F2-F\hl6#-F]z6%F=FbpF[w>F2-F=Fi`m>F2F_ju>FP-Fibl6%FiqFf\vF\ehl>F
P-F^[y6$-Fjw6$FPF`]lF\q-F^a[lFfrF$F$F$Fc[elRFecs6&F_yF\bqF^fuF/F$F$C(>F2Fgy>F=F
gbq>FP7#-F@6$,&*&-F_io6#&FbpF^\lF:F=F:F:FJF:/FJFbcq>FP-FecqF\r>FP-F]z6%RF_\mF$F
$F$FjcqF$F$F$FPF=7#-F@6$&FbpF]i^lFbghlF$F$F$FdjrR6%%(vector1G%(vector2GFe]l6$%%
vec1G%%vec2GF$F$C%>F2-F`jmFew>F=-F`jmF`bp-Fc]n6#,&-Fe_l6$/-Ff[nFj`yF=F2F:-Fe_l6
$/F_ihlF2F=F\qF$F$F$F]`^lRF^^^l6'F`vF,F\dmFhdu%'lderivGF$F$C%>FJFdq?&F=FfwFQC%>
FPFQ?&F2-Fibl6%FiqF[wF=FQC$>FX-Fdgx6%FbpF2F=@$5/FXF%/F`esF^oC$>FPF%F^]p@$FP>FJ6
$FJF=F^hzF$F$F$Fh[lR6'F^hn%&pder1GF_hn%&pder2GFfo6'FhoF/F`dalFadalFbdalF$F$C$@'
/FbpFfw>F2,&F[wF:F`wF\q3Fhifl/FfeclFNC&-FT6%F\jlFbiuFc_[l-FT6%-FL6$F:FfwFbiuFg\
y>F2Fafal@$Fhew?(F=F:F:F_ealFhew>F2,&FfealF:-FL6$F=F\xF\qC$>F2,&FehnF:FiwF\q@$F
hewC%>FX-Fb[q6$F_eal-Fdt6#F\x?(F=F:F:FXFhew>F2Fb\il@$FhewC%Fh[ilFj[il>F2FafalF2
F$FiiflF$%)map2DerAGRFghz6)F``lFgdtFihzF+Fe]lFjhzF/F$F$@)Fb`l-F`hmF[[[lFc`lC$>F
2Fbp@'-F_q6$-FL6$F^oFgalFeizC&>F=7#-FL6#F^^il>FDFfhm>FP7#-F@6$,&FDF:-Fdt6#-Fe_l
6$Fc]alF=F\qF[][l-F`hm6%-FL6#FgalFPF_\[l-F_q6$F^^ilFhz-F`hm6%Fa_il7#-FL6$FNF^^i
lF_\[l-F`hm6%Fa_ilFhjzF_\[lFj`lC*>FJFdq>F2Fbp?(F$F:F:F$F\iyC$>FJFc[[l>F2Fcal>FJ
F^hz>FDF\`o>FXF_\[l>F=Fa\[l-F`hmF]][l-F`hmF_][lF$F$F$F`cuRF\y6/%)sim_leqnG%'l1t
ermGFics%)zero_eleGFeoF/F]bl%&oldfnGF__dlF^`lFjoFd]l%%lfnDGF$F$C(>Ffbr-F]z6$Feh
zF[w>F2Fbp>F=Fcam>FPF^o?(F$F:F:F$0FPFctC'>FPFct>FDF^o?&FXF2FQC$>FDF[^w@$4FedsC%
@$-F_q6$FXFaq>FX-F`[\l6$FXFfbr@%5Fcbil-F_q6$FXF_rC$>F2-Fen6$/FDFdqF2>F=Fh_lC%>F
2-F]z6%FddtF2FX>F=-F]z6%FddtF=FX>F=Fh_lF^]p>F2-Fibl6$RF\_nF$FbzF$@%F^sF%FQF$F$F
$F2>F=-Fibl6$RF\_nF$FbzF$@%F^sF%FQF$F$F$F=FefsF$F$F$F\^nRFe_m6+F]]nF\]n%&lnargG
%(numfuncG%%lnumGFcipFhoFdvF/F$F$C,>F2Fggs>F=Fh`m>FP-F]z6$RF`zF$FbzF$@%FfzFgyF^
oF$F$F$F=@$/FctFbb`lF\]t>FJRF^`pF$F$F$-Fjw6$-F]z6%RFgidlF$FbzF$@%Ff[zF:F^oF$F$F
$F[wFbpF^xF$F$F$>FX-F]z6%FJF2F=>FDFgagl>Fb[lFdq?(FgjlF:F:F`jqFQ@$/&FXFd_tFD>Fb[
l6$Fb[l&F2Fd_t7$Fef\lFDF$F$F$F_bnRFe_m6.%)ode_funcG%(ode_ordG%(ode_varG%(ode_eq
nGFhoFioF`vF_^rF__dlF[h]l%'HDTodeGF/F$F$C1>F2-Fc``l6%FbpFh`mFfw@$/F2F%F\]t>Fb[l
Ffam>FPFe\hl>F=Fidfl>F2Ffhw@$4-FdjwFheoF\]t>FJ-Fcdw6&FbpFb[lF2FP>FX-%'dsolveG6$
FJ-F[fwF\r@$5/FXFdq-Fiq6$FX<%F[d^l%&DESolG%$sumGF\]t>FX-Fip6#-Fe_l6$/FchilF2,&-
FaenFjnF:-F``sFjnF\q>F_elFdq?(F[foF:F:F=FQ@$-Fiq6$FX(%#_CGF[foC%>Fgjl-Faeq6$-Fe
_l6$Fijq7#F]`zF[w>F_el6$F_el-F^imFd_t>FX-Fe_l6$/FjiilFgjlFX>Fi[s-FdjyFheo@%4-Fi
q6$Fi[sF4C$-F]pFb^t7$FX7$FibmF_elC%@%2F^oF=>Fhio-F\al6$F2-F]dnFf`p>FhioF2>FD-Fc
]n6#,&FbpF:*&Fi[sF:FhioF:F\q7%FXFc[jl7$Fi[sFDF$F$F$F$
