MVR5
I'Desolv=6"%&falseGE\[lby%%icdeGR6%%%leqnG%$lfnG%%lvarG6%%&tleqnG%*LSwitchesG%#
iiG6#%QCopyright~(c)~1998~by~Khai~Vu~and~John~CarminatiGF$C(@%2""$9#>8%7#-%$seq
G6$&9"6#8&/FA;""%F7>F97"-%6desolv/EvaluateSwitchG6#F9>8$-%$mapG6%%2desolv/conve
rtm2DG9$9%>FK-FM6$R6#%#exGF$6$%)operatorG%&arrowGF$@%-%%typeG6$FP%"=G,&-%$lhsG6
#FP"""-%$rhsGF]o!""FPF$F$F$FK>FK-%/desolv/icdeDerG6%FKFQ9&-%2desolv/convertD2mG
6#FKF$F$F$%<desolv/DifferentialVariableGR6#%%FuncGF$6#%>Copyright~(c)~1997~by~K
hai~VuGF$-%$zipG6%%"$G-%#opG6$F6FP-Fep6$""#FPF$F$F$%<desolv/RemoveVanishUnknown
sGR6'%%LEqnG%&LIneqG%%LSolG%$LFnG%)LastFlagG6.%+LVanishUnkG%'EqnFnsG%'SolFnsG%*
NotSolFnsG%&TUnknG%%NLFnG%%TFnsG%&TFns2G%%TEqnG%&NLEqnG%'NLIneqGF06#%HCopyright
~(c)~September~1999~by~Khai~VuGF$C(>F9-%(convertG6$-%/desolv/fmemberG6$FP9'%$se
tG>FA-Fcr6$-Ffr6$FfoFhrFir@%39(/-%%nopsG6#-%'selectG6%FgnFA%)functionG""!C$>FK-
Fgs6%FgnF9Fis>F9-%&minusG6$F9FK>FK<">8'-Fat6$F9FA?(8/F^oF^o-FdsF]o%%trueGC&>8,-
Fep6$FjtFP>8*-Fcr6$-Ffr6$F_uFhrFir>8+-Fat6$FcuFft@$/-Fds6#FiuFjsC$>Fiu-Ffr6$-%'
subsopG6$/Fjt%%NULLGFPFcu@%F]v>FK-%&unionG6$FKFcu>FK-F\w6$FK-%>desolv/CombinedV
anishUnknownsG6&F_uFdvFfoFhr@%/-FdsFioFjs7%FPFQFhrC+>8)Fhr?&8(FKF\u>Fjw-%%subsG
6$/F\xFhvFjw>FK-FM6$%/desolv/map2DerGFK>8-FP?&F\xFKF\u>Fgx-FM6%%0desolv/subsfun
DGFgx/F\xFjs>Fgx-FM6$RFVF$FXF$@%/FPFjsFhvFPF$F$F$Fgx>8.FQ?&F\xFKF\u>Fey-FM6%F\y
FeyF]y7%FgxFeyFjwF$F$F$%Ddesolv/ConditionalHighestDerivativeGR6%%$EqnGF]p%*Orde
rListG6&%(ListDerG%&HiDerG%%TDerG%(HiPowerGF^pF$C(>FK-%5desolv/GetDerivativeG6$
FQFP@$2FfwF^o-%'RETURNG6#FF>F9-Fep6$F^oFK?&FAFKF\u@$2Fjs-%3desolv/HigherOrderG6
'FAF^oF9F^oFfo>F9FA>Fft-%'degreeG6$FPF97$F9FftF$F$F$%2desolv/NormalFormGR6&%&ol
eqnGF+F,F_z6.F0%#jjG%#kkG%$eqnGF[r%%LderG%&nLderGF*%&eqnumG%%lHDTG%%narbG%'nume
qnGF1F$C)@$3-%)assignedG6#&%*infolevelG6#F'2FjsFb]l-%&printG6#%EReduce~the~syst
em~into~a~normal~formG>Fiu-%/desolv/SortHDTGFgr>FKF^o>FjtF^v?(F$F^oF^oF$2FKF^vC
/>FK,&FKF^oF^oF^o@$2F^vFjt@$F^]lC$>FjtF^v-Fg]l6#-%$catG6$%/#~equations~:~GFjt>F
_u&&FiuFio6#F^o>Fjw&Fa_l6#Fip>Fft-%.desolv/ReduceG6%&Fa_l6#F67#-F<6$&FiuF@/FA;F
^o,&FKF^oFaoF^oFhr@$/FftFjsC%>Fiu-Fev6$/FKFhvFiu>FKFb`l%%nextG>Fft-%-desolv/rmf
acG6$FftFQ>Fcu-%1desolv/diffparseG6$FftFhr@$2%/DesolvSubLimitG-%'lengthG6#&FcuF
[`lC'-Fg]l6$%3before~replacementG-Fhal6#Fft>Fey-%.desolv/newarbG6#-Ffr6$FjalFfo
>%/DesolvSubsListG7$/Fey&%(eqnHeadGF[`l-Fep6#Fhbl>Fft-%'expandG6#,&FftF^o*&,&Fe
yF^oFjalFaoF^o)&FcuFb_l&FcuFe_lF^oF^o@$3F_]l2FipFb]l-Fg]l6#-F\_l6$%UArbitrary~r
eplacement~;~The~coefficient's~length~is~GFgal>Fft-%.desolv/mysimpGF_al>Fcu-Fev
6%/F6Fhv/FDFhvFcu>Fft7%F_uFcuFft@%/FcuFjw>Fiu-Fev6$/FKFftFiuC&>FiuFg`l?(F9F^oF^
oF^vF\u@$1-Ff[l6%-Fep6#Fcu-Fep6#&&FiuFIFe_lFhrFjs%&breakG@%2F^vF9>Fiu7$-FepF_vF
ft>Fiu-Fev6$/F96$FftF]flFiu>FKF9@$3F_]l2F^oFb]l-Fg]l6#-F\_l6&F^_lF^v%3~;~total~
length~:~G,&-FhalF_vF^o!"$F^oFiuF$6$FhblFfalF$%Gdesolv/OrderVarsOfFunctionByAut
oWeightGR6#%1FunctionNWeightsG6'%-ListVariableG%+ListWeightG%'VDigitGF0%,MainWe
ightsG6#%QCopyright~(c)~1998~by~John~Carminati~and~Khai~VuGF$C+@$/-Fds6#&FPFe_l
Fjs-F][lF]o>FKFfhl>F\x&FPF[`l>F\x-FM6$R6#%"aGF$FXF$-%%evalG6#,&"%****F^oFPFaoF$
F$F$F\x>FA-Fhal6#Ffw>F97#-F<6$,&*&&F\xF`blF^o)"#5FAF^oF^oFftF^o/Fft;F^oFfw>F9-%
%sortGFI>F9-FM6%%$modGF9F`jl7$&FPFb_l-FM6%FepF9FKF$F$F$%;desolv/ExtendInfinites
imalGR6%%%funcG%%ivarG%%dvarG6*F0%%var_G%'resultG%)zeta_km1G%&fnameGF,%&ldvarG%
$posG6#%>Copyright~(c)~1998~by~Khai~VuGF$C$>Fjw-Fep6#-%-desolv/gendlG6%FQFfoFjs
@'4-Fgn6$FP%(indexedG@'-%'memberG6%FPFQFiu-&%9Desolv_Infinitesimal_IndGF_v6#Fjw
-F\]m6%FPFfoFiu-&%9Desolv_Infinitesimal_DepGF_vFa]mF^o/F[uF^oC$>F\x-Fep6$FjsFP@
%-F\]m6%F\xFfoFiuC%>FA7#-F<6$*&-%%diffG6$-F\x6#-Fep6#FQ&FQFioF^o-Ff^m6$-&F`]mFi
oFa]m-FepF]oF^o/FK;F^o-FdsF[_m>FA-Fcr6$FA%"+G,&-Ff^m6$Fd]mFa_mF^oFAFaoF^oC(>F\x
Fj]m>Fcu7#Fa_m>Fft-F_[m6%-Fev6$/F^oFhvFPFQFfo>FA7#-F<6$*&-Ff^m6$Fh^m-Fep6#-Fev6
$/F^oF\_mFcuF^o-Ff^m6$F__m-Fep6$F^oFcuF^oFb_m>FAFf_m,&-Ff^m6$FftFdamF^oFAFaoF$6
$F`]mFf]mF$%2desolv/intlintestGR6&Fe\lFb[m%$ldfG%%xvarG6#F[\mF^pF$@+4-%$hasG6$F
P-Fep6$FjsFQ%&nofctG4-Fdbm6$FQFhr%'noxdepG0-F[\l6$FPFfoF^o%'nonlinG2Fjs-Fcr6$-F
ep6$FipFfoFh_mC$-F\]m6%Fhr-Fep6$F6FfoFK@%2Fjs-Fep6$FKFdcmF%%'noxderG%&noderGF$F
$F$%-desolv/separGR6%Fe\lF+F,6'%$eqlG%'eqlistG%&DIfacG%'tDIfacG%'mulfacGF^pF$@%
43-Fgn6$FP%)fractionG2Fjs-Fds6#-Ffr6$-Ffr6$FPFQFfoC&>FK-%-desolv/eqsepG6%7$F^oF
PFQFfo>F9Fhv?(F$F^oF^oF$0FKFFC&>FA&FKFb_l>FK-Fev6$Fe`mFK?&FftFKF\uC&>F\x*&&FftF
e_lF^o&FAFe_lFao@$2Fjs-Fds6#-Ffr6$F\x7$-Fep6#FfoFj^mF[al>FK-F_x6$/FftFhvFK>FA7$
,&&FAFb_lF^o*&F\xF^o&FftFb_lF^oF^oF[gm>F9Fht7#F97#F\fmF$F$F$%7desolv/FindInclus
iveICGR6&F*F+F,F_z6,F0Fc\l%'neweqnG%(num_eqnG%&newicG%%HDTiG%%HDTjGF\clFj\l%&Nu
mICGF1F$C*@$F^]l-Fg]l6#%ICalculate~integrability~conditions~(ICs)G>FftF[u>F\xFh
v>FgxFjs?(FKFipF^oFft2Fgx%2DesolvIterationICGC$>Fjw&&FPFioFe_l?(F9F^oF^oFb`lFei
mC2>Fcu&&FPFIFe_l@$0-Fep6$F^o&FjwFb_l-Fep6$F^oFgclF[al>FA-%:desolv/CrossDiffere
ntiateG6$FjimF_jm@$/FAFjsF[al@$%-DesolvReduceGC$@$Fjcl-Fg]l6#-F\_l6$%+Reduce~IC
~G7$FKF9>FA-Fh_l6%FAFPFhr@$F\[nF[al@$Fjcl-Fg]l6$%8New~IC~whose~length~is~G-Fhal
F@>Fgx,&FgxF^oF^oF^o>FA-F^al6$FAFQ>%/DesolvEqnCountG,&Ff\nF^oF^oF^o>Fiu-Fbal6$F
AFhr@$2Ffal-Fhal6#&FiuF[`lC&>F_u-Fcbl6#-Ffr6$FAFfo>Fhbl7$/F_uF_]nF]cl>FA-Facl6#
,&FAF^o*&,&F_uF^oF_]nFaoF^o)&FiuFb_l&FiuFe_lF^oF^o@$Fjcl-Fg]l6#-F\_l6$F`dlF]]n>
FA-FcdlFd\n>Fiu-Fev6%FgdlFhdlFiu>FA7%Ff\nFiuFA>F\x6$F\xFA>F\x7#F\x@$F[glC$-Fg]l
6#-F\_l6&%4~#~new~equations~:~G-Fds6#F\xFagl,&-FhalFj_nF^oFdglF^o@$F[dl-Fg]l6#-
F\_l6$%+~length~:~G-FM6$FhalF\xF\xF$6(%0DesolvICHasDoneGFf\nF^[nFhblFfalFfimF$%
<desolv/ConvertNonDerivativeGR6#F^z6$%%LDerG%%LSubGF\\mF$C&>FK-Fcr6$-%8desolv/G
etAllDerivativeGF]o%%listG>FK-Fgs6$R6#%"fGF$FXF$-%&evalbG6#/-%9desolv/Different
ialOrderGF]oFjsF$F$F$FK>F9-FM6$RF_ilF$FXF$/FP-Fep6$F^oFPF$F$F$FK-Facl6#-Fbil6#-
F_x6$F9FPF$F$F$%)separateGRFcdmF$F^pF$-FM6$RFVF$F$F$FfhlF$F$F$-Fadm6%FPFQFfoF$F
$F$%;desolv/SelectLHSDerivativeGR6$F]q%$LfnG6$F^z%&SLEqnGF\\mF$C'>F9Fhv?&FKFPF\
u@$3-Fgn6$FKFin-Fgn6$-F\oFioFis>F9Fbt>F9F`hm@$0-FdsFIF[u>F9-%9desolv/PickLHSDer
ivativeGFfem-FM6%FOF9FQF$F$F$%-desolv/RmfacGR6$Fe\lF+6$%%lfacG%$facGF^pF$@%-Fgn
6$FP%"*GC%>FKFhv?&F9FPF\uC&@$/-Fds6#-Ffr6$F9FQFjsF[al@$-Fgn6$F9%"^G>F9-Fep6$F^o
F9@$3-Fgn6$F9Fis-F\]mFafnF[al>FK6$FKF9-Fcr6$7#FKFgenFPF$F$F$%5desolv/CollectToL
istGR6$F^z%*CVariableG6#F[rF^pF$C&>FK-FaclF]o@%-Fgn6$FKFh_m>FK-Fcr6$FKFcan>FKFb
gn>FK-FM6%%6desolv/SeparateToListGFKFQ-%3desolv/CollectListGFioF$F$F$%/desolv/e
qsortDGR6#F*6'%'sum_ntG%'ot_ordGF0FghmFe\lF^pF$C.>FftF[u@$2FftFipFP>FKFhv>F9Fhv
>FAFjs?&F\xFPF\uC$>FA,&FAF^oF^oF^o@%-Fgn6$F\xFh_m>FK6$FK,&*&Fi_nF^o)FajlF_blF^o
F^oFAF^o>F96$F9,&*&-Fcr6$-Fep6$FipF\xFh_mF^oF`jnF^oF^oFAF^o>FK-Ffjl6#Fbgn>FK-FM
6%R6$F`il%"bGF$F$F$-%%iremG6$FP)Fajl-FhalF[_mF$F$F$FKFft>F9-Ffjl6#F`hm>F9-FM6%R
F`[oF$F$F$Fb[oF$F$F$F9Fft>FK7$-FepFI-FepFio7#-F<6$&FP6#&FKF@/FAFcjlF$F$F$FHR6#F
/6#%(TSwitchGF1F$C0>FK-Fgs6%FdbmFP%)simplifyG@%32FjsFfw-Fgn6$-F`o6#Fbfm%(numeri
cG>%/DesolvSimplifyGFg]o>F[^oF^o>FK-Fgs6%FdbmFP%)excludedG@%Fc]o>%2DesolvExclus
iveICGFg]o>Fc^oFjs>FK-Fgs6%FdbmFP%*iterationG@%Fc]o>%0DesolvIterationGFg]o>F[_o
%)infinityG>FK-Fgs6%FdbmFP%*redundantG@%Fc]o>%0DesolvRedundantGFg]o>Fd_oFjs>FK-
Fgs6%FdbmFP%)orderingG@%Fd]o@)-Fgn6$Fg]oFcanC$>%/DesolvOrderingG""&>%0DesolvOrd
erListGFg]o-F\]m6$Fg]o7&F^oFipF6FD>F``oFg]o/Fg]o%%autoG>F``oFjs>F``oF^o>F``oF^o
>FK-Fgs6%FdbmFP<$%,IterationICG%,iterationicG@%33Fd]o-Fgn6$FbfmFinFe]o>FfimFg]o
>Ffim"$+"@$4-F`]l6#F^[n>F^[nF\u@$4-F`]l6#Ffal>Ffal"'++5F$6+Fc^oF[_oF^[nFc`oFfal
F``oF[^oFd_oFfimF$%0desolv/diffcofsGR6$%&Lder1G%&Lder2G6%%'maxindG%&lind1G%&lin
d2GF1F$C&>FK-Fap6%%$maxGFgp-Fep6$FipFQ>F9-Fap6%R6$%"xG%"yGF$FXF$-Fbil6#,&FPF^oF
QFaoF$F$F$FKFgp>FA-Fap6%RF[doF$FXF$F^doF$F$F$FKFeco7$F9FAF$F$F$%4desolv/Functio
nNameGRF\pF$F^pF$FebnF$F$F$%3desolv/intminderivGR6%%%expnGF+F,6'Fb[m%$varG%%fla
gG%&livarG%'lderivGF^pF$C%>FftFhv?&F9FfoF\uC%>FAF\u?&FK-Fgs6%FdbmFQF9F\uC$>F\x-
%/desolv/lderivxG6%FPFKF9@$5/F\xF%/FejnFjsC$>FAF%F^fl@$FA>Fft6$FftF97#FftF$F$F$
%.desolv/rmfac2GRF`en6%%%teqnG%%termGFg[mF^pF$C%-%1desolv/storecaseG6#-%&denomG
F]o>FK*&-%'factorG6#-%&numerGF]oF^o-%)icontentGF]oFao@'-Fgn6$FKFgenC'>FAFhv?&F9
FKF\u@$2FjsF^fn>FA6$FAF9>FA-Fcr6$7#FAFgen-F_go6#*&FKF^oFAFao@%/FAF^o-%&ERRORG6$
%.contradictionGFcyFA-Fgn6$FKFefn@%2Fjs-Fds6#-Ffr6$F`[lFQC$-F_go6#-Fep6$FipFKF`
[lF_ioFKF$F$F$%0desolv/isderivDGR6$%%der1G%%der2G6$%'lindexG%'linnegGF^pF$C%>FK
-Fap6%RF[doF$FXF$F`doF$F$F$FPFQ>FK-F_x6$/FjsFhvFK@%Few.%&equalGC$>F9-Fgs6%FgnFK
%)negativeG@'/FhdnFjsF\u/FhdnFfw.%'negderGF%F$F$F$%/desolv/map2indGR6#%$derG6'F
j[mFb[mFfjoF,F0F^pF$C$@%/Fj]mFf^mC&>FKFhv>F9FP?(F$F^oF^oF$/-Fep6$FjsF9Ff^mC$>FK
6$FK-Fep6$FipF9>F9Fgfn>F9F[]pC&>FA7#-Fep6#-Fep6$FjsFj]m>FftF_`m>FK-F<6$&Fft6#&F
AFj_n/F\x;F^o-FdsF@>F9-Fep6#Fj]m&F9FioF$F$F$%0desolv/condoperGRF`\p6&%%fvarGFi[
m%#l1G%(ord_derGF^pF$@%-Fdbm6$Fj]m%"DGC$>FAFi]p@%-Fdbm6$FAFb_pC'>FKFa_m@&-Fgn6$
FKFis>FK-Fep6$FjsFK4-Fgn6$FK%%nameGFghl>F9Ff^p@%2F^oFd^p>Fft-Fep6$FipFA>FftF^o-
Ff^m6$-F9Fio-Fcp6$FKFftFPFPF$F$F$%2desolv/transform1GR6'Fe\lF+%%lstrGFc[mFd[m6+
Fg[mF0%%lderGFa\p%%tderG%'lfdvarG%$funG%'lsubfnG%$lfdGF^pF$C->Fiu-Fap6%RF[doF$F
XF$/FP-Fcr6$FQ%'stringGF$F$F$FQFfo>FK-F_x6$FiuFP>F_u-%'indetsGF\`p>Fjw-Fgs6%R6$
Fian%#ldGF$FXF$-F\]m6$Fj]mFQF$F$F$F_uFas>F_u-Fat6$F_uFjw>Fiu-FM6$RFhanF$FXF$/FP
Fj]mF$F$F$Fjw>FK-F_x6$FiuFK>F_u-F_x6$FiuF_u>F_u-Fgs6$RFhanF$FXF$33F`_p-F[bn6#Fg
]m-Fgn6$Fa_mFc`pF$F$F$F_u>Fiu-FM6$RF`\p6$%'ordderGFb[mF$F$C%@%-Fgn6$Fi]pFis>FK-
Fep6$FipFi]p>FKF^o>F9-Ff^p6#Fa_m/FP-Ff^m6$F9-Fcp6$F`\oFKF$F$F$F_uF[dpF$F$F$%7de
solv/SubsDerFunctionGR6$Fe\l%'subeqnG6'Fj[mFg[mFa\pFeap%%sfunGF^pF$C(>FftFaan>F
ft-Fgs6%FdbmFft-Fep6$F^o-F\oF[_m>F\xF]gp>F9FP?(F$F^oF^oF$2Fjs-FdsF`blC*>FA-Fep6
$F^oFft>Fft-Fev6$Fe`mFft>FK-Fap6%RF[doF$FXF$F^doF$F$F$Fh`pFgjn@$-Fdbm6$-FM6$RF_
ilF$FXF$-F[bn6#2FPFjsF$F$F$FKF\uF[al>FK-Fap6%Fcp-Fep6$F6F\xFK@%Few>F9-Fa]o6#-F_
x6$FQF9>F9-Fa]o6#-F_x6$-Ff^m6$FQFa\oF9>Fft-FbanFI>FftFifpF9F$F$F$%.desolv/gcdDe
rGR6$%'Hterm1G%'Hterm2G6+%&acof1G%&acof2G%$ex1G%$ex2G%&sfac1G%&sfac2GFcen%%tfac
G%&scindGF1F$C0>FA-Ffgo6#-Fa]oF]o>Fft-Ffgo6#-Fa]oF[_m@%-Fgn6$FAFgen>F\x<#-FepF@
>F\x<#FA@%-Fgn6$FftFgen>Fjw<#-FepF`bl>Fjw<#Fft>F\x-FM6$RF_ilF$FXF$@%-Fgn6$FPFef
n7$FebnFgp7$FPF^oF$F$F$F\x>Fjw-FM6$RF_ilF$FXF$@%Fh\qFj\qF[]qF$F$F$Fjw>F_u-%*int
ersectG6$-FM6$RF_ilF$FXF$F\[mF$F$F$F\x-FM6$RF_ilF$FXF$F\[mF$F$F$Fjw>FKFhv>F9Fhv
?&FcuF\xF\u@%-F\]m6$FgclF_uC&>Fiu-Fgs6%RF`[oF$FXF$@%/F\[mFQF\uF%F$F$F$FjwFgcl>F
iu-Fep6$F^oFiu@%33-Fgn6$FhclFi]o-Fgn6$Fa^nFi]o2Fa^nFhcl>F96$F9)Fgcl,&FhclF^oFa^
nFao>FK6$FK)Fgcl,&Fa^nF^oFhclFao>Fjw-Fat6$Fjw<#Fiu>F96$F9Ffcl>Fjw-FM6$RF_ilF$FX
F$)F\[mFfhlF$F$F$Fjw>FK*&-Fcr6$<#FKFgenF^o-Fcr6$FjwFgenF^o>F9-Fcr6$<#F9FgenFf[n
F$F$F$FehnR6$%%TermG%*SVariableG6%%+ProductSetG%1DependentFactorsG%3Independent
FactorsGF^pF$C&@%Feen>FK-Fcr6$FPFir>FK<#FP>F9-Fgs6%FdbmFKFQ>FA-FatF_gn7$-FcrFe[
q-Fcr6$F9FgenF$F$F$%:desolv/EquateCoefficientsGR6$F^z%%LVarG6#%*ListCoeffGF^pF$
C%>FK-%.desolv/sumsepGFfem>FK-FM6$RFVF$FXF$FfhlF$F$F$FKFKF$F$F$%-desolv/splitGR
F`en6%Fg[mF[goFcenF^pF$C'F^go>F9Fdgo@'-FgnFhbqC&>FKFhv?&FAF9F\u@'4-Fdbm6$FA%$De
rG-F_goF@-Fgn6$FAFefnC$-F_go6#Fh`p>FK6$FKFA>FKF]eq@$Fi`lF_io>FKFbgnFcfn@%-Fdbm6
$FgfnFedqC$-F_go6#F`]p>FK7#FgfnF_io>FKF`hm@$2F^oFfw@$3-F`]l6#&Fc]l6#%(pdesolvG2
FjsF`fqC$-Fg]l6#-F\_l6%%6Braching~occurs~with~GFfw%+~factor(s)G-Fg]l6#FgoFKF$F$
F$%4desolv/SortDerPowerGR6$FeapF_z6'%)SortListG%+RemainListG%&TListG%'numderG%)
DerOrderGF^pF$C)>FftF[u@$FainFghl>FKFhv>F9-Fcr6$FPFcan>F\xFjs?(F$F^oF^oF$2FjsFh
dnC&>FA-Fgs6%R6$FianF`ilF$FXF$-F[bn6#/-F_bn6#F\[mFQF$F$F$F9F\x>F9-Fgs6%RFehqF$F
XF$-F[bn6#0FihqFQF$F$F$F9F\x@$2Fjs-Fds6#%&TlistGC$>FA-%/desolv/SortDerGFd\n>FK6
$FKFh[q>F\x,&F\xF^oF^oF^oFbgnF$F$F$%:desolv/ConvertAngleTodiffGR6$%&VectrGF\cq6
#%(VecFuncG6#%QCopyright~(c)~1997~by~Khai~Vu~and~Colin~McIntoshGF$C$>FK-%$_PPGF
i^m@%-Fgn6$FPFh_m-FM6%R6$%#vtG%#vfGF$FXF$*&FebnF^o-Ff^m6$FQFgpF^oF$F$F$FPFK*&Fe
bnF^o-Ff^m6$FKFgpF^oF$F$F$FbwR6&F^zF]qF_qF`q6*FdqFeqFfqF[rFiqFjq%&VFlagGF0F^rF$
C)>FK-Fcr6$-Ffr6$<$FPFQFhrFir>F9F[s>FAFdbq>FKFbr>FcuF\u?(FiuF^oF^oFd_mFcuC&>Fft
-Fep6$FiuFQ>F\x-Fcr6$-Ffr6$<$FPFftFhrFir>Fjw-FatF__n@%/-FdsFa]mFjsC$>Fjw-Ffr6$-
Fev6$/FiuFhvFQF\x@$2FjsFh]r>FcuF%>FcuF%@%FcuFKFdtF$F$F$%2desolv/lineartestGR6%F
[eo%&derivGF_bm6#%(var_numGF^pF$@%33-Fdbm6$FQFfo/-F[\lFfemF^o2Fjs-Fcr6$FecoFh_m
C$-F\]m6%Ffo-Fep6$F6FQFK@%2Fjs-Fep6$FKFecoF\uF%F%F$F$F$%-integratepdeGRFcdm6#F[
goF^pF$C%>FK-FOFfem>FK-%.integratepdeDG6&FKFQFfoF\u@%/FK%%FAILGFj`r-Fep6#-FhoFh
]oF$F$F$%&genncGR6%F*Fd[mFc[m61%'linvarG%(ldfeqnsGFe\lF0%(inf_fnsGF,Fb[m%%liscG
%&gliscG%'geniscG%&lleqnG%%lconG%&tlconG%&lcaseG%%caseGF^pF$C.>FgxR6&%%UvarGFd[
mFc[mF,6$F0%'isceqnGF$F$C$>F9-F<6$*&-&%#xiG6#&FfoFio6#-Fep6#FhrF^o-Ff^m6$-FP6#F
cgmF]crF^o/FK;F^o-FdsFdgm,&-Fcr6$F`hmFh_mF^o-&%$etaGF]oF^crFaoF$F$F$@%F5C$>817#
-Fep6#&F?6#FD>Fadr-FM6%R6$F\do%"nGF$FXF$-Fbil6#,(FQF^oFPFaoF^oF^oF$F$F$FadrFgcr
>Fadr7#-F<6$,&FgcrF^oFftFao/Fft;Fjs,&FgcrF^oFaoF^o-%2desolv/checkinputGFccn>%)s
pe_caseGFdt>FK-Fc\m6%FfoFQFjs>FjwFbgm>F\x-F<6$-&F[cr6#&FfoF`bl6#-FepFa]m/FftFfc
r>F\x7$F\x-F<6$-&F]dr6#&FQF`blFgfr/FftFc_m>F_u-FM6'FgxFQFQFfoFK>FeyFhv?&82FadrF
\uC2>Fjt/-&F[cr6#&Ffo6#Fhgr6#Fa\oF^o?(Fft,&FhgrF^oF^oF^oF^oFgcrF\u>Fjt6$Fjt/-Fd
frFahrFjs>Fjt7#Fjt>Fiu-F_x6$FjtF_u?(FftF^oF^oFd_mF\uC$>Fcu-Ff^m6$-FagrFdcrF_hr>
Fiu-Fev6$/Fft/Fcu-%&solveG6$&FiuF`blFcuFiu>Fjt-F_x6$/Fa\oFhfrFjt>F9-%1desolv/gd
efiningG6&FPFiuFQFfo>F9-FM6&FOF9F\xFF>80-FM6%R6$FWF+F$FXF$7$-FO6%F[oFQFFF_oF$F$
F$FjtF\x>Fjt-FM6$RFVF$F$F$FjnF$F$F$Fjt>F9-FM6&%/desolv/reducelGF9FhjrF\x>F9-%1d
esolv/reductionG6$F9F\x>F9-FM6$FhoF9>F9-FM6%F^alF9F\x>F9-%.desolv/eqsortG6#7$-F
ep6#FjtF`\o>Fey6$FeyF97%7#FeyF\xFjwF$6#F[frF$%2desolv/LieBracketGR6%%'Vectr1G%'
Vectr2GF\cq6%%&Part1G%&Part2GFdjqFejqF$C&>FA-FjjqFdcr>FK-F_jqF_cm>F9-F_jqF__r-F
a]o6#,&-F_x6$/FAF9FKF^o-F_x6$/FAFKF9FaoF$F$F$%/desolv/partintGR6&Fe\lF+F,Fc[m6,
FbenFcen%)new_condGF0Fc\l%%nlfnGFg[m%&newfnG%&hipowG%#PxGF^pF$C(>FK-%1desolv/gi
ntorderG6%FjgnFhrFfo>FjwFhv>FAFhv>FcuFhv?&F9FKF\uC&>Fiu-%.desolv/newunkG6$-Ffr6
$&F9Fb_lFfoFQ>Fjw6$Fjw-%/desolv/Der2mapGF_v>F_u-F[\l6$&F9Fe_lFhr@'/F_uFj`r-F][l
6#Fjs/F_uFjsC$>FA6$FA,&F_`sF^o-Ff^m6$FiuFhrFao>Fcu6$Fcu*&Fg`sF^oFiuF^oC&>FA6$FA
,&F_`sF^o-Ff^m6$Fiu-Fcp6$Fhr,&F_uF^oF^oF^oFao>Fcu6$Fcu*&Fg`sF^o-Ff^m6$Fiu-Fcp6$
FhrF_uF^o>Fcu6$Fcu-F<6$*()FaoF\xF^o-Ff^m6$Fg`s-Fcp6$FhrF\xF^o-Ff^m6$Fiu-Fcp6$Fh
r,&F_uF^oF\xFaoF^o/F\x;F^o,&F_uF^oFaoF^o>Fcu6$Fcu*()FaoF_uF^oFiuF^o-Ff^m6$Fg`sF
dbsF^o@%/FcuFhvFjs7$7$-Facl6#-Fcr6$7#FcuFh_mFA7$Fj^mFjwF$F$F$%.desolv/subsOTGRF
bfp6'Fi[mFa\pFeap%%lindG%%lsubGF^pF$C)>FK-Fep6$F^oFQ>FA-Fgs6$R6#%#fnGF$FXF$-F[b
n6#/Fj]m.FedqF$F$F$-Fibp6$FPFis>FA-Fgs6%FdbmFAFK>F\xFhv?&F9FAF\uC%>Fft-Fap6%RF[
doF$FXF$F`doF$F$F$F`]pFeco>Fft-FM6%FgnFftFh[p@$4-Fdbm6$FftF\u>F\x6$F\x/F9Fjs>F\
xFa_n@%/Fi_nFjsFP-Facl6#-F_x6$F\xFPF$F$F$%5desolv/FindExcludeICGR6'F]qFgcn%%Lva
rGF_z%&LevelG6,%&NewICG%)OrderLfnGF0Fc\lF]p%.ListExcludeICG%'NoFuncG%&TLEqnG%&E
qniiG%&EqnjjGF1F$C)@$2F[uFipF\[l>FcuFd_m>F9-FM6%R6$Fian%#lfGF$FXF$-Fep6#-Fgs6%F
dbmFQFPF$F$F$FhrFQ>F97#-F<6$&F96#,(FcuF^oFAFaoF^oF^o/FA;F^oFcu>FjwFhv?(FAF^oF^o
F[uF\u?(FftFhinF^oF[uF\uC%@$-F\]m6$<$&&FPF@Fb_l&&FPF`blFb_lFf`nF[al@$5/Fc^oF^o/
FasF^o@$/-Fep6$F^o&&FejsFe_lFb_l-Fep6$F^o&&FgjsFe_lFb_lF[al?&F\xF9F\uC$>FK-%Ade
solv/ExcludeCrossDifferentiateG6&FejsFgjsF\xFhr@$0FKFjsC%>Fiu-Fev6%/FAFhvFhgmFP
>FK-Fh_l6%FKFiuFhr@$F]\tC%>FK-F^al6$FKFQ>Fjw6$FjwFKF^fl7#FjwF$6#Ff`nF$%.desolv/
conderGR6%Fe\lFd[mFc[m6)F+%(lfn_strGFg[mF0Fiap%(selDderGFeapF^pF$C->FjwRFhanF$F
$F$33-FgnF\fsF`_pFg]mF$F$F$>FKF[fs>Fcu-Fcr6$-Fgs6%%/desolv/isderivGFKFQFcan>F9-
FM6$%/desolv/der2strGFcu>F\x-Fap6%R6$Ffes%$fnsGF$F$F$/FPFQF$F$F$FcuF9>F\x7$-Fep
Fj_n-F<6$/FbirFagrFbgr>FAFfgs>FK-F_x6$F\xFK>FK-FgsF\]t>F\x-FM6$RFhanF$F$F$/FP-F
i^pF]oF$F$F$F_hn7$-F_xF__nF9F$F$F$FfrRF`enF$F^pF$-Fgs6%R6$FianFWF$F$F$@%Fj]t-Fd
bm6$FQFj]m-FdbmFizF$F$F$FQFPF$F$F$%/desolv/isfunofGR6$Fb[m%&lnameGFeesF^pF$@'-F
dbm6$Fj]mFf^mC%>FKFebn?(F$F^oF^oF$-Fdbm6$F^`pFf^m>FKF`[l-F\]m6$F^`pFQF`_p-F\]m6
$Ff^pFQF`cpF$F$F$%.desolv/intpdeGR6&F[eoF+F,Fc[m6-F^eo%&iflagG%%tlfnGFe\lFb[mFe
ap%$lcoG%)int_lderGFg[m%*temp_lderG%(ex_partGF^pF$C)>FKF%>FftFjgn>F_uFhv>FAFeem
>FA-Fgs6%FdbmFQFhr?(F$F^oF^oF$32FjsFd^p4FKC'>F\x-Fep6$F^oFA>Fjw-F\fo6%FftF\xFhr
>F9-F[bm6&FftF\xFjwFhr?(F$F^oF^oF$3F`ct/F9F%C%>Fcu-%'lcoeffG6$FftFjw@%-%4desolv
/intcoefftestG6&-%.desolv/lderivG6%FcuF\xFQFjwFhrFQC)>Fiu-%$intG6$FjwFhr>Fcu-F]
et6$-F_x6$/Fiu%%DUM_GFcuFeet>Fcu-F_x6$/FeetFiuFcu>F_u6$F_uFcu>Fft-Fa]o6#,&FftF^
o-Ff^m6$FcuFhrFao>Fgx-Fhdt6%FftF\xFQ@%/FgxFjw>FKF\u>FjwFgx>FKF\u>F9Fict@)4-Fdbm
6$Fft-Fep6$FjsF\x>FA-Ffr6$FftFA2Fi_nFgcrC$>FA-Fev6$Fe`mFA>FKF%50F9F%FK>FKF\u>FK
F%@%3F`ct0F_uFhvC&>Fey-F]et6$-%4desolv/explicitpartG6%FftFQFhrFhr>F_u6$F_uFey>F
ft,&FftF^o-Ff^m6$FeyFhrFao7$-Fcr6$7#F_uFh_m-Fa]oF`blFjsF$F$F$%Adesolv/Different
ialOrderVariableGR6$F]p%%XVarG6#%%XPosGF^pF$@%-F\]m6%FQFdp.FK-FepFi[rFjsF$F$F$F
\^lR6$F*F_z6'F[]lF0Fc\l%&nleqnG%)lastswapGF^pF$C'>FKF[u@$2FKFipFghl>FftF[hq?(F9
FKFaoFipF\u?(FAF^oF^o,&F9F^oFaoF^oF\u@$2Fjs-Ff[l6%-Fep6#&&FftF@Fe_l-Fep6#&&Fft6
#FhinFe_lFQ>Fft-Fev6%/FAFf[u/FhinFb[uFftFftF$F$F$%0desolv/odecheckGRFcdm6.%&fun
cxGFeap%%lordG%'lowderG%&hiderG%'lowordG%&hiordGF[\m%*order_ODEGF_bmF\_p%'finde
xGF^pF$C$?&FgxFfoF\uC7>FK-Fgs6%FdbmFQFgx@$0FfwF^oF[al>FKF`[l>F9-Fgs6%FdbmF[fsF^
`p>F9-Fcr6$F9Fcan>FA-FM6$RF`\pF$F$F$-Fcr6$FgpFh_mF$F$F$F9>Fjw-%$minG6#Fh[q>Fcu-
FdcoF`^u-F\]m6%FcuFA.Fiu>F\x-Fep6$FiuF9-F\]m6%FjwFAFe^u>FftFg^u>F9-Fev6$F_^rF9>
FA-Fev6$F_^rFA>Fey-Fep6$F6Fft>Fjt-Fep6$FipFft-F\]m6%FgxFeyFe^u>Fjt-Fev6$/FiuFjs
Fjt>Fft-Fev6$/FipFjtFft>Fjw-Fcr6$FjtFh_m@&Fj[p@%F\elF[alC$>F_u,&FcuF^oFjwFao-F]
[l6#7&FKFftFgxF_u4-Fdbm6$-FM6&%0desolv/IsXDerivGF9FftFgxF%C$>F_u,&FejnF^o-Fcr6$
Ff_uFh_mFaoFj`uF%F$F$F$%5desolv/InfinitesimalGR6%%'IndVarG%'DepVarGF\cq6#%'InfF
unGF\\mF$C$>FK-FM6%R6$FianFa[oF$FXF$-FPFi^mF$F$F$F`]mFfo7$Fa\o-Fep6#-FM6%RFdbuF
$FXF$FebuF$F$F$Ff]mFfoF$FjamF$F_bnRF\pF$F^pF$F[^uF$F$F$FghnR6#%,SumListListG6%%
(SumListG%-CollectedSumG%-CommonFactorGF^pF$C&>FKFP>F9Fhv?(F$F^oF^oF$Fd]oC&>FA-
Fgs6%RF`[oF$FXF$-F[bn6#/Ffhl&FQFe_lF$F$F$FKFbfm>FK-Fgs6%RF`[oF$FXF$-F[bn6#0Ffhl
F`duF$F$F$FKFbfm@%Ff`p>FA7$-Fa]o6#-Fcr6$-FM6$RF_ilF$FXF$F\[mF$F$F$FAFh_m&F\hmFe
_l>FAF\hm>F9FhtF`hmF$F$F$%3desolv/collecttermGR6#Fben6&Fg[m%&lxfunGFcen%&tlfacG
F^pF$C&>FKFhv>F9Fhv?&FAFPF\u@$4-F\]m6$F\hmF`hmC&>F96$F9F\hm>Fft-Fgs6%R6$FWFianF
$F$F$@%Fg^qF\uF%F$F$F$FPF\hm>Fft-Fcr6$-FM6$RFVF$F$F$FfhlF$F$F$FftFh_m>FK6$FK7$F
\hmFftFbgnF$F$F$%?desolv/RemoveRedundantUnknown2GR6'F]q%'LInEqnG%&LSolnG%&LUnkn
GF\cq6,%'LConstG%'LFunctG%)LSolnFunG%)LSolnConG%(LRedConGF\an%&TSolnGFgqF[r%'TI
nEqnGF^pF$C.>FK-Fgs6%FgnFhrFc`p>F9-Fgs6%FgnFhrFis>FA-FM6%%7desolv/RemoveConstan
t2GFfoFK>Fft-Fap6%RF`[oF$FXF$/F[o-Facl6#,&F_oF^o-F`oF[_mFaoF$F$F$FfoFA>F\x-%>de
solv/FindRedundantConstant2G6%FftFKFas>Fjw-FM6$RF_ilF$FXF$FcyF$F$F$F\x>F_u-Fbil
6#-F_x6$FjwFP>Fgx-Fbil6#-F_x6$FjwFQ>Fcu-Fbil6#-F_x6$FjwFfo>Fjw-FM6$RF_ilF$FXF$/
FPFhvF$F$F$F\x>Fiu-Fbil6#-F_xF^et7&F_uFgxFcuFiuF$F$F$%3desolv/rank_factorGR6%Fb
enF+F,6'%(DmoduleGFjeu%"iGFe\lFg[mF^pF$C(@$FghsFghl>FK7'%(onetermG%'gensepG%(so
lveDEGF^`r%)odesolveG>F\xFhv>F9FP?(FAF^oF^oFa`oF\u?&FftF9F\uC%@$/-Fg\o6%FftFQFf
oFj`rF[al>F\x6$F\xFft>F9-F_x6$FhgmF97$F\xF`\oF$F$F$FbfqR6%F*%%lofnGF,6.FajtF+%&
nlofnGFg[m%$sysG%%lsysG%%lsolG%&tlinqG%%linqG%&largsG%*starttimeGF0F^pF$C6@$4-F
gn6$FP<$FcanFir-F`io6#%TINPUT~ERROR~!!~first~argument~must~be~a~set~or~listG@$4
-Fgn6$FQFf^v-F`io6#%UINPUT~ERROR~!!~second~argument~must~be~a~set~or~listG@$4-F
gn6$FfoFf^v-F`io6#%TINPUT~ERROR~!!~third~argument~must~be~a~set~or~listG>F_uFF>
Fgx7#-F<6$&F?Fh\s/FjtFC>Fey-%%timeGF$-%6desolv/CheckParameterG6#Fgx>FA-FM6&FOFQ
FQFF>FA-FM6$R6#FhapF$F$F$/FPFPF$F$F$FA>F9FQ>FK-FM6&FOFPF9FF>FK-FM6$RFVF$FXF$@%F
fn-Facl6#-Figo6#FjnFPF$F$F$FK>FKF_hn>Fiu-FM6&FOF_uF9FF>Fjw7#7&FKFiuFAF9>FcuFhv?
(F$F^oF^oF$Fa^rC'@$F]fq-Fg]l6#-F\_l6$%,#~systems~:GFh]r>F\xFdjm>Fjw-Fev6$Fe`mFj
w>Fft-%0desolv/pdesolveG6$F__tFfo@'/FftFj`rF[al/FbgpF^o>Fcu6$FcuF^hm>Fjw7$F`\qF
hfr>Fcu-Fho6#Ffds@$F]fqC$-Fg]l6$%6Modules'statistics~:~G%/desolv_modstatG-Fg]l6
$%6Computational~time~:~G,&F``vF^oFeyFao@%/-FdsFielF^oFhelFcuF$F$F$%2desolv/tra
nsform2GR6%Fe\lF+Fcap6-Fg[mFeapFa\p%$memGF0%&func1G%&func2GFiap%&fvarsGF\_pF^eo
F^pF$C*>Fiu-Fap6%RF[doF$FXF$/-Fcr6$FPFcbpFQF$F$F$FfoFQ>FKFebp>F9Fhbp>F9-Fgs6$RF
hanF$FXF$F_atF$F$F$F9>F9-%/desolv/sortderGFI>Fiu-Fap6%RF[doF$FXF$F\_tF$F$F$FfoF
Q?&FAF9F\uC)>FjwFA?(F$F^oF^oF$-Fdbm6$-Fep6$FjsFjwFf^m>Fjw-Fep6$F^oFjw>F_u7#Fhfr
>FcuFA>FeyF%?&FgxF_uF\u@$-F\]m6%FgxFQ.FftC$>Fcu-F_x6$/FgxFffrFcu>FeyF\u@$FeyC%>
Fcu-Fcr6$FcuFb_p>Fcu-F_x6$FiuFcu>FK-F_x6$/FAFcuFKFKF$F$F$%0desolv/eqsortD2GR6&F
*%%lhilG%&NleqnG%&NlhilG6,%%ssysG%'unssysGF0Fc\l%&nlhilG%(luninfoG%'lsinfoG%&in
fo1G%&info2G%%infoGF^pF$C.>FgxR6#Fe\lF$F$F$@%F\[r7%F[u-FhalF]oFjs7%F^oFbjvF[^uF
$F$F$>FKFP>F9Fhv?(FAFd_mFaoF^oF\u@$/&FQF@FjsC$@$0Fg\oFjs>F96$F9Fg\o>FK-Fev6$Fb\
tFK>F9F`hm>Fjw-FM6$FgxF9>Fcu-FM6$FgxFK>F\x-F_x6$F`[pFQ>FftFjs?&FiuFjwF\uC'>Fft,
&FftF^oF^oF^o@$5-F\]m6$&F9F`blFK-F\]m6$,$Fe\wFaoFKF[al>FAFjs?&F_uFcuF\uC$@&2F`^
n&F_uFb_lC&>Fcu-Fev6$/Fhin-Fep6#7$FiuF_uFcu>FK-Fev6$/Fhin-Fep6#7$Fe\w&FKFg[uFK>
F\x-Fev6$/Fhin-Fep6#7$Fjs&F\xFg[uF\xF^fl/F`^nF^]w@&/F`^nF^o@$2F_]n&F_uF[`lC&>Fc
uFa]w>FKFh]w>F\xF`^wF^fl2Fa^n&F_uFe_lC&>FcuFa]w>FKFh]w>F\xF`^wF^fl>FAFhin@$/FAF
fwC%>Fcu7$FhelFiu>FK7$Fa\oFe\w>F\x7$F__tFjs>FfoFK>FhrF\xF$F$F$FbjrR6&F*FearFd[m
Fc[m6.FbarFcarFe\lF0%%exprG%%lfndGF,Fb[m%&sleqnGFdar%'tsleqnG%'symconGF^pF$C+>F
K-Fc\m6%FhrFfoFjs>Fcu7$F_crFcgm>F9Fhv>Fjw-FM6%R6$Fian%"lGF$FXF$FebuF$F$F$FfoFhr
>F_u-FM6$RFVF$F$F$FgavF$F$F$FP>F_u-FM6&%3desolv/clearderdepGF_uFQFjw>F_u-FM6%%/
desolv/LDsolveGF_uFjw?&FAF_uF\uC'>F\x-%/desolv/prolongG6%-Facl6#-Figo6#,&-F\oF@
F^o-F`oF@FaoFhrFfo>F\x-F_x6$/Fa\oFhelF\x>F\x-F_bw6%F\x7$Fj^m-Fep6#F_uFjw>F\x-F_
]t6%F\xFfoFhr@%/-Fds6#&F\xFe_lFjs>F96$F9&F\xFb_lC%>Fjt-Fbcq6#F__t>Fjt-FM6$RF_il
F$FXF$FfhlF$F$F$Fjt>F96$F9Fg\sF`hmF$F$F$FhdtR6%Fe\lFb[mF+6&Fi[mFf\lFa\pFeapF^pF
$C)@$43-Fgn6$FQFisFcbmFj`s>FKFfbm>Fft-Fgs6%FdbmF[fsFK>F9Fegp>FftFhgp?&FAFftF\u@
$2Fjs-%.desolv/comDerG6%FAF9Ffo>F9FAF9F$F$F$Fc\mR6%Fc[mFd[mF_ep6(F_bm%(deplistG
%&lder1G%&lder2G%&tlderGF0F^pF$C&>F9Fhv>FA-FM6%R6$Fi[mF,F$F$F$FebuF$F$F$FQFP?(F
jwF^oF^oFfoF\uC&>FftFhv?&FKFPF\uC%>F\x-FM6%Ff^mFAFK>F\x-Fgs6%R6$F`ilFgawF$F$F$@
%-F\]mFfemF%F\uF$F$F$F\xFgfo>Fft6$FftF__t>F96$F9Fh[q>FAFgfo7%Fa_mF9Fh[qF$F$F$%=
desolv/FunctionWeightForAutoGR6$%,ListFWeightGFaz6%%-NListFWeightG%(TWeightG%)T
FWeightGF1F$C%>FKFhv?&FAFPF\uC%>F9-Fds6#-Fgs6%FdbmFQF\hm>FA-Fev6$/FipF9FA>FKF]e
qFbgnF$F$F$%/desolv/PVectorGR6$Fe\l%*parameterG6#%&PnameGF^pF$C$@%Fgew>FKFfbm>F
KFQ@'F\[r-Fgs6%FdbmFPFK-Fdbm6$FPFKFPFjsF$F$F$Fg[sR6%Fe\lF*F+6,Fa\pFeapF[gwFg[mF
0Fc\lF\es%%seqnG%$resG%*reduceDerGF^pF$C(>FgxR6$Fa\pF*6$Fa[xFb[xF$F$C$?&FKFQF\u
C$>F9-%/desolv/subsDerG6%FPFKF^o@$0F9Fjs-F][lFIFhvF$F$F$@$/Fd_mFjsFghl>FftFjgn>
F9-Fgs6$RFeesF$FXF$FgesF$F$F$-Fibp6$FftFis?(F$F^oF^oF$F_hqC%>Fcu-FM6%FgxF9FQ@%F
hdv>F9-Fgs6$RFeesF$FXF$FgesF$F$F$-Fibp6$-F`o6#FdamFis>F9-FM6$RFVF$FXF$-Fep6#-Fg
s6$RFeesF$FXF$FgesF$F$F$-Fibp6$F_oFisF$F$F$Fcu>Fft-F_xF[v-Facl6#-FigoF`blF$F$F$
Ff[lR6'%'DerFn1G%'Power1G%'DerFn2G%'Power2G%+ListFunctnG6'Fg[m%'fnord1G%'fnord2
G%'LIndexGF0F1F$C$@%/FPFfo>FK,&FQF^oFhrFaoC$>FK,&F^bnF^o-F_bnFdgmFao@$/FKFjs@%0
Febn-Fep6$F^oFfoC%-F\]m6%FebnFas.F9-F\]m6%Fb`xFas.FA>FK,&F9F^oFAFaoC%Fe`x>Fft-F
ap6%RF`[oF$FXF$F`doF$F$F$FgpFdcm?(F\xFbgpFaoF^oF_`x>FK&FftFj_nFKF$F$F$%*genfact
orGRF`ar6/F+Fg`wF.%&deropG%'lhsderGFifw%)deplist2G%'lopeqnGF0%%operGFe\lFg[mFi[
mF^pF$C4>FA-FM6$RFVF$FXF$FgavF$F$F$FP>Fft-%)adjointDG6$-%(frechetG6%FAFQFfoFfo>
Fft-Fcr6$Fft%)listlistG>FjwF]fr>FK7#-F<6$-&%'lambdaGFjcwFgfr/F_uFc_m>Fiu-FM6%RF
`[oF$FXF$-Fcr6$-Fap6%%/desolv/appoperGFPFQFh_mF$F$F$FftFK>Fcu-FM6$Fe^tFjw>Fiu-F
_x6$/FhfrFhelFiu>F9-FM6%RFfawF$F$F$FebuF$F$F$FQFfo>Fiu-FM6&F_bwFiuFPF9>Fiu-FM6&
F_]tFiuFQFfo>Fjt-FM6$RF_ilF$FXF$-Fep6#-%3desolv/separation2G6%F\[mFFFfhlF$F$F$F
iu>FK-F_x6$FedxFK>Fjt-FM6%F^alFjtFK>Fjt-FM6&FOFjtFKFF>Fjt-Fj[s6$FjtFK>Fjt-FM6$F
hoFjt7%-Fd\sFh\sFKFcuF$F$F$%1desolv/splithnihGRFcdm6)%&lovarGFg]v%'maxvarG%%lmf
nG%(hompartG%*inhompartG%"sGF^pF$C+>FK-FfrF_cm>F9Feem>FAFfw>Fft-FM6%R6$FianF\er
F$FXF$@%3Fj]t/F[uFQFj]mFhvF$F$F$F9FA@$/FbgpFjs-F][l6#7$FjsFP>F\xFhv>FjwFhv@'F\[
r?&FcuFPF\u@%-FdbmF[v>F\x6$F\xFcu>Fjw6$FjwFcu-Fdbm6$FPFft>F\xFcu>FjwFcu7$-Fcr6$
Fa_nFh_m-Fcr6$F]]tFh_mF$F$F$%4desolv/SolveOneTermGRFcdm6(%)SolutionG%.SolveFunc
tionG%-NewFunctionsG%&PolynG%3IntVariableNOrdersG%*VarNOrderGF^pF$C$@$-Fgn6$Fjg
nFh_m-F][l6#Fj`r@%F]bn7$FcyFQC(>F\x-Fap6%RF[erF$FXF$7$FPFQF$F$F$FdpFgp>FAFhv>FK
Fjs?&FjwF\xF\u@$2Fjs&FjwFe_lC%>Fft-%2desolv/PolynomialG6%Fhfr-%9desolv/Function
VariablesGF]oFQ>FA6$FA-Fep6#Fjfm>FK,&FKF^oF^hmF^o>F9-Fev6$/Fip7#-Fcp6$Fjs-Fds6#
Fj[yFP7$/F9FK7$Fj^mFAF$F$F$%7desolv/WeighVarsInListGR6$F\cqFaz6'F[hl%0WeighOfVa
riableG%*TVariableG%*TotalVarsGFfitF1F$C'>FK-FM6$RFhanF$FXF$-Fep6#-FjoF]oF$F$F$
-Fcr6$FQFcan>FftFfw>F9Fhv?&F\xFPF\uC$>FA-F_x6$FaxFK>F96$F9-Fbil6#,&FftF^oFd^pFa
oF`hmF$F$F$%'comtabGR6$%%LVecGF\cq63F0Fc\l%*NumVectorG%(CommtorGF[r%'GenEqnG%*L
CoeffEqnG%*ConstSolnG%*StructEqnG%)LVecDiffG%(LLambdaGFdcx%$psiG%'LGammaG%&gamm
aG%+GammaCountG%*ComMatrixGF^pF$C1>FAF[u>Fgx-FM6$%:desolv/ConvertAngleToDiffGFP
@$3-F`]l6#&Fc]l6#F[_y2FjsFg`yC$-Fg]l6#%2The~vectors~are~:G?(FKF^oF^oFAF\u-Fg]l6
#/&FhjrFio&FgxFio>Fey7#-F<6$&FjtFio/FK;F^oFA>847#-F<6$7#-Fcp6$FjsFAFiay>F\by-Fc
r6$F\by%'matrixG>FadrFhv>83Fjs>Fgx-FM6%F_jqFPFQ>Fjw-Fcr6$-Fap6%RF[doF$FXF$*&FPF
^oFQF^oF$F$F$FeyFgxFh_m>F_u-Fcr6$7#-F<6$*&FhayF^oFbayF^oFiayFh_m?(FKF^oF^o,&FAF
^oFaoF^oF\u?(F9Fc^lF^oFAF\uC$>Fft-F^]s6%FjimF_jmFQ@$0FftFjsC'>F\x-Fa]o6#-Figo6#
,&FjwF^oFftFao>Fcu-Fcr6$-Fep6$F^o-%.desolv/sepsumG6$F\xFQFir@$/FcuF%F[al>Fiu-Fi
ir6$Fcu-Fcr6$FeyFir@%F_^rC&>Fiby,&FibyF^oF^oF^o>Fadr6$Fadr/&Fhgr6#FibyFft>&F\by
F_gnFafy>&F\byFbt,$FafyFaoC$>FdfyF^dp>Fffy,$FdfyFao>Fadr7#Fadr@$Fd`yC%-Fg]l6#%9
The~Commutator~table~is:G-Fg]l6#F\by@$2FjsFibyC$-Fg]l6#%&whereG-Fg]l6#Fadr@%Ffg
y7$-%%copyGFdgyFadrF_hyF$F$F$FacvR6'F*%&lineqGF\^vF+F,6A%(lmoduleG%*flag_mod7GF
\\v%,startmoduleGFajtF[_sF`bt%&nlsolG%'nlineqG%'tlineqGF`ivFg[mFe\l%%eqnsGF[goF
cfpFbenF[^v%'nleqn2GFi[m%%sys1G%%sys2GF0Fc\l%'eqn_noG%%llenG%'maxlenG%'minlenG%
)totallenG%+flag_delayG%)NumUnknsGF^pF$C<>FK7)%)onetermDG%/directseparateGFd\v%
1indirectseparateGFf`r%*odesolveDG%*decoupleDG>FjwFhr>FiuFfo>F_uFQ>F\xFP?(8:F^o
F^oFi_nF\uC$>F\by-Fgcq6$&F\x6#FcjyFjw@%/-FdsFdgyF^o>F\x-Fev6$/Fcjy-FepFdgyF\xC(
>F\x-Fev6$/FcjyFhvF\x>85Fhv?(8;F^oF^oF\[zF\uC$>Fgx7$-Fep6#-Fev6$/Fj[zFhvF\byFic
w>Fh[z6$Fh[z7&7$&F\by6#Fj[zF__tFgxFiuFjw>Fh[z7#Fh[z@$F]fqC$-Fg]l6#%8The~splitti
ng~systems~:G@$2F6F`fq-Fg]l6#-Fho6#Fh[z-F][lFe]z>Fey7#-Fcp6$FaoFi_n@$%0desolv_s
implifyG>F\x-FM6$FcdlF\x>F\x-FhhnFj_n>F9Fjs>FftF^o?(F$F^oF^oF$2Fjs-Fds6#FeyC)?(
FcjyF^oF^oFf^zF\u@$2&FeyFijyF^oC$>F\byFfjy@%F[[z>F\xF^[zC(>F\xFd[z>Fh[zFhv?(Fj[
zF^oF^oF\[zF\uC$>FgxF]\z>Fh[zFd\z>Fh[zFj\z@$F]fqC$F]]z@$Fa]zFb]zFf]z-F]iv6&F\xF
ey.F\x.Fey-%@desolv/SaveInterMediateSolutionG6'F\xF_uFiuFjwFas>8<F^o>FjtFj`r?(F
cjyF^oF^oFf^zF\uC$@%54-Fgn6$FhjyFh_m2-Fds6#Fhjy%1DesolvScanFactorGC%>Fd`zFcjy>F
hjrFhjy?(FAF^oF^o%-DesolvModuleG2FA""(@$2F\_zFAC%@$/FAFa`o@$/-Fds6#-Ffr6$FhjrFj
wF^oC$>Fjt-F\jy6%FhjrFjwFas@$0FjtFj`rC$>FA""'F^fl>Fjt-Fg\oFfbz@%/FjtFj`r>Fey-Fe
v6$/FcjyFAFeyF^fl?(FAF^oF^oFeazFfazC$?(Fj[zF^oF^oFf^zF\u@$2&FeyFh\zFAC'>Fd`zFj[
z>Fhjr&F\xFh\z@$F\bz@$F^bzC$>FjtFebz@$FhbzC$>FAF[czF^fl>FjtF]cz@%F_cz>Fey-Fev6$
/Fj[zFAFeyF^fl@$FhbzF^fl@$FhbzF^fl@'FhbzC%@.F^ioC->Fbdv-Fev6$/F^o,&&FbdvFb_lF^o
F^oF^oFbdv>8AF%>87-Fep6$F^o-F\o6#&FjtFb_l>F\x-Fev6$/Fd`zFhvF\x>Fey-Fev6$FcfzFey
>F_u-FM6%F\yF_uF_fz@$-F\]m6$FjsF_uFejx>F_u-Fgs6$RFVF$FXF$4-Fgn6$FPFi]oF$F$F$F_u
>Fiu-FM6%R6$FW%#seGF$FXF$/F[o-F\y6$F_oFQF$F$F$FiuF_fz?(Fj[zFf^zFaoF^oF\u@$-Fdbm
6$F]dzFjezC$>Fhgr-F\y6$F]dzF_fz@'/FhgrFjsC$>F\x-Fev6$Fb\zF\x>Fey-Fev6$Fb\zFey2*
&-Fhal6#F]dzF^o%3DesolvExpandFactorGF^o-FhalF`hr>FhezF\uC$>F\x-Fev6$/Fj[zFhgrF\
x>Fey-Fev6$/Fj[zFjsFey@%FhezC%>F\x7$F__t,&F]fzF^o-F`oF^fzFao>Fey7$-FepFg^zFajl>
Fjw&FjtFe_l>Fjw-F_x6$/-Fc`s6#F]fzFhvFhjz/FAFipC&>Fbdv-Fev6$/Fip,&&FbdvFe_lF^oF^
oF^oFbdv@$3-F`]l6#F\^zF\^z>Fjt-Fev6$/F^o-FM6$FcdlF_fzFjt>F\x7$-FepF^fz-Fep6#Faf
z>Fey7$-Fcp6$Fjs-FdsF^fz-Fep6#Fefz/FAF6C%>Fbdv-Fev6$/F6,&&FbdvF[`lF^oF^oF^oFbdv
>FhezF%@%2F^o-FdsFh\sC->86-Fev6$/Fd`z-Fep6#FhjzF\x>897&F[^[lF_uFiuFjw>F_uFhfz@$
F[gz-F][l6#7#Fb^[l>F_u7$Ficw&FhjzFb_l>F\xFafz>F\x-FM6%F\yF\xF_fz>Fiu-FM6%RFhgzF
$FXF$FjgzF$F$F$FiuF_fz>Fjw-F_x6$F\[[lFjw>887&F\xF_uFiuFjw-F][l6#7$Fh_[lFb^[lC+>
FjezF[fz>F\xFafz>FeyFefz>F_uFhfz@$F[gzFejx>F_u-Fgs6$RFVF$FXF$FagzF$F$F$F_u>Fiu-
FM6%RFhgzF$FXF$FjgzF$F$F$FiuF_fz?(Fj[zFf^zFaoF^oF\u@$F_hzC$>FhgrFchz@'FfhzC$>F\
xFihz>FeyF\izF^iz>FhezF\uC$>F\xFgiz>FeyF[jz@%FhezC$>F\x7$F__t-Facl6#-Figo6#Fbjz
>FeyFejz>FjwFe_[l/FAFDC'>Fbdv-Fev6$/FD,&&FbdvFfdrF^oF^oF^oFbdv@%Fh[[l>FadrF_\[l
>FadrF_fz>F\x7$-FepF\hyFd\[l>Fey7$-Fcp6$Fjs-FdsF\hyF[][l@&/&FjtF[`lFip>FjwFhjz/
Ffc[lF6C$>F\x7$F__tFhjr>Fey7$Ffjz"#6F\bzC'>Fbdv-Fev6$/Fa`o,&&Fbdv6#Fa`oF^oF^oF^
oFbdv@%Fh[[l>FadrF_\[l>FadrF_fz>F\xF]c[l>FeyF`c[l>FjwFhjz/FAF[czC$>Fbdv-Fev6$/F
[cz,&&Fbdv6#F[czF^oF^oF^oFbdv@%2F6Fh][lC)>F[^[l-Fev6$/Fd`z-Fep6#Ffc[lF\x>Fb^[lF
c^[l>F\x-Fev6$/Fd`zF_fzF\x>F_u7$Ficw&Ffc[lFb_l>FjwFhjz>Fh_[lFi_[lFj_[lC&@%Fh[[l
>Fhgr-FcdlF^fz>FhgrF_fz>F\x7$FhgrFd\[l>Fey7$FjsF[][l>FjwFhjz-%9desolv/PrintModu
leOutputG6)FAFd`z7$FhjrF_fzF\xFiuFjwFey@$F[gzFejx-Fdbm6$FeyFajlC&-F\]m6%FajlFey
.Fd`z>Fhjr&F\x6#Fd`z>Fjt-Fd\vFfbz@$FhbzC'-Fg]l6#%6Do~delay~substitutionG>F\xFaf
z>FeyFefz>FjezF[fz?(Fj[zFf^zFaoF^oF\u@$F_hzC$>FhgrFchz@%FfhzC$>F\xFihz>FeyF\izC
$>F\xFgiz>FeyF[jzC->Fcu-Ffr6$7%F\xF_uFiuFjw>Fiu-%>desolv/RemoveRedundantUnknown
G6'F\xF_uFiuFcuFas>F\xF`^n>F_uFa^n>FiuF_]n>Fjw-Fjp6'F\xF_uFiuFcuF%>F\xFdjm>F_uF
d[y>Fjw&FjwF[`l>FjwFhi[l@%/FidvFh]rC-@$52Ff^zFip2FeazFgazF^fl>Fbdv-Fev6$/Fgaz,&
&Fbdv6#FgazF^oF^oF^oFbdv@$3F^fq2F^oF`fq-Fg]l6#%CModule~(~7~)~applies~to~the~sys
temG@%2F9Fip>%.desolv_torderGF^o>Fa\\lFip>Fjt-Fdo6%F\xFjwFas>F\xFjt>FeyFh]z-%:d
esolv/PrintModule7OutputG6%F\xFeyFjw>FftF^o>F9,&F9F^oF^oF^o@$1%/DesolvIcdeCallG
F9F^fl>Fey7#-Fcp6$FjsFi_n>Fiu-F]j[lFb`z>F\xF`^n>F_uFa^n>FiuF_]n>Fjw-Fjp6'F\xF_u
FiuFjwF\u>F\xFdjm>F_uFd[y>FjwFhj[l>F\x-FhoFj_n>F_u-FhoFjcw>Fiu-FhoF_v>FjwFhi[l>
Fjw-%/desolv/sortfunGFa]m7#Fi_[lF$6*FeazFbizFbdvF`azFa\\lF\^zF`]\l%+DesolvSaveG
F$%'genvecGR6%F\^v%%lparGF,6'%)argumentG%$parG%*GeneratorG%(LievecsG%'symvecGFe
jqF$C*@$2F7F6-F`io6#%2Arguments~missingG?&FKF?F\u@$45-FgnF`hn-Fgn6$FKFir-F`io6#
%@Arguments~must~be~lists~or~setsG?&FK&F?Fb_lF\u@$4F_dn-F`io6#%QElements~of~the
~first~argument~must~be~equationsG>FA-Fap6%RF`[oF$FXF$*&F_oF^o-%%DiffG6$FjjqFQF
^oF$F$F$FPFfo>FA-Facl6#Ff_m>FftFhv?&F9FQF\uC(>F\x-F_jwFeho@$F]yF[al>F\x-%:desol
v/ConvertDiffToAngleGFj_n@$-Fgn6$F9Fc`p>F\x-Fbil6#-F_x6$/F9F^oF\x>Fft6$FftF\x@$
3-F`]l6#&Fc]l6#F]_\l2FjsFac\lC$-Fg]l6$FbjwF9-Fg]l6#-Fb`yFj_nFgfoF$F$F$%/desolv/
backintGR6&Fe\l%&lvlcoGF+F,65F^eoF[goF[_sFabtF]eo%#coGFg]v%%fnewG%%ieqnG%&fnew1
GFhap%#tfGF\_pFg[m%%ilcoG%%var1G%$tcoG%(reverseGF`btF^pF$C*>Fh[zR6#Fgaw6#F]\vF$
F$7#-F<6$&FP6#,(F[uF^oFKFaoF^oF^o/FK;F^oF[uF$F$F$>FKF\u>F9FP>FAFhv>F[^[lFfo?&Ff
tFQFKC%>F\xF^hm>FftFjfm?&FjwFftF\uC%>Fcu-Ffr6$F9F[^[l>F_u-%,integrateDEG6%F9Fcu
Fa_n@%0F_uFj`rC)>F[^[lFb_w>F_u&F^]wFb_l>Fiu-Fgs6%RF]isF$F$F$4FbhwF$F$F$F[^[lFcu
>FKF\u>Fadr-Facl6#*&FjwF^oF_uF^o?&FeyFiuF\uC&>FjtFey>Fhjr-Fat6$<#F_cr<#Ffjz?&Fh
gr-Fh[zF[_mF\uC$>Fiby&FhgrFb_l@$-Fdbm6$FeyFibyC$>Fhgr&FhgrFe_l?&F\byFhgrF\u@$2F
js-Fds6#-Ffr6$F\byFhjr>Fjt-Ff^m6$*&FjtF^oF\byFaoFiby@$0FeyFjt>Fadr-F_x6$/FeyFjt
Fadr>FA6$FAFcfl>FKF%>FAFiho@%FK7$FadrFA7$FPFAF$F$F$%0desolv/islinearGR6$Fe\lFb[
m6$F\go%&nopsfGF^pF$C$>F9R6$F\goFfesF$F$F$-Fds6#-Fgs6%FdbmF[fsFQF$F$F$@'F\[rC$?
&FKFPF\u@$2F^o-F96$FKFfbm-F][l6#F%F\u2F^o-F9FebmF%F\uF$F$F$%1desolv/depallvarGR
6%F[eoF_bmF,6$%(set_facGF_eoF^pF$C%@%Feen>FK<#Fa_m>FKF_bq>F9-Ffr6$FabqFfo@%/F9F
foF\uF%F$F$F$Fj[sR6$F*F+6/%)sim_leqnG%'l1termGFghm%)zero_eleGFe\lF0Fi[m%&oldfnG
F\_sFd[mFcenFben%%lfnDGF^pF$C(>Fhjr-FM6$FexFQ>FKFP>F9FF>FAFjs?(F$F^oF^oF$0FAFfw
C'>FAFfw>FjwFjs?&F\xFKF\uC$>Fjw,&FjwF^oF^oF^o@$4FjinC%@$-Fgn6$F\xFgen>F\x-F^en6
$F\xFhjr@%5Fd^]l-Fgn6$F\xFefnC$>FK-Fev6$/FjwFhvFK>F9Fd]vC%>FK-FM6%FhdsFKF\x>F9-
FM6%FhdsF9F\x>F9Fd]vF^fl>FK-Fgs6$RFVF$FXF$@%FcyF%F\uF$F$F$FK>F9-Fgs6$RFVF$FXF$@
%FcyF%F\uF$F$F$F9F_\oF$F$F$%0desolv/Der2DiffGR6#Ff`w6$FeapF\esF^pF$C%>FK-Fgs6$R
FhanF$FXF$/Fj]mFedqF$F$F$F[fs>F9-FM6$RFhan6$Fi[mFj[mF$F$C%>FKFebn>F9F`p@%Fj[p/F
PFK/FP-Fda\l6$FKF`\oF$F$F$FKF[cnF$F$F$%0desolv/reducecnGRF_[x6.FeapFg[mFcfpFa\p
F`bt%&lLderGF0%&ltrigG%(indetfnG%+reducediffGF]atF\esF^pF$C,>FgxR6%Fa\pF*F]at6$
Fa[xFf\lF$F$C$?&FKFQF\uC$>F9-FexFh]o@&0FebnFgfnF[al-F\]m6$-F`jo6$FgpF`]p<$F\uFc
[p-F][l6#/-Fc`sF]o-%.desolv/mindifG6%F9&FKFe_lFgpFhvF$F$F$>F_uR6$%%expxG%$lneG6
#F+F$F$C$@%Ffn>FK-Fgs6%Fa^tFc^xFQ>FK-Fgs6%Fa^tF[fsFQ-Fep6#FcxF$F$F$@$Fd\xFghl>F
ey-FM6$RFhanF$FXF$@%Fj]tFj]mFPF$F$F$Ffo>F9Fjgn>Fjw-FM6$RFVF$FXF$-FexFjhqF$F$F$F
Q>FK<#-F_u6$F9Fey?(F$F^oF^oF$Fd]oC%>Fjt-FM6&FgxFKFQFey@%/Fh][lF^o>FK<#-F_u6$-Fe
p6$F^oFjtFey>FK-FM6%F_uFjtFey>F9-F_x6$FjtF9>Fiu7:%$sinG%$cosG%$tanG%$secG%$cscG
%$cotG%%sinhG%%coshG%%tanhG%%sechG%%cschG%%cothG%'arcsinG%'arccosG%'arctanG%'ar
csecG%'arccscG%'arccotG%(arcsinhG%(arccoshG%(arctanhG%(arcsechG%(arccschG%(arcc
othG@%-Fdbm6$F9Fiu-Fa]o6$-Facl6#-FigoFI%%trigGF^i]lF$F$F$F_goR6#F[eo6%Fg[mFcenF
[goF^pF$@$FagzC%>FA-FfgoF]o@%Fd[qC$>FKFhv?&F9FAF\u@$4-Fgn6$F9Fi]o>FKF_gn>FKFA@%
-F`]l6#%0desolv_nzassumeG>Ffj]l-F\w6$Ffj]lFj`q>Ffj]lFj`qF$Fej]lF$%8desolv/Print
FinalOutputGR6&F]qFhguF_q%%LFunGF$FejqF$@$F]fqC$-Fg]l6#%6Complete~one~solutionG
@$2FipF`fqC(@%/F[uFjs-Fg]l6#%@The~system~is~completely~solvedGC$-Fg]l6#%7Unsolv
ed~equation(s)~:G-Fg]l6#-FhoF]o@$2FjsFd_mC$-Fg]l6#%1Inequation(s)~:~G-Fg]l6#-Fh
oF[_m-Fg]l6#%.Solution(s)~:G-Fg]l6#-FhoFdgm-Fg]l6#%GUnknown~function(s)~and~con
stant(s)~:~G-Fg]lF`crF$F$F$FgexRFcdm6%%)new_eqnsG%)fn_occurG%&lsvarGF^pF$C&>FKF
%>F9Feem>FA-Fgs6%R6$F\doF^isF$F$F$4Fi`tF$F$F$FfoF9@%3F_ct-Fdbm6$FPFA-Fep6$F^o-F
aey6%FjgnFAFfo7#FPF$F$F$%;desolv/AccumulateVarWeightGR6%%,ListVWeightG%(VWeight
G%*MFunctionG6%F0%)TVWeightG%+NewVWeightGF1F$C%>FAFP?(FKF^oF^oFd_mF\uC$>F9F\_m@
%/FfoF_`sC$>F9-Fev6$/F6Fg`s&FAFio>FA-Fev6$/FKF9FAC%>F9-Fap6%RF`[oF$FXF$,&FPF^oF
QF^oF$F$F$Fg`s&Fb`^lFfdr>F9-Fev6$/FDF9Fb`^l>FAFd`^lFAF$F$F$%-desolv/felimGR6&Fe
\lF+F,%%elfnG6-Fd[m%&lfnv1G%&lfnv2G%%lsepGFg[m%%ldcoGFabtF[go%&oelfnGF]gxF_d\lF
^pF$C(>FgxFegx>F_uFhr>F\xFP>FjwFhv?(F$F^oF^oF$2Fjs-FdsFjcwC/>FK-Fep6$F^o-%0deso
lv/stardep2G6$F_uFgx>F9-Fgs6%R6$Fian%"vGF$FXF$4-FdbmFfemF$F$F$FQFK>FA-Fgs6%RF`c
^lF$FXF$Fbc^lF$F$F$F_uFK>FA-FM6$RFhanF$FXF$Fj]mF$F$F$FA@$2Fjs-Fds6#-Ffr6$-FbgoF
j_nFAFejx>Fft-Fadm6%F\xF9Fbgn>Fcu-FM6%RFjfuF$FXF$@%-Fdbm6$FfhlFQF\[mFhvF$F$F$Ff
tFA>Fcu-%/desolv/sortlenGFiel?(F$F^oF^oF$2FjsFidvC(>FeyFdam>Fcu-Fev6$Fe`mFcu@$/
FeyFjsF[al>Fjw6$FeyFjw>Fcu-FM6&R6%FW%$co1G%#v1GF$FXF$-Fa]o6#-Ff^m6$*&FPF^oFQFao
FfoF$F$F$FcuFeyFK>Fft-FM6&RF_f^l6#%$texGF$F$C$>FK-Fa]o6#-Ff^m6$*&F\[mF^oFQFaoFf
o@%F_`xFhv-Fev6$/F^oFKFPF$F$F$FftFeyFK>Fft-FM6%FcrFftFgen>F\x-Fcr6$FftFh_m>Fgx-
Ffr6$F\xFgx>F_u-Ffr6$F\xF_u7$F\xF]]tF$F$F$F^alRF`en6'F[goF\goFbenFcenFg[mF^pF$C
%F^go>FKFdgo@'F]hoC'>FAFhv?&FftFKF\uC%@$/-Fds6#-FfrF_alFjsF[al@$-Fgn6$FftFefn>F
ftFegp>FA6$FAFft>F\xFgho-F_go6#*&FKF^oF\xFao@%/F\xF^oF_ioF\xFcio@%FfioC$F\joF`[
lF_ioFKF$F$F$FexRF\av6*Fa\pF[esF_coFj[mF,%&orderGF0%&D2DerGF^pF$C$>FiuRF\av6&F^
coF_co%%fordGF0F$F$@'-Fgn6$Fi]pFi\mC&>FKFf]p>FAFfw>F97#-F<6$,&FAF^o-Fds6#FfgmFa
o/FftFde\l-Fedq6%Ff^pF9F_`mFbep-Fedq6%Ff^p7#FeepF_`m-Fedq6%Ff^p7#F^oF_`mF$F$F$@
)4Fj]t-Fedq6%FPFFFFF`_p-FiuF]oF_atC&>FftFhv>FKFP?(F$F^oF^oF$FdatC$>Fft6$FftF^jo
>FKF`[l@%-Fdbm6$F^`pFb_p-Ff^m6$-FiuFioFftC'>FftFgfo>FjwFbgp>F\x7#Fa\o>F97#-F<6$
,&FjwF^o-Fds6#-F_x6$/&F\xFielFhvFftFao/FcuFcjl-Fedq6%F^`pF9F\x-Fedq6%Fj]m7#-Fcp
6$FjsF[uF_`mF$F$F$%.desolv/simpscGRFci]lFdi]lF^pF$C%>FKFhv@$FagzC$>FA-Ffgo6#-Fa
cl6#*&FPF^oFjgoFao@'Fd[q?&F9FAF\u@$F^j]l>FKF_gn3Fgdq-Fgn6$Fh`pFi]o>FKFcct>FKFAF
KF$F$F$%7desolv/DimensionVectorGR6%F_q%%LUnkGF\cq6(%)TUnknownG%&TLUnkGFd_\lFe_\
l%'SymVecG%*IndFactorG6#%DCopyright~(c)~March~1999~by~Khai~VuGF$C(>FA-Fap6%RF`[
oF$FXF$Fba\lF$F$F$FPFfo>FA-Facl6#-Fa]oFha\l>FftFhv>F9-Fi^\lF[_m?&FKF9F\uC'>F\x-
F_jw6$FAFK@$F]yF[al>F\xF`b\l>F\x7$FKF\x>FftF\c\lFgfoF$F$F$%=desolv/RemoveVanish
Unknowns2GR6&F]qF^qF_qF`q6)Fcq%$LT1G%$LT2GFgqF\rF]rFhqFi`_lF$C'>F9Fbr>FAF[s>FKF
gt>FK-Fgs6$RFhanF$FXF$Fa\_lF$F$F$FK@%FewFgwC+>FcuFhr?&FftFKF\u>Fcu-F_x6$FhgmFcu
>FKFcx>F\xFP?&FftFKF\u>F\x-FM6%F\yF\xFd`l>F\x-FM6$RFVF$FXF$@%FcyFhvFPF$F$F$F\x>
FjwFQ?&FftFKF\u>Fjw-FM6%F\yFjwFd`l7%F\xFjwFcuF$F$F$%>desolv/FunctionWeightInSys
temGR6$F[anF]q6%FczF_iw%(LWeightGF^pF$C%>FAFhv?&FKFPF\uC$>F9-FM6%RF`[oF$FXF$@%F
cc^lF^oFjsF$F$F$FQFK>FA6$FA-Fcr6$F9Fh_mFihoF$F$F$Ff`rRFjdo63%(lfn_occGF_eoFc[m%
%lintG%(num_intGF0FajtF^eo%)int_expnG%'varnumG%&expn1GF*%)int_consGF[_s%(genfla
gGF]gxF_gxF^pF$C1>FKFeem>FibyFegx>F\by-FM6%RF\hxF$FXF$@%F^hxFPFhvF$F$F$FK-FdsFb
fy@$/F\[zFjsFejx@%F5>FhgrFedr>FhgrF\u@%Fhgr>F9Ffo>F9-Fhdo6%FPFKFfo>F_uFjgn>FiuF
\u>FftFhv>FcuFhv>F9-Fgs6%R6$Fac^lFWF$F$F$Fi`tF$F$F$F9FP?&FAF9F\u@%-FdbmF]eqC&>F
\xFjs>FeyF_u?(F$F^oF^oF$0FeyFjsC$>Fey-F[bt6&F_uFKFibyFA@$Ffh_l@'/&FeyFe_lFjsC$>
F_u&FeyFb_l>F\xF^jq4Fhgr>FeyFjsC%>Fjt-Ff^s6&F^i_lFKFibyFA@$0FjtFjsC'>F_u,&Fai_l
F^o&F_fzFb_lF^o>F\xF^jq>FKFhjz>Fcu6$Fcu-Fep6#-Fev6$Fe`mF_fz>FiuF%>FeyFjs>Fft6$F
ft7$FAF\xC$>Fey-Facl6#-F]et6$F_uFA@$4-Fdbm6$FeyF]etC$>F_uFey>Fft6$Fft7$FAF^o>Ff
tFgfo>F\x-Fcr6$-FM6$RFVF$FXF$FfhlF$F$F$FftFh_m@'F]yFj`r-Fdbm6$F_uF]etC$@$3-F`]l
6#&Fc]l6#%*integrateG2FjsFj\`l-Fg]l6$%2cannot~integrate:GF_uFj`rC&>FadrFhv?&Fgx
FftF\u@$2Fjs&FgxFe_lC%>Fhjr-%&polyDG6%-FepFd`vFibyFQ>Fadr6$Fadr-Fep6#&FhjrFe_l>
F_u,&F_uF^o&FhjrFb_lF^o>FK-Fcr6$<$Fa\oFj^mFcan7$7$F_uFcu7$Fa\oFadrF$F$F$Fb\vRFc
dmF``rF^pF$C%>FKFc`r>FK-FiiyFeo@%Fi`rFj`rF]arF$F$F$%2desolv/InitialiseGRFi\]l6#
Fb[mF1F$C'>Ff\nF[u>FhblFF>Ff`nFdt>%0DesolvMinimalICG-%&tableGF$?&FKFQF\u@%F[`p>
&Fj_`l6#F^`pFF>&Fj_`lFioFFF$6&Ff\nFhblFf`nFj_`lF$F^\xR6%F_evFcfp%'statusG6#F^eo
F^pF$C%@$0Febn-Fep6$F^o&FQFb_lFj`s>FK-F`jo6$Fgp-Fep6$FipF^a`l@)/FKF\u/FP-Fic]l6
$Fj^mFgp/FKFc[p/F^a`lF`du/FfoF^oFjsFPF$F$F$%3desolv/splitimpexpGR6%Fe\lF+F_bm6%
%%lxfnG%(imppartG%(exppartGF^pF$C'>FK-Fgs6%FdbmFQFfo>FK-FM6$RFhanF$FXF$@%Fj]tFj
]mFPF$F$F$FK>F9-Fgs6%FdbmFjgnFK>FA-Facl6#,&FPF^oF9FaoFedoF$F$F$Fb`yR6#FbjqF$Fej
qF$@%F\[r-FM6$R6#Fb[rF$FXF$*&FebnF^o-Fda\l6$FjjqFgpF^oF$F$F$FPFjc`lF$F$F$%)diff
/DerGRF$6(FhapFj[mF]eoFfjo%"pGF0F^pF$C'>FK-Fedq6#&F?6#;F^oF6>F97#-F<6$&F?Fa]m/F
jwFC>FftF^jo?&FAF9F\u@%-F\]m6%FA-Fep6$F6FKF\x>Fft-Fev6$/F\x,&FdaxF^oF^oF^oFftFj
`s-Fev6$/FipFftFKF$F$F$%.desolv/approlGR6$%&prolnGFe\lFggnF^pF$C%@%-Fgn6$FQFin>
FK,&F]gpF^oFfiuFao>FKFQ>FK-FM6%R6$F\bxF[eoF$F$F$*&F\[mF^o-Ff^m6$FQFfhlF^oF$F$F$
FPFK-Facl6#-FcrF]hnF$F$F$FaeyR6$%(sep_eqnG%(sep_varG6&%(indlistGFifw%(temlistGF
\goF^pF$@%F\[rC&>FKFF>F9FF?&FftFPF\uC%>FA-%0desolv/colterm2G6%FKF9-%1desolv/sep
factorGF_al>FKF\hm>F9F[gmFf[nC$>FA-Feh`lFfem7$7#F\hm7#F[gmF$F$F$FcblR6#F,6#%(ar
bnameGF^pF$C&@$4-F`]l6#%1desolv_arbitraryG>Fgi`lFjs>Fgi`l,&Fgi`lF^oF^oF^o>FK-F\
_l6$%$AAAGFgi`l@%Fh[^lFK-FKFjepF$Ffi`lF$FbauR6%FcjoFdjoF_bm6#%%dordGF^pF$C$>FK,
&F[^uF^oFc_rFao@%32FjsFK/FP-Ff^m6$FQ-Fcp6$FfoFKF\uF%F$F$F$%6desolv/ModdingProce
ssGR6%Fe\lF*%%ljetG6-%)PostLDerG%(PreLDerGFg[m%+ReduceFlagG%&CountGFbhsF[r%&TLD
erG%&TLJetG%(LAllJetGF\anF\\mF$C0@$FcyFj`s@$Fd\xFghl>FAF^[q>FftF\u>F9Ffo>F\xFjs
>Fiu-Fcr6$-FbanF[_mFcan>F_u-FM6$%7desolv/DerToJetBundle2GFiu>Fey-Fap6%RF`[oF$FX
F$F\_tF$F$F$FiuF_u>Fjw-Fbil6#-F_x6$FeyFQ>Fgx<$FcgmFicw?(F$F^oF^oF$3Fft2F\xFajlC
'>FftF%?&FcuFjwF\u>FA-Fbil6#-F_x6$FcuFA>FA-Fa]oF@>FK-Ffr6$FAFgx@$0F9FKC%>F9FK>F
\xF^jq>FftF\u@$1FajlF\x-F`io6#%NFinite~loop~occurs~during~the~Modding~processG7
$FAFKF$F$F$%)classifyGRF$6%%%lunkGF]atF0F^pF$C(@%/F7Fjs>FKFF>FK7#-F<6$F>/FA;F^o
F7@$4-F`]l6#%0desolv_liassumeG>Fd`alFdt@$4Fdj]l>Ffj]lFdt>F9-FM6$RFVF$FXF$@%Fj]t
Fj]mFPF$F$F$FK>Ffj]l-FM6$F[__lFfj]l@%F_hqC&-Fg]l6#%7non-zero~assumptions~:G?(FA
F^oF^o-FdsFej]lF\u@$-Fdbm6$-Fep6$FAFfj]lF9-Fg]l6#F[bal-Fg]l6#%Clinearly~indepen
dent~assumptions~:G?(FAF^oF^o-FdsFc`alF\u@$-Fdbm6$-Fep6$FAFd`alF9-Fg]l6#FgbalC&
Fcaal?(FAF^oF^oFgaalF\uF]balF_bal?(FAF^oF^oFcbalF\uFibalF$6$Fd`alFfj]lF$%>desol
v/GetNonLinearDerivativeGR6#F]q6)F[anFcz%&NLDerGFbhs%+LinearFlagGF^z%(TDegreeGF
\\mF$C'>FKFaan>Fft-FM6$RFVF$FXF$@%FfnFgavFPF$F$F$FP>FAFhv?&F9FKF\uC%>F\xF\u?&Fj
wFftF\uC$>Fcu-F[\l6$FjwF9@$5/FcuFj`r2F^oFcu>F\xF%@$F\x>FAFehoFihoF$F$F$FiiyRFcd
m6(Fg[mFb[mF[_s%(int_conGF_eo%$lvnGF^pF$C$@$FcjxFejx@%/F[^uFjsFhjxC(>F\x-Fap6%R
6$F\erF\doF$F$F$7$FQFPF$F$F$FgpFdp>FAFhv>FKFhv?&FjwF\xF\u@$Fc[yC%>Fft-Fj]`l6%Fh
frFdpFQ>FAF]\y>FK6$FKF^hm>F9-Fev6$/Fip7#-Fcp6$Fjs-Fds6#FdpFP7$/F9-Fcr6$FbgnFh_m
F]]yF$F$F$FiiuR6%FiguF\huF\cq6%F]q%%TSolG%%LRHSGF^pF$C)>FK-FM6$F`oFP>FK-FM6%R6$
FW%#lvGF$FXF$-Fep6#-FibqFfemF$F$F$FKFfo>FK-FcrFc`\l>F9-Fiir6$FK-Fcr6$FQFir>F9-F
gs6$RFVF$FXF$4-F[bnF]oF$F$F$F9>FA-FM6$F`oF9-FfrFd\nF$F$F$%5desolv/OrderFunction
GR6$F\iwF_z6%F\hl%'FDigitGF0F1F$C'>F9-Fhal6#F[u@%5/F``oF^o/F``oFip>FK7#-F<6$,(*
&Fa[tF^o)Fajl,&FDF^oF9F^oF^oF^o*&&FejsF[`lF^o)FajlF9F^oF^oFAF^o/FAFde\l>FK7#-F<
6$,(*&,&FeilF^oFa[tFaoF^oFjjalF^oF^o*&,&FeilF^oF][blFaoF^oF^[blF^oF^oFAF^oF_[bl
>FK-FfjlFio>FK-FM6%FjjlFKF^[bl-FM6%FepFKFQF$6#F``oF$Feh`lR6$%(sep_facGFfg`l6#%(
dep_facGF^pF$@'FeenC$>FK-Fgs6%FdbmFPFQ7$*&FPF^oFKFaoFKFcc^lF\fmF[]qF$F$F$FdbxR6
$%*op_matrixGFc[m6%%)admatrixGF0Fc\lF^pF$C%>FK-&%'linalgG6#%*transposeGF]o?(F9F
^oF^o-&Ff]bl6#%'rowdimGFioF\u?(FAF^oF^o-&Ff]bl6#%'coldimGFioF\u>&FKFht-%1desolv
/adjointD2G6$Fd^blFQ-F`hyFioF$F$F$Fj]`lR6&F]eo%%degpGF,F+6&Fe_alFg[mF0%(unknown
GF^pF$C&>FKFhv>F9Fhv?(FAF^oF^oFQF\uC%>Fft-F[`s6$-F_x6$Fb[vFfoFhr>F96$F9*&)FPF\d
yF^oFftF^o@%/-Fds6#Fc_uFjs>FK6$FKFegp>FK6$FK-Fegp6#-FepF_`bl7$FicrFbgnF$F$F$Fjb
^lR6$F+F,6)%%freqG%%min_GFg[mFb[m%%ldepG%%varxGF0F^pF$C+>F\xFhv?&FftFPF\u@$-Fgn
F[]x>F\x6$F\xF`\q>F\xFa_n>FK7#-F<6$Fjs/FcuFc_m?&FjwF\xF\uC$-F\]m6%FjwFQ.Fcu>FK-
Fev6$/Fcu,&&FKFielF^oF^oF^oFK>F9-F_^uFahr>FAFhv?(FcuF^oF^oFd_mF\u@$/FfbblF9>FA6
$FA&FQFielFihoF$F$F$FfhtRF^b`lF$F^pF$@'F\[r-Fgs6&R6%FWF^is%#xvGF$FXF$4-FdbmFdem
F$F$F$FPFQFfoFgcblFPFjsF$F$F$%.desolv/SEPAR2GR6$Fe\l%%svarG6'%&lsfacGFjeuFg[m%%
dfacG%%ifacGF^pF$@%F\[rC&>FK-FM6%Feh`lF[hqFQ>FAFhv?(F$F^oF^oF$F_fmC&>FftF`[l>F9
-Fgs6%RF`[oF$F$F$@%/FfhlFQF\uF%F$F$F$FKFjfm>FK-Fgs6%RF`[oF$F$F$@%F_eblF%F\uF$F$
F$FKFjfm@%/FhdnF^o>FAFhi^lC$>F\x-Fcr6$-FM6$RF_ilF$F$F$F\[mF$F$F$F9Fh_m>FA6$FA7$
F\xFjfmFiho7#Fjh`lF$F$F$%/desolv/diffvarGR6#Fh^r6'F,Fj[mFfjo%%ind_GFb[mF^pF$@'-
Fdbm6$FPFb_pC'>FKF_`m>FAFf]p>F9Fhv?&FftFAF\u>F96$F9&FKF`blF9-Fdbm6$FPFf^mC&>F\x
FP>F9Fhv?(F$F^oF^oF$-Fdbm6$F\xFf^mC$>F96$F9Fgjn>F\x-Fep6$F^oF\xF9F^oF$F$F$Fc\vR
FcdmF``rF^pF$C%>FKFc`r>FK-%(gensepDGFeo@%Fi`rFj`rF[arF$F$F$Fd\sRFjhnF[inF^pF$C.
>FftF[u@$FainFP>FKFhv>F9Fhv>FAFjs?&F\xFPF\uC$>FAFhin@%Fjin>FKF]jn>F96$F9,&*&-%.
desolv/orderdGFj_nF^oF`jnF^oF^oFAF^o>FKFjjn>FK-FM6%RF`[oF$F$F$Fb[oF$F$F$FKFft>F
9Fh[o>F9-FM6%RF`[oF$F$F$Fb[oF$F$F$F9Fft>FKF_\oFb\oF$F$F$FhzR6$F]p%%ExpnGF$F^pF$
@%Fj]t-Fgs6%FdbmFf\alFj]m-Fgs6%FdbmFf\alFPF$F$F$Ff^blR6$%&doperGFc[m6%F[eo%&Adj
opGF+F^pF$C'>F9R6$F\bxF,F$F$F$@%/FfhlF^o*&F\[mF^o-%$QQ_GFi^mF^o*&)FaoFdhlF^o-Ff
^m6$Fh[cl-FepFehlF^oF$F$F$@%F\[r>FK-Fa]o6#-FM6%F9FPFQ>FK-F9Ffem>FA-Fgs6%FdbmFhb
pFj[cl>FA-Fbfv6$-Fcr6$FAFcanFao-%/desolv/conm2opGF]eqF$F$F$F[\yRF\pF$F^pF$FdpF$
F$F$%9desolv/RearrangeVariableGR6$%&DerivG%(FuncVarG6%Fd_xF\cqFfitF1F$C'>FKFhv>
F9F`du?&FAF9F\u>FK6$FK-FcitFg^^l>FKFbgn-Fev6%/FipFK/F6F9FPF$F$F$FjiyRFcdm6(F^zF
g[mFg]vF]gxFj]^l%&lDfacGF^pF$C'>FK-FaclFggo>FA-FfrFj\t>Fft-Ffr6$FKFfo>F\x-%/des
olv/stardepG6%FKFAFft@%5/F\xFj`r0FbdwFjsFj`rC(>F9-FadmFc_cl>Fjw-Fcr6$-FM6$RFVF$
FXF$F\[mF$F$F$F9Fir>F9-FM6$RFVF$FXF$FfhlF$F$F$F9@$3-F`]l6#&Fc]l6#%+separationG2
FjsFi`clC)-Fg]lF$-Fg]l6#%JDirect~separation~applied~to~the~equationGF`\^l-Fg]l6
#%Zwith~assumption~that~expressions~are~linearly~independentG-Fg]lFa]m-Fg]l6#%&
givesG-Fg]l6#-FhoFI@%Fb`al>Fd`al<$Fjw-FepFc`al>Fd`al<#Fjw7%F9FFF^oF$Fc`alF$%)ma
p2DerAGRF\av6)Fa\pF[esF_coFj[mF,Fcj^lF0F^pF$@)Fa\_l-%%DerAGFc\_lF`_pC$>FKFP@'-F
gn6$-Fep6$FjsF^`pFi\mC&>F97#-Fep6#F^ccl>FjwFhdn>FA7#-F<6$,&FjwF^o-Fds6#-F_x6$F_
dsF9FaoFc^_l-Fhbcl6%-FepFa``lFAFg]_l-Fgn6$F^cclFis-Fhbcl6%Fadcl7#-Fep6$FipF^ccl
Fg]_l-Fhbcl6%FadclF_\_lFg]_lF_atC*>FftFhv>FKFP?(F$F^oF^oF$FdatC$>FftF[]_l>FKF`[
l>FftFgfo>FjwFbgp>F\xFg]_l>F9Fi]_l-FhbclFe^_l-FhbclFg^_lF$F$F$%<desolv/ApplyExt
endOperator2GR6&Fe\lFc[mFd[mFicn6.F,F^[^lF[an%&lfindG%&lfjetG%%tjetGF\esFbar%'e
qnindG%%prolGFg[mF0F^pF$C:>FiuFb\m>FK7$Fj^mFcgm>F9-FM6%RF`[oF$FXF$FebuF$F$F$Ffo
FQ>FAF_an>Fft-FM6$%2desolv/DerToIndexGFA>F\x-FM6$%6desolv/DerToJetBundleGFA>Fcu
-Fap6%RF`[oF$FXF$F\_tF$F$F$FAF\x>F_u-Facl6#-F_x6$FcuFP>Fgx-F<6$7$-%&ZETA_G6%&FQ
Fh\sFQFfoFbhcl/FjtFc_m>Fgx6$Fgx-F<6$7$-F`hcl6%&FfoFh\sFQFfoF[icl/FjtFfcr>Fgx7$F
gx-F<6$7$-F`hcl6%&FftFh\sFQFfo&F\xFh\s/FjtFc^p>Fgx-Fcr6$-Fcr6$FgxFirFcan>Fey-FM
6%R6$F\bxFWF$F$F$F]g`lF$F$F$FgxF_u>Fey-Fcr6$FeyFh_m>Fey-Facl6#-Fbil6#-F_x6$/F`h
clF_[mFey>Fey-F_x6$/FcflFa\oFey>Fcu-Fap6%RF`[oF$FXF$F\_tF$F$F$F\xFA>Fey-Facl6#-
Figo6#-Fbil6#-F_x6$FcuFey@$3-F`]l6#&Fc]l6#%'gendefG2FjsFa\dl-Fg]l6#-F\_l6$%B#~T
erms~before~Modding~process~:~GFf^z>Fey-FOFj\s>Fey-%/desolv/ModdingGF]it>Fey-Fg
`nFg^z@$F^\dl-Fg]l6#-F\_l6$%A#~Terms~after~Modding~process~:~GFf^zFeyF$F$F$%:de
solv/WeighFunctionOnLHSGR6$FazF_z6%F\iw%)FuncNVarG%)TListDerGF1F$C%>FKFhv?&F9FQ
F\uC$>FA-Fgs6%FdbmFPF_`s>FK6$FKFd^pFbgnF$F$F$Fdf\lRFcdm6$Fad\l%%oeqnGF^pF$C$>FK
-F^`r6&FPFQFfoF%@%Fi`rC$>F9-Fe\vFccn@%/F9Fj`rFj`r7$7#F_`sFg`sC$>F9-Fe\v6%&FbfmF
b_lF[d]lFfo@%Fb_dlFKFc_dlF$F$F$%3desolv/SetUpSystemGRFdhm6'F0F\clFj\lFe\lFajtF1
F$C&>F\xFhv>FKFjs?&FftFPF\uC'>FKFc^l>F9Faal@$2Ffal-Fhal6#&F9F[`lC&@$Fjcl-Fg]l6#
-F\_l6$F`dlFg`dl>FA-Fcbl6#-Ffr6$Fi`dlFfo>Fhbl7$/FAFi`dlF]cl>Fft-F_x6$/Fi`dlFAFf
t>F9-Fev6%FgdlFhdlF9>F\x6$F\x7%FKF9FftFa_nF$FeglF$%4desolv/splitimpexp2GRF^b`l6
&F`b`lFab`lFbb`l%%xfacGF^pF$C)>FKFeb`l>FK-FM6$RFhanF$FXF$@%Fj]tFj]mFPF$F$F$FK>F
9F]c`l>FAF`c`l>F9-FfgoFI@%F]dq>Fft-Fgs6&R6%F`ilF\doFianF$FXF$3Fcc^l4-FdbmF_cmF$
F$F$F9FfoFK>FftF^o@%/FftF^oFedo7$-Facl6#*&F9F^oFftFao-Facl6#*&FAF^oFftFaoF$F$F$
FbfvR6#Feap6&%(num_derGFa\uF0%%num_GF^pF$C)>FKF[u@$FfjtFghl>Fft)Fajl-FhalFio>F9
7#-F<6$,&*&-Fiibl6#FejsF^oFftF^oF^oFAF^o/FA;F^oFK>F9Fejl>F9-FM6%RF`[oF$F$F$-Fc[
oFfemF$F$F$F9Fft@%32F^oF7/&F?Fe_lFao7#-F<6$&FP6#&F96#,&FKF^oFAFao/FA;FjsFb`l7#-
F<6$&FP6#&F9F@FgedlF$F$F$Fc`vR6#F_^v6$F\]oF06#%>Copyright~(c)~1999~by~Khai~VuGF
$CA>%/DesolvConstantGFjs>Fd`alFdt@%Fdj]lC$>%1desolv_pnzassumeGFfj]l>Ffj]lFdt>Ff
j]lFdt>Fbdv7#-Fcp6$FjsFgaz@$4Fi[[l>F\^zF%>FK-Fgs6%FdbmFP%'moduleG@%Fc]oC$>FeazF
g]o@$2FgazFeaz>FeazFgaz>FeazFgaz>FK-Fgs6%FdbmFP<$%-expandfactorG%-ExpandFactorG
@%Fc]o>FbizFg]o>Fbiz"%+5>FK-Fgs6%FdbmFP<$%+scanfactorG%+ScanFactorG@%Fc]o>F`azF
g]o>F`azF^o>FK-Fgs6%FdbmFP<$%(access7G%(Access7G@%Fc]o>F`]\lFg]o>F`]\lFD>FK-Fgs
6%FdbmFPFis@%Fc]o>%/DesolvFunctionGFg]o>Ff[elFjs>FK-Fgs6%FdbmFP%)constantG@%Fc]
o>F[hdlFg]o>F[hdlFjs>FKF_]o@%Fc]o>F[^oFg]o>F[^oF^o>FKFf^o@%Fc]o>F[_oFg]o>F[_oF]
_o>FKFg_o@%Fd]o@)F\`oC$>F``oFa`o>Fc`oFg]oFd`o>F``oFg]oFh`o>F``oFjs>F``oF^o>F``o
Fjs>FK-Fgs6%FdbmFP%%saveG@%3Feao/Fg]oF\u>F\_\lF\u>F\_\lF%>FKF^ao@%Fdao>FfimFg]o
>FfimFjao>Fc^oFjs>Fd_oFjs@$F\bo>F^[nF\u@$Fabo>FfalFeboF$67Ff[elF[hdlF`hdlFfj]lF
bizFeazFd`alFbdvF`azF\^zF`]\lFc^oF[_oF^[nFc`oFfalF``oF[^oFd_oF\_\lFfimF$Fh[yR6%
%*PVariableG%'POrderG%*VariablesG6'%*FunctionsG%+NewUnknownGF]jxF^jxF0F^pF$C(@%
F5>FKFedr>FKFF@$2FQFjsFejx>FAFhv>FftFhv?(F\xFjsF^oFQF\uC%>F9-F[`s6$Ff_blFK>Fft6
$Fft*&)FPF\xF^oF9F^o>FA6$FA-Fc`sFI7$-Fcr6$FgfoFh_mFihoF$F$F$FbcqR6$Fe\lF]eo6$%%
numeG%%denoGF^pF$C%>F9-%.desolv/facsepG6$FagoFQ>FKFj^cl@%F\hnC%>FK-FM6%F\aelF_h
nFQ>FK-%(coltermGFio-FM6%R6$%$numG%$denGF$FXF$7$*&F\[mF^oF^a`lFao*&FfhlF^oF`duF
aoF$F$F$FKF9C$>FK-F\aelFj\t7#7$*&FbfmF^oF_`sFao*&F[d]lF^oFg`sFaoF$F$F$%9desolv/
HighestDerivativeGR6$F^zF_zF`zF^pF$C(>FKFaan@$F[[lF\[l>F9F`[l?&FAFKF\u@$2Fjs-Ff
[l6'FAF^oF9F^oFQ>F9FA>FftFj[lF]\lF$F$F$F[jyRFcdm64Fg]vF]gxF*Fj]^l%&tlvarGF[_sF`
btF[b^lFg^cl%$vl1G%#L1GFfa^lF^zF[goF]eoFh`elF\dblFg[mF^pF$C(>FhjrFj^cl>FK-Ffr6$
FhjrFQ>F9-Ffr6$FhjrFfo>Fft-Fb_cl6%FhjrFKF9>FiuFhv@)FdcvFj`r/FjfmFjs-F][l6#-Fjiy
Fddel/FjfmF^oC$?&F\byF^hmF\uC$>Fh[z-%3desolv/indsepTypeIG6&FhjrFKF9F\by@$0Fh[zF
j`r-F][l6#7%&Fh[zFb_l7$Fj^m-Fep6#&Fh[zFe_lFipFejxC->FAFF>FeyFj`r>FftF^hm>FadrFh
jr?(F$F^oF^oF$FagpC/>FhgrF^hm>FftFhgp>Fjt-Fgs6%FdbmFKFhgr>Fjt-Ffr6$FadrFjt>FeyF
j`r>F\xFj`r>Fey-Fca^l6&FadrFKF9Fjt@$/FeyFj`rFejx>Fiu6$Fiu-Fep6#F^i_l>FadrFai_l@
$/FadrFjsFejx>F\x-Fb_cl6%Fadr-Ffr6$FadrFKF9@$5Ff_cl/FbdwFjs>FftFF@%/-FbgoF\hyF^
o>FcuFFC%>FibyFehel>Fadr-Facl6#-FigoF\hy>FcuF]hel@%/FcuFF>FgxF9>Fgx-Fgs6%R6$F\d
oFcd\lF$F$F$Fc^^lF$F$F$F9Fcu@%330F\xFj`rFahel0FadrFjs>Fey-F]cn6%FadrFKFgx>FeyF]
gy@$-F`]lFbfy>Fey-FM6%R6$FW%%DenoGF$F$F$Fff^lF$F$F$FeyFiby@$Ff`clC'-Fg]l6#%<Ind
irect~separation~type~IIG-Fg]l6#%3Separated~equationG-Fg]lF]o-Fg]l6#%4Resultant
~equationsG-Fg]lFg^z7%FeyFFF6F$Fc`alF$%0desolv/SEPARATEGRFcdm6(F*Fh]^lF0%'occfn
sG%(occvarsGFa[xF^pF$C&>FK7#Fjgn>F9Fhv?(F$F^oF^oF$Fd]oC$>FAFjs?&FjwFKF\uC&>FAFh
in>Fft-FfrFhju>F\x-Fgs6&R6%F`ilFa[o%"cGF$F$F$3Fi`t4-Fdbm6$FfoFPF$F$F$FfoFjwFft@
%FcgsC$>F96$F9Fjw>FKF`[wC$>FK7$-Fep6#F`[w-Fep6#-Fep6$F^o-%/desolv/sepsum1G6$FPF
\xF^flF`aqF$F$F$Fi\\lR6%F]q%%LHisGF^[^l6'F`iyFbiyFaiyFciyFeiyFejqF$@$Fi[\lC)>FK
-FM6$FhalFP>F9Fhbbl>FA-FdcoFahr>FftFbg`l>F\x-Fds6#Fegx-Fg]l6#-F\_l6&%3No.~equat
ion(s)~:~GFd_m%4,~No.~unknown(s)~:~GF\x@$Fa]z-Fg]l6#-F\_l6(%1System~length~:~GF
ft%3,~~Smallest~eqn~:~GF9%2,~~Largest~eqn~:~GFAF$F$F$Fic]lR6%Fh^rFe\l%'maxderG6
$Fa\uFj[mF^pF$C%>FK-Fap6%RF`[oF$FXF$F^doF$F$F$FfoFgp>F9-Fap6%R6$F^abl%%ordxGF$F
$F$-FcpFfemF$F$F$FdpFK@%Fj[pFQ-Ff^m6$FQF`\oF$F$F$%(int/DerGRF$6(FfjoF`d`lFb[m%(
DerpartG%(nonDerpGFg[mF^pF$C%>FAF?-F\]m6%.%#_XG-Fep6$F6FAFg`x@%2Fjs-Fep6$F9Fh`p
C%>FKFh`p>FK-Fev6$/F9,&&FKFIF^oFaoF^oFK-Fev6$Fc^clFAFj`rF$F$F$%9desolv/AllLevel
ExcludeICGR6&F]qFgcnF[hsF_z6)F\hs%)LevelSysG%.AllExcludeICsGF[r%'HDTermGF0%(TSy
stemGF1F$C'>FAFhv>F9FP>FKFjs?(F$F^oF^oF$32FKFajlF_hqC'>FKFc^l>F9-Fhgs6'F9FQFfoF
hrFK@$F_hqC$?(FjwF^oF^oFhdnF\uC'>Fft&F9Fa]m>F\x-FgbelFcal>Ff\nFg\n>Fft7%Ff\nF\x
Fft>F9-Fev6$/FjwFftF9>F9-F^\l6&F9FQFfoFhr@$/Fd_oFjs>F9-%7desolv/RemoveRedundant
GF\\l@$F_hqC%>Fcu-FM6$RFVF$FXF$FjhlF$F$F$F9>FA6$FAFcu@$F[glC$-Fg]l6#-F\_l6$%:Th
e~excluded~ICs~-~level~GFK-Fg]l6#-Fe`]lFielFihoF$6$Ff\nFd_oF$F^]dlR6$Fe\lF*6(Fe
[alFf[alFg[mFg[alFh[alF[rF\\mF$C+@$FcyFj`s@$Fd\xFghl>FAF^[q>FftF\u>F9-FbanF@>F\
xFjs?(F$F^oF^oF$Fg]alC'>FftF%?&FjwFQF\u>FA-Fbil6#-F_x6$FjwFA>FAFb^al>FKFdgfl@$F
g^alC%>F9FK>F\xF^jq>FftF\u@$F]_alF^_alFAF$F$F$%2desolv/initializeGRF$6#F0F^pF$C
+>Ff[elFjs>F[hdlFjs@%Fdj]lC$>F`hdlFfj]l>Ffj]lFdt>Ffj]lFdt@%-F`]l6#Feaz@$Fbidl>F
eazFgaz>FeazFgaz>FbdvFdhdl>Fd`alFdt@$4-F`]l6#F`az>F`azF^o@$4-F`]l6#Fbiz>FbizF^j
dl@$Fhhdl>F\^zF%F$6,Ff[elF[hdlF`hdlFfj]lFbizFeazFd`alFbdvF`azF\^zF$%2desolv/che
ckICminGR6#%)ICvertexG6%Fi[m%$pICGF^eoF1F$C'>FKFebn@$4-F`]l6#Fc``l>Fc``lFF?&F9F
c``lF\uC$>FAFac]l@&5/FAF\u/FAFb[pFe[]l/FAF\\p>Fc``l-F_x6$/F9FhvFc``l>Fc``l7$-Fe
pFb[glFPF\uF$6#Fj_`lF$Fb_sR6%Fe\lF\dblF,6(%'deneqnGF[]lF^dblFieu%&lpfacG%&lffac
GF^pF$C(>FKFago>F9Fj^cl>Fft-Fgs6%Fdbm-FibpF\gnFQ@%/FftFdt>FA7$F^oF9C%>F\x-%'coe
ffsG6%F9Fft.Fjw@$-Fdbm6$-FM6&Fi[]lF]]tFQFfoF\uF\[l>FA-Fap6%R6$%%fac1G%%fac2GF$F
$F$F^[yF$F$F$F]]tFa_n@$0FKF^o@%-FdbmFj\t>FA-FM6%R6$FWF\belF$F$F$7$F\[m*&FfhlF^o
FQFaoF$F$F$FAFK>FA-FM6%RF__glF$F$F$7$Fcg^lFfhlF$F$F$FAFK-FM6%R6$FcenF]eo6#%&lop
owGF$F$@%-Fgn6$FfhlFh_mC$>FK-%(ldegreeGF\e^l@%Fjj`l7$*&F\[mF^o)FQFKF^o*&FfhlF^o
Fg`glF^oFPFPF$F$F$FAFQF$F$F$FcdlRFci]l6'FeapFg[mFa\p%&stermGF\goF^pF$C)>FK-Fgs6
$RFhanF$FXF$-F[bn6#F^a]lF$F$F$F[fs@%Fafdl>FK-%0desolv/sortDer2G6$FKFcfdl>FK-%0d
esolv/sortDer3GFio>FftFjgn@%-FgnF]h^l>Fft-Fcr6$FftFir>FftFb\q>F9Fhv?&FAFKF\uC%>
F\x-Fgs6%FdbmFftFA>Fft-FatF\c\l>F96$F9-Ffgo6#-Fa]o6#-FcrF[jn-Fcr6$F`aqFh_mF$F$F
$FgeflR6$F]q%&LExICG6'%'NewICsG%&ICEqnGF^zF[r%.RedundantFlagGF1F$C%>FKFhv?&F9FQ
F\uC%>F\xF%?&FAFP4F\xC$>Fft-Fa]o6#,&&FAF[`lF^oFi`dlFao@%Fd`l>F\xF\uC$>Fft-Fa]o6
#,&FfdglF^oFi`dlF^o@$Fd`l>F\xF\u@$F`dgl>FKF_gnFbgnF$F$F$%*diff/DerAGRF$F_d`lF^p
F$C'>FK-FhbclFdd`l>F9Fid`l>FftF^jo?&FAF9F\u@%Fae`l>FftFfe`lFj`sFje`lF$F$F$%=des
olv/VariableWeightForAutoGR6$F`_^lF[an6&Fd_^l%5DerivativeOfFunctionG%,WeighOfVa
rsG%-NListVWeightGF1F$C%>FftFhv?&FKFPF\uC&>F9-Fgs6%FdbmFQFbfm>FA-F^]y6$F[d]lF9>
FK-Fev6$/F6FAFK>Fft6$FftFKFgfoF$F$F$FjiqRF_gq6(Fb[xF0Fc\lFbjtFdgq%*AcsendingGF^
pF$C(>F\xF[u@$2F\xFipFghl@%2FipF7>Fjw&F?F[`l>FjwF^o>FKF[hq?(F9F\xFaoFipF\uC%>Ff
tF^o?(FAF^oF^oFjjtF\u@$2Fjs*&-Ff[l6%-FepFf\o-Fep6#F^^wFQF^oFjwF^oC$>FK-Fev6%/FA
F^^w/FhinFg\oFK>FftFA@$2FftFjjt>F9FftFKF$F$F$FfaelR6#%%lgrpG6(F0%$lDIG%%tgrpG%&
grpDIG%$temG%%tsumGF^pF$C(@$Fh[^lFghl>FAF\bq>FKFjs>F9F[``l?(F$F^oF^oF$0FAFdtC'>
FKFc^l>FftFcct>FAFggt>F\x-Fgs6%R6$%%gdi1G%%gdi2GF$F$F$@%/F\[mF^a`lF\uF%F$F$F$FA
Fft@%/F\xFdt>Fh^pFftC%>Fjw,&-Fcr6$-FM6$R6#%$gdiGF$F$F$FfhlF$F$F$F\xFh_mF^oFjfmF
^o>Fh^p7$F^hmFjw>FA-Fat6$FAF\xFe]uF$F$F$%2desolv/CheckInputGRF`ar6+F.Fe\lFa\pFe
ap%%HderG%&lHderGF[es%(degHderGF06#%9Copyright~(c)~by~Khai~VuGF$C,@$4-FgnF\hq-F
`io6#%?The~1st~argument~is~not~a~listG@$4-FgnF^^y-F`io6#%?The~2nd~argument~is~n
ot~a~listG@$4-Fgn6$FfoFcan-F`io6#%?The~3rd~argument~is~not~a~listG@$-Fdbm6$-FM6
$R6#%"uGF$FXF$-Fgn6$FPFc`pF$F$F$FQF%-F`io6#%OElements~of~the~2nd~argument~are~n
ot~variablesG@$-Fdbm6$-FM6$R6#F\doF$FXF$Fi^hlF$F$F$FQF%-F`io6#%OElements~of~the
~3rd~argument~are~not~variablesG>FKFc`r?(F_uF^oF^oFfwF\uC&>F9&FKFjcw>FftFjip@$F
ahx-F`io6%F9%A~does~not~contain~any~variables~GFQ@$333-Fgn6$F9Fin4-Fgn6$-F\oFI<
%FgenFefnFh_m-Fgn6$F[ahlFis/-Fep6$FjsF[ahlFedq@$-Fdbm6$-Fev6$/F_uFhvFKF[ahl-Fg]
l6#%IWarning:~CIRCULAR~substitution~may~occurG>F`]m-FM6$RF_ilF$FXF$&F[crF]oF$F$
F$Ffo>Ff]m-FM6$RF_ilF$FXF$F\drF$F$F$FQF\uF$FjamF$%0desolv/subsDer3GRFbfpFdfpF^p
F$C(>Fft-Fgs6$RFhanF$FXF$FaaglF$F$F$F[fs>FftFifp>F\xF]gp>F9FP?(F$F^oF^oF$FagpC*
>FAFegp>FftFhgp>FK-Fap6%RF[doF$FXF$F^doF$F$F$Fh`pFgjn@$-Fdbm6$-FM6$RF_ilF$FXF$F
dhpF$F$F$FKF\uF[al>FKFhhp@%Few>F9-FaclF_ip>F9-FaclFdip>Fft-Fgs6$RFhanF$FXF$Faag
lF$F$F$F`]gl>FftFifp@%Fh[[l-FcdlFIF9F$Fj[[lF$%+liebracketGR6%%(vector1G%(vector
2GF,6'Fg[m%&lvec1G%&lvec2G%%vec1G%%vec2GF^pF$C%@%F\[r>Fft-FM6%RF`[oF$FXF$*&F\[m
F^o-Ff^m6$FijqF_\clF^oF$F$F$FPFfo>Fft*&F\[mF^o-Ff^m6$Fh]sF_\clF^o@%-Fgn6$FQFh_m
>F\x-FM6%RF`[oF$FXF$FiehlF$F$F$FQFfo>F\x*&F^a`lF^o-Ff^m6$Fh]s-Fep6#F`duF^o>FK-F
acl6#,&-F_x6$/Fh]sF\xFftF^o-F_x6$/Fh]sFftF\xFaoF$F$F$Fa^tR6$Fb[mFd[mF`\pF^pF$@)
Fa\_lF%F`_pFiatF_atC%>FKFP?(F$F^oF^oF$-Fdbm6$FKFf^m>FKF`[lFgatF%F$F$F$F[enR6$F]
qF`q6)F[anF\an%%LJetG%)LSortJetGFbhs%)PickLEqnGF0F\\mF$C&>F\xFL>FK-F_calFj_n>FK
-Fgs6$RFhanF$FXF$-F[bn6#2FjsF^bnF$F$F$FK@%1Fi_nFfwC*>FA-FM6$FagclFK>F9-Fap6%RF`
[oF$FXF$F\_tF$F$F$FKFA>F\x-Fbil6#-F_xF[\s>Fft-%<desolv/SortByEquationWeightGFd\
hl>F\x-Fcr6$F\xFir>FjwFhv?(FcuFi_nF^oFfwF\uC$>Fjw-Fiir6$F\x-Fcr6$&Fft6#F[jsFir@
$0FjwFhvF^fl@%F`[ilC&>Fjw-Fgs6$RFVF$FXF$FaialF$F$F$Fjw>F9-Fap6%RF`[oF$FXF$F\_tF
$F$F$FAFK>Fjw-Fcr6$-Fbil6#-F_xFc]flFcan-FM6$FhoFjw-F`io6#%OCannot~pick~out~the~
left~hand~side~derivativesGFc\ilF$F$F$FijmR6$%)UnitEqn1G%)UnitEqn2G6-%'EqnNo1G%
%Eqn1G%(HDTerm1G%'EqnNo2G%%Eqn2G%(HDTerm2G%(Hcoeff1G%(Hcoeff2G%'ACof12G%+CrossI
ndexG%0CrossDerivativeGF1F$C6>FKF\[m>FAFfhl>F9Fjhl>FftF^a`l>FjwF`du>F\x&FQF[`l@
$-F\]m6$<$FKFftFf`nFj`s@$0-Fep6$F^oF\hmFbjmFj`s>Fgx-Fap6%Fdco-Fep6$FipF\hm-Fep6
$FipFdjm>Fey-Fev6$/FipFgxF\hm@$4-FgjflFg^zFj`s>Ff`n-F\w6$Ff`nFa^il>F9-Fic]l6%F\
hmF9Fgx>F\x-Fic]l6%FdjmF\xFgx>Fcu-%3desolv/coefficientGF[f]l>Fiu-F_`il6$F\xFey>
F_u-F\jp6$FcuFiu-F_go6#F^]w-F_go6#Fb_w-Facl6#,&*&F^]wF^oF9F^oF^o*&Fb_wF^oF\xF^o
FaoF$F^]tF$FierRF`ar6#Ff[mF^pF$C&@$F_]hl-F`io6#%BThe~first~argument~must~be~a~l
istG@$Fe]hl-F`io6#%CThe~second~argument~must~be~a~listG?&FKFQF\u@$F[`p-F`io6#%Y
Each~element~of~the~second~argument~cannot~be~a~functionG@$F[^hl-F`io6#%BThe~th
ird~argument~must~be~a~listGF$F$F$F\jyRFcdm60%)ode_funcG%(ode_ordG%(ode_varG%(o
de_eqnG%*LConstantG%'TConstGFg[mF[goFb[mF[_sF\_s%'HcoeffG%'HDTodeGF0F^pF$C2>FK-
F]\u6%FPFeemFfo@$/FKF%Fejx>F_uFbfm>FA&FKF[`l>F9&FKFfdr>FKF[d]l@$4-Faj\l6$FPF_uF
ejx>Fft-%2desolv/ConvertODEG6&FPF_uFKFA>Fcu-%'dsolveG6$Fft-%&ODE_YGF@@$5F_ds-Fd
bm6$Fcu<%F]et%$sumG%&DESolGFejx>Fcu-Figo6#-F_x6$/FedilFK,&-F\oFielF^o-F`oFielFa
o>F\x-Fgs6$%6desolv/DSolveConstantG-FibpFiel>FgxFhv?&FjwF\xF\uC%>Fey-F[`s6$-F_x
6$Fb\t7#FicwFQ>Fgx6$Fgx-Fc`sFg^z>Fcu-F_x6$/FjwFeyFcu>Fjt-%1desolv/HighCoeffGF\d
il@%4-Fdbm6$FjtFedqC$-F_goFh\s7$Fcu7$Fj^mFgxC%@%2FjsF9>Fhjr-Ff^m6$FK-FcpFeho>Fh
jrFK>Fiu-Facl6#,&FPF^o*&FjtF^oFhjrF^oFao7%FcuFfgil7$FjtFiuF$F$F$FajhlRFid_l6&Fc
zF[e_l%'NumDerGF0F^pF$C)>FKF[hq>F9-Fgd_lFj\t>FAFgil>F97#-F<6$,&*&Fe\wF^oF`jlF^o
F^oFftF^oFbjl>F9Fejl>F9-FM6%Fc[oF9F`jl7#-F<6$&FK6#Fe\w/Fft;F^oFhdnF$F$F$%?desol
v/FindRedundantExpressionGR6$%*LVecGroupGF\cq6(%'TGroupG%'RGroupG%,RemainGroupG
Fh`_lF^[^l%+LRedundantGFi`_lF$C'>FAFP>FjwFhv?(F$F^oF^oF$F_ctC'>FKF\hm>FAFggt?&F
9FAF\u@$-%3desolv/IsSubVectorG6%Fg`sF[d]lFQC$>Fjw6$Fjw/F_`sFjs>FA-F_x6$F_\glFA>
Fft-%=desolv/FindIndependentFactorG6$F[d]lFQ@$0FftF^o>Fjw6$Fjw/Fbfm*&FbfmF^oFft
Fao>F\x-FM6$R6#%#tgGF$FXF$FebnF$F$F$FP-FO6$F]]tF\xF$F$F$FcfwR6%FcjoFdjoF+6(Fg[m
F0Fb_xFc_x%&countGF]atF^pF$C&>Fjw-FM6$RFhanF$FXF$@%Fj]tFj]mFPF$F$F$Ffo>FKFgj`l@
$F_`xC%>F\x-F_^u6$-Fds6#Fgp-Fds6#Feco?(F9F^oF^oF\xF_`x>FK,&-Fep6$F9FgpF^o-Fep6$
F9FecoFao@$F_`xC%-F\]m6%FebnFjwFj`x-F\]m6%F_esFjwFjgv>FK,&FftF^oFAFaoFKF$F$F$%1
desolv/diffSpolyGR6%%%eqn1G%%eqn2GF_z6,%'lHead1G%'lHead2G%&Head1G%&Head2G%%acof
G%&SpolyGFj[m%%dcofG%&deqn1G%&deqn2GF1F$C)>FK-FbalF_cm>F9-FbalF__r>FA*&)FbfmF[d
]lF^oFecilF^o>Fft*&)F_`sFg`sF^oFi`dlF^o>F\x-F\jpFhi^l@%/FbfmF_`s@%3/FedwFA/Fbdw
Fft>FjwFjs>Fjw-Facl6#,&*&FedwF^oFQF^oF^o*&FbdwF^oFPF^oFaoC->Fiu-Fgbo6$FbfmF_`s>
Fcu-Fap6%Fcp-Fep6$F6FbfmF`^n@%/FidvFjs>F_uFP>F_u-Ff^m6$FPFhel>Fcu-Fap6%Fcp-Fep6
$F6F_`sFa^n@%Fcbjl>FgxFQ>Fgx-Ff^m6$FQFhel>FK-Fbal6$F_uFfo>F9-Fbal6$FgxFfo>FAFe`
jl>FftFh`jl>F\xF[ajl>Fjw-Facl6#,&*&FedwF^oFgxF^oF^o*&FbdwF^oF_uF^oFao-FcdlFa]mF
$F$F$F]\jlR6$%$VecGF\cq6%Fh`_l%)CoeffEqnG%*IntFactorGFi`_lF$C'@%F\[r>F9-FM6$R6#
%#vcGF$FXF$FebnF$F$F$FP>F9Febn>FA-F[hoFI>F9-Fbil6#*&F`cdlF^oFAFao@'F]dq>FK-Fgs6
%R6$%#ftGFbhalF$FXF$Fbc^lF$F$F$F97$Fj^mFedq-Fdbm6$F9F^fjl>FKF^o>FKF9*&FKF^oFAF^
oF$F$F$F_dilR6&Fe\lFb[mFh^rF_bm6$FeapFd\uF^pF$C'>FK-Fgs6$RFeesF$FXF$FgesF$F$F$F
[fs>FK-Fgs6%FdbmFKFfbm>F9Fbcm>FK-FM6&R6%Fa\pF]eoFa\u6#F_epF$F$C$>FK,&F[^uF^oFfo
Fao@%F_`x/FP-FfdilF[_m/FP-Ff^m6$F[hjl-Fcp6$FQFKF$F$F$FKFhrF9-F_x6$FKFPF$F$F$%0d
esolv/SortFuncGRFad]l6'%'num_fnGFa\uF0Fhddl%&lndepGF^pF$C*>FKF[u@$FfjtFghl>F\x-
FM6$RFhanF$FXF$@%Fj]tF[uFjsF$F$F$FP>FftF]edl>F97#-F<6$,&*&&F\xF@F^oFftF^oF^oFAF
^oFgedl>F9Fejl>F9-FM6%RF`[oF$F$F$F^fdlF$F$F$F9Fft@%F`fdlF^gdlFdfdlF$F$F$%;desol
v/ApplyExtendOperatorGRFjecl6/F,F^[^lF[anF\fclF]fcl%'lfjet2G%(ljetModGF\esFbarF
_fclF`fclFg[mF0F\\mF$CD>F_uFb\m>FKFdfcl>F9-FM6%RF`[oF$FXF$FebuF$F$F$FfoFQ>FAF_a
n>FftF[gcl>F\xF_gcl>Fjw-FM6$Fj\alFA>FcuFjw>Fiu-Fap6%RF`[oF$FXF$F\_tF$F$F$FAF\x>
Fgx-Facl6#Febp>Fey-F<6$7$-F`hcl6%&FQ6#FhjrFQFfoF^\[m/FhjrFc_m>Fey6$Fey-F<6$7$-F
`hcl6%&FfoF_\[mFQFfoFh\[m/FhjrFfcr>Fey7$Fey-F<6$7$-F`hcl6%&FftF_\[mFQFfo&F\xF_\
[m/FhjrFc^p>Fey-Fcr6$FheyFcan>Fjt-FM6%RF`jclF$F$F$F]g`lF$F$F$FeyFgx>FjtFc`u>Fjt
-Facl6#-Fbil6#-F_x6$F[[dlFjt>Fjt-F_x6$/FicwFa\oFjt>Fiu-Fap6%RF`[oF$FXF$F\_tF$F$
F$F\xFjw>Fjt-F_x6$FiuFjt>Fjt-FOF^g]l>FA-Fcr6$-FbanFh\sFcan>FjwF^[[m>Fiu-Fap6%RF
`[oF$FXF$F\_tF$F$F$FAFjw>Fjt-Facl6#-Figo6#-Fbil6#F\_[m>Fcu<$FhfrFhel>Fcu-Fat6$F
cu-Fcr6$FfoFir@$F^\dl-Fg]l6#-F\_l6$Fi\dlFh][l>Fjt-F`[al6%FjtFhrFcu>FcuFhjz>FjtF
_fz>FcuFc`[m@$F^\dl-Fg]l6#-F\_l6$Ff]dlFh][l7$FjtFcuF$F$F$FaeelR6&Fe\lF+F,F\dbl6
1Fg]v%%lnfnG%%nunkG%%xfunGF]gx%%fdepG%'fnodepG%(nunkdepG%'divfacG%,implicitsumG
%,explicitsumG%&exsumGFeg`lFh]^lF[goF^pF$C:>FKFeem>F\xFegx>Fft-Fgs6%FdbmFKFhr@$
0FbgpF^oFejx>FftFegp>Fjw7#F`\q>Fcu-Fgs6%R6$F\do%#lxGF$FXF$F`g\lF$F$F$F\xFjw>Fhg
rFjgn>Fgx-Fgs6%FdbmFhgr-Fep6$FjsFft>Fey-Fbil6#,&FhgrF^oFgxFao>Fgx-FfgoFd`v>F_uF
^o@$-Fgn6$FgxFgen>F_u-Fgs6%FdbmFgxFcu@$0F_uF^oC&@$2Fjs-Fds6#-Ffr6$F_uFKFejx-F_g
oFjcw>Fgx*&FgxF^oF_uFao>Fey-Fa]o6#-Facl6#*&FeyF^oF_uFao@$-Fdbm6$FgxFcuFejx>Fey-
FbcqF]it>FhjrFhv>FadrFhv>F9Fhv?&FjtFeyF\u@'2Fjs-Fds6#-Ffr6$F_fzFcu>Fadr6$FadrFh
jz-Fdbm6$FhjzFcuC(>Fiu-Ffr6$FhjzFjw>FA-F[`sF]]r>Fadr6$Fadr,&FAF^oFhjzFao>Fhjr6$
Fhjr*&FAF^oF_fzF^o>FA-Fc`sF@>F9Fht>Fhjr6$Fhjr*&F_fzF^oFhjzF^o>Fhjr-Facl6#,&FgxF
^o-Fcr6$7#FhjrFh_mF^o>Fadr-Fa]o6#7$FhjrFadr@$Ff`clC*F^acl-Fg]l6#%FIndirect~sepa
ration~type~I~applied~toG-Fg]l6#-F\_l6$%5the~equation~w.r.t.~GFhrF`\^lFfacl-Fg]
l6#-FhoF\hy-Fg]l6#%=with~new~unknown~function(s)G-Fg]lFi[o7$FadrF`hmF$F$F$%<des
olv/AccumulateFuncWeightGR6%F\iw%(FWeightGFb_^l6%F0F`iw%+NewFWeightGF1F$C%>FAFP
?(FKF^oF^oFd_mF\uC$>F9Fb`^l@%F\`^lC$>F9-Fev6$/FipF\_mF9>FAFd`^lC$>F9-Fev6$/F6,&
F\_mF^oFi`dlF^oF9>FAFd`^lFAF$F$F$FagclRFd_`l6#Fj[mF^pF$C$>FKF\^y-F\_l6%Febn%"_G
Fa\oF$F$F$F_e^lRFjhn6&F[]lFhddlF`iyF0F^pF$C(>FKF[u>F9F]edl>FA7#-F<6$,&*&-Fhal6#
FgjsF^oF9F^oF^oFftF^o/FftFhedl>FA-FfjlF@>FA-FM6%RF`[oF$F$F$F^fdlF$F$F$FAF97#-F<
6$&FP6#&FAF`blFj\\mF$F$F$%'ordervGR6$F[eoFc[m6'%*deriv_varG%'fn_varGFj[mF0Fb[mF
^pF$C(>F9F_`m@$0FhdnFd_m-F`io6#%TINPUT~ERROR~!!!~Check~input~functions~and~vari
ablesG>FAFf]p>F\x-(%#D_GFf^pFi^m>FK-F<6$&F9Fe]\m/FftFc^p-Ff^mFf_tF$F$F$Fj[tR6&F
h\ilFi\ilF]pF_zFj\ilF1F$C9>FKF\[m>FAFfhl>F9Fjhl>FftF^a`l>FjwF`du>F\xF]^il@$F_^i
lFj`s>FA-F[z6%F9FfoFhr@$/Fd^pFjsFj`s>Fjw-F[zF]dw@$Fg]rFj`s>FgxFg^il>FeyF^_il@$F
b_ilFj`s>Ff`nFe_il>F9Fh_il>F\xF[`il>FcuF^`il>FiuFa`il>F_uFd`ilFf`ilFh`ilFj`ilF$
F^]tF$%0desolv/isderiv2GRFbjo6#FfjoF^pF$@%Fb[hl@%2FfhlF`duF%C$>FK-Fap6%R6$%#a1G
%#a2G6#%"dGF$F$C$>FKF`do@%2FKFjsF%FKF$F$F$FjhlF]^il@%-Fdbm6$FKF%F%F\uF%F$F$F$Fc
\dlRF`ar6.FbarFcarFe\lF0Ff`wFg`wFdarF,Fb[m%&lzetaGFh`wFj`wF\\mF$C5-Fe\hlFccn>FK
F]fr>FiuFbgm>Fjw-FM6%RFfawF$FXF$FebuF$F$F$FQFfo>Fcu-Fhau6%FfoFQFiu>Fey-Fdcn6$FP
Fjw>F9Fhv@$F^\dl-Fg]l6#%=Calculate~defining~equationsG?&FAFeyF\uC(@$F^\dl-Fg]l6
#-F\_l6$%7processing~equation~:~G-Fe`]lF@>F\x-F^jjl6&F^cwFfoFQFey>Fjt-FibqFidw>
Fjt-FM6%F^alFjtFcu@$F^\dl-Fg]l6#-F\_l6$%8#~defining~equations~:~GFh][l>F9F_ew>F
9F`aq@$F^\dl-Fg]l6#-F\_l6$%A#~defining~equations~in~total~:~GFhdn>F9-FM6&FOF9Fc
uFF>F9-Fj[s6$F9Fcu>F9F]\s>F9-FM6%F^alF9Fcu>F9-Fd\s6#Fe]u>F9-Fcr6$F9Ff^m@$F^\dl-
Fg]l6#-F\_l6$%H#~defining~equations~after~reduction~:~GFhdn7%F9FcuFiuF$F$F$%1de
solv/conDif2opGRFci]l6%Fja^lFg[m%%itemGF^pF$C(@$FcyFj`s>FK-Faey6$FPFjjq>FK-Fap6
%RF`[oF$FXF$7$FP7#FecoF$F$F$FbfmF[d]l>F9Fjs?&FAFKF\u>F9,&F9F^oFAF^oF9F$F$F$%6de
solv/IsSubComponentGR6%%&Comp1G%&Comp2GF\cq6&F[r%%TUnkG%%TVarG%%FVarGFi`_lF$@%/
FgpFecoC%>FK/F_esFebn>F9-Fban6#F_es@%F_hqC(>F9Fgfn>FK-FiirF_gn>FA-Fcr6$-Fgs6%R6
$%#vrGFWF$FXF$Fi`tF$F$F$FfoFKFir>Fft-Fcr6$-Fep6$F6F9Fir>FA-FatFhi^l@%F]`\mF\uF%
F%F%F$F$F$Fd\vRFcdm6.Feap%%lfunGFghjl%'maxdepG%'maxfnsGFb[m%%fsolGF0F[cilF[go%%
ldegGF[\mF^pF$C7>FK-Fcr6$-Fgs6$RFeesF$FXF$FgesF$F$F$F[fsFcan>FA-FM6$RFhanF$F$F$
F^galF$F$F$FK>FftFb^u@$2FftF__flFejx>F\xFhv?(FiuF^oF^oFd^pF\u@$/-Fep6$FiuFAFft>
F\x6$F\x-FepF\dp>F\xFa_n>FKFhv>F9Fhv?&FjwF\xF\u@%2Fjs-Fcr6$-Fep6$FipFjwFh_m>FK6
$FKFjw>F9Fc]fl>FKFbgn>F9F`hm@$F[fqFejx>F9-FM6%RF`[oF$FXF$@%-Fdbm6$FQFebnFhvFPF$
F$F$F9FK>Fey-FM6%RF`[oF$FXF$-F[\lFizF$F$F$F9FP-F\]m6%-F_^u6#FfjzFeyFjt>Fjw-FepF
^g]l>Fcu-FiirFac\m@$5F_ds2F^o-Fds6#<#FcuFejx>F_u-FadtFac\m@%4-Fdbm6$F_uFedqC$@$
Fjd[mFbe[m7#/FjwFcuC$>Fgx-Facl6#,&FPF^o*&F_uF^o)Fjw-F[\lFac\mF^oFao7$Fg_]m7$F_u
FgxF$F$F$FjemR6%FedmF+F,6)%$EqlG%'sepgrpGFg[mF\dblF*%%olfnG%&olvarGF^pF$C&>FKF\
_^l>FAFhv?(F$F^oF^oF$F_fmC(>F9Fbfm>FKFdfm>Fjw-Ffr6$Fg`sFQ>Fcu-Ffr6$Fg`sFfo>Fft-
Fgs6%R6$F]eoF^isF$F$F$Fc^^lF$F$F$FcuFjw@%Fahx>FAFehoC%>F\x-Fbcq6$Fg`sFft>F\x-FM
6%RF`[oF$FXF$7$*&F\[mF^oFQF^oFfhlF$F$F$F\xF_`s>FK7$F__tFa\oFihoF$F$F$F_`ilR6$Fe
\lF\goF``rF1F$C$>FKFjgn@'F\hnC$>FKFabq-Facl6#*&FKF^oFQFaoFj^glF`c]mFjsF$F$F$Fh_
lR6%Fe\lF*F_z64FeapFg[m%(result2GF0Fd\l%$DTkG%%pDTkGFihm%&pHDTiGF`]jl%#EiG%'coe
ffiG%'coeffrG%'acofriG%)exponentGFi\lFa\pFg[alF1F$C)@$FcyFj`s@$Fd\xFghl@%F5>Fgx
Fedr>FgxFd_m>F9Fjgn>Fh[zF\u?(F$F^oF^oF$Fh[zC&>Fh[zF%>FKFjip>FK-FM6$RFVF$FXF$F[]
qF$F$F$FK?(FftF^oF^oFgx4Fh[zC&>Fiu&FagrFe_l>F_uFa^n>FiuF`^n?&F\byFKF`e]mC/@$4-F
dbm6$F9)&F\byFb_l&F\byFe_lF[al>FjwF]f]m>FcuF^f]m@$0F`gvFi^qF[al@'/FjwFiuC%@$2Fc
uF_uF[al>Fey&FagrF[`l>Fjt-F_`il6$Fey)FiuF_u/-F`jo6$F\]]m-Fep6$FipFiuF\uC$>Fey-F
ic]l6%FiuFif]mF\]]m>Fjt-F_`ilFje^lF[al>Fhjr-F_`il6$F9)FjwFcu>Fadr-F\jp6$FhjrFjt
@%Fdf]m>Fhgr,&FcuF^oF_uFao>Fhgr,&FcuF^oFaoF^o@%/F[^oF^o>F9-Fa]o6#,&*&&FadrFb_lF
^oF9F^oF^o*(&FadrFe_lF^o)FjwFhgrF^oFeyF^oFao>F9-FaclFih]m-F_go6#F\i]m-F_go6#F^i
]m@$F`gsFj`s>Fh[zF\u-Fa]o6#*&F9F^oFbejlFaoF$6#F[^oF$%2desolv/gintorder1GR6$F[eo
F\dbl6&%*psvar_facG%*fsvar_facGFcen%(fac_eqnGF^pF$C'@%4Feen>FftF_bq>FftFP>FKFhv
>F9Fhv?&FAFftF\u@%554-FdbmFd\n/FQFA33Fgdq/FcctFQF\`_l>FKF]eq>F9Fht7$-Fcr6$F`hmF
genF`gnF$F$F$Fe^tRFefbl6'Fi[mF,Fj[mFfjoF0F^pF$@)Fa\_lFPF`_pC'>F9F_`m>FKFf^p>Fft
Ff]p>FA7#-F<6$&F96#Fdax/F\x;F^oFbgp-F\_l6%FKF\\\mFh[qF_atC'>FKFP>FAFhv?(F$F^oF^
oF$F]hhlC$>FA6$FAF^jo>FKF`[l>FKF^`p-F\_l6%FKF\\\mFAFj]mF$F$F$%;desolv/OrderVars
OfFunctionGRFhgl6(F[hlF\hlF]hlF0F^hl%+SubWeightsGF1F$C,@$FchlFghl>FKFfhl>F\xFjh
l>Fjw&FPFfdr@$5Fbjal/F``oF6C$>F\x-FM6$RF_ilF$FXF$FailF$F$F$F\x>Fjw-FM6$RF_ilF$F
XF$FailF$F$F$Fjw>FAFgil>F97#-F<6$,(*&F_jlF^o)Fajl,&FDF^oFAF^oF^oF^o*&&FjwF`blF^
oF`jlF^oF^oFftF^oFbjl>F9Fejl>F9FhjlF[[mF$F`\blF$F]gclRFd_`lFg[\mF^pF$C$>FKF\^y@
%Fd]o&FebnFahrFebnF$F$F$%6desolv/SortDerivativeGRF_gq6(%+HigherListG%*LowerList
G%)FirstDerGFcz%(DerListG%%backGF^pF$C,@%F[hgl>FjwF]hgl>FjwF^o@$FghsFghl>FAFebn
>F\xFc`m>F9Fhv>FKFhv?&FftF\xF\u@%2Fjs-Ff[l6%Fh[qF`\qFQ>F96$F9Fft>FKF_ap>FK-F]`^
m6$FbgnFQ>F9-F]`^m6$F`hmFQ@%/FjwF^o7%Fa\oFAF`\o7%F`\oFAFa\oF$F$F$Fab\lRF_jv6#Fj
a^lF^pF$C&@$FcyFj`s>FK-FcgnFeg\m>FK-FM6$R6#%"VGF$FXF$-%-anglebracketG6$Febn-Fep
6$FipFgpF$F$F$FKFbg`lF$F$F$FiiblRFd_`l6$F^_pFa\pF^pF$@'FcgblC&>FKFjs>F9FP?(F$F^
oF^oF$-FdbmFef\mC$>FKFc^l>F9FgfnFKF`_p@'F[[_l-FdsFh]pFbepFeepF^oFjsF$F$F$F\foR6
%Fe\lFb[mF_bm6(FeapFa\pFf\l%)ord_derxGF^_p%(var_posGF^pF$C+@$FeewFe[]l>FK-Fgs6%
FdbmF[fsFfbm>FAF`[l>FKFdfm-F\]m6%FfoF]bflFj]gl>F\x-Fcr6$Fh`pFh_m>Fft-Fep6$FjwFh
`p?&F9FKF\u@$52Fft-Fep6$FjwF`]p3/Ffe^mFft2Fft-Fcr6$F`]pFh_mC$>F\xF[f^m>FftFfe^m
FAF$F$F$FdoR6%F*Fg`]mF,64F0FajtFhhm%%loopGF`iyF^eoFe\lFb[mF+Fchy%'hcoeffGF\clFj
\l%'lanameG%+TOrderListGF_zF`hs%,ProcessTimeGF1F$C6>Fh[zF``v-Fb_`lFfem@$2Ff\nFi
pFghl>F_uFfa_l>F9-FM6%F^alFPF_u>FibyFF@'/F``oFa`oC'>FibyFc`o>F9-%4desolv/Rearra
ngeAllG6$F9Fiby>FhgrFiby>Fiby-FM6$RFhanF$FXF$F\[mF$F$F$Fiby>F9-F[`dl6&F9F_uFfoF
iby/F``oFjsC(>Fiby-FM6$RFhanF$FXF$@%Fj]t7$Fj]mF_`m7$FPFFF$F$F$F_u>Fiby-%Kdesolv
/Order1Function2VariableByAutoWeightGFig^m>F9Fgg^m>FhgrFiby>Fiby-FM6$RFhanF$FXF
$F\[mF$F$F$Fiby>F9F`h^mC,>Fiby-FM6$RFhanF$FXF$@%Fj]tFj]mFPF$F$F$F_u>F9F`h^m>Fib
y-FM6$RFhanF$FXF$@%Fj]tFih^mFjh^mF$F$F$F_u>FibyF\i^m>Fiby-%?desolv/Order1Functi
on2VariableG6%F9F_uFiby>F9Fgg^m>FhgrFiby>Fiby-FM6$RFhanF$FXF$F\[mF$F$F$Fiby>F9-
FM6$RFVF$FXF$FjhlF$F$F$F9>F9F`h^m@$F^]lC$-Fg]l6#%DOrdering~of~functions~and~var
iablesG-Fg]lF`hr?&FcuF9F\uC$>Fey-F_`il6$Fjal)&FhclFb_l&FhclFe_l@%-Fdbm6$FeyFedq
>Fgx<$FeyF\^`l-F_goFg^z>FftFjs>FjwF\u?(F$F^oF^oF$3Fjw2FftF[_oC*>FjwF%>FftF`\w@$
F^]lC%-Fg]l6#-F\_l6$%-Iteration~:~GFft>F\x7#-F<6$-Fhal6#Fh^p/FKF^jil@%/F\xFF-Fg
]l6#-F\_l6%%0~#~equations~:~GFi_n%5~;~lengths~:~[0,0,0]G-Fg]l6#-F\_l6+Fa^_mFi_n
%/~;~lengths~:~[G-F_^uFidw%",G-FdcoFidwFi^_mF`cgl%"]G>F9-F^\lFah^m@$2Fjs-FdsF^c
lC%>Fadr-FM6$F\oFhbl>Fadr-FM6$FfdoFadr?(FKF^oF^oFhdnF\u@$-Fdbm6$Fh^pFadrC%>Ff\n
Fg\n>Fcu-Fev6$/F^oFf\nFh^p>F9-Fev6$/FKFcuF9>FA-FbhmFah^m@$F_ctC$>F97$Fh[qF`\o>F
jwF\u@$F___mC'>FadrFc__m>FadrFf__m?(FKF^oF^oFhdnF\u@$-Fdbm6$F9FadrC%>Ff\nFg\n>F
cuF_`_m>F9Fc`_m>F9-%?desolv/SubstituteIntroFunctionGFI>FhblFF@$F^]l-Fg]l6#-F\_l
6&F^_lFhdnFagl,&-FhalFIF^oFdglF^o@$F___m>F9Fja_m@$F[glC$-Fg]l6#%1The~final~syst
emG-Fg]l6#-Fe`]l6#-FM6$RFVF$FXF$FjhlF$F$F$F9@%/Fc^oFjs>F9-FM6$RFVF$FXF$FjhlF$F$
F$F9C%>F\by-F]cflFah^m>F9-FM6$RFVF$FXF$FjhlF$F$F$F9>F97$F9F\by>Fh[z,&F``vF^oFh[
zFao@$F[gl-Fg]l6#-F\_l6$%<Total~computational~time~:~GFh[z>Fiby-FM6$RFhanF$FXF$
Fih^mF$F$F$F_uFgg^mF$6(Ff\nFc^oFhblFfalF[_oFc`oF$FjeilR6#%&ConstG6'Fg[m%(String
1G%)LIntegerGFjbilF0FejqF$@%Fi^hlC$>FftFhev@%32FipF_bl/-%*substringG6$Fft;F^oFi
pQ#_CF$C&>FKF\u>FA<,Q"0F$Q"1F$Q"2F$Q"4F$Q"3F$Q"5F$Q"6F$Q"7F$Q"8F$Q"9F$?(F\xF6F^
oF_blFKC$>F9-F\f_mF\c\l@$4-FdbmFeho>FKF%FKF%F%F$F$F$%0desolv/subsDer2GRFbfpFdfp
F^pF$C+>%(des_resGFjs>Fft-Fgs6$RFhanF$FXF$FaaglF$F$F$F[fs>FftFifp>F\xF]gp?&FAFf
tF\uC&>FK-Fap6%RF[doF$FXF$F^doF$F$F$Fh`pFgjn@$-Fdbm6$-FM6$RF_ilF$FXF$FdhpF$F$F$
FKF\uF[al>FKFhhp@%Few>Fft-F_x6$/FAFQFft>Fft-F_x6$/FAFgipFft>Fft-Fgs6%FgnFftFin>
Fjg_m-F_x6$FftFP>F9-Facl6#Fjg_m@%Fh[[lFhdhlF9F$6$F\^zFjg_mF$%7desolv/SortDerFun
ctionGRFeddl6'FgddlFa\uF0FhddlFi[alF^pF$C+>F\xFP>FKF[u@$FfjtFghl>FftF]edl>F9-FM
6$F_bnF\x>F97#-F<6$,&*&FcgdlF^oFftF^oF^oFAF^oFgedl>F9Fejl>F9-FM6%RF`[oF$F$F$F^f
dlF$F$F$F9Fft@%F`fdl7#-F<6$&F\xFhfdlF\gdl7#-F<6$&F\xFbgdlFgedlF$F$F$%/desolv/de
ptermGR6$F[eoFb[mF$F^pF$@'F\[r-Fgs6%FdbmFPFfbmFcbmFPFjsF$F$F$F]iuR6$%%SolnGF\hu
6#%&LSubsGF^pF$C$>FK-FM6$RF_ilF$FXF$FcyF$F$F$FQ-Facl6#FahjlF$F$F$Fhg^mR6$F]q%)L
istFuncG6&%/NewDerivativesGF]pFcgq%/ListDerivativeGF1F$C'>FftFaan>FKFhv?&F9FQF\
uC%>FA-Fgs6%FdbmFftF_`s>FA-FM6%R6$%#f1G%#f2GF$FXF$/FP-Fc]clFfemF$F$F$FAF9>FKF\j
q>FKFbgn-FbilFb]`mF$F$F$F]j[lR6'F]q%%LInqGF_qF^[^lF\cq6(Fbhs%&TLInqG%&TLSolG%&R
ExprGFbjilFgjilFi`_lF$C%>F\x-F``_l6%FfoFhrFas>Fjw-F_jil6$F\xFas@%Fa^rC'>FKFP>F9
FQ>FAFfo?&FftFjwF\uC%>FK-F\yF_ap>F9-F\yFda^m>FA-F\yFhi^l7%FKF9FA7%FPFQFfoF$F$F$
Fa`zR6'F]qF^qF_qF`qF\cq6(%-pdesolv_eqnsG%.pdesolv_ineqsG%,pdesolv_solG%,pdesolv
_lfnG%-pdesolv_lvarG%3pdesolv_parametersGF^pF$C)>FKFb\^l>F9F[]^l>FAFa]^l>FftFhr
>F\xFas>Fjw7'FKF9FAFftF\xB6)FKF9FAFftF\xFjw%*pdesolv.mGF$F$F$Fbh`lR6%Fhg`lFifw%
(add_facG6#%(ele_numGF^pF$@'Fh[^l7$7#&FfoFb_l7#&FfoFe_l-F\]m6%Fhb`mFQF\jt7$-Fev
6$/FK,&FjimF^oFfb`mF^oFPFQ7$7$Fa_mFfb`m7$Fj^mFhb`mF$F$F$%Adesolv/OrderFunctionB
yAutoWeightGRFiialFjialF1F$C'>F9F^jal>FK7#-F<6$,&*&Fa[tF^oF^[blF^oF^oFAF^oF_[bl
>FKFj[bl>FKF\\blF^\blF$F$F$F^gilRFcj\l6)Fi[mFeap%(HighDerG%*HighPowerG%)MaxOrde
rGFg[m%'LOrderGF^pF$C.@%Fgew>FKFfbm>FKFQ>F9-Fgs6$RFeesF$FXF$FgesF$F$F$F[fs>F9-F
gs6%FdbmF9FK@$Fj[pFj`s>F9Fe]u>Fcu-FM6$F_bnF9>F\x-Fdco6#Fhel-F\]m6%F\xFcuFj`x>FA
-FepFeho>Fft-F[\lFg^^l>Fjw-Fgs6%FdbmFjgn)FAFft-Facl6#*&FjwF^oFaf`mFaoF$F$F$Fdj^
mR6%F*F+F_z6+F\iwF`_^lF0%*FuncOrderGF.Fff^mFd_^lF`iw%&LHDerGF1F$C*>F9Fhv>FKFhv?
&FftFfoF\uC%>FA-FdsF_\y>F96$F97%F`\qF`byF`by>FK6$FK7%F^hmFjsFjs>F9F`hm>FKFbgn?&
FftFfoF\uC.>Fjw-F_x6$FhgmFfo>Fjw7$FhfrFft>Fjw-FM6$RFhanF$FXF$F\[mF$F$F$Fjw>F\x-
FM6$RFVF$FXF$FjhlF$F$F$FP>F\x-F[`dl6&F\xFQFFFjw>F\x-F^\lFjh`m>F_u-FM6$RFVF$FXF$
&FfhlFb_lF$F$F$F\x>F\x-FM6$RFVF$FXF$FjhlF$F$F$F\x>Fiu-Fg]dlFdcjl>FK-F]j[m6%FKFi
uF^hm>Fcu-%<desolv/WeighVarsOfFunctionsG6$FfoF\x>F9-F]_^l6%F9FcuF^hm>Fjw-FgialF
_gn-FM6$Fe]^mFjwF$F$F$FcbwRF`en6(Fg[mFeap%&LDderG%(new_eqnGFa\pF]atF^pF$C'>Fjw-
FM6$RFhanF$FXF$@%Fj]tFj]mFPF$F$F$FQ>FftFjgn>Fft-FOF_al>F9-Fgs6$RFhanF$FXF$Faagl
F$F$F$Fj\x@%Fj[p-F`io6#-F\_l6&%.The~equation~GFP%C~does~not~contain~any~functio
n(s)~GFj^mC(>FAFgfn>F9-Fev6$Fe`mF9?&F\xF9F\u@$2-Fcfw6%FAF\xFjwFjs>FAF\x>FK-Fiir
Fcgt@$Fi`l-F`io6#-F\_l6&F\\amFP%N~cannot~be~solved~for~the~highest~derivative~G
F[h[m-Fho6#/FA-FaclFioF$F$F$Fi^\lRFad]lFehjlF^pF$C*>FKF[u@$FfjtFghl>F\x-FM6$RFh
anF$FXF$@%Fj]tF[uFjsF$F$F$FP>FftF]edl>F9Fbijl>F9Fejl>F9-FM6%RF`[oF$F$F$F^fdlF$F
$F$F9Fft@%F`fdlF^gdlFdfdlF$F$F$F]i^mRF_jt6(F\iwF`_^lF0Fhf`mFff^mF[anF1F$C,>F9Fh
v>FKFhv?&FftFQF\uC%>FAF`g`m>F96$F97$F`\qF`by>FK6$FK7$F^hmFjs>F9F`hm>FKFbgn>FjwF
aan>FK-FihwF_]]m>F9-F\fglFc]fl>F\x-Fcc`mF_gn-FM6$FfglF\xF$F$F$FORF`en6(F]at%'ln
ameAGF\fcl%'lfindAG%%larbG%'ConEqnGF^pF$C3@%F[hgl>F\xF]hgl>F\xFF>FK-Fgs6%FgnFQF
is>FK-FM6$RFhanF$FXF$@%Fj]tFj]mFPF$F$F$FK>FAF[fw>FA-FM6$RF_ilF$FXF$/FP-FexF]oF$
F$F$FA>FA-Fgs6%RF`[oF$FXF$-F\]m6$-Fep6$F^oF_oFQF$F$F$FAFK>Fjw-F_x6$FAFP>FK-Fgs6
%FgnFQFc`p>FA-Fgs6%Fdbm-FibpF]oFK>FA-FM6$RF_ilF$FXF$FfaamF$F$F$FA>FA-Fgs6%RF`[o
F$FXF$F\bamF$F$F$FAFK>Fjw-F_x6$FAFjw>F9-FM6$RFhanF$FXF$@%Fj]tFj]mFPF$F$F$F\x>Ff
t-Fgs6%FdbmFibamF9>Fft-FM6$RF_ilF$FXF$FfaamF$F$F$Fft>Fft-Fgs6%RF`[oF$FXF$F\bamF
$F$F$FftF9-Fbil6#-F_xFbdtF$F$F$%.desolv/findscGRF_jv6&F+FhapF_d\lF[goF^pF$C$>FK
-Fgs6$RFeesF$FXF$FgesF$F$F$F[fs?&F9FKF\uC%>FA-%&coeffGF\\l@$5F^ioFb\tF[alFfdqF$
F$F$FhoRF_jv6&%%DeqnGFeap%&larbfGFg[mF^pF$C&>F9F[fs>FA-Fgs6$RFhanF$FXF$-F[bn6#/
Fj]mFhbclF$F$F$F9>F9-Fgs6$RFhanF$FXF$FaaglF$F$F$F9@%Fj[pFPC&>F9-FM6$RF`\pF$F$F$
/FPFgc]lF$F$F$F9>FftF[cn@$F_ctC$>FA-FM6$RF`\pF$F$F$/FP-%0desolv/Der2mapAGF]oF$F
$F$FA>Fft-F_xFhi^lFftF$F$F$%)printDerGRFci]l6#%%lDerGF^pF$C&>FK-Fgs6$RFeesF$FXF
$FgesF$F$F$F[fs>FK-FM6$RFVFg[\mF$F$C$>FK-Fap6%RF[erF$F$F$F^aflF$F$F$FdpFgp@%/FK
FFFdbn/FP&Febn6#-F\_lFahrF$F$F$FK-Fg]lFb]`mFhvF$F$F$Fj\alRFd_`lFg[\mF^pF$C$>FKF
\^y@%FewFebnFj[\mF$F$F$%/restartpdesolvGRFbhy6,FajtFi]vFg[mFj]vF[^vFhhyFihyF_^v
F`^vF0F^pF$C0>Fiu7#-F<6$&F?Fd`v/Fgx;F[czF7>F_uF``v-Fc`vF_v>Fjw-FM6&FOFfoFhrFF>F
K-FM6&FOFPFhrFF>FK-FM6$RFVF$FXF$@%FfnFgavFPF$F$F$FK>FKF_hn>Fcu-FM6&FOFQFhrFF>F\
x7#7&FKFcuFjwFhr>FjwFhv?(F$F^oF^oF$2FjsFi_nC'@$F]fq-Fg]l6#-F\_l6$FjbvFi_n>FftFe
dw>F\x-Fev6$Fe`mF\x>FA-Facv6$F`\qFas@'/FAFj`rF[al/Fd^pF^o>Fjw6$FjwF\hm>F\x7$Fh[
qF__t>Fjw-Fho6#F]]t@$F]fqC$F_dv-Fg]l6$Fedv,&F``vF^oF_uFao@%/Fh]rF^oFhfrFjwF$F$F
$%/desolv/subleqnGR6)F*%$ofnGF`ivFcfpFajt%%nofnGFfiv6(F0Fi[mFghmFe\l%&tlhilGF.F
^pF$C->F9-Fep6$F^o-F\oF`cr>FAFgcr>F\xFfo>FjwFP>F\x-FM6%R6$%#eqGFi[mF$FXF$@%Fcc^
lFjsF^oF$F$F$FPF9>F\x-Fap6%RF`[oF$FXF$@%/FQFjsFjsFPF$F$F$FfoF\x>Fjw-F\yFgr?(FKF
AFaoF^oF\u@$/&FjwFioFjsC$>Fjw-Fev6$Fi`lFjw>F\x-Fev6$Fi`lF\x>FasFjw>9)-FM6&R6%Fe
\lFa[xFianF$F$F$@%FicdlFjgzFPF$F$F$FQFhrF9>9*F\xF$F$F$FbalR6$Fe\lF_z6)F[goFeapF
h\hlFa\p%'HpowerGF[cilF^eoF1F$C)@%F[hgl>FcuF]hgl>FcuFjs>FKFjgn>F9-Fgs6$RFeesF$F
XF$FgesF$F$F$F[fs>FAFgfn?&FftF9F\u@$2Fjs-Ff[l6'FftF^oFAF^oFQ>FAFft>F\x-F[\lF]eq
@%/FcuFjsC'>Fjw-F_`il6$FK)FAF\x@%F\hn>F9Fghal>F9Fj`q>F9-FM6$RFVF$FXF$@%FeenFa_m
FPF$F$F$F9>F9-Fcr6$-Fgs6$RFeesF$FXF$FgesF$F$F$F9Fcan7&FAF\xFjwF9F[[`lF$F$F$%)sy
mmetryGR6%F]q%%DVarG%%IVarG6'%'DeteqsG%'ICEqnsG%&LArgsG%%TArgGF06#%QCopyright~(
c)~1999~by~Khai~Vu~and~John~CarminatiGF$C'>FA7#-F<6$&F?Fj_n/F\xFC>FK-Fc\dlFccn>
Fft-Fgs6%FdbmFA%%mordG@$33Fagp-Fgn6$F^hmFin/-F`o6#F^hmF^oC&>Fft-Fgs6%FdbmFAFi_o
@$Fahx>FA7$Fh[q/Fi_oFi`o>F9-F'6$Fa\oFh[q>FK-Fev6$/F^oF9FK-FbfqFbfbmF$F$F$%0deso
lv/exactintGRF]bt6/F^eoF_btF`btFe\lFb[mFeapFabtFbbtFg[mFcbtFdbtFab`lFbb`lF^pF$C
+>FKF%>Fjt-F\b`l6%FPFQFhr>Fhjr-Fep6$FipFjt>FjtFgf]l>F_uFhv>FAFeem>FAF[ct?(F$F^o
F^oF$F^ctC'>F\xFcct>Fjw-F\fo6%FjtF\xFhr>F9-F[bm6&FjtF\xFjwFhr?(F$F^oF^oF$F\dtC%
>Fcu-Fadt6$FjtFjw@%FddtC)>FiuF\et>FcuF`et>FcuFget>F_uF[ft>Fft-Fa]o6#,&FjtF^oF`f
tFao>Fgx-Fhdt6%FjtF\xFQ@%Ffft>FKF\u>FjwFgx>FKF\u>F9F^hbm@)4-Fdbm6$FjtF_gt>FA-Ff
r6$FjtFAFdgtC$>FAFggt>FKF%Fjgt>FKF\u>FKF%@%F_htC'>Fhjr-F]et6$FhjrFhr>Fey-F]et6$
-Ffht6%FjtFQFhrFhr>F_u6%F_uFeyFhjr>Fjt,&FjtF^oF\itFao7$F_it-Fa]oFh\sFjsF$F$F$Fe
dtR6&FcjoFdjoF_bmF+F$F^pF$@)FcyF\uFc_bmF%2-Fcfw6%-F]etF__rFPFhrFjsF%F\uF$F$F$Fd
g[lR6(%%snumG%(eqn_numG%%solnGF*F[^bmF+6'F`iyFaiyFbiyFciyFeiyF^pF$@$Fi[\lC(-Fg]
l6#-F\_l6'%*Module~(~GFP%:~)~applies~to~equation~(~GFQ%#~)G@$Fe[^lC%-Fg]l6#-Fho
6#Ffb`m-Fg]l6$%'gives~G-Fho6#Fhb`m@$2FDF`fqC)-Fg]l6#%6The~current~system~:~G-Fg
]l6#-FhoF`cr-Fg]l6#%:The~current~unknown(s)~:~G-Fg]l6#Fc`bm-Fg]l6#%;The~current
~solution(s)~:~G-Fg]l6#-Fho6#Fas@$32Fa`oF`fq2F[czF7-Fg]l6#-F\_l6$%.history~list
:G&F?Fg[\l@%2Fjs-FdsF`crC&>FK-FM6$FhalFhr>FAFhbbl>F9F\_fl>FftFbg`lC%>FAFjs>F9Fj
s>FftFjs>F\x-Fds6#-Ffr6$FhrFc`bm-Fg]l6#-F\_l6&Fe_flFe_cm%5,~~No.~unknown(s)~:~G
F\x@$Fa]z-Fg]l6#-F\_l6(F\`flFftF]`flFAF^`flF9F$F$F$FfhblRFcdm65F^zFg]vF]gxF*F]e
oF[goF]d\lF`btFj]^lFicelFgcelFh`elFg^clFfa^lF[b^lFhcelF[_sF\dblFg[mF^pF$C(>FKFj
^cl>F9F\_cl>FAF^_cl>F_u-Fb_cl6%FKF9FA>FhgrFhv@'Fi`sFj`r/Fb_wFjsC(>F[^[l-FadmFha
cm>Fhjr-FM6$RFVF$FXF$F\[mF$F$F$F[^[l>F[^[l-FM6$RFVF$FXF$FfhlF$F$F$F[^[l@$Ff`clC
)F^aclF_aclF`\^lFbacl-Fg]lF_\[mFfacl-Fg]l6#-Fho6#F[^[l>Fd`al<$F_bcl-Fcr6$FhjrFi
r7%F[^[lFFF^o@%/Fb_wF^oC$?&Fh[zF^]wF\uC$>F[^[l-Faeel6&FKF9FAFh[z@$0F[^[lFj`r-F]
[l6#7%&F[^[lFb_l7$Fj^m-Fep6#&F[^[lFe_lFipFejxC->FftFF>FgxFj`r>F_uF^]w>FjwFK?(F$
F^oF^oF$Fcb^lC.>F\xF^]w>F_u-Fev6$Fe`mF_u>Fadr-Fgs6%FdbmF9F\x>Fadr-Ffr6$FjwFadr>
FgxFj`r>FeyFj`r>Fgx-Fca^l6&FjwF9FAFadr@$/FgxFj`rFejx>Fhgr6$Fhgr-Fep6#Ff]`l>Fjw&
FgxFb_l>Fey-Fb_cl6%Fjw-FfrFedalFA@$5FbgelF]i_l>F_uFF@%/-FbgoFa]mF^o>FiuFFC%>Fjt
F\gcm>Fjw-Facl6#-FigoFa]m>FiuFffcm@%/FiuFF>FibyFA>Fiby-Fgs6%RFeielF$F$F$Fc^^lF$
F$F$FAFiu@%330FeyFj`rF]i_l0FjwFjs>Fgx-F]cn6%FjwF9Fiby>FgxF]]t@$-F`]lFh\s>Fgx-FM
6%RFejelF$F$F$Fff^lF$F$F$FgxFjt@$Ff`clC'FijelF\[flF_[flF`[fl-Fg]lFd`v7%FgxFFF6F
$Fc`alF$F\aelR6$F\goF]eo6#F_dblF^pF$@'FeenC$>FKFi\bl7$FKF\]blFcc^lF[]qF\fmF$F$F
$Fb_clRFcdm6+F]gxFg]v%&lnargG%(numfuncG%%lnumGF\ablFg[mF[\mF0F^pF$C,>FKFegx>F9F
eem>FA-FM6$RFhanF$FXF$@%Fj]tF[uFjsF$F$F$F9@$/FfwFb^uFejx>FftRFja]mF$F$F$-Fcr6$-
FM6%RF`c^lF$FXF$@%Fcc^lF^oFjsF$F$F$FQFPFh_mF$F$F$>F\x-FM6%FftFKF9>FjwFh^_m>FcuF
hv?(F_uF^oF^oFi_nF\u@$/&F\xFjcwFjw>Fcu6$FcuF\`hl7$FfdsFjwF$F$F$Fc`sRF`\p6$Fj[mF
,F^pF$@%/F]^jlFjsFebnC%>F9Fdp>FK-Fap6%RF\aflF$F$F$F^aflF$F$F$F9Fgp@%Few-Febn6#F
`\o-Ff^m6$Fd\dmFa\oF$F$F$%-desolv/infinGRFa[m6)F0Ff[mFg[mFh[mFi[mF,Fj[mF^pF$C$>
FjwF`\m@'Ff\m@%Fbhw-F_bhlFa]m-F\drFa]mFg]mC&>F\xFj]m>FcuF_`m>FA7#-F<6$*&Fe^mF^o
-Ff^m6$-&F[cr6#F\_mFa]mFa_mF^oFb_m,&-Ff^m6$-&F]drFj_nFa]mFa_mF^oFga\lFaoC'>F\xF
j]m>FcuF_`m>Fft-Fh\dmFb`m>FA7#-F<6$*&F[amF^o-Ff^m6$F[^dmFdamF^oFb_m,&FhamF^oFf_
mFaoF$F$F$F^dxR6$F\bxFb[mF$F^pF$@'F\[r-FMF]dxFg[clFfb]m*&F\[mF^o-Ff^m6$FQF_\clF
^oF$F$F$%(findunkGRF`en6'Feap%&lnderG%&ltermGF\goF]atF^pF$C(>F\x-FM6$RFhanF$F$F
$@%Fj]tFj]mFPF$F$F$FQ>FKFhv>F9Fhv>FAF_bq?(F$F^oF^oF$FejglC%>FftFcct>FAFggt@*-F\
]mF_al>F9Fda^mFcabl@$3-Fdbm6$Fjc[m<$Ff^mFb_p-FdbmF\c\l>FKF_ap5F]bglF\\q>FA<$Fh[
qF`\q3Fdi^l-Fgn6$Ff_uFi]o>FA<$Fh[qFegp7$-Fcr6$Fj`qFcan-Fcr6$F`aqFcanF$F$F$F]j`m
R6$F_zF]q6'%/AllDerivativesGF]pF`fglFafgl%,ListWeightsGF1F$C&>FK-FM6$RFVF$FXF$-
Fep6#FaanF$F$F$FQ>F\xFhv?&F9FPF\uC&>FA-Fgs6%FdbmFKF_`s>Fft-F^]y6$Fg`sFA>F9-Fev6
$F\f`lF9>F\x6$F\xF9Fa_nF$F$F$F\yRFbfpFdfpF^pF$C+>Fjg_mFjs>Fft-Fgs6$RFhanF$FXF$F
aaglF$F$F$F[fs>FftFifp>F\xF]gp?&FAFftF\uC&>FK-Fap6%RF[doF$FXF$F`doF$F$F$Fh`pFgj
n@$-Fdbm6$-FM6$RF_ilF$FXF$FdhpF$F$F$FKF\uF[al>FKFhhp@%Few>FftF`i_m>FftFdi_m>Fft
Fhi_m>Fjg_mF[j_m>F9F^j_m@%Fh[[lFhdhlF9F$Faj_mF$%/desolv/subsfunGR6$F[eoFcfp6(F,
Fj[mFg[mF0Fa\pFeapF^pF$@%4-Fgn6$F]gpFis-Fa]o6#-Facl6#-F_xFizC)>FjwF[fs>Fjw-Fgs6
%FdbmFjw-Fep6$FjsF]gp>Fjw-Fat6$Fjw-Fgs6%FdbmFjwF]gp>Fjw-Fgs6$RFhanF$F$F$F`_pF$F
$F$Fjw>FAFP?&F\xFjwF\uC&>FK7#F__t>F97#-Fep6#-Fep6$FjsF_gt>F9Fiiil>FA-F_x6$/F\x-
Ff^m6$FfiuF9FA-Fa]o6#-Facl6#-F_x6$FQFAF$F$F$F[b_mRFacal6'%&NLeqnG%&LNameG%%SEqn
GFe\lF0F1F$C'>F9Fc__m>F9-FM6$FfdoF9>FKFP?(F\xF^oF^oFfwF\uC$>Fft&FKFj_n@$-FdbmFf
foC&?&FAFhblF\u>Fft-Fa]o6#Fjgam>Ff\nFg\n>Fft-Fev6$Fa`_mFft>FK-Fev6$/F\xFftFKFKF
$6#Ff\nF$%(comtab2GR6$%%lvecGF,63F0Fc\l%)gamcountG%%commGFe\lF[go%&leqnsGFdcxFh
_yF^\cm%'numvecG%%ceqnG%(llambdaGFj_y%'lgammaG%*commatrixG%(con2mapGF^pF$C0>F\b
yR6$%$vecGF,F$F$F$@%F\[r-FM6%RF`[oF$FXF$FiehlF$F$F$FPFQFiehlF$F$F$>FeyF[u@$Fd`y
C$F[ay?(FKF^oF^oF[uF\u-Fg]l6#/&F_uFio-FepFbhjl>Fhjr7#-F<6$Fa_l/FK;F^oFey>Fjt-Fc
r6$7#-F<6$*&Fa_lF^oFi\emF^oF_]emFh_m>Fiby7#-F<6$7#-Fcp6$FjsFeyF_]em>Fiby-Fcr6$F
ibyFfby>FhgrFhv>FAFjs>Fjw-Fcr6$-Fap6%RF[doF$FXF$FccyF$F$F$Fhjr-FM6%F\byFPFQFh_m
?(FKF^oF^o,&FeyF^oFaoF^oF\u?(F9Fc^lF^oFeyF\uC$>Fft-FidhlFady@$FcdyC'>F\x-Fa]o6#
-FaclFgdy>Fcu-Fe[fl6%F\xFFFQ@$FdeyF[al>Fgx-Fiir6$FcuF`ccm@%/FgxFhvC'>FAFhin>Fft
-F]g\mF`bl>Fhgr6$Fhgr/&FadrF@Fft>&FibyF_gnFh`em>&FibyFbt,$Fh`emFaoC$>Fj`em-F_x6
$FgxFjt>F\aem,$Fj`emFao>Fhgr7#Fhgr@$Fd`yC%-Fg]l6#%9The~commutator~table~is:G-Fg
]lFbfy@$2FjsFAC$FhgyFf[_m@%F]bem7$-F`hyFbfyFhgrFabemF$F$F$FhbwR6%Fe\lFc[mFd[m6)
F+F0F\fclF_fclF`fclFf`wFg[mF^pF$C->FK-Fgs6%Fa^tF[fsFfo>FK-Fbfv6$F_hnFao>FA-FM6$
F^\pFK>Fft-F`ap6'FPFKFAFQFfo>F\x-F<6$7$-F`hcl6%&FQFIFQFfoFhcem/F9Fc_m>F\x6$F\x-
F<6$7$-F`hcl6%&FfoFIFQFfoFadem/F9Ffcr>F\x7$F\x-F<6$7$-F`hcl6%&FAFIFQFfo-Fcr6$Fj
demFcbp/F9Fc^p>Fcu-FM6%RF`jclF$F$F$F]g`lF$F$F$F\xFft>Fcu-Fcr6$FcuFh_m>Fcu-Facl6
#-Fbil6#-F_x6$/F`hclFh\dmFcu>Fcu-Fjdv6'FcuFKFAFQFfoF$F$F$FgbxRF`ar6.Ff[mFe\lF.%
(fmatrixGF0Fc\lF+Fa\p%(ope_eleG%%cof_G%(subexpnGFj[mF^pF$C)@%-Fgn6$F\[mFin>FA-F
M6$RFVF$F$F$FjnF$F$F$FP>FAFP?&FKFQF\uC$>Fey/-FKFdcr,&FcgemF^o*&%%eps_GF^o-(%"QG
FKFdcrF^oF^o>FA-Fgedm6$FAFey>FA-Ff^m6$FAFfgem>FA-F_x6$/FfgemFjsFA>Fft-%&arrayG6
$Fde\lFc_m?(F\xF^oF^oFd^pF\uC$>F9Fa^p?(FjwF^oF^oFd_mF\uC(>FcuF`]gl>Fcu-Fgs6%Fdb
mFcu(Figem&FQFa]m>Fcu-Fbfv6$-Fcr6$FcuFcanFao>F_uFhv?&FiuFcuFa\xC'>Fgx-Fdeam6%F9
FiuF^o@$/FgxFjsF[al>F9,&F9F^o*&FgxF^oFiuF^oFao>Fjt-FcfblF_v@%/FjtF^o>F_u6$F_u7$
FgxF^o>F_u6$F_u7$FgxFihr>&Fft6$F\xFjwF_it-F`hyF`blF$F$F$Fa]clR6$F[eoF+6'Fg[mF[g
oFa\pFdfemFj[mF^pF$C'>FKFhv>F9FP?&FAFQFa\xC&>Fft-Fdeam6%F9FAF^o>F9,&F9F^o*&FftF
^oFAF^oFao>F\x-FcfblF@@%F^j^l>FK6$FK7$FftF^o>FK6$FK7$Fft7#F`\fm@$Fa\x>FK6$FK7$F
9F^oFbgalF$F$F$FbanR6#FgjblF$F^pF$-Fgs6$RFhanF$FXF$FaaglF$F$F$F[fsF$F$F$Fb[jlR6
%%%Vec1G%%Vec2GF\cq6'%'LComp1G%'LComp2G%'TComp1G%'TComp2GFicelFi`_lF$C%@%F\[r>F
KF_`m>FKF\_^l@%Fafhl>F97#Fj^m>F97#FQ@%/FfwFhdnC$?&FftF9F\uC$>F\x-FM6&F`h\mFKFft
Ffo@$4-Fdbm6$F\xF\uFe[]lF\uF%F$F$F$F_bwRF_[xF$F^pF$-F\b]l6%FP-FM6$RFVF$FXF$F_`m
F$F$F$FQFfoF$F$F$F[`sR6$F,F+F$F^pF$@%Fh[^lC&@$4-F`]l6#F[hdl>F[hdlFjs>F[hdl,&F[h
dlF^oF^oF^o?(F$F^oF^oF$-Fdbm6$FQ-F\_l6$%#C_GF[hdl>F[hdlF``fm-Fedq6%Fd`fmFFFFC&@
$4-F`]l6#Ff[el>Ff[elFjs>Ff[el,&Ff[elF^oF^oF^o?(F$F^oF^oF$-Fdbm6$FQ-F\_l6$%#F_GF
f[el>Ff[elFaafm-Fedq6%FeafmFh^_lFPF$6$Ff[elF[hdlF$%(solvedeGRFcdmF``rF^pF$C%>FK
Fc`r>FK-Fd\vFeo@%Fi`rFj`rF]arF$F$F$F$
