M7R0
I8Multiply_Row_Col_Matrixf*6$%#X1G%#X2G6#%$auxG6"F)C$>8$-%(ROW_COLG6$-%#opG6#9$
-F16#9%-%'RETURNG6#F,F)6#F.F)F),
I0Multiply_Matrixf*F$6)F(%#aiG%#ajG%#akG%%aux1G%%aux2G%%aux3GF)F)C$?(8%"""FG%&N
BVARG%%trueGC$>8(7#-%$seqG6$&F36$FF8&/FS;FGFH?(8'FGFGFHFIC$>8)7#-FO6$&F66$FSFWF
T>&%(MATPRODG6$FFFW-F"6$FLFZ-F86#-%%evalG6#F\oF)6%F.FHF\oF)F),
I7Multiply_Matrix_Vectorf*F$F=F)F)C%>FZ7#-FO6$&F66#FF/FFFU?(FFFGFGFHFIC$>FLFM>&
%(VECPRODGF^pF^o-F86#7#-FO6$FdpF_pF)6%F.FHFepF)F),
I2Krawczyk_Analyzerf*6$%"XG%%nbeqG60F>F?F(FA%$au1G%$au2G%$au3G%"WG%"YG%"UG%$as1
G%#asG%$AU1G%$AU2GF)F)C4?(F,FGFGFHFIC)>FS&F3F9>FL-F16$FGFS>FZ-F16$""#FS@$/FLFZC
%>8*FL>&%#AYGF9F\s\>F\s-%&evalfG6#*&,&FLFGFZFGFG$Fgr""!!""@$51F\sFL1FZF\s?(F)FG
FGF)F[tC$>%'DigitsG,&FatFGFgrFG>F\s-Fcs6#*&FfsFGFgsFis>F^sF\s?(F,FGFGF6FIC$Z%>&
%#FYGF9-&%$FEQGF96#-FO6$&F_sF^pF_pF)-F86#Fhs>&%#fyGF9F\u?(F,FGFGF6FI?(FFFGFGFHF
IC%Z%>&%#YPG6$F,FF-&%$GEQGFav6#-FO6$&F_s6#%#alG/FjvFUF)Feu>FSF_v@%/-%%nopsG6#FS
FG>&%#AWGFavF_v>Fcw*&-Fcs6#,&FbrFGFerFGFGFgsFisZ%>8,-%(inverseG6#FdwF)Feu?(F,FG
FGF6FI?(FFFGFGFHFIZ%>&%#FPGFav-Fcv6#F0F)Feu>Fdw-F;6$F\xFex?(F,FGFGFHFI?(FFFGFGF
HFIC%>FS$FhsFhs@$/F,FF>FS$FGFhs>&%$AAUGFav-_%(intpakXG%2Interval_subtractG6$FSF
cw?(F,FGFGFHFIC$>8/F_y?(FFFGFGFHFIC%>FSFey@%F^wC$>8.-_Fiy%-Interval_addG6$F_z-%
$absG6#Fey>F_zFfzC)>FLFbr>FZFer@$31FLFhs1FhsFZ@%1-F\[l6#FL-F\[l6#FZ>FS7$F_y-Fcs
6#Fj[l>FS7$F_y-Fcs6#Fh[l@$2FhsFL>FS7$Fb\lF^\l@$2FZFhs>FS7$F^\lFb\l>Ffz-Fhz6$F_z
FS>F_zFfz@$2FG-F16$FgrF_zC$@$/%,ALIAS/debugGFgr-%'printfG6$Q+Kraw~s:%a|+F)F_zFe
u?(F,FGFGFHFI>&%#AUGF9-Fhy6$F`rF^s>80-Ffo6$FfyF_^l>81-Ffo6$F\xFiu?(F,FGFGFHFIC$
>F^^l-Fhy6$&Fc^lF9&Fg^lF9>&%#AKGF9-Fhz6$F^sF^^l@$Ff]l-Fi]l6$Q+Kraw~K:%a|+F)-Fco
6#Fc_l?(F,FGFGFHFI@$52-F16$FgrFb_l-F16$FGF`r2-F16$FgrF`r-F16$FGFb_l[@$0F,,&FHFG
FGFGC$@$Ff]l-Fi]l6%Q2No~inter:~%a,~%a|+F)Fb_lF`r-F86#Fis?(F,FGFGFHFI@$52Fe`lF``
l2Fg`lFb`lFi`l@$/F,F\al-F86#FGFeuF)6.FHF`uFdvF_sF]uFiuF`vFc_lFdwFexFfyF_^lF)F),
I0Krawczyk_Solverf*F^q65F>F?%$ar0G%$ar1GF(FAFbqFcqFdqFeqFfqFgqFhqFiqFjqF[r%$AY1
G%%flagG%)flag_solGF)F)C8>%,RESULT_KRAWG7"?(F,FGFGFHFIC)>FLF`r>F\s-F16$FGFL>8+-
F16$FgrFL@$/F\sFaclC%>F\xF\s>F^sF\xF`s>F\x-Fcs6#*&,&F\sFGFaclFGFGFgsFis@$51F\xF
\s1FaclF\x?(F)FGFGF)F_dlC$>FatFbt>F\x-Fcs6#*&F]dlFGFgsFis>F^sF\x>86FG>Fat"$+"?(
F,FGFGF6FIC%Z%>F\uF^uF)Feu>FhuF\u@%3/-F`w6#FhuFG2%/ALIAS/fepsilonG-F\[lFgel>F[e
lFhsC%>F\s-F16$FGFhu>Facl-F16$FgrFhu@$552Fiel-F\[l6#F\s2Fiel-F\[l6#Facl2Fiel,&F
aclFGF\sFis>F[elFhs@$/F[elFGC$?(F,FGFGFHFI>Fhbl7$-F16#FhblF^sF[bl?(F,FGFGF6FI?(
FFFGFGFHFIC%Z%>F_vFbvF)Feu>FLF_v@%/-F`wFi[lFG>FcwF_v>Fcw*&-Fcs6#,&F^clFGFbclFGF
GFgsFisZ%>FfzF]xF)Feu?(F,FGFGF6FI?(FFFGFGFHFIZ%>FdxFfxF)Feu>Fdw-F;6$FfzFex?(F,F
GFGFHFI?(FFFGFGFHFIC%>FLF_y@$Fay>FLFcy>Fey-Fhy6$FLFcw?(F,FGFGFHFIC%>Fg^lF_y?(FF
FGFGFHFIC%>FLFey@%F^hlC$>Fc^l-Fhz6$Fg^lF[[l>Fg^lFc^lC)>F\sF^cl>FaclFbcl@$31F\sF
hs1FhsFacl@%1FgflFjfl>FL7$F_y-Fcs6#Fjfl>FL7$F_y-Fcs6#Fgfl@$2FhsF\s>FL7$Fc[mF_[m
@$2FaclFhs>FL7$F_[mFc[m>Fc^l-Fhz6$Fg^lFL>Fg^lFc^l@$2FG-F16$FgrFg^lC$@$Ff]l-Fi]l
6$F[^lFg^lFeu@%/F,FG>FSFg^l@$2FerFc\m>FSFg^l?(F,FGFGFHFI>F^^lF`^l>82Fd^l>83-Ffo
6$FfzFiu?(F,FGFGFHFIC$>F^^l-Fhy6$&Fb]mF9&Fd]mF9>Fb_lFd_l@$Ff]lFg_l?(F,FGFGFHFI@
$5FhalFgalFi`l@$F[alFeu>85FG?(F)FGFGF)/Fe^mFGC.>Fe^mFhs?(F,FGFGFHFIC*>&%#APGF9F
b_l>FLF]_m>F\sF^cl>FaclFbcl@$FeclC%>F\xF\s>F^sF\xF`s>F\x-Fcs6#*&F]dlFGFgsFis@$F
_dl?(F)FGFGF)F_dlC$>FatFbt>F\x-Fcs6#*&F]dlFGFgsFis>F^sF\x>F[elFG?(F,FGFGF6FIC%Z
%>F\uF^uF)Feu>FhuF\u@%Fdel>F[elFhsC%>F\sF^fl>FaclFafl@$Fdfl>F[elFhs@$F`glC$?(F,
FGFGFHFI>FhblFdglF[bl?(F,FGFGF6FI?(FFFGFGFHFIC%Z%>F_vFbvF)Feu>FLF_v@%F^hl>FcwF_
v>Fcw*&FchlFGFgsFisZ%>84F]xF)Feu?(F,FGFGF6FI?(FFFGFGFHFIZ%>Fdx-Fcv6#-FO6$&F^_mF
ivF[wF)Feu>FdwF]il?(F,FGFGFHFI?(FFFGFGFHFIC%>FLF_y@$Fay>FLFcy>FeyFfil?(F,FGFGFH
FIC%>Fg^lF_y?(FFFGFGFHFIC%>FLFey@%F^hlC$>Fc^lFajl>Fg^lFc^lC)>F\sF^cl>FaclFbcl@$
Fhjl@%F\[m>FLF^[m>FLFb[m@$Ff[m>FLFh[m@$Fj[m>FLF\\m>Fc^lF^\m>Fg^lFc^l@$Fb\mC$@$F
f]lFg\mFeu@%Fj\m>FWFg^l@$2-F16$FgrFWFc\m>FWFg^l@$52F`emFer/-F`wFfglFhsC,>FatFbt
?(F,FGFGFHFI>F^^l-Fhy6$F]_mF^s>Fb]mFd^l>Fd]m-Ffo6$F`bmFiu?(F,FGFGFHFIC%>F^^lFj]
m>Fb_lFd_l?(FFFGFGFHFI>&FfzFav&F`bmFav>FSFW?(F,FGFGFHFI@$52Fg`l-F16$FGF]_m2-F16
$FgrF]_mF``lFi`l@$F[alFeu>Fe^mFGF`sFeuF)60FHF`uFdvFhblF^_mF_sF]uFiuF`vFc_lFdwFe
xFfyF_^lF)F),
I0Newton_Intervalf*F^q65F>F?F@FjvF(FAFBFC%%aux4GFbqFcqFdq%#abGFdblF]uFfq%#MMGFi
q%$au4GF)F)C.>FatF]el>Fg^lFhs?(F,FGFGFHFIC*>F]_mF`r>FLF`r>8-F^cl>FfzFbcl@$/Ffhm
FfzC%>F_zFfhm>F^sF_zF`s>F_z-Fcs6#*&,&FfhmFGFfzFGFGFgsFis@$51F_zFfhm1FfzF_z?(F)F
GFGF)FcimC$>FatFbt>F_z-Fcs6#*&FaimFGFgsFis>F^sF_z?(F,FGFGF6FIZ%>F\^mF^uF)Feu>%%
FDIMG7#-FO6$F\^m/F,;FGF6?(F,FGFGF6FI?(FFFGFGFHFIC&Z%>F_v-Fcv6#-FO6$&F_s6#FW/FWF
UF)Feu>FLF_v@%F^hl>FcwF_v>Fcw*&FchlFGFgsFisZ%>F_v-Fcv6#-FO6$&F3Fb[nFc[nF)FeuZ%>
Fd]mF]xF)Feu>F`bm-F;6$Fd]mF`v>Fc^l-Ffo6$Fd]mFbjm?(F,FGFGFHFIC4>FL&F`bm6$F,F,>Ff
hmF^cl>FfzFbcl@$Ff]l-Fi]l6&Q.MM[%d,%d]:%a|+F)F,F,F[]n@$31FfhmFhs1FhsFfzF`s>Fe^m
F_y?(FWFGFGFHFIC'@$/FWF,F`s>FLFe^m>&F_^lFb[n-Fhy6$Fa[nF_\n>FZ-_Fiy%/Interval_ti
mesG6$&F`bm6$F,FWF^^n>Fe^m-FhzF_o>FZ-Fhy6$Fe^mF__l>F\s-Fc^n6$FZ-_Fiy%4Interval_
reciprocalG6#F[]n>Facl-Fhz6$F^sF\s@$Ff]lC%-Fi]l6$Q8Newton~on~component~%d|+F)F,
-Fi]l6%Q.old:%a,~Z:%a|+F)F`rFacl-Fi]l6$Q&U:%a|+F)Fe^m@$52-F16$FgrFaclFb`l2Fe`l-
F16$FGFaclC$@$Ff]l-Fi]l6#Q;Box~eliminated~par~Newton|+F)Fbal>FfhmFh`n>FfzFe`n>F
_zFb`l>F[elFe`l@$52F[elFfhm2FfzF_zFbal@$52F_zFfhm2FfzF[elC'@$Fian>F_zFfhm@$Fjan
>F[elFfz>F]_m7$F_zF[el>Fg^lFG@$Ff]l-Fi]l6%Q/Improv~%a->%a|+F)F`rF]_m@$/Fg^lFG-F
86#F^_mFeuF)64%%ABOXG%&NBBOXGFH%$SOLG%%NBEQGF`uFdv%(TYPE_EQG%&NBSOLG%&TIME1G%&L
ABOXG%/ALIAS_SIMPLIFYGF_^lF^_mF_sF`vFdwFbjmF)F),
I'Taylorf*6$%"iGF_q6,F(FAFbqFcqFdqF]hm%#u1G%#u2GF>F?F)F)C(?(F\xFGFGFHFIC*>F,&F6
6#F\x>FS-F16$FGF,>FW-F16$FgrF,@$/FSFW-F86#7$%%FAILGF]en>FL-Fcs6#*&,&FWFGFSFGFGF
gsFis@$51FLFS1FWFL?(F)FGFGF)FdenC$>FatFbt>FL-Fcs6#*&FbenFGFgsFis>&F_sFadnFL>&F_
^lFadn-Fhy6$F`dn7$FLFLZ%>&F]uF2-&F`uF26#-FO6$&F_s6#Ffhm/FfhmFUF)Fjdn?(FfhmFGFGF
HFIZ%>&F`v6$F3Ffhm-&FdvFdgn6#F4F)Fjdn>F,-F"6$7#-FO6$FcgnF_gn7#-FO6$&F_^lF^gnF_g
n>FF-Fhz6$FgfnF,-F8F^pF)6*%#EQGF`uFdvFHF_sF]uF`vF_^lF)F),
I,Method3B_GRf*6#%$anbG62%$epsGF>F?F@F(FjcnF[dn%#u3GFbqFcq%+old_DigitsGFdbl%-fl
ag_successG%$afpG%$av1G%$av2GF)F)C&?(FFFGFGFHFIC->FL&F\cn6$F3FF>FZF^cl>F\sFbcl>
Facl,&F\sFGFZFis@$52%,ALIAS/Max3BGFacl2Facl%.ALIAS/Delta3BGF`s>Fccn7#-FO6$&F\cn
6$F3FWFc[n>F,Fcjn>F_zFG>Fc^lFhs?(F)FGFGF)/F_zFGC1>F_zFhs>FccnFejn>FLFhin>FZF^cl
>F\sFbcl>F\xFZ>Ffhm,&FZFGF,FG?(F)FGFGF)1FfhmF\xC$>FatFbt>FfhmFg[o>&FccnF^p7$F\x
Ffhm>FLFhs@$3/FdcnFG0%2ALIAS/Use_Simp_3BGFhs>FL-%)SimplifyG6#7#-FO6$&FccnFawFT@
$30FLFhs0FLFis?(FSFGFGFHFI>F]]o-F16$FSFL@$/FLFisC,>Fhin7$FfhmF\s>FLFhs@$Fb\o>FL
-Fh\o6#7#-FO6$&F\cn6$F3FSFT@$Fg]oFbal@$F_]oC%?(FSFGFGFHFIC$>Fc^oFd]o>F]]oFd]o>F
Z-F16$FGFhin>F\s-F16$FgrFhin>Fc^lFG@$2F\s,&FZFGF,""$Fi`l>F,,$F,Fgr>F_zFGF`s?(FW
FGFGF_cnFIC+>FfzFat>FLF\en>FL-%(Eval_GRG6$FWFj\o>Fb]mF^cl>Fd]mFbcl@$5/Fb]mF]en/
Fd]mF]enF`s@$3/&F`cnFb[nFhs52FhsFb]m2Fd]mFhsFi`l@$3/F[ao!"#1FhsFb]mFi`l@$3/F[ao
FisF]aoFi`l@$0FW,&F_cnFGFGFGC'>FhinFj]o>Fc^lFG@$Fd_oFi`l>F,Fh_o>F_zFG@$/Fc^lFG?
(FWFGFGF_cnFIC+>FfzFat>FLF\en>FL-F``o6$FWF`^o>F\xF^cl>FfhmFbcl@$5/F\xF]en/FfhmF
]enF`s@$3Fj`o52FhsF\x2FfhmFhsFbal@$3Faao1FhsF\xFbal@$3FfaoFbcoFbal>Fc^lFhs?(FFF
GFGFHFIC,>FLFhin>FZF^cl>F\sFbcl>FaclF]jn@$F_jnF`s>FccnFejn>F,Fcjn>F_zFG?(F)FGFG
F)F^[oC1>F_zFhs>FccnFejn>FLFhin>FZF^cl>F\sFbcl>F\x,&F\sFGF,Fis>FfhmF\s?(F)FGFGF
)Fi[oC$>FatFbt>F\xF\eo>F^\oF_\o>FLFhs@$Fb\o>FLFg\o@$F_]o?(FSFGFGFHFI>F]]oFd]o@$
Fg]oC,>Fhin7$FZF\x>FLFhs@$Fb\o>FLF^^o@$Fg]oFbal@$F_]oC%?(FSFGFGFHFIC$>Fc^oFd]o>
F]]oFd]o>FZF]_o>F\sF`_o>Fc^lFG@$2,&F\sFGF,!"$FZFi`l>F,Fh_o>F_zFGF`s?(FWFGFGF_cn
FIC+>FfzFat/FLF\en>FLF_`o>Fb]mF^cl>Fd]mFbcl@$Fe`oF`s@$Fi`oFi`l@$F`aoFi`l@$FeaoF
i`l@$FhaoC'>FhinF\fo>Fc^lFG@$F[goFi`l>F,Fh_o>F_zFG@$Fabo?(FWFGFGF_cnFIC+>FfzFat
>FLF\en>FLFgbo>F\xF^cl>FfhmFbcl@$F\coF`s@$F`coFbal@$FecoFbal@$FhcoFbalFeuF)6-F\
cnF]cnFHF^cnF_cnF`uF`cnFacnFbcnFccnFdcnF)F),
I)Method3Bf*Fjhn61F]inF>F?F@F(FjcnF[dnF^inFbqFcqF_inFdblF`inFbinFcinF)F)C&?(FFF
GFGFHFIC->FLFhin>FZF^cl>F\sFbcl>FaclF]jn@$F_jnF`s>FccnFejn>F,Fcjn>F_zFG>Fc^lFhs
?(F)FGFGF)F^[oC1>F_zFhs>FccnFejn>FLFhin>FZF^cl>F\sFbcl>F\xFZ>FfhmFg[o?(F)FGFGF)
Fi[oC$>FatFbt>FfhmFg[o>F^\oF_\o>FLFhs@$Fb\o>FLFg\o@$F_]o?(FSFGFGFHFI>F]]oFd]o@$
Fg]oC,>FhinFj]o>FLFhs@$Fb\o>FLF^^o@$Fg]oFbal@$F_]oC%?(FSFGFGFHFIC$>Fc^oFd]o>F]]
oFd]o>FZF]_o>F\sF`_o>Fc^lFG@$Fd_oFi`l>F,Fh_o>F_zFGF`s?(FWFGFGF_cnFIC+>FfzFat>FL
F\en>FL-&F`uFb[n6#F[]o>Fg^lF^cl>Fb]mFbcl@$5/Fg^lF]enFf`oF`s@$3Fj`o52FhsFg^l2Fb]
mFhsFi`l@$3Faao1FhsFg^lFi`l@$3FfaoFe]pFi`l@$FhaoC'>FhinFj]o>Fc^lFG@$Fd_oFi`l>F,
Fh_o>F_zFG@$Fabo?(FWFGFGF_cnFIC+>FfzFat>FLF\en>FL-F[]p6#Fa^o>F\xF^cl>FfhmFbcl@$
F\coF`s@$F`coFbal@$FecoFbal@$FhcoFbal>Fc^lFhs?(FFFGFGFHFIC,>FLFhin>FZF^cl>F\sFb
cl>FaclF]jn@$F_jnF`s>FccnFejn>F,Fcjn>F_zFG?(F)FGFGF)F^[oC1>F_zFhs>FccnFejn>FLFh
in>FZF^cl>F\sFbcl>F\xF\eo>FfhmF\s?(F)FGFGF)Fi[oC$>FatFbt>F\xF\eo>F^\oF_\o>FLFhs
@$Fb\o>FLFg\o@$F_]o?(FSFGFGFHFI>F]]oFd]o@$Fg]oC,>FhinF\fo>FLFhs@$Fb\o>FLF^^o@$F
g]oFbal@$F_]oC%?(FSFGFGFHFIC$>Fc^oFd]o>F]]oFd]o>FZF]_o>F\sF`_o>Fc^lFG@$F[goFi`l
>F,Fh_o>F_zFGF`s?(FWFGFGF_cnFIC+>FfzFatFcgo>FLFj\p>Fg^lF^cl>Fb]mFbcl@$F`]pF`s@$
Fc]pFi`l@$Fh]pFi`l@$F[^pFi`l@$FhaoC'>FhinF\fo>Fc^lFG@$F[goFi`l>F,Fh_o>F_zFG@$Fa
bo?(FWFGFGF_cnFIC+>FfzFat>FLF\en>FLFi^p>F\xF^cl>FfhmFbcl@$F\coF`s@$F`coFbal@$Fe
coFbal@$FhcoFbalFeuF)F^ioF)F),
I(Eval_GRf*Fgcn62F?Fdbl%$LOWG%%HIGHGF(F`inFA%$lowG%%highGFbqFcq%&flag1G%&flag2G
FBFdqF]hmF)F)C6>FFFG>FSFibl>FWFibl?(F,FGFGFHFIC%@$/-F`w6#&F6F9FgrC%>Ffhm-F16$FG
Fcep>Ffz-F16$FgrFcep@$2%2ALIAS/maxgradientG,&FfzFGFfhmFisFi`l>FS7$-F1FawFcep>FW
7$-F1Fb[nFcep@$F[alC%>FLF\enZ%>F\s-FifnFggnF)-F8Fi[l-F8Fhfl>FFFG>FZFhs?(F)FGFGF
)/FFFGC$>FFFhs?(F,FGFGFHFIC(>FL&FSF9@$/F^clFbclF`sZ%>FL-&Fdv6$F3F,6#FafpF)F`s@$
F^hlC%@$Fe\lC&>FFFG>FZFG>FS-%'subsopG6$/F,7$-F16$FGFfgpF[ipFSF`s@$2FLFhsC&>FFFG
>FZFG>FS-Fghp6$/F,7$-F16$FgrFfgpFgipFSF`sF`s@$2FhsF^clC&>FFFG>FZFG>FSFfhpF`s@$2
FbclFhsC%>FFFG>FZFG>FSFcip>FFFG?(F)FGFGF)F`gpC$>FFFhs?(F,FGFGFHFIC(>FL&FWF9@$Fh
gpF`sZ%>FL-F\hp6#FdfpF)F`s@$F^hlC%@$Fe\lC&>FFFG>FZFG>FW-Fghp6$/F,7$-F16$FgrF\[q
F]\qFWF`s@$F^ipC&>FFFG>FZFG>FW-Fghp6$/F,7$-F16$FGF\[qFh\qFWF`sF`s@$FjipC&>FFFG>
FZFG>FWFi[qF`s@$F`jpC%>FFFG>FZFG>FWFd\q@$/FZFGC$Z%C$>FL-FifnF^hp>F\s-FifnFa[qF)
Fjdn-F86#7$F^cl-F16$FgrF\s>FLF\en>F_zFG>Fc^lFGZ%>F\sFjfpF)>F_zFhsZ%>Fg^l-Fecn6$
F3F6F)>Fc^lFhs@$3F^[o/Fc^lFhs>FLF\s@$3/F_zFhsFabo>FLFg^l@$333F^[oFabo0-F16$FGFg
^lF]en0Fc\mF]enC*@$Ff]l-Fi]l6&QEEvaluation~normal~%d:~%a,~Taylor~%a|+F)F3F\sFg^
l>Ffhm-F16$FGF\s>FfzF`^q>Fb]mFj_q>Fd]mFc\m@$2FfhmFb]m>FfhmFb]m@$2Fd]mFfz>FfzFd]
m>FL7$FfhmFfzF[gpF)6%FdvFHF`uF)F),
I1Smear_Bissectionf*Fjhn6+F>F?FbqFcq%%amaxG%#ASGF(FA%&ajbonGF)F)C$?(F,FGFGF_cnF
IC%>FLFhs>FZ-FO6$-Fcs6#,&F`_oFGF]_oFisF_p?(FFFGFGFHFIC,Z%>F\s-Fcv6#-FO6$FhinF_p
F)Fbal@$/-F`wFhflFGF`s>FSFc`q>FWF`^q@$31FSFhs1FhsFW@%1-F\[lFaw-F\[lFb[n>F\s7$F_
y-Fcs6#Fecq>F\s7$F_y-Fcs6#Fdcq@$2FhsFS>F\s7$F\dqFhcq@$2FWFhs>F\s7$FhcqF\dq>Facl
-Fc^n6$F\s&FZF^p@$3/-F`wF[glFG2FLFaclC%>F\xFF>FLFaclF`s@$2FLFe`nC%>F\xFF>FLFe`n
F`s-F8FadnF)6$FdvF\cnF)F),
I6GradientSolve_intpakXf*6%%#eqG%$varG%%initG61F>F?F@FbqFcqFdqF]hmF[inF_qF(FA%$
AUXG%'apointG%%INITGFainF)F)CV>%)NAME_SOLGQ&_Sol_F)@$0-%'lengthG6#%)ALIAS/IDGFh
s>Fefq-%$catG6$Fefq-%(convertG6$F\gq%'stringG>Fb]m-%&fopenG6$Fefq%&WRITEG-%'fcl
oseG6#Fb]m>FdcnFhs@$0%/ALIAS/time_outGFhs>Fbcn-%%timeGF)>F_cn-F`wF2>FH-F`wF5@$0
-F`w6#9&FH-F86#Q=Each~var~should~have~a~rangeF)>%)NBBOXMAXG"%+5>F\cn-%&arrayG6$
;FGFaiqFU>F^_m-Feiq6#FU>F_^lFiiq>F_sFiiq>F]uFiiq>FiuFiiq>F`v-Feiq6$;FGF_cnFU>Fd
wF`jq>FfyF`jq>F\hmF`jq>Fc_lFiiq>Fex-Feiq6$FUFU>F\oFhjq>FepFiiq>FccnFiiq>F^cnFib
l>%%FSOLGFibl>FacnFhs>%2SOL_KRAW_INTERVALGFibl>%)SOL_KRAWGFibl>Fg^lFibl>Fat"$+#
>Fc^lFG?(F,FGFGFjhqFIC'>Ffhm&F\iqF9>FW-Fcs6#-F16$FGFfhm>FL-Fcs6#-F16$FgrFfhm@$0
FLFW>Fc^lFhs>Fg^l7$-F16#Fg^l7$FWFL>F`uFibl>F`cnFibl?(F,FGFGF_cnFIC'>FfhmF`r@$/-
%%typeG6$Ffhm%"<GFIC&>F`cn7$-F16#F`cnFbao>FfhmF`\r>F`u7$-F16#F`u-%(inapplyG6$Ff
hmF4F`s@$Fe]rC&>F`cn7$F]^rFis>FfhmF`\r>F`u7$Fb^r-Fe^r6$FfhmF6F`s>F`cn7$F]^rFhs>
F`uF]_r>FdvF`jq?(F,FGFGF_cnFI?(FFFGFGFHFIC%>Ffhm-%%diffG6$F`rF]p>FatFg[r>FcvF^_
r>FfhmFhs?(F,FGFGFHFI>Ffhm,&FfhmFG*&(%#axGF,FG(%#ayGF,FGFG>FWFibl?(F,FGFGFHFI>F
W7$FdfpFb`r?(F,FGFGFHFI>FW7$FdfpFd`r>F.-Fe^r6$FfhmFdfp?(F,FGFGFHFI>&F\cnFddnF`_
l?(F,FGFGF_cnFIC$Z%>Ffhm-F_u6#-FO6$&F\cn6$FGFFF_pF)C$>Ffhm-%(sprintfG6%QepAn~in
terval~evaluation~error~was~caught|+|+|*|*when~evaluating~equation~%d|+|+|*|*th
e~error~was:%aF)F,%*lasterrorG-F8F^gn@$/-Fg]r6$Ffhm%)intervalG%&falseGC$>Ffhm-F
`br6$Q<Cannot~evaluate~equation~%dF)F,Fdbr@$FaboC$?(F,FGFGF_cnFIC1>F_inFatZ%>Ff
hmFgarF)C$>FfhmF_brFdbr>FWF`\r>FLFe\r@$5533/&F`cnF9FhsFcdqFe\l3/FadrFisFe\l3/Fa
drFbao1FhsFLC$>FdblFG?(F)FGFGF)/FdblFGC*>FdblFhs>Fat,&FatFG"#5FGZ%>Ffhm-F_u6#-F
O6$&F\cn6$FaclFFF_pF)C$>FfhmF_brFdbr>%%au11GF`\r>%%au22GFe\r@$2FhsF[frFi`l@$2F]
frFhsFi`l@$552Fiel-F\[l6#,&FWFGF[frFis2Fiel-F\[l6#,&FLFGF]frFis2Fiel-F\[l6#,&-F
\[l6#,&FLFGFWFisFG-F\[l6#,&F]frFGF[frFisFisC%>FdblFG>FWF[fr>FLF]fr>FWF`\r>FLFe\
r>FatF_in@$Ff]l-Fi]l6%Q*eq~%d:%a|+F)F,Ffhm@$3FedrFacqC$>FacnFhs-F86#F^cn@$3Fcdr
2FhsFWC$>FacnFhsFfhr@$3F`dr5FjhrF^ipC$Ffhr>FacnFhs@$3FedrF^ip>%(NBVERIFG,&FeirF
GFGFG@$3FcdrFd[l>FeirFfir@$333F`dr1FWFhsFfdr2FcgrFielC$>FeirFfir>F_z7$-F16#F_zF
fhm@$/FeirF_cnC*>FfhmFibl?(F,FGFGFHFI>Ffhm7$-F1F^gnFbar>F^cn7$-F1FghrFfhm>F_[r7
$-F16#F_[rF_z>Facn,&FacnFGFGFG>Fb]m-Fggq6$Fefq%'APPENDG-%(fprintfG6%Fb]mQ$%a|+F
)F^cnFjgq>FaclFG>F]cnFG?(F)FGFGF)1FaclF]cnC@@$F_hq@$2*&F`hqFG$"#gFhsFG,&FbhqFGF
bcnFis-F86#7#Q)TIME~OUTF)@$2FhsFg]lC&-Fi]l6'QGCurrent~box:%d/%d,~remain~%d~(Sol
:%d)|+F)FaclF]cn,&F]cnFGFaclFisFacn-Fi]l6#Q(Box~is|+F)?(F,FGFGFHFI-Fi]l6%Q-[%10
f~%10f]~F)-F16$FG&F\cn6$FaclF,-F16$FgrF_^s-Fi]l6#Q"|+F)@$Ff]lC%Ff]s?(F,FGFGFHFI
-Fi]l6$Q$%a~F)F_^sFc^s@$/9#""%C'>FdcnFG>%(AREPEATG&9"6#F__s>Ffhm-Fh\o6#7#-FO6$F
_^s/F,FU@$/FfhmFisC$>Facl,&FaclFGFGFGF`s@$30FfhmFhs0FfhmFisC$?(F,FGFGFHFI>F_^s-
F16$F,Ffhm@$Ff]lC%-Fi]l6#Q=After~simplification~Box~is|+F)?(F,FGFGFHFIFi^sFc^s@
$/%)ALIAS/3BGFGC&@$Ff]l-Fi]l6#Q-Entering~3B|+F)@$/-FhhnF[glFisC$>FaclFb`sF`s@$F
f]l-Fi]l6#Q,Exiting~3B|+F)@$Ff]lC%-Fi]l6#Q,Box~is~now|+F)?(F,FGFGFHFIFi^sFc^s@$
2Fhs-F`w6#Fb[r?(F,FGFGF[csFIC$?(FFFGFGFHFIC*>Ffhm-F16$FF&Fb[rF9>FWF`\r>FLFe\r>F
fzFfer>FZ-F16$FGFfz>F\s-F16$FgrFfz@$31FWFZF\tF`sFi`l@$/FFF\alC$>FaclFb`sF`s?(F,
FGFGFHFIC$>Ffhm,&Fa^sFGF]^sFis@$1%0ALIAS/maxnewtonGFfhmFi`l@$3/FHF_cn2FfhmF[esC
&@$Ff]l-Fi]l6#Q1Entering~Newton|+F)>Ffhm-Fggm6$Fj_sF_cn@$3/-F`wF^gnFGF_`sC$>Fac
lFb`sF`s@$3Fe`s0FfhmFG?(F,FGFGFHFI>F_^s&FfhmF9?(F,FGFGFHFIC$>FfhmFhds@$1%.ALIAS
/maxkrawGFfhmFi`l@$3F^es2FfhmFifsC&@$Ff]l-Fi]l6#Q3Entering~analyzer|+F)>Ffhm-F\
qFges@$Ff]l-Fi]l6$Q:Exiting~analyzer~with~%a|+F)Ffhm@$/FfhmFGC%@$Ff]l-Fi]l6#Q9K
raw~DONNE~UNE~SOLUTION|+F)>Ffhm-F^blFges@%F`fs@$Ff]l-Fi]l6#QAMais~le~solver~n'a
~pas~converge|+F)C.>FfhmFibl?(F,FGFGFHFI>Ffhm7$F][sF_^s>Fb[r7$-F1F\csFfhm>FfhmF
ibl?(F,FGFGFHFI>Ffhm7$F][s&FhblF9>F^cnF_[s>FacnFf[s>Fb]mFh[sF[\sFjgq>FaclFb`sF`
s>FeirFhs>%'NBVIOLGFhs>F_zFibl?(F,FGFGF_cnFIC2>F_inFat>Ffhm-F``o6$F,7#Fder>FWF`
\r>FLFe\r@$5/FWF]en/FLF]enF`s@$F\drC$>FdblFG?(F)FGFGF)FjdrC+>FdblFhs>FatF^er>Ff
hmF_js>F[frF`\r>F]frFe\r@$F_frFi`l@$FafrFi`l@$2FGFfgrFi`l@$FcfrC%>FdblFG>FWF[fr
>FLF]fr>FWF`\r>FLFe\r>FatF_in@$Ff]lF_hr@$FchrC%>FaclFb`s>FiisFGFi`l@$FihrC%>Fac
lFb`s>FiisFGFi`l@$F^irC%>FaclFb`s>FiisFGFi`l@$Fcir>FeirFfir@$Fhir>FeirFfir@$F[j
rC$>FeirFfir>F_zFcjr@$/FiisFGF`s@$FgjrC,>FfhmFj_s>F^cnF_[s>F_[rFb[s>FacnFf[s>Fb
]mFh[sF[\sFjgq@$3/%5ALIAS/stop_first_solGFG1%-ALIAS/maxsolGFacnFfhr>FaclFb`sF`s
@$32,&FaiqFGF]cnFisF_er2Facl"#?-F86#7#Q'FAILEDF)@$Fe^tC&>FdblFhs?(F,FaclFGF]cnF
IC$?(FFFGFGFHFI>&F\cn6$,(F,FGFaclFisFGFGFF&F\cnFav>Fdbl,&FdblFGFGFG>FaclFG>F]cn
Fdbl?(F,F]cnFisFb`sFI?(FFFGFGFHFI>&F\cn6$,&F,FGFGFGFFFg_t@$5/%7ALIAS/single_bis
ectionGFG/FHFGC(>FfaqFG>Ffhm&F\cn6$FaclFG>FWF`\r>FLFe\r>FdaqFcgr?(F,FgrFGFHFIC'
>FfhmF_^s>FWF`\r>FLFe\r>%&amax1GFcgr@$2FdaqFeatC$>FfaqF,>FdaqFeat@$3/Fe`tFgr2FG
FHC$>Ffaq-FaaqF[gl@$/FfaqFisC(>FfaqFG>FfhmFj`t>FWF`\r>FLFe\r>FdaqFcgr?(F,FgrFGF
HFIC'>FfhmF_^s>FWF`\r>FLFe\r>FeatFcgr@$FgatC$>FfaqF,>FdaqFeat?(F,FGFGFHFI>&F\cn
6$Fb`sF,F_^s>Ffhm&F\cn6$FaclFfaq>FWF`\r>FLFe\r>FZ-Fcs6#*&,&FWFGFLFGFGFgsFis@$51
FZFW1FLFZ?(F)FGFGF)FcdtC$>FatFbt>FZ-Fcs6#*&FadtFGFgsFis>&F\cn6$Fb`sFfaq7$FZFL>F
ict7$FWFZ>F]cn,&F]cnFGFGFGFfhrF)6PF\cnF]cnFHF^cnF_[rF]inF_cnF`uF`cnFacnFbcnFdvF
_inFaiqFdblF[frF]frFeirFiisFfaqFdaqFeatFh\oFdcnFc_sF.Fd[rFb[rFhblF^_mF_sF]uFiuF
`vFc_lF\oFepFdwFexF_^lFfyFccnFefq%'TAYLORGF\hmFbjmF)F),
I)Simplifyf*6#F_q63%%AauxGFeqF>F?F@FbqFcqFbinFcinF(FAFBFCFdbl%&flaggG%#XXG%(rea
peatGF)F)C+>Fd]mF3>Ffz-F16$FGFd]m>Fb]mFhs>F`bmFG>FFFibl?(FLFGFGFHFI>FF7$-F1F^p-
Fcs6#,&&Fd]m6$FLFgrFG&Fd]m6$FLFGFis?(F)FGFGF)/F`bmFGC(>F`bmFhs?(FSFGFG-F`w6#%$C
S2GFIC1>Fg^lFhs>Fc^l-F16$FSFigt>FL-F16$FGFc^l>FfhmF\enZ%>Ffhm--F16$FgrFc^l6#-FO
6$&Fd]mFi[l/FLFUF)F`s>FaclF`gt>F\xF^gt>F\s-Fcs6#-%%sqrtGFd\r>FZF^\r@$30F\sF]en2
Fc\rFhsFbal@$3Fgit2F\sFaclFbal@$3Fgit2F\x,$F\sFisFbal@$2FhsFZC(>FZ-Fcs6#-FcitF[
\l@$32,$FZFisFacl2F\xFZFbal@$2FaclF_jtC%>Fb]mFG>Fg^lFG>FaclF_jt@$2F\sF\xC%>Fb]m
FG>Fg^lFG>F\xF\s@$32FjjtF\xF[[uC%>Fb]mFG>Fg^lFG>F\xFjjt@$32FaclFZFijtC%>Fb]mFG>
Fg^lFG>FaclFZ@$1FZFhsC&>FZ,$F`itFis@$5Fjit30FZF]enF[[uFbal@$3F^]uFa\uC%>FaclFZ>
Fg^lFG>Fb]mFG@$3FgitFc[uC%>F\xF\s>Fg^lFG>Fb]mFG@$Fhbn>Fd]m-Fghp6$/FL7$FaclF\xFd
]m?(FSFGFG-F`w6#%#CSGFIC/>Fg^lFhs>Fc^l-F16$FSFd^u>FLF`ht>FfhmF\enZ%>FfhmFehtF)F
`s>FZF`\r>F\sFe\r>FaclF`gt>F\xF^gt@$F\]uFbal@$F`]uC%>FaclFZ>Fg^lFG>Fb]mFG@$Ff]u
C%>F\xF\s>Fg^lFG>Fb]mFG@$Fhbn>Fd]mF]^u?(FSFGFG-F`w6#%$CS4GFIC/>Fg^lFhs>Fc^l-F16
$FSFb`u>FLF`ht>FfhmF\enZ%>FfhmFehtF)F`s>FZF`\r>F\sFe\r>FaclF`gt>F\xF^gt@$F\]uFb
al@$F`]uC%>FaclFZ>Fg^lFG>Fb]mFG@$Ff]uC%>F\xF\s>Fg^lFG>Fb]mFG@$Fhbn>Fd]mF]^u?(FS
FGFG-F`w6#%$CS3GFIC,>Fg^lFhs>Fc^l-F16$FSF`bu>FfhmF\en>F,F\enZ%C$>Ffhm-F`htFhht>
F,FehtF)F`s>FZF`\r>F\sFe\r>FaclFcdn>F\xFfdn@$533F^]u0F\xF]enF[[u33Fgit0FaclF]en
F[jtFbal?(FLFGFGFHFIC%>FfhmF[gt@$2Fc_s-F\[l6#,&&FFFi[lFGFfhmFis>F`bmFG>FF-Fghp6
$/FLF[gtFF@$/Fb]mFG-F86#Fd]mFeuF)6)Fd^uFigtF`buFb`u%$VARGFHFc_sF)F),
I*filter_2bf*6#%"xG6&F(%%CONSG%#swGF>F)F)C'>FFFibl>F,-%*decomposeGF2>F,7#-F1F9?
(FWFGFG-F`wF9FIC+>%1_EnvAllSolutionsGFI>FS-%&solveG6$-Fbgq6%F3%)rationalG%&exac
tG-F16$FWF,@$0-F`w6#7#FSFGF`s@$0-%+searchtextG6$Q#_NF)-Fbgq6$FSFdgqFhsF`s@$0-Fc
gu6$Q#_BF)FfguFhsF`s@$0-Fcgu6$Q#_ZF)FfguFhsF`s@$/-Fg]r6$-FcsFaw%(numericGFIF`s@
$0-Fcgu6$Q'RootOfF)FfguFhsF`s>FF7$Fjft7$FifuFSFehnF)F)F)F),
I5GeneralSolve_intpakXf*F[fq6/F>F?F@FbqFcqFdqF]hmF[inF_qF(F`fqFafqFainF)F)C>>Fe
fqFffq@$Fhfq>FefqF^gq>Fc^lFfgq-F[hq6#Fc^l>FdcnFhs@$F_hq>FbcnFbhq>F_cnFehq>FHFgh
q@$FihqF]iq>FaiqFbiq>F\cnFdiq>F^cnFibl>F_[rFibl>FacnFhs>FbfqFibl>FatFg[r>F_zFG?
(F,FGFGFjhqFIC'>FfhmF\\r>FWF^\r>FLFc\r@$Fh\r>F_zFhs>Fbfq7$-F16#FbfqF^]r>F`uFibl
>F`cnFibl?(F,FGFGF_cnFIC'>FfhmF`r@$Fe]rC&>F`cnF\^r>FfhmF`\r>F`uFa^rF`s@$Fe]rC&>
F`cnFj^r>FfhmF`\r>F`uF]_rF`s>F`cnFa_r>F`uF]_r?(F,FGFGFHFI>Fbar&FbfqF9?(F,FGFGF_
cnFIC$Z%>FfhmFgarF)C$>FfhmF_brFdbr@$FfbrC$>FfhmF]crFdbr@$F^[oC$?(F,FGFGF_cnFIC1
>F_inFat>FfhmFgar>FWF`\r>FLFe\r@$F\drC$>FdblFG?(F)FGFGF)FjdrC*>FdblFhs>FatF^er>
FfhmFber>F[frF`\r>F]frFe\r@$F_frFi`l@$FafrFi`l@$FcfrC%>FdblFG>FWF[fr>FLF]fr>FWF
`\r>FLFe\r>FatF_in@$Ff]lF_hr@$FchrC$>FacnFhsFfhr@$FihrC$>FacnFhsFfhr@$F^irC$Ffh
r>FacnFhs@$Fcir>FeirFfir@$Fhir>FeirFfir@$F[jrC$>FeirFfir>Ffz7$-F16#FfzFfhm@$Fgj
rC*>FfhmFibl?(F,FGFGFHFI>FfhmF\[s>F^cnF_[s>F_[r7$Fc[sFfz>FacnFf[s>Fc^lFh[s-F\\s
6%Fc^lF^\sF^cnFhiu>FaclFG>F]cnFG?(F)FGFGF)Fb\sC1@$F_hq@$2*&F`hqFGFh\sFGFj\sF[]s
@$F`]sC&Fb]sFf]s?(F,FGFGFHFIFj]sFc^s@$Ff]lC%Ff]s?(F,FGFGFHFIFi^sFc^s@$F]_sC'>Fd
cnFG>Fc_sFd_s>FfhmFh_s@$F_`sC$>FaclFb`sF`s@$Fd`sC$?(F,FGFGFHFI>F_^sFj`s@$Ff]lC%
F^as?(F,FGFGFHFIFi^sFc^s@$FcasC$@$/-F_ioF[glFisC$>FaclFb`sF`s@$Ff]lC%Febs?(F,FG
FGFHFIFi^sFc^s>FeirFhs>FiisFhs>FfzFibl?(F,FGFGF_cnFIC2>F_inFatZ%>FfhmFberF)F`s>
FWF`\r>FLFe\r@$FejsF`s@$F\drC$>FdblFG?(F)FGFGF)FjdrC+>FdblFhs>FatF^er>FfhmFber>
F[frF`\r>F]frFe\r@$F_frFi`l@$FafrFi`l@$Fe[tFi`l@$FcfrC%>FdblFG>FWF[fr>FLF]fr>FW
F`\r>FLFe\r>FatF_in@$Ff]lF_hr@$FchrC%>FaclFb`s>FiisFGFi`l@$FihrC%>FaclFb`s>Fiis
FGFi`l@$F^irC%>FaclFb`s>FiisFGFi`l@$Fcir>FeirFfir@$Fhir>FeirFfir@$F[jrC$>FeirFf
ir>FfzF^`v@$Fd]tF`s@$FgjrC->FfhmFibl?(F,FGFGFHFI>FfhmFjhs>F^cnF_[s>F_[rFh`v>Fac
nFf[s>Fc^lFh[sF[avFhiu@$F]^tFfhr>FaclFb`sF`s@$Fd^tFi^t@$Fe^tC&>FdblFhs?(F,FaclF
GF]cnFIC$?(FFFGFGFHFI>Fd_tFg_t>FdblFi_t>FaclFG>F]cnFdbl?(F,F]cnFisFb`sFI?(FFFGF
GFHFI>F_`tFg_t@$5Fd`tF]btC1>FfaqFG>FfhmFj`t>FWF`\r>FLFe\r>FdaqFcgr?(F,FgrFGFHFI
C'>FfhmF_^s>FWF`\r>FLFe\r>FeatFcgr@$FgatC$>FfaqF,>FdaqFeat?(F,FGFGFHFI>FfctF_^s
>FfhmFict>FWF`\r>FLFe\r>FZ-Fcs6#*&FadtFGFgsFis@$Fcdt?(F)FGFGF)FcdtC$>FatFbt>FZ-
Fcs6#*&FadtFGFgsFis>F^etF`et>FictFbet>F]cnFdetFfhrF)6<FbfqF\cnF]cnFHF^cnF_[rF]i
nF_cnF`uF`cnFacnFbcnF_inFaiqFdblF[frF]frFeirFiisFfaqFdaqFeatFh\oFdcnFc_sFefqF)F
)