MVR4
I(amoeba1:6'%"pG%"yG%$dimG%*abserrmaxG%,iter_returnG60%%nptsG%"jG%%inhiG%$iloG%
$ihiG%"iG%%yprrG%$yprG%(abserr2G%#prG%$prrG%%pbarG%%iterG%+abserrmax2G6"F9C+>8-
-%&arrayG6#;"""9&>8.F=>8/F=>8$,&FBFAFAFA>80""!>81*$9'""#?(F9FAFAF9%%trueGC/>8'F
A@%2&9%6#FQ&FZ6#FAC$>8(FA>8&FQC$>FjnFQ>F\oFA?(8)FAFAFHFSC$@$2&FZ6#Fao&FZ6#FV>FV
Fao@&2&FZ6#FjnFeoC$>F\oFjn>FjnFao2&FZ6#F\oFeo@$0FaoFjn>F\oFao>8,FL?(8%FAFAFBFS>
Fhp,&FhpFA*$,&&9$6$FVFjpFA&F`q6$FjnFjp!""FQFA@$2FhpFN%&breakG@$/FK"%+?C$>&9(Fgn
FK-%'RETURNG6#Fdq>FK,&FKFAFAFA?(FjpFAFAFBFS>&FF6#FjpFL?(FaoFAFAFHFS@$Fep?(FjpFA
FAFBFS>Ffr,&FfrFA&F`q6$FaoFjpFA?(FjpFAFAFBFSC$>Ffr*&FfrFAFBFdq>&F<Fgr,&Ffr$"#?F
dqFbq$!#5Fdq>8+-%,amoeba~funcG6#F<@'1F[tFgoC%?(FjpFAFAFBFS>&FDFgr,&FdsFfsFfrFhs
>8*-F]t6#FD@%2FgtFgoC$?(FjpFAFAFBFS>FbqFdt>F\pFgtC$?(FjpFAFAFBFS>FbqFds>F\pF[t1
FbpF[tC&@$2F[tF\pC$?(FjpFAFAFBFS>FbqFds>F\pF[t?(FjpFAFAFBFS>Fdt,&Fbq$""&FdqFfrF
_v>FgtFht@%2FgtF\pC$?(FjpFAFAFBFS>FbqFdt>F\pFgt?(FaoFAFAFHFS@$0FaoFVC$?(FjpFAFA
FBFSC$>Fds,&F]sF_vF_qF_v>F]sFds>FeoF\tC$?(FjpFAFAFBFS>FbqFds>F\pF[t>F]rFKFVF9F9
F9,
I-amoeba~debug%&falseGF9,
I%lsqd:6)'%&xlistG-%%listG6#%(numericG'%&ylistGF^x'%#eqG/%%nameG%*algebraicG'%)
indepvarGFgx'%+firstguessG-F_x6#/FgxFax'%)simpsizeGFax'%(simpaccG%(integerG62%"
GG%&ypredG%%devvG%%ssqqG%"kG%"QG%$ansG%'depvarG%$fcnG%%consG%#npG%$outG%&ndataG
%&paramG%(diffsetGF06"FezC5>8+-%$lhsG6#9&>8,-%$rhsGF[[l>8--%$mapG6$Fjz9(>8.-%%n
opsG6#Fb[l>80-Fj[l6#9$@$0-Fj[l6#9%F]\l-%&ERRORG6#%AThe~lists~must~have~equal~le
ngthG>81-%#opG6#-%'selectG6%%%typeG-%&minusG6$-%'indetsG6#F^[l<#9'Fgx>82-Fd]l6$
<#F[]l<#-F]]lF[\l@$0F\^l<"-Fg\l6#%dqInconsistent~input~-~check~independent~vari
able~name~and~that~initial~guesses~provided~for~all~parametersG>8$-%(unapplyG6$
F^[lFj]l>8%-Fd[l6$Fi^lF`\l>8&-%$zipG6%:6$%"sG%"tGFez6$%)operatorG%&arrowGFez*$,
&F`\l"""Fe\l!""""#FezFezFe\lF^_l>8'""!?&8(Fb_lFS>Fc`l,&Fc`lF_`lFf`lF_`l>8)-F[_l
6$Fc`lFa^l>8*-%'amoebaG6&Fj`l-Fd[l6$F`[lFf[l9))"#5,$9*F``l>8/7#-%$seqG6$/&Fb[l6
#83-%&evalfG6#&F^alFabl/Fbbl;F_`lFh[l>%$SSQG-Fj`l6#-F]]l6#F^al/Fhz-%%subsG6$Fja
l-Fi^l6#Fj]lFez6#FjblFez,
I'amoeba:6''%"fG%*procedureG'%&startG-F_x6#Fax'%&scaleGFax'%'abserrGFax%)noeval
hfG6*F'%'EVALHFG%%fargGF%F&F7F.F,F9F9C2>%,amoeba~iterGFL>FH-Fj[l6#FZ@%/-Fj[l6#7
#-F]]l6$FA-%%evalG6#F`qFA>F]tFeel>F]t-F[_l6$-.F`q6#-F]bl6$&F\o6#F[t/F[t;FAFHF\o
@$Fgw-%&printG6#F]t>Fjp29#F`v@$FjpC%>FV-F>F]el-%*traperrorG6#-%'evalhfG6#-F]tFh
o@$/%""G%*lasterrorGC$>FjpFhw-Fgfl6#%MWarning:~cannot~use~hardware~floating~poi
nt.G>FV7$FZ-F]bl6$-%'subsopG6$/F[t,&&FZFbflFAFBFAFZFcfl>Fjn-%'vectorG6#-Fd[l6$F
]tFV>FV-%'matrixGFho>Fao-F>6#;FAFA@%Fjp-Fdgl6#-F"6'-%$varGFho-F]jlF]pFHFP-F]jlF
foC%>Fjn-Fd[l6$FdblFjn>FV-Fd[l6$FdblFV-F"6'FVFjnFH-Fdbl6#FPFao>Fgt-%&roundG6#Fi
gl>Fjdl-F][m6#&FaoFgn@$/FgtFdq-Fg\l6#%5too~many~iterations.G@$Fgw-Fgfl6&FgtFVFj
nFjdl7#-F]bl6$&FV6$FgtF[tFcflF96$F]tFjdlF9
