https://github.com/pdieterle/diffWavePropAndInit
Tip revision: f8d9feffd57d05f47c8c14c6d9850643b2858d0a authored by pdieterle on 18 December 2020, 14:53:17 UTC
Add files via upload
Add files via upload
Tip revision: f8d9fef
figureCode.nb
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 0, 0]
NotebookDataLength[ 123615, 2762]
NotebookOptionsPosition[ 119733, 2692]
NotebookOutlinePosition[ 120125, 2708]
CellTagsIndexPosition[ 120082, 2705]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["\<\
In this notebook, we calculate the dynamics of the wave front with a \
Heaviside activation function; we explicitly plot the concentration profiles \
from figures 1 and 2 and the concentration gradients a la figure 4 of \
Dieterle et al. (2020).\
\>", "Subsubsection",
CellChangeTimes->{{3.812824043669132*^9, 3.812824164196375*^9}, {
3.812893622589428*^9, 3.812893623529263*^9}, {3.814016904415378*^9,
3.814016905516777*^9}, {3.8141863983235483`*^9,
3.8141863991134987`*^9}},ExpressionUUID->"9f649eb5-553d-479c-a38a-\
504bc2532627"],
Cell["\<\
First, we define the system parameters; vBase is the desired wave speed, Dc \
is the diffusion constant, cellDim is the number of dimensions in which the \
cells live, diffDim is the number of dimensions in which diffusion takes \
place, Cth is the threshold concentration (implicitly normalized by the \
emission rate, a, and cell density, rho) which gives the desired wave speed \
vBase given the other system parameters, and rInv is the initial signaling \
radius.\
\>", "Text",
CellChangeTimes->{{3.8128200535738564`*^9, 3.812820245902184*^9}, {
3.812820298392951*^9, 3.812820403549923*^9}, {3.812820472733918*^9,
3.812820600792253*^9}, {3.812820650546049*^9, 3.812820677309263*^9}, {
3.8128210027545433`*^9,
3.812821003567994*^9}},ExpressionUUID->"8110de1f-e0f6-4e12-82cf-\
66fe630a4d28"],
Cell[BoxData[{
RowBox[{
RowBox[{"vBase", " ", "=", " ",
RowBox[{"2.", " ",
SuperscriptBox["10",
RowBox[{"-", "6"}]]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Dc", "=", " ",
RowBox[{"1.", " ",
SuperscriptBox["10",
RowBox[{"-", "10"}]]}]}], ";"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"cellDim", " ", "=", " ", "2"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"diffDim", " ", "=", " ", "2"}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{
RowBox[{"If", " ", "diffDim"}], " ", "\[Equal]", " ", "cellDim"}], ",",
" ",
RowBox[{
RowBox[{
RowBox[{"this", " ", "gives", " ", "Cth"}], " ", "=", " ",
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}]}], ";", " ",
RowBox[{
RowBox[{"if", " ", "diffDim"}], " ", "\[Equal]", " ",
RowBox[{"cellDim", "+", "1"}]}]}], ",", " ",
RowBox[{
RowBox[{"it", " ", "gives", " ", "Cth"}], " ", "=", " ",
RowBox[{
RowBox[{"2", "/", "Pi"}], "/", "vBase"}]}]}],
"*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Cth", " ", "=", " ",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"diffDim", "+", "1", "-", "cellDim"}], ")"}],
FractionBox["Dc",
SuperscriptBox["vBase", "2"]]}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"diffDim", "-", "cellDim"}], ")"}],
FractionBox["2",
RowBox[{"Pi", " ", "vBase"}]]}]}]}], ";"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"rInv", " ", "=",
RowBox[{"3", " ",
RowBox[{"Dc", "/", "vBase"}]}]}], ";"}]}], "Input",
CellGroupingRules->{GroupTogetherGrouping, 10000.},
CellChangeTimes->{{3.756040794283849*^9, 3.7560409737156277`*^9}, {
3.75604106891358*^9, 3.756041077856443*^9}, {3.756041118891117*^9,
3.7560411311359787`*^9}, {3.756041625359209*^9, 3.7560416281275873`*^9}, {
3.7560433596173677`*^9, 3.7560433894801188`*^9}, {3.7560439116459417`*^9,
3.756043968491229*^9}, {3.7560523295764217`*^9, 3.756052329715748*^9}, {
3.756052391065074*^9, 3.75605239429084*^9}, {3.7560524918373213`*^9,
3.756052492851796*^9}, {3.7560529547784967`*^9, 3.756052967544126*^9}, {
3.756053131028384*^9, 3.7560531359480658`*^9}, {3.756054050356715*^9,
3.756054057909041*^9}, {3.756054158357193*^9, 3.7560541595083923`*^9}, {
3.7560544723714314`*^9, 3.756054472487576*^9}, {3.756054509041484*^9,
3.7560545091296597`*^9}, {3.756054541721381*^9, 3.7560545422732363`*^9}, {
3.75605457420815*^9, 3.756054574288047*^9}, {3.75605461527072*^9,
3.756054615461893*^9}, {3.756054650617159*^9, 3.756054650728394*^9}, {
3.756054691068261*^9, 3.756054691859223*^9}, {3.756054836257433*^9,
3.756054846600268*^9}, {3.756055364428463*^9, 3.7560553653563766`*^9},
3.756056041355762*^9, {3.7560562541949387`*^9, 3.756056254531012*^9}, {
3.756056501928475*^9, 3.756056502207961*^9}, {3.756056584761302*^9,
3.756056585417697*^9}, {3.7560566316870213`*^9, 3.7560566325669613`*^9}, {
3.756057374537209*^9, 3.75605737801688*^9}, 3.756058053543765*^9, {
3.756058174253119*^9, 3.7560581818555527`*^9}, {3.756062017703886*^9,
3.756062017870302*^9}, {3.756077700097196*^9, 3.756077707824307*^9}, {
3.756129800455793*^9, 3.7561298014220543`*^9}, {3.7561301288431187`*^9,
3.756130132403393*^9}, {3.7561302004613857`*^9, 3.756130235950314*^9}, {
3.756130336113912*^9, 3.756130337192069*^9}, {3.756131697429245*^9,
3.756131709521144*^9}, {3.756133289263569*^9, 3.756133315655239*^9}, {
3.75613545873929*^9, 3.756135459795003*^9}, {3.756136034547943*^9,
3.756136052125804*^9}, {3.756138470964869*^9, 3.7561384719709463`*^9}, {
3.7561385510140343`*^9, 3.7561385516375628`*^9}, {3.756139531021203*^9,
3.756139531317412*^9}, {3.756140720671722*^9, 3.756140720780244*^9}, {
3.756140797703512*^9, 3.756140823970233*^9}, {3.756140885828114*^9,
3.75614088668432*^9}, {3.75614111265945*^9, 3.7561411134728327`*^9}, {
3.756141154142771*^9, 3.7561411553858852`*^9}, {3.756141283642742*^9,
3.756141283840839*^9}, {3.7561414168529453`*^9, 3.7561414174134407`*^9}, {
3.756143269378447*^9, 3.756143312105802*^9}, {3.75614352401438*^9,
3.756143604793655*^9}, {3.7561441925726223`*^9, 3.7561441927551527`*^9}, {
3.7561442823231916`*^9, 3.7561442824713573`*^9}, {3.756144877101622*^9,
3.7561449253353567`*^9}, {3.7561449893570642`*^9, 3.756145101672133*^9}, {
3.756145282971572*^9, 3.756145287343123*^9}, {3.756145525624873*^9,
3.756145549957613*^9}, {3.7561458827337933`*^9, 3.756145887315836*^9}, {
3.756146038737795*^9, 3.756146046181807*^9}, {3.756147030546771*^9,
3.756147031178623*^9}, {3.756147332640277*^9, 3.7561473329688683`*^9}, {
3.756148536161846*^9, 3.756148536945239*^9}, {3.756148689713982*^9,
3.756148690855893*^9}, {3.7561487312103643`*^9, 3.75614873236971*^9}, {
3.756148825568096*^9, 3.7561488267585993`*^9}, {3.756148921095241*^9,
3.756148921317053*^9}, {3.75614938009613*^9, 3.756149380294298*^9}, {
3.756149474545126*^9, 3.756149477432151*^9}, {3.7561495783231697`*^9,
3.756149578643935*^9}, {3.7561496697457952`*^9, 3.7561496698943567`*^9}, {
3.756150228291881*^9, 3.756150230801107*^9}, {3.7561503310093403`*^9,
3.756150333471835*^9}, {3.756150453489944*^9, 3.7561504546419897`*^9}, {
3.75615049204533*^9, 3.756150492164599*^9}, {3.75615054336139*^9,
3.756150543632666*^9}, {3.7561505975490303`*^9, 3.75615059776423*^9}, {
3.7561511506690083`*^9, 3.756151150803212*^9}, {3.756151488570496*^9,
3.756151492394701*^9}, {3.7561545099136763`*^9, 3.756154518270939*^9}, {
3.756154957620084*^9, 3.7561549577684813`*^9}, {3.756155005039806*^9,
3.756155007967*^9}, {3.756156402734787*^9, 3.756156406614202*^9}, {
3.756157009377481*^9, 3.75615704667699*^9}, {3.756157084897894*^9,
3.756157146872017*^9}, 3.756157198258122*^9, {3.7561572468841057`*^9,
3.7561572854120092`*^9}, {3.756157350934924*^9, 3.756157351292574*^9}, {
3.756157407226466*^9, 3.756157411198298*^9}, {3.756157508339491*^9,
3.756157508681473*^9}, {3.756157619912163*^9, 3.756157627280142*^9}, {
3.756225604850662*^9, 3.7562256057676*^9}, {3.756226169436708*^9,
3.756226173034774*^9}, {3.756226872294127*^9, 3.756226872980144*^9}, {
3.756226982384701*^9, 3.756226983696072*^9}, {3.756227067923295*^9,
3.756227068569716*^9}, {3.756227145070764*^9, 3.756227145590602*^9}, {
3.756227206004986*^9, 3.756227206644208*^9}, {3.75623193018738*^9,
3.7562319347713947`*^9}, {3.756232287059037*^9, 3.756232305886351*^9}, {
3.75623234003465*^9, 3.756232340190064*^9}, {3.7562362670251513`*^9,
3.756236383052437*^9}, {3.7562400522381077`*^9, 3.756240053379797*^9},
3.756240090003368*^9, {3.756243508506473*^9, 3.7562435736628447`*^9}, {
3.7562436429445972`*^9, 3.7562436430247383`*^9}, {3.756243678116371*^9,
3.7562437111141033`*^9}, {3.756243746122505*^9, 3.756243746234112*^9}, {
3.7562439470694437`*^9, 3.756243985330662*^9}, {3.756255728600194*^9,
3.756255760457664*^9}, {3.756255908961433*^9, 3.756255909815662*^9}, {
3.756255956328353*^9, 3.75625596257421*^9}, {3.756255997934292*^9,
3.756256004277091*^9}, {3.756256106186247*^9, 3.756256126182205*^9}, {
3.756256271020405*^9, 3.75625631255354*^9}, {3.756258665022428*^9,
3.7562587623777523`*^9}, {3.7562588182789*^9, 3.756258893038398*^9}, {
3.756258929697414*^9, 3.756258978078143*^9}, {3.75625908773316*^9,
3.7562591502087107`*^9}, {3.756305013958003*^9, 3.7563050342124853`*^9}, {
3.756309377701272*^9, 3.756309383868926*^9}, {3.756327113985896*^9,
3.7563271399046183`*^9}, {3.756327169917973*^9, 3.75632717005973*^9}, {
3.7565623309094353`*^9, 3.756562358135408*^9}, {3.7565705899457493`*^9,
3.756570590314167*^9}, {3.7565754282311172`*^9, 3.756575431388774*^9},
3.756577570086112*^9, {3.756582890822172*^9, 3.756582891793272*^9}, {
3.7565851475372257`*^9, 3.75658514816711*^9}, {3.756588375421011*^9,
3.756588375593259*^9}, {3.7565884616333847`*^9, 3.756588462992519*^9}, {
3.7566054736868563`*^9, 3.756605473994266*^9}, {3.7566062285836163`*^9,
3.7566062510611143`*^9}, {3.756644135538597*^9, 3.756644165033538*^9}, {
3.756647362147204*^9, 3.756647446673822*^9}, {3.756647655500361*^9,
3.756647686162757*^9}, {3.756647825792848*^9, 3.756647825917273*^9}, {
3.7566479774605093`*^9, 3.7566479781880836`*^9}, {3.756649211056983*^9,
3.7566492352283697`*^9}, {3.756650912841012*^9, 3.756650948528545*^9}, {
3.756660251464552*^9, 3.756660252646566*^9}, {3.756662518002542*^9,
3.7566625288461313`*^9}, {3.756663417501663*^9, 3.756663418355077*^9}, {
3.756664174590687*^9, 3.756664214142672*^9}, {3.7566883060751057`*^9,
3.756688311417356*^9}, {3.756711794572661*^9, 3.756711812139995*^9}, {
3.756711855513734*^9, 3.7567118565851803`*^9}, {3.7567275301170473`*^9,
3.756727540037142*^9}, {3.756727570318012*^9, 3.75672757111659*^9}, {
3.756731608090773*^9, 3.7567316282893267`*^9}, {3.7567334289044952`*^9,
3.756733429716608*^9}, {3.7567337693452787`*^9, 3.756733769461317*^9}, {
3.756734057198391*^9, 3.7567340572703333`*^9}, {3.756810300859964*^9,
3.756810320209433*^9}, {3.756810356015952*^9, 3.75681036109608*^9}, {
3.7568133106316023`*^9, 3.756813342178834*^9}, {3.756813503944121*^9,
3.7568135107634993`*^9}, {3.756815167990375*^9, 3.756815170391296*^9}, {
3.756815258259582*^9, 3.756815293004778*^9}, {3.7568313421480227`*^9,
3.756831376874262*^9}, {3.756831431933427*^9, 3.756831433593974*^9}, {
3.7568357924817133`*^9, 3.756835792952022*^9}, {3.756835902459509*^9,
3.756835968142323*^9}, {3.757614447230679*^9, 3.757614468804822*^9}, {
3.7576162872473173`*^9, 3.7576162912146683`*^9}, {3.757616769618321*^9,
3.75761677051409*^9}, {3.757616996498233*^9, 3.757616997314636*^9}, {
3.757688941673254*^9, 3.757688941894353*^9}, {3.757693357947563*^9,
3.757693368537472*^9}, {3.7577163318016033`*^9, 3.757716360831294*^9}, {
3.757716398075338*^9, 3.757716493209999*^9}, {3.7577165522055492`*^9,
3.757716553356564*^9}, {3.7577165925009813`*^9, 3.7577165932526817`*^9}, {
3.7577166355876837`*^9, 3.757716635762971*^9}, 3.75771674579595*^9, {
3.7577618165947027`*^9, 3.7577618171286077`*^9}, {3.75795092926453*^9,
3.757950931590828*^9}, {3.7580324893461*^9, 3.758032545703188*^9}, {
3.7580326448631573`*^9, 3.758032726528969*^9}, 3.758033589275426*^9, {
3.758034474721138*^9, 3.758034497015956*^9}, {3.758034530801717*^9,
3.75803454190239*^9}, {3.758039142235528*^9, 3.7580391463050137`*^9}, {
3.758039621582275*^9, 3.758039622020924*^9}, {3.758051314188912*^9,
3.758051325120873*^9}, {3.758055534883383*^9, 3.7580555356410933`*^9}, {
3.758057440784362*^9, 3.758057449046328*^9}, {3.7580594693287582`*^9,
3.758059469413857*^9}, {3.758065152182499*^9, 3.758065152325204*^9}, {
3.758069848696683*^9, 3.758069856238051*^9}, {3.758070102731283*^9,
3.758070103111908*^9}, {3.758070193489493*^9, 3.7580701955039473`*^9}, {
3.758072526535928*^9, 3.758072527567719*^9}, {3.7580741825930433`*^9,
3.758074183391226*^9}, {3.758116061326193*^9, 3.7581160622687473`*^9}, {
3.758117440070446*^9, 3.758117475003985*^9}, {3.758117511283813*^9,
3.758117511457959*^9}, {3.758117541690556*^9, 3.7581175418333483`*^9}, {
3.758117573001584*^9, 3.7581176076867933`*^9}, {3.758129690828404*^9,
3.758129725896001*^9}, {3.758129862252152*^9, 3.758129863905121*^9}, {
3.758130416969091*^9, 3.758130447775424*^9}, {3.758130610350943*^9,
3.758130612511983*^9}, {3.758130705513077*^9, 3.7581307522555637`*^9}, {
3.758151178354039*^9, 3.758151179041065*^9}, {3.758153219766374*^9,
3.758153225171309*^9}, {3.75820089675028*^9, 3.758200907972134*^9}, {
3.758293022387354*^9, 3.7582931038199177`*^9}, {3.7582935585037622`*^9,
3.75829355855153*^9}, {3.7583060340359373`*^9, 3.758306040313528*^9}, {
3.758306105596333*^9, 3.758306117503619*^9}, {3.759755953014497*^9,
3.759755953525093*^9}, {3.759756583806087*^9, 3.759756583900646*^9}, {
3.759757533879621*^9, 3.7597575339822073`*^9}, {3.759759637953779*^9,
3.759759638110388*^9}, {3.759760526663843*^9, 3.7597605308629227`*^9}, {
3.759761054142308*^9, 3.759761071910428*^9}, {3.759761123124078*^9,
3.759761125289895*^9}, {3.7597624430786467`*^9, 3.759762473043676*^9}, {
3.759767809506373*^9, 3.7597678114464197`*^9}, {3.7597694935207376`*^9,
3.759769493788946*^9}, {3.7597695309348583`*^9, 3.759769531123754*^9}, {
3.759769562625852*^9, 3.759769562769313*^9}, {3.759771090021666*^9,
3.759771090596549*^9}, {3.759775602138888*^9, 3.7597756023064013`*^9}, {
3.759775931864544*^9, 3.759775932047456*^9}, {3.759784590960981*^9,
3.7597845910714207`*^9}, {3.759796993540801*^9, 3.759796994145352*^9}, {
3.759831611172058*^9, 3.759831611948435*^9}, {3.759837328060796*^9,
3.759837356305894*^9}, {3.759837428903674*^9, 3.759837465891781*^9}, {
3.7598375826726437`*^9, 3.7598375827727633`*^9}, {3.759842982504743*^9,
3.759842993910335*^9}, {3.759843076410368*^9, 3.759843076671542*^9}, {
3.759843114352209*^9, 3.759843114805949*^9}, 3.7598536363670597`*^9, {
3.759854076811966*^9, 3.759854077328336*^9}, {3.759854115715261*^9,
3.759854155417315*^9}, {3.75985419861312*^9, 3.759854206051838*^9}, {
3.759854467838683*^9, 3.7598544787779713`*^9}, {3.759858530112516*^9,
3.7598585302381907`*^9}, {3.759858562081483*^9, 3.759858562173325*^9}, {
3.7598607377311697`*^9, 3.7598607754454813`*^9}, {3.759861276322631*^9,
3.7598612899767*^9}, {3.759862102141767*^9, 3.759862226794729*^9}, {
3.7598622836245813`*^9, 3.759862292038004*^9}, {3.759882009692123*^9,
3.759882064740757*^9}, 3.759923623059926*^9, {3.7605412996072187`*^9,
3.7605413002844152`*^9}, {3.760541371175109*^9, 3.760541374452176*^9}, {
3.7605498117384157`*^9, 3.760549814232473*^9}, {3.760550114545183*^9,
3.760550126749456*^9}, {3.760551823385467*^9, 3.760551849511896*^9}, {
3.760552362181348*^9, 3.7605523662827806`*^9}, {3.760552435565304*^9,
3.760552438122716*^9}, {3.760552478101305*^9, 3.760552509736437*^9}, {
3.7605525912717733`*^9, 3.760552616640378*^9}, {3.760552697629859*^9,
3.760552744919639*^9}, {3.760552779576543*^9, 3.760552816871423*^9}, {
3.760553136713372*^9, 3.760553172236228*^9}, {3.760553794793456*^9,
3.760553796352722*^9}, {3.760553843913649*^9, 3.760553862198566*^9}, {
3.760553892925988*^9, 3.7605539131894283`*^9}, {3.760553953268116*^9,
3.76055397069937*^9}, {3.760555533635358*^9, 3.760555535480258*^9}, {
3.760555868330042*^9, 3.760555891443494*^9}, {3.7605559459626007`*^9,
3.7605559488348103`*^9}, {3.760556019205399*^9, 3.7605560224811983`*^9},
3.760558262830044*^9, {3.760981974913828*^9, 3.7609820063102283`*^9}, {
3.760982057998498*^9, 3.760982058284637*^9}, {3.760982206579865*^9,
3.760982256006054*^9}, {3.760982808205669*^9, 3.760982844074758*^9}, {
3.7609828752321043`*^9, 3.760982969131261*^9}, {3.760988875859577*^9,
3.7609889037133017`*^9}, {3.7609890358434477`*^9, 3.760989036706109*^9}, {
3.760989069474533*^9, 3.760989194904731*^9}, {3.7609892305019617`*^9,
3.76098924699638*^9}, {3.7609904096161537`*^9, 3.760990413844576*^9}, {
3.760990450052904*^9, 3.760990606886175*^9}, {3.760990769070304*^9,
3.760990778610754*^9}, {3.7609917696035624`*^9, 3.760991796241848*^9}, {
3.7609989641520977`*^9, 3.76099898489468*^9}, {3.760999037222323*^9,
3.7609990389170647`*^9}, {3.760999088025362*^9, 3.760999102357555*^9}, {
3.761000613563283*^9, 3.761000877460992*^9}, {3.761001488841176*^9,
3.761001520702045*^9}, {3.7610422310801077`*^9, 3.761042232030057*^9}, {
3.761042657116989*^9, 3.761042683130018*^9}, {3.761055036230589*^9,
3.761055049719014*^9}, {3.761388153124906*^9, 3.761388188145022*^9}, {
3.761388291933524*^9, 3.76138830843546*^9}, {3.761388347556218*^9,
3.761388356711541*^9}, {3.7614015039805593`*^9, 3.7614015319517736`*^9}, {
3.761401563457027*^9, 3.7614015640789557`*^9}, {3.761495595812365*^9,
3.761495595930252*^9}, {3.76149563963737*^9, 3.7614956638762407`*^9}, {
3.761495873960781*^9, 3.761495875716236*^9}, {3.761496027934783*^9,
3.761496028406333*^9}, 3.761496086903048*^9, {3.761496202827729*^9,
3.761496210473881*^9}, {3.7614963335937967`*^9, 3.761496339346119*^9}, {
3.761496394045225*^9, 3.7614963958357487`*^9}, 3.7614965487006617`*^9, {
3.761496780596045*^9, 3.761496780960857*^9}, {3.761501199152636*^9,
3.761501242849309*^9}, {3.761501763620661*^9, 3.761501810342886*^9}, {
3.761502096914987*^9, 3.7615021335277576`*^9}, {3.761502743710353*^9,
3.7615027502373047`*^9}, 3.7615028920049334`*^9, {3.7615066822403927`*^9,
3.761506729415237*^9}, {3.76150694746014*^9, 3.76150699772266*^9}, {
3.7615777670497093`*^9, 3.761577767374855*^9}, {3.761577945719439*^9,
3.76157796745548*^9}, {3.761648991561228*^9, 3.76164905230686*^9}, {
3.7616490920748177`*^9, 3.7616491231428127`*^9}, {3.761649262453775*^9,
3.76164937414085*^9}, {3.761649420400898*^9, 3.7616494774589357`*^9}, {
3.7616495079083557`*^9, 3.761649561031261*^9}, {3.76164960147747*^9,
3.761649623488641*^9}, {3.761649708937586*^9, 3.761649760878421*^9}, {
3.7616498544567423`*^9, 3.761649858575508*^9}, {3.761649903639536*^9,
3.76164992470897*^9}, {3.761650052860979*^9, 3.76165009027686*^9}, {
3.761650135001959*^9, 3.761650140119583*^9}, 3.761650186492283*^9, {
3.761650282328125*^9, 3.761650290501809*^9}, {3.761650336300702*^9,
3.761650340948407*^9}, {3.7616505220062838`*^9, 3.761650534705616*^9}, {
3.761650584797076*^9, 3.7616505891947527`*^9}, 3.761650619335676*^9, {
3.761650695980661*^9, 3.761650756403591*^9}, {3.761650821118183*^9,
3.7616509705141363`*^9}, {3.761651263808712*^9, 3.761651284163604*^9}, {
3.7616513438348713`*^9, 3.7616514283517036`*^9}, {3.761651459792959*^9,
3.761651601189658*^9}, {3.761651638182486*^9, 3.761651662557208*^9}, {
3.7632280812835913`*^9, 3.7632282065657578`*^9}, {3.763228301481269*^9,
3.763228301527849*^9}, {3.763229642680606*^9, 3.763229642869248*^9}, {
3.763230689681563*^9, 3.763230690359836*^9}, {3.763231932242154*^9,
3.763231949359858*^9}, {3.763244423058126*^9, 3.7632444504527607`*^9}, {
3.76324493246815*^9, 3.7632449650181713`*^9}, {3.7632931179183493`*^9,
3.763293118012298*^9}, {3.763323681864105*^9, 3.763323682311089*^9}, {
3.763326528907772*^9, 3.763326529819607*^9}, {3.763326592180969*^9,
3.763326605137082*^9}, {3.763326647396179*^9, 3.763326648298625*^9}, {
3.76332668035071*^9, 3.7633266936490602`*^9}, {3.763326810938946*^9,
3.763326845378707*^9}, {3.764328990600479*^9, 3.764329008575946*^9}, {
3.764329070263544*^9, 3.764329080075526*^9}, {3.764329173975051*^9,
3.764329260266923*^9}, {3.764331412349525*^9, 3.764331412427043*^9}, {
3.764334919777192*^9, 3.764334944679041*^9}, 3.7643349876602993`*^9, {
3.764763733966786*^9, 3.7647637357665167`*^9}, {3.764763969198152*^9,
3.764763970267178*^9}, {3.764765026509206*^9, 3.764765064249237*^9}, {
3.7648700476102123`*^9, 3.7648702241726227`*^9}, {3.764870494352401*^9,
3.764870521523547*^9}, {3.764951701214465*^9, 3.764951729051528*^9},
3.76527737344252*^9, {3.7652774346944447`*^9, 3.765277436331203*^9}, {
3.765279235695506*^9, 3.765279235999269*^9}, {3.765545148305408*^9,
3.765545162853837*^9}, {3.765545767617886*^9, 3.765545769609212*^9}, {
3.766241687956572*^9, 3.766241703328446*^9}, {3.766242802737154*^9,
3.766242805584771*^9}, {3.7662429991237383`*^9, 3.766243066807025*^9}, {
3.766243452901267*^9, 3.766243542166916*^9}, {3.766243708938624*^9,
3.766243740094718*^9}, {3.766243903413422*^9, 3.7662439087750998`*^9}, {
3.766243945925503*^9, 3.766243967442655*^9}, {3.77187108303666*^9,
3.771871095793943*^9}, 3.771871920661705*^9, {3.771873260593762*^9,
3.771873261991804*^9}, {3.7721903173287287`*^9, 3.772190331804928*^9}, {
3.772190376357546*^9, 3.772190425857958*^9}, {3.772478127290346*^9,
3.772478143216392*^9}, 3.7724798308037853`*^9, {3.772480362892712*^9,
3.772480363042886*^9}, {3.772820226588163*^9, 3.772820226842*^9}, {
3.772820270787107*^9, 3.7728202920646687`*^9}, {3.7728806956332417`*^9,
3.7728807983413887`*^9}, {3.772880895741919*^9, 3.7728808995796423`*^9}, {
3.772881089474718*^9, 3.772881089616796*^9}, {3.772885181480646*^9,
3.7728851816298*^9}, {3.77288523482024*^9, 3.772885237477182*^9}, {
3.7728853050123863`*^9, 3.7728853051788883`*^9}, {3.772894362137125*^9,
3.772894362336488*^9}, {3.7728944498886003`*^9, 3.77289445005425*^9}, {
3.772894506837974*^9, 3.7728945083327723`*^9}, {3.7728977013022947`*^9,
3.772897701484643*^9}, {3.7729006807319736`*^9, 3.772900684530842*^9}, {
3.7729011626245613`*^9, 3.772901243792779*^9}, {3.772901383037828*^9,
3.772901408382078*^9}, {3.7729024675914516`*^9, 3.772902491096961*^9}, {
3.77290252588194*^9, 3.772902536277769*^9}, {3.7729026451903753`*^9,
3.77290265894727*^9}, {3.772902701468636*^9, 3.772902720238291*^9}, {
3.772902878106348*^9, 3.772902878247653*^9}, {3.772906747128436*^9,
3.7729067474866247`*^9}, {3.772909090098812*^9, 3.7729090902571583`*^9}, {
3.772909755027151*^9, 3.772909817079022*^9}, 3.7729099096756363`*^9,
3.772913300950192*^9, {3.772913699416568*^9, 3.772913700055944*^9}, {
3.772919146570347*^9, 3.7729191681032124`*^9}, 3.7729234337867403`*^9, {
3.8128202519311323`*^9, 3.8128202921832113`*^9}, {3.8128204066764402`*^9,
3.812820460232698*^9}, {3.81282051280145*^9, 3.8128205131293087`*^9}, {
3.81282056839869*^9, 3.812820572965818*^9}, {3.8128206036888227`*^9,
3.812820626544346*^9}, {3.8128206653704023`*^9, 3.812820670793482*^9}, {
3.812820708040132*^9, 3.8128207898104343`*^9}, {3.812820858326275*^9,
3.812820968811678*^9}, {3.812826130091071*^9, 3.8128261554157953`*^9}, {
3.812881136400981*^9, 3.812881177846438*^9}, {3.8128812201695757`*^9,
3.8128812348364353`*^9}, {3.812882567308715*^9, 3.812882568728355*^9}, {
3.812916088127894*^9, 3.812916088581066*^9}, {3.812927565973686*^9,
3.812927566368517*^9}, {3.8129306505014057`*^9, 3.812930652538053*^9}, {
3.813060446973151*^9, 3.813060447752293*^9}, {3.813930719355957*^9,
3.813930720162648*^9}, {3.814017150987935*^9, 3.8140171535413322`*^9}},
CellLabel->"In[1]:=",ExpressionUUID->"aa0f6726-1167-4f75-a1b6-d3e71c46e215"],
Cell["\<\
Next, we write down the Green's function for each dimensionality; this \
calculation assumes a spherically symmetric environment. The Green's function \
describes the concentration generated at radius r by a spherical (circular \
ring for sources in 2D, point for sources in 1D) shell at radius R; the shell \
is assumed to have made a momentary emission at time tt and the concentration \
is measured at a later time t = tt+T. The diffusion constant is Dc. GF is the \
desired Green\[CloseCurlyQuote]s function given diffDim and cellDim. gradGF \
is the radial gradient of these functions\
\>", "Text",
CellChangeTimes->{{3.812819861190261*^9, 3.812819974453225*^9}, {
3.812821016655105*^9, 3.812821018303409*^9}, {3.812823689764305*^9,
3.81282369873234*^9}, {3.812826184405582*^9, 3.81282619445477*^9}, {
3.812927913060802*^9, 3.8129279245019703`*^9}, {3.814095272888217*^9,
3.814095287951631*^9}},ExpressionUUID->"f9c5f87b-a4fe-4512-b47f-\
8a0fd2619f52"],
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"GF11", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=",
RowBox[{
FractionBox["1",
RowBox[{"Sqrt", "[",
RowBox[{"4", " ", "Pi", " ", "Dc", " ", "T"}], "]"}]],
RowBox[{"Exp", "[",
RowBox[{"-",
FractionBox[
SuperscriptBox[
RowBox[{"(",
RowBox[{"r", "-", "R"}], ")"}], "2"],
RowBox[{"4", " ", "Dc", " ", "T"}]]}], "]"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"GF12", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox["1",
RowBox[{"2", " ", "Pi", " ", "Dc", " ", "T"}]],
RowBox[{"Exp", "[",
RowBox[{"-",
FractionBox[
SuperscriptBox[
RowBox[{"(",
RowBox[{"r", "-", "R"}], ")"}], "2"],
RowBox[{"4", " ", "Dc", " ", "T"}]]}], "]"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"GF22", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox["R",
RowBox[{"2", "Dc", " ", "T"}]],
RowBox[{"Exp", "[",
RowBox[{"-",
FractionBox[
RowBox[{
SuperscriptBox["r", "2"], "+",
SuperscriptBox["R", "2"]}],
RowBox[{"4", " ", "Dc", " ", "T"}]]}], "]"}],
RowBox[{"BesselI", "[",
RowBox[{"0", ",",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]]}], "]"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"GF23", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox["R",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{"4", "Pi"}], "]"}], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{"Dc", " ", "T"}], ")"}],
RowBox[{"3", "/", "2"}]]}]],
RowBox[{"Exp", "[",
RowBox[{"-",
FractionBox[
RowBox[{
SuperscriptBox["r", "2"], "+",
SuperscriptBox["R", "2"]}],
RowBox[{"4", " ", "Dc", " ", "T"}]]}], "]"}],
RowBox[{"BesselI", "[",
RowBox[{"0", ",",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]]}], "]"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"GF33", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox[
RowBox[{"R", " "}],
RowBox[{"r", " ",
RowBox[{"Sqrt", "[",
RowBox[{"Pi", " ", "Dc", " ", "T"}], "]"}]}]],
RowBox[{"Sinh", "[",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]], "]"}],
RowBox[{"Exp", "[",
RowBox[{"-",
FractionBox[
RowBox[{
SuperscriptBox["r", "2"], "+",
SuperscriptBox["R", "2"]}],
RowBox[{"4", " ", "Dc", " ", "T"}]]}], "]"}]}]}], ";"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"GF", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], ":=",
RowBox[{
RowBox[{"ToExpression", "[",
RowBox[{"\"\<GF\>\"", "<>",
RowBox[{"ToString", "[", "cellDim", "]"}], "<>",
RowBox[{"ToString", "[", "diffDim", "]"}]}], "]"}], "[",
RowBox[{"r", ",", "R", ",", " ", "T"}], "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"gradGF11", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=",
RowBox[{
FractionBox[
RowBox[{"(",
RowBox[{"R", "-", "r"}], ")"}],
RowBox[{"4", " ",
SqrtBox["\[Pi]"], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{"Dc", " ", "T"}], ")"}],
RowBox[{"3", "/", "2"}]]}]],
SuperscriptBox["\[ExponentialE]",
RowBox[{"-",
FractionBox[
SuperscriptBox[
RowBox[{"(",
RowBox[{"r", "-", "R"}], ")"}], "2"],
RowBox[{"4", " ", "Dc", " ", "T"}]]}]]}]}], " ",
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"gradGF12", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox[
RowBox[{"(",
RowBox[{"R", "-", "r"}], ")"}],
RowBox[{"4", " ",
SuperscriptBox["Dc", "2"], " ", "\[Pi]", " ",
SuperscriptBox["T", "2"]}]],
SuperscriptBox["\[ExponentialE]",
RowBox[{"-",
FractionBox[
SuperscriptBox[
RowBox[{"(",
RowBox[{"r", "-", "R"}], ")"}], "2"],
RowBox[{"4", " ", "Dc", " ", "T"}]]}]]}]}], " ",
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"gradGF22", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox["R",
RowBox[{"4", " ",
SuperscriptBox["Dc", "2"], " ",
SuperscriptBox["T", "2"]}]],
SuperscriptBox["\[ExponentialE]",
RowBox[{"-",
FractionBox[
RowBox[{
SuperscriptBox["r", "2"], "+",
SuperscriptBox["R", "2"]}],
RowBox[{"4", " ", "Dc", " ", "T"}]]}]], " ",
RowBox[{"(",
RowBox[{
RowBox[{"R", " ",
RowBox[{"BesselI", "[",
RowBox[{"1", ",",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]]}], "]"}]}], "-",
RowBox[{"r", " ",
RowBox[{"BesselI", "[",
RowBox[{"0", ",",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]]}], "]"}]}]}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"gradGF23", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox["R",
RowBox[{"4", " ",
SqrtBox["\[Pi]"], " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{"Dc", " ", "T"}], ")"}],
RowBox[{"5", "/", "2"}]]}]],
SuperscriptBox["\[ExponentialE]",
RowBox[{"-",
FractionBox[
RowBox[{
SuperscriptBox["r", "2"], "+",
SuperscriptBox["R", "2"]}],
RowBox[{"4", " ", "Dc", " ", "T"}]]}]],
RowBox[{"(",
RowBox[{
RowBox[{"R", " ",
RowBox[{"BesselI", "[",
RowBox[{"1", ",",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]]}], "]"}]}], "-",
RowBox[{"r", " ",
RowBox[{"BesselI", "[",
RowBox[{"0", ",",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]]}], "]"}]}]}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"gradGF33", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], " ", ":=", " ",
RowBox[{
FractionBox[
RowBox[{
SuperscriptBox["R", "2"], "/", "r"}],
RowBox[{"2", " ",
SqrtBox["\[Pi]"],
SuperscriptBox[
RowBox[{"(",
RowBox[{"Dc", " ", "T"}], ")"}],
RowBox[{"3", "/", "2"}]]}]],
SuperscriptBox["\[ExponentialE]",
RowBox[{"-",
FractionBox[
RowBox[{
SuperscriptBox["r", "2"], "+",
SuperscriptBox["R", "2"]}],
RowBox[{"4", " ", "Dc", " ", "T"}]]}]],
RowBox[{"(",
RowBox[{
RowBox[{"Cosh", "[",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]], "]"}], "-",
RowBox[{
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{
SuperscriptBox["r", "2"], "+",
RowBox[{"2", " ", "Dc", " ", "T"}]}], ")"}], " "}],
RowBox[{"r", " ", "R"}]],
RowBox[{"Sinh", "[",
FractionBox[
RowBox[{"r", " ", "R"}],
RowBox[{"2", " ", "Dc", " ", "T"}]], "]"}]}]}], ")"}]}]}], ";"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"gradGF", "[",
RowBox[{"r_", ",", "R_", ",", " ", "T_"}], "]"}], ":=",
RowBox[{
RowBox[{"ToExpression", "[",
RowBox[{"\"\<gradGF\>\"", "<>",
RowBox[{"ToString", "[", "cellDim", "]"}], "<>",
RowBox[{"ToString", "[", "diffDim", "]"}]}], "]"}], "[",
RowBox[{"r", ",", "R", ",", " ", "T"}], "]"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.812819529994217*^9, 3.8128195956339893`*^9}, {
3.812819731398355*^9, 3.8128198508177643`*^9}, {3.8128199774568644`*^9,
3.8128200328102493`*^9}, {3.81282101977234*^9, 3.8128210311758423`*^9}, {
3.8128235174479094`*^9, 3.812823563828326*^9}, {3.812823603565918*^9,
3.812823686663846*^9}, {3.812825943849143*^9, 3.812825985601365*^9}, {
3.812916912412298*^9, 3.812916916907324*^9}, 3.812916948341258*^9,
3.812917024770651*^9, {3.8129175276207037`*^9, 3.812917530652295*^9}, {
3.81291757681463*^9, 3.812917580638914*^9}, {3.812927926414612*^9,
3.812927944437067*^9}, {3.8129281955438747`*^9, 3.812928499887857*^9}, {
3.814015830179119*^9, 3.814016020709448*^9}},
CellLabel->"In[7]:=",ExpressionUUID->"1d5411fc-c466-4794-844c-f005c3e23b2f"],
Cell["\<\
Since our goal is to calculate the wave front dynamics in our system, we must \
compute the concentration profile as it evolves in time. For Heaviside \
activation, computing the wave front position will directly allow for the \
computation of the concentration profile, since the latter involves an \
integral of the Green\[CloseCurlyQuote]s function against the former. To \
start, we define the numerical precision and desired duration of this \
computation.\
\>", "Text",
CellChangeTimes->{{3.8128210529240923`*^9, 3.812821053805209*^9}, {
3.812821103756871*^9, 3.812821115811273*^9}, {3.8128231724050713`*^9,
3.812823176034902*^9}, {3.8128232075594263`*^9,
3.812823360164195*^9}},ExpressionUUID->"28ba037a-8215-4606-9435-\
94131c549906"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
RowBox[{
RowBox[{"Duration", " ", "of", " ", "the", " ", "computation"}], ";",
" ",
RowBox[{
"this", " ", "should", " ", "ideally", " ", "be", " ", "many", " ",
"multiples", " ", "of", " ", "the", " ", "characteristic", " ",
"time"}]}], ",", " ",
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"tMax", " ", "=",
RowBox[{"10", " ",
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{
RowBox[{
"Time", " ", "step", " ", "of", " ", "the", " ", "computation"}], ";",
" ",
RowBox[{
"this", " ", "should", " ", "ideally", " ", "be", " ", "a", " ",
"small", " ", "fraction", " ", "of", " ", "the", " ", "characteristic",
" ", "time"}]}], ",", " ",
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[CapitalDelta]t", " ", "=", " ",
RowBox[{
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}], "/", "5"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"Print", "[",
RowBox[{"\"\<Number of simulation steps: \>\"", "<>",
RowBox[{"ToString", "[",
RowBox[{"tMax", "/", "\[CapitalDelta]t"}], "]"}]}], "]"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{
"Spatial", " ", "precision", " ", "of", " ", "the", " ", "wave", " ",
"front", " ", "finder"}], ";", " ",
RowBox[{
"this", " ", "should", " ", "be", " ", "a", " ", "small", " ",
"fraction", " ", "of", " ", "the", " ", "distance", " ", "travelled",
" ", "by", " ", "the", " ", "wave", " ", "front", " ", "in", " ", "one",
" ", "time", " ", "step"}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Delta]r", " ", "=", " ",
RowBox[{"vBase", " ",
RowBox[{"\[CapitalDelta]t", "/", "100"}]}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{
"Define", " ", "a", " ", "maximum", " ", "step", " ", "size", " ", "of",
" ", "the", " ", "wave", " ", "front", " ", "finder"}], ";", " ",
RowBox[{
"this", " ", "needs", " ", "only", " ", "to", " ", "be", " ", "a", " ",
"few", " ", "multiples", " ", "of", " ", "the", " ", "distance", " ",
"travelled", " ", "by", " ", "the", " ", "wave", " ", "front", " ",
"in", " ", "one", " ", "time", " ", "step"}]}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"\[CapitalDelta]r", " ", "=", " ",
RowBox[{"N", "[",
RowBox[{"2", "vBase", " ", "\[CapitalDelta]t"}], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"Define", " ", "the", " ", "number", " ", "of", " ", "steps", " ",
"needed", " ", "to", " ", "achieve", " ", "a", " ", "maximum", " ",
"step", " ", "size", " ", "of", " ", "\[CapitalDelta]r", " ", "and", " ",
"a", " ", "precision", " ", "of", " ", "\[Delta]r"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"searchNum", " ", "=", " ",
RowBox[{"Ceiling", "[",
RowBox[{"Log2", "[",
RowBox[{"\[CapitalDelta]r", "/", "\[Delta]r"}], "]"}], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"Define", " ", "the", " ", "step", " ", "sizes"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"rSteps", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{
FractionBox["\[CapitalDelta]r",
SuperscriptBox["2", "i"]], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "searchNum"}], "}"}]}], "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"How", " ", "many", " ", "time", " ", "steps", " ", "do", " ", "you", " ",
"want", " ", "to", " ", "take", " ", "before", " ", "the", " ",
"solver", " ", "updates", " ", "its", " ",
RowBox[{"progress", "?"}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"updateNum", " ", "=", " ", "10"}], ";"}]}]}]], "Input",
CellChangeTimes->{{3.812823774531157*^9, 3.812823802344368*^9}, {
3.812824273665534*^9, 3.812824704172463*^9}, {3.81282478153769*^9,
3.812824789064204*^9}, {3.812824838505707*^9, 3.812824937616825*^9}, {
3.812881694326172*^9, 3.8128816961829576`*^9}, {3.8128825040811377`*^9,
3.8128825247548323`*^9}, 3.8128828490533123`*^9, 3.812892123254572*^9, {
3.812893836667419*^9, 3.812893836710829*^9}, 3.812909332668779*^9, {
3.812910245364442*^9, 3.812910248537593*^9}, 3.8129110114150267`*^9, {
3.812911042666991*^9, 3.812911044024642*^9}, {3.8129123477333107`*^9,
3.812912348317111*^9}, 3.812913416887574*^9, {3.8129275326519337`*^9,
3.812927538963936*^9}, {3.81306802704496*^9, 3.8130680276457453`*^9}},
CellLabel->"In[19]:=",ExpressionUUID->"d9a87681-af07-467e-9402-8a19d89d974c"],
Cell[BoxData["\<\"Number of simulation steps: 50.\"\>"], "Print",
CellChangeTimes->{{3.812824697653908*^9, 3.8128247046384153`*^9},
3.812824941847777*^9, {3.812881142295257*^9, 3.812881239400709*^9},
3.812881640527275*^9, {3.812882506306921*^9, 3.8128825269781857`*^9},
3.812882849426075*^9, {3.812892118739249*^9, 3.8128921238891077`*^9},
3.8128923212312117`*^9, 3.812893628207087*^9, 3.812893837651287*^9, {
3.8129090473978977`*^9, 3.81290906390023*^9}, 3.81290914121111*^9,
3.812909185436738*^9, 3.812909370058434*^9, 3.812910249190072*^9, {
3.8129123420022078`*^9, 3.8129123575971603`*^9}, 3.812912414954829*^9,
3.812913424664153*^9, 3.812915685898703*^9, 3.812916094266293*^9,
3.8129165620567837`*^9, 3.8129169535260067`*^9, 3.812917033432378*^9,
3.812917582871915*^9, {3.812927548278817*^9, 3.8129275716833563`*^9},
3.8129286945830603`*^9, 3.8129306655024633`*^9, 3.81306019180785*^9,
3.813060452593286*^9, 3.813060571424655*^9, 3.813060726516947*^9,
3.81306094283323*^9, 3.813061082757299*^9, 3.813068028111765*^9,
3.8130697784020767`*^9, 3.813930725335512*^9, 3.814016038201579*^9,
3.814016910338563*^9, 3.8140172648658333`*^9, 3.8140961298748207`*^9,
3.814186404220089*^9},
CellLabel->
"During evaluation of \
In[19]:=",ExpressionUUID->"bbddf319-3c30-43f2-9362-f3ee0699935d"]
}, Open ]],
Cell["\<\
Next, we calculate the initiation time of the wave by finding the time at \
which the concentration at rInv exceeds Cth. To do so, we step halfway to \
tMax and check if the concentration exceeds Cth. If it does, we step a \
quarter way; if it doesn\[CloseCurlyQuote]t we step three quarters of the \
way. Repeatedly dividing the time interval in half, we eventually converge to \
the initiation time tInit up to a precision of \[CapitalDelta]t. In our \
numerical integration, we use the LocalAdaptive option to handle the \
singularity around T=0. Thusly, we have defined a function, initTime, which \
takes as arguments rInv, tMax, and \[CapitalDelta]t.\
\>", "Text",
CellChangeTimes->{{3.8128247314659243`*^9, 3.8128247399456987`*^9}, {
3.812880373037813*^9, 3.8128804392167263`*^9}, {3.812880496952327*^9,
3.812880497200509*^9}, {3.81288137916747*^9, 3.8128814902534313`*^9}, {
3.812881701316517*^9, 3.812881723885068*^9}, {3.812893290687099*^9,
3.812893291768963*^9}, {3.8140953074093723`*^9, 3.8140953083072166`*^9}, {
3.814095874148307*^9,
3.81409589341568*^9}},ExpressionUUID->"e17f4ff8-50db-4563-944d-\
7608557cef97"],
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
"How", " ", "many", " ", "time", " ", "steps", " ", "are", " ", "needed",
" ", "to", " ", "divide", " ", "the", " ", "interval", " ", "of", " ",
"length", " ", "tMax", " ", "to", " ", "resolution", " ",
RowBox[{"\[CapitalDelta]t", "?"}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"nTimeSteps", " ", "=", " ",
RowBox[{"Ceiling", "[",
RowBox[{"Log2", "[",
FractionBox["tMax", "\[CapitalDelta]t"], "]"}], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"Calculate", " ", "the", " ", "time", " ", "steps"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"timeSteps", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{
FractionBox["tMax",
SuperscriptBox["2", "i"]], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "nTimeSteps"}], "}"}]}], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"An", " ", "array", " ", "that", " ", "keeps", " ", "track", " ", "of",
" ", "how", " ", "far", " ", "to", " ", "walk"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"timeBools", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{"1", ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "nTimeSteps"}], "}"}]}], "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{"recMax", " ", "=", " ", "15"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"nIntOp", " ", "=", " ",
RowBox[{"{",
RowBox[{"\"\<LocalAdaptive\>\"", ",",
RowBox[{"\"\<MaxRecursion\>\"", "\[Rule]", "recMax"}]}], "}"}]}],
";"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"initTime", "[",
RowBox[{"rInv_", ",", "tMax_", ",", "\[CapitalDelta]t_"}], "]"}], ":=",
RowBox[{"Do", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"numSteps", " ", "=", " ",
RowBox[{"Ceiling", "[",
RowBox[{"Log2", "[",
RowBox[{"tMax", "/", "\[CapitalDelta]t"}], "]"}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"tSteps", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{
FractionBox["tMax",
SuperscriptBox["2", "i"]], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "numSteps"}], "}"}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"tBools", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{"1", ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "numSteps"}], "}"}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"i", "=", "1"}], ",",
RowBox[{"i", "<",
RowBox[{"numSteps", "+", "1"}]}], ",",
RowBox[{"i", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"test", " ", "=", " ",
RowBox[{"\[CapitalDelta]t", " ",
RowBox[{"Ceiling", "[",
RowBox[{
FractionBox["1", "\[CapitalDelta]t"],
RowBox[{"Total", "[",
RowBox[{
RowBox[{"(",
RowBox[{"tSteps", " ", "tBools"}], ")"}], "[",
RowBox[{"[",
RowBox[{"1", ";;", "i"}], "]"}], "]"}], "]"}]}], "]"}]}]}],
";", "\[IndentingNewLine]",
RowBox[{"cTest", " ", "=", " ",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{"GF", "[",
RowBox[{"rInv", ",", "R", ",", "tt"}], "]"}], ",",
RowBox[{"{",
RowBox[{"R", ",",
RowBox[{
RowBox[{"-", "rInv"}], "+",
RowBox[{"rInv", " ",
RowBox[{"HeavisideTheta", "[",
RowBox[{"cellDim", "-", "1.5"}], "]"}]}]}], ",", "rInv"}],
"}"}], ",",
RowBox[{"{",
RowBox[{"tt", ",", "0", ",", "test"}], "}"}], ",",
RowBox[{"Method", "\[Rule]", "nIntOp"}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"If", "[",
RowBox[{
RowBox[{"cTest", ">", "Cth"}], ",",
RowBox[{
RowBox[{"tBools", "[",
RowBox[{"[", "i", "]"}], "]"}], "=", "0"}]}], "]"}], ";",
"\[IndentingNewLine]",
RowBox[{"If", "[",
RowBox[{
RowBox[{"i", "\[Equal]", "numSteps"}], ",",
RowBox[{"Return", "[",
RowBox[{"\[CapitalDelta]t", " ",
RowBox[{"Ceiling", "[",
RowBox[{
FractionBox["1", "\[CapitalDelta]t"],
RowBox[{"Total", "[",
RowBox[{"tSteps", " ",
RowBox[{"tBools", "[",
RowBox[{"[",
RowBox[{"1", ";;", "i"}], "]"}], "]"}]}], "]"}]}],
"]"}]}], "]"}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}],
";"}], ",", "\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"n", ",", "1"}], "}"}]}], "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"tInit", " ", "=", " ",
RowBox[{"initTime", "[",
RowBox[{"rInv", ",", " ", "tMax", ",", " ", "\[CapitalDelta]t"}],
"]"}]}]}]}]], "Input",
CellGroupingRules->{GroupTogetherGrouping, 10000.},
CellChangeTimes->{{3.756582667444396*^9, 3.756582723920148*^9}, {
3.7565827689724283`*^9, 3.756582772539909*^9}, {3.756582827076426*^9,
3.75658282735647*^9}, {3.756583550440131*^9, 3.75658358359859*^9},
3.756588407356176*^9, {3.7565966697881413`*^9, 3.7565966709780207`*^9}, {
3.756644116333273*^9, 3.756644129839808*^9}, {3.7566475025949087`*^9,
3.756647584798318*^9}, {3.7566476333877974`*^9, 3.75664768062906*^9}, {
3.756647840470337*^9, 3.7566478658014803`*^9}, 3.756676413234355*^9, {
3.757950991352098*^9, 3.75795100607218*^9}, {3.757951111992612*^9,
3.757951113220964*^9}, {3.758069830220339*^9, 3.758069840234666*^9}, {
3.758072243913925*^9, 3.758072250720022*^9}, {3.7581160706559887`*^9,
3.758116071247487*^9}, {3.760988919450873*^9, 3.7609889199320507`*^9}, {
3.761502102491959*^9, 3.761502103055112*^9}, {3.772880908031684*^9,
3.772880927952507*^9}, {3.772880958319522*^9, 3.772880964197537*^9}, {
3.7728810465165043`*^9, 3.7728810820841208`*^9}, {3.772881134324651*^9,
3.772881220432993*^9}, {3.7728813457018547`*^9, 3.772881351028648*^9}, {
3.77288142565316*^9, 3.772881464222458*^9}, {3.77289549624722*^9,
3.772895498212614*^9}, 3.7729234337876463`*^9, {3.81288040324541*^9,
3.8128804051293163`*^9}, {3.812880460142144*^9, 3.812880461478045*^9}, {
3.8128807928098583`*^9, 3.8128808100866327`*^9}, {3.812880954097247*^9,
3.812880981476653*^9}, {3.812881049133359*^9, 3.812881100927861*^9}, {
3.812881336248263*^9, 3.812881368184155*^9}, {3.8128815111694517`*^9,
3.812881676631277*^9}, {3.8128817349626083`*^9, 3.812881740809799*^9}, {
3.812881875630233*^9, 3.8128818760127068`*^9}, {3.812882314816566*^9,
3.812882469306391*^9}, {3.812892139353924*^9, 3.812892184810541*^9}, {
3.8128922913355713`*^9, 3.8128923090744543`*^9}, {3.812909169708407*^9,
3.812909172839779*^9}, {3.8130604339578133`*^9, 3.813060473645567*^9}, {
3.8130605289269*^9, 3.813060554184684*^9}, {3.813060720597538*^9,
3.813060721273993*^9}, {3.813060816380041*^9, 3.813060817576763*^9},
3.813060936138639*^9, {3.813060996010458*^9, 3.813060996145543*^9}, {
3.813061053695283*^9, 3.813061101932407*^9}, {3.813068003592874*^9,
3.813068013685808*^9}, {3.813069781810883*^9, 3.813069782863886*^9}, {
3.814016044590507*^9, 3.814016045450076*^9}, {3.814095849899873*^9,
3.81409586515805*^9}, {3.814096149246613*^9, 3.814096149413217*^9}, {
3.8141852128232317`*^9, 3.814185218244156*^9}},
CellLabel->"In[27]:=",ExpressionUUID->"1f0cee4e-a7c2-46d6-bede-6f85096a6f5c"],
Cell[BoxData["60.`"], "Output",
CellChangeTimes->{3.814096152638624*^9, 3.814186407802071*^9},
CellLabel->"Out[33]=",ExpressionUUID->"ad77fcab-354a-4dc8-ae23-ea34b9556d85"],
Cell["\<\
We then compute the location of the wave front at every subsequent time point \
up to tMax. To do so, we keep track of the wave front location and construct \
a linear interpolation of the wave front position, then integrate this \
against the Green\[CloseCurlyQuote]s function.\
\>", "Text",
CellChangeTimes->{{3.81288174649824*^9, 3.812881764990966*^9}, {
3.812881800761746*^9,
3.81288187245322*^9}},ExpressionUUID->"b026b330-3d2e-4447-b0e2-\
fcae81df7333"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
"An", " ", "array", " ", "that", " ", "keeps", " ", "track", " ", "of",
" ", "the", " ", "wave", " ", "front", " ", "location"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"nSteps", " ", "=", " ",
RowBox[{"Ceiling", "[",
FractionBox[
RowBox[{"tMax", "-", "tInit"}], "\[CapitalDelta]t"], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"rCrit", " ", "=", " ",
RowBox[{"ConstantArray", "[",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{"nSteps", "+", "2"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"rCrit", "[",
RowBox[{"[", "1", "]"}], "]"}], " ", "=", " ",
RowBox[{"{",
RowBox[{"0", ",", "rInv"}], "}"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"rCrit", "[",
RowBox[{"[", "2", "]"}], "]"}], "=",
RowBox[{"{",
RowBox[{"tInit", ",", "rInv"}], "}"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"Compute", " ", "the", " ", "wave", " ", "front", " ", "locations", " ",
"in", " ", "time"}], "*)"}], "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"i", "=", "1"}], ",",
RowBox[{"i", "<",
RowBox[{"nSteps", "+", "1"}]}], ",",
RowBox[{"i", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"t", " ", "=", " ",
RowBox[{"tInit", "+",
RowBox[{"i", " ", "\[CapitalDelta]t"}]}]}], ";",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"rStart", " ", "=", " ",
RowBox[{"rCrit", "[",
RowBox[{"[",
RowBox[{
RowBox[{"i", "+", "1"}], ",", "2"}], "]"}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"rBools", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{"1", ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "searchNum"}], "}"}]}], "]"}]}], ";",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"j", " ", "=", " ", "1"}], ",",
RowBox[{"j", "<",
RowBox[{"searchNum", "+", "1"}]}], ",",
RowBox[{"j", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"rTest", " ", "=", " ",
RowBox[{"rStart", "+",
RowBox[{"Total", "[",
RowBox[{
RowBox[{"(",
RowBox[{"rBools", " ", "rSteps"}], ")"}], "[",
RowBox[{"[",
RowBox[{"1", ";;", "j"}], "]"}], "]"}], "]"}]}]}], ";",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"fTest", " ", "=", " ",
RowBox[{"Interpolation", "[",
RowBox[{
RowBox[{"Append", "[",
RowBox[{
RowBox[{"rCrit", "[",
RowBox[{"[",
RowBox[{"1", ";;",
RowBox[{"i", "+", "1"}]}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "rTest"}], "}"}]}], "]"}], ",",
RowBox[{"InterpolationOrder", "\[Rule]", "1"}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"concTest", " ", "=", " ",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{"GF", "[",
RowBox[{"rTest", ",", "R", ",",
RowBox[{"t", "-", "tt"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"tt", ",", "0", ",", "t"}], "}"}], ",",
RowBox[{"{",
RowBox[{"R", ",", "0", ",",
RowBox[{"fTest", "[", "tt", "]"}]}], "}"}], ",",
RowBox[{"Method", "\[Rule]", "nIntOp"}]}], "]"}]}], ";",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"If", "[",
RowBox[{
RowBox[{"concTest", "<", "Cth"}], ",",
RowBox[{
RowBox[{"rBools", "[",
RowBox[{"[", "j", "]"}], "]"}], "=", "0"}]}], "]"}], ";",
"\[IndentingNewLine]",
RowBox[{"If", "[",
RowBox[{
RowBox[{"j", "\[Equal]", "searchNum"}], ",",
RowBox[{
RowBox[{"rCrit", "[",
RowBox[{"[",
RowBox[{"i", "+", "2"}], "]"}], "]"}], "=",
RowBox[{"{",
RowBox[{"t", ",",
RowBox[{"rStart", "+",
RowBox[{"Total", "[",
RowBox[{
RowBox[{"(",
RowBox[{"rBools", " ", "rSteps"}], ")"}], "[",
RowBox[{"[",
RowBox[{"1", ";;", "j"}], "]"}], "]"}], "]"}]}]}],
"}"}]}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], ";",
"\[IndentingNewLine]",
RowBox[{"If", "[",
RowBox[{
RowBox[{
RowBox[{"Mod", "[",
RowBox[{"i", ",", "updateNum"}], "]"}], "\[Equal]", "0"}], ",",
RowBox[{"Print", "[",
RowBox[{
RowBox[{"ToString", "[",
RowBox[{"100.",
RowBox[{"i", "/", "nSteps"}]}], "]"}], "<>",
"\"\<% complete\>\""}], "]"}]}], "]"}]}]}], "\[IndentingNewLine]",
"]"}]}]}]], "Input",
CellChangeTimes->{{3.812881009877967*^9, 3.812881025681052*^9}, {
3.8128817776451073`*^9, 3.812881788174006*^9}, {3.812881879311839*^9,
3.812881908079247*^9}, {3.812881946700988*^9, 3.812882010154969*^9}, {
3.81288205378146*^9, 3.81288209650896*^9}, {3.8128916427285767`*^9,
3.812891680945418*^9}, {3.812891726992648*^9, 3.8128918739994717`*^9}, {
3.812891926323234*^9, 3.812892097641461*^9}, {3.812893100661244*^9,
3.8128931613900747`*^9}, {3.8128934683188972`*^9,
3.8128934686710587`*^9}, {3.812893610881176*^9, 3.812893613238139*^9}, {
3.812893727853513*^9, 3.812893776077338*^9}, {3.812893810973349*^9,
3.812893817817876*^9}, 3.812909029122549*^9, {3.812909128539427*^9,
3.812909130287129*^9}, {3.8129107806618423`*^9, 3.812910783734931*^9}, {
3.8129108455440493`*^9, 3.81291084773763*^9}, {3.812910998646398*^9,
3.812911006352943*^9}, {3.81291246717484*^9, 3.8129124982029123`*^9}, {
3.812913438018888*^9, 3.81291343984125*^9}, 3.813060565358432*^9, {
3.81306084516887*^9, 3.813060917684245*^9}, {3.81306099261185*^9,
3.813061062205534*^9}, 3.813068036149791*^9, {3.8140160629364843`*^9,
3.8140160632910557`*^9}, {3.814185227669751*^9, 3.814185284190486*^9}},
CellLabel->"In[34]:=",ExpressionUUID->"84ac632a-74dc-4374-87e8-ba9f2924b487"],
Cell[CellGroupData[{
Cell[BoxData["\<\"26.3158% complete\"\>"], "Print",
CellChangeTimes->{3.8130609205719767`*^9, 3.813060950692954*^9,
3.813061001451462*^9, 3.813061070092955*^9, 3.813061154366481*^9,
3.8130680885773993`*^9, 3.81306983753084*^9, 3.813930739180245*^9,
3.8140160837797127`*^9, 3.814016923992962*^9, 3.814017324359872*^9,
3.814096217116522*^9, 3.814186454144249*^9},
CellLabel->
"During evaluation of \
In[34]:=",ExpressionUUID->"b655968e-f828-45b9-9f09-a3ad4e132d7a"],
Cell[BoxData["\<\"52.6316% complete\"\>"], "Print",
CellChangeTimes->{3.8130609205719767`*^9, 3.813060950692954*^9,
3.813061001451462*^9, 3.813061070092955*^9, 3.813061154366481*^9,
3.8130680885773993`*^9, 3.81306983753084*^9, 3.813930739180245*^9,
3.8140160837797127`*^9, 3.814016923992962*^9, 3.814017324359872*^9,
3.814096217116522*^9, 3.8141865133184443`*^9},
CellLabel->
"During evaluation of \
In[34]:=",ExpressionUUID->"64fb176c-beac-46d7-a4ff-09b6e44d8fc7"],
Cell[BoxData["\<\"78.9474% complete\"\>"], "Print",
CellChangeTimes->{3.8130609205719767`*^9, 3.813060950692954*^9,
3.813061001451462*^9, 3.813061070092955*^9, 3.813061154366481*^9,
3.8130680885773993`*^9, 3.81306983753084*^9, 3.813930739180245*^9,
3.8140160837797127`*^9, 3.814016923992962*^9, 3.814017324359872*^9,
3.814096217116522*^9, 3.814186571869341*^9},
CellLabel->
"During evaluation of \
In[34]:=",ExpressionUUID->"86dcb18b-968d-4fd2-a27a-143c9230089c"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
"Plot", " ", "the", " ", "wave", " ", "front", " ", "radius", " ", "in",
" ", "time"}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"ListPlot", "[",
RowBox[{"rCrit", ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",", " ",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<time\>\"", ",", " ", "\"\<wave front radius\>\""}],
"}"}]}]}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"Also", " ", "plot", " ", "the", " ", "rate", " ", "of", " ", "change",
" ", "of", " ", "the", " ", "wave", " ", "front", " ", "location"}],
"*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"wfTimes", " ", "=", " ",
RowBox[{
RowBox[{"Transpose", "[", "rCrit", "]"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"wfRad", " ", "=", " ",
RowBox[{
RowBox[{"Transpose", "[", "rCrit", "]"}], "[",
RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"wfSpeeds", " ", "=", " ",
FractionBox[
RowBox[{
RowBox[{"wfRad", "[",
RowBox[{"[",
RowBox[{"2", ";;"}], "]"}], "]"}], "-",
RowBox[{"Drop", "[",
RowBox[{"wfRad", ",",
RowBox[{"-", "1"}]}], "]"}]}], "\[CapitalDelta]t"]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"Show", "[",
RowBox[{
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{
RowBox[{"wfTimes", "[",
RowBox[{"[",
RowBox[{"2", ";;"}], "]"}], "]"}], ",", "wfSpeeds"}], "}"}],
"]"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",", " ",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<time\>\"", ",", " ", "\"\<d(wave front radius)/dt\>\""}],
"}"}]}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{"0", ",",
RowBox[{"1.1", "vBase"}]}], "}"}]}]}], "]"}], ",",
"\[IndentingNewLine]",
RowBox[{"Plot", "[",
RowBox[{"vBase", ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",",
RowBox[{"Last", "[", "wfTimes", "]"}]}], "}"}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Dashed"}]}], "]"}]}],
"]"}]}]}]], "Input",
CellChangeTimes->{{3.8128923509241858`*^9, 3.8128923549415903`*^9}, {
3.812893867391596*^9, 3.812893907120316*^9}, 3.8128939600509253`*^9, {
3.812907588388957*^9, 3.812907640772477*^9}, {3.8129080256356373`*^9,
3.8129081974160213`*^9}, {3.812926972474286*^9, 3.812927005932507*^9}, {
3.8129308465054283`*^9, 3.8129308517469797`*^9}, {3.814186425538015*^9,
3.814186441350005*^9}},
CellLabel->"In[44]:=",ExpressionUUID->"36cd9677-5a82-481e-ba8a-cd3be4e7b70d"],
Cell[BoxData[
GraphicsBox[{{}, {{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJxTTMoPSmViYGDQAGIQDQNJRqEGWiuV7SE8PwcUvkOAw0JRiU27rkL5DYEO
Lwzi4xpqVSD8A0EO2sJ1K+ebq0L1hzhMaBZKu/UPyncIdbDx3JX26awaVH+Y
g0FK6ernS9Sh+sMd/m2fJCvdqAHVH+mwzEVob0+SJlR/lIMG15dOG08tqP5o
h1YOixYtU22o/hiHPfb3Voep6UD1xzkkTDZ9pSilC9Uf71DIbGffI6AH5Sc4
xKxZcdCERx/Cf5DgcMsiSdeO1QBqfqLD3pUerRtboXyFJIcFUvp3VdgMofYl
OXxYfszbvBPKT0h2KDu8zXUmrxHU/hSHP/eW2f+eCOUvSHHY2fH+yBpxY6j9
qQ5rD0y4unYelP8g1eF6W1fhPzUTqP1pDoy+Lbz+66F8hXSHpIXXTN5bmELt
T3c4/FUjpv8QlJ+Q4RB14zJfva8Z1P5MBw3P6z77bkL5CzId9ms3/V+bYg61
P8vhsNa9EzafoPwHWQ4TZz9kq6mzgNqf7XBGY1p0FK8l1P4chwvqn1RezoLy
D+Q4cG9dVnlc0wpqf66Dm1P0uaU7oHyGPIftsRnrLntYQ+3Pc/j0SoZ503Uo
3yHfYdaUpjz9TBt7AO/Yr5c=
"]]}}, {{}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{
FormBox["\"time\"", TraditionalForm],
FormBox["\"wave front radius\"", TraditionalForm]},
AxesOrigin->{0, 0},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
Method->{
"OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
Identity[
Part[#, 1]],
Identity[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
Identity[
Part[#, 1]],
Identity[
Part[#, 2]]}& )}},
PlotRange->{{0, 250.}, {0, 0.000433515625}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.812893423242381*^9, 3.812894735230872*^9, {3.812907616075445*^9,
3.8129076426467333`*^9}, 3.812908197766658*^9, 3.812909280335478*^9,
3.8129101914759197`*^9, 3.812910473873308*^9, 3.812913022330679*^9,
3.8129139676777077`*^9, 3.8129280536432056`*^9, 3.812929178297533*^9, {
3.812930825834597*^9, 3.812930852930791*^9}, 3.813060351675219*^9,
3.813060590133526*^9, 3.813060963504059*^9, 3.813061012268002*^9,
3.813061359829157*^9, 3.81306829715047*^9, 3.813070033486019*^9,
3.813930773103009*^9, 3.814016118031665*^9, 3.814016956899363*^9,
3.814017490823669*^9, 3.8140963857867823`*^9, {3.814186629169458*^9,
3.814186629688572*^9}},
CellLabel->"Out[44]=",ExpressionUUID->"8ed4f077-da0f-4883-b769-6e6bb29e548f"],
Cell[BoxData[
GraphicsBox[{{{}, {{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[
0.012833333333333334`], AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJxdzUsoRGEYBuA/skFiQZzEUIhi4yRSemuilOvkkvsxjOMyRrIQSWdDYiGX
kmnKJEljM2HcJ8eGyFmJbERSkstCSixcz6f8vvr7e3q/ty/G3GGy+DDGEr7e
969PAdjfQRHkvF1BPh7L+rFSDOXB43zZc+tWTVCOc1c3HUu6WQkE/7OgnPtl
3SiFOCa5pj0r1C+DHH6b/dDvoX45tP1Hnzi/VepXQOx+3qkeIaMS9vi3vsmQ
NepXQTj9yDiaIqvVUAb8Xnyj1qlfC7Fm6CS6now63IgBK5lzv5aQfx3c2ZW0
oftSAtueKBxtJyv10CbCUhbdZIOZt2qGaIy8u0rb1C01QBNmDt57yKyRt7MR
9qeYBcFLhgXy4dxgGtui+xaw2QSLyUhWmngbZGg9LqNtkKz+s9QMuTg5dviA
zFogJrrZfOA23W/hc7RCY6kXu4WUX7by+0ob7Gce7/k42WDl91Urn0vtvJkN
ojvd8Xrye9/G5+gAG9rqDY3wZn0CjQS+2g==
"]]}}, {{}, {}}}, {{{}, {},
TagBox[
{RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6],
Opacity[1.], Dashing[{Small, Small}], LineBox[CompressedData["
1:eJxTTMoPSmViYGAwAWIQ/frTHx+3tKt2vW+3Lvh+7ICdnLWm+8cFzA4wvt+1
Ix/eLxCG8x0mNnO7lqjA+YKfwzZMTNaD8+12VwTWCpjA+QuK9629E2EJ59/z
vHvpt5sdnJ/qXfU26bwjnO+46zvL+TwXOJ97x0MG6YducL72R/c8mw5PON9q
p9QBWyYfOH+9Z7zqsXw/OH8vl8A/wc4AOF9qx7LbGw4HwvlKzb8Ld0kEw/lx
RpyT2tNDEP7T9F1yZH4onJ9v65+f8ykMzneZ2a7TohcB5x8RqlvE1hoJ55cX
6C7yuRoF5y8UFZ7CxB+DkL+xzU45ORah/9GE25sXxCHMv6jMxHEjHs7fkber
4PGUBDi/U6LuZpVRIpxvcHHBPPe7CP6H3/M1vrUkwfk5XrtWxygmI8JnuXH3
j1MI/imNE8k+mSlw/kdpS8YmllQ4/60s5+83axD8pc5f3vJ7pcH5RhV2thlf
EPzqC4tbKyanw/l3Zz2dlKWfAeeLlrT9mXgDwecsdIovLMuE8y/VxMq1KWbB
+Vt/hSyqOYngB89Jl9dOy4bzj/vc/MEnkAPnf1NIO71qHYJvn34+ScMrFxHf
e6aH3f6E4FuYVju7TchDuIczRUb5BYJ/647P1//2+XA+ANhM2e0=
"]]},
Annotation[#, "Charting`Private`Tag$20613#1"]& ]}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{
FormBox["\"time\"", TraditionalForm],
FormBox["\"d(wave front radius)/dt\"", TraditionalForm]},
AxesOrigin->{56.041666666666664`, 0},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
Method->{
"OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
Identity[
Part[#, 1]],
Identity[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
Identity[
Part[#, 1]],
Identity[
Part[#, 2]]}& )}},
PlotRange->{{56.041666666666664`, 250.}, {0, 2.2*^-6}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {0, 0}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.812893423242381*^9, 3.812894735230872*^9, {3.812907616075445*^9,
3.8129076426467333`*^9}, 3.812908197766658*^9, 3.812909280335478*^9,
3.8129101914759197`*^9, 3.812910473873308*^9, 3.812913022330679*^9,
3.8129139676777077`*^9, 3.8129280536432056`*^9, 3.812929178297533*^9, {
3.812930825834597*^9, 3.812930852930791*^9}, 3.813060351675219*^9,
3.813060590133526*^9, 3.813060963504059*^9, 3.813061012268002*^9,
3.813061359829157*^9, 3.81306829715047*^9, 3.813070033486019*^9,
3.813930773103009*^9, 3.814016118031665*^9, 3.814016956899363*^9,
3.814017490823669*^9, 3.8140963857867823`*^9, {3.814186629169458*^9,
3.8141866297211447`*^9}},
CellLabel->"Out[48]=",ExpressionUUID->"dbc7163f-6eb9-4dd9-aa56-adc784cf03c1"]
}, Open ]],
Cell["\<\
Next, we wish to visualize the concentration profile at the end of our \
numerical simulation, so we look at it from r = 0 to 1.5 times the largest \
wave front radius achieved.\
\>", "Text",
CellChangeTimes->{{3.81288174649824*^9, 3.812881764990966*^9}, {
3.812881800761746*^9, 3.81288187245322*^9}, {3.812907671618938*^9,
3.812907733150786*^9}, {3.812907875227005*^9,
3.812907875368112*^9}},ExpressionUUID->"bd73eded-ac7d-4683-b811-\
b6086b68897c"],
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
"Number", " ", "of", " ", "points", " ", "at", " ", "which", " ", "to",
" ", "evaluate", " ", "the", " ", "concentration"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"nConcPts", " ", "=", " ", "50"}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"Maximum", " ", "radius", " ", "to", " ", "observe"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"rMax", " ", "=", " ",
RowBox[{"1.5",
RowBox[{
RowBox[{"Last", "[", "rCrit", "]"}], "[",
RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"Radii", " ", "at", " ", "which", " ", "to", " ", "calculate", " ", "the",
" ", "concentration", " ", "profile"}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"rPts", " ", "=", " ",
RowBox[{"Subdivide", "[",
RowBox[{"rInv", ",", "rMax", ",", "nConcPts"}], "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"Time", " ", "at", " ", "which", " ", "we", " ", "wish", " ", "to", " ",
"observe", " ", "the", " ", "concentration"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"tEval", " ", "=", " ",
RowBox[{
RowBox[{"Last", "[", "rCrit", "]"}], "[",
RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"Calculate", " ", "the", " ", "concentration", " ", "using", " ", "the",
" ", "previously", " ", "attained", " ", "wave", " ", "front", " ",
"locations"}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"interpFunc", " ", "=", " ",
RowBox[{"Interpolation", "[",
RowBox[{"rCrit", ",",
RowBox[{"InterpolationOrder", "\[Rule]", "1"}]}], "]"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"conc", " ", "=", " ",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{"GF", "[",
RowBox[{"rPts", ",", "R", ",",
RowBox[{"tEval", "-", "tt"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"tt", ",", "0", ",", "tEval"}], "}"}], ",",
RowBox[{"{",
RowBox[{"R", ",", "0", ",",
RowBox[{"interpFunc", "[", "tt", "]"}]}], "}"}], ",",
RowBox[{"Method", "\[Rule]", "nIntOp"}]}], "]"}]}], ";"}]}]}]], "Input",\
CellChangeTimes->{{3.812893972651964*^9, 3.812894041380836*^9}, {
3.8128940808386497`*^9, 3.812894200853098*^9}, {3.812895090006173*^9,
3.812895130056717*^9}, {3.812895182901517*^9, 3.812895188768941*^9}, {
3.8128952940600986`*^9, 3.812895294240314*^9}, {3.812895481670691*^9,
3.812895603938406*^9}, {3.812907736516466*^9, 3.8129078158870068`*^9},
3.81306058574338*^9, {3.81306125415239*^9, 3.813061255044582*^9},
3.813068054341132*^9, {3.813068794578361*^9, 3.8130688045261297`*^9}, {
3.8140160913552523`*^9, 3.81401609172374*^9}},
CellLabel->"In[49]:=",ExpressionUUID->"d749fb9b-107e-4226-bd85-8077c148babb"],
Cell["\<\
Next, we calculate the concentration profile according to our asymptotic \
theory. This varies depending on system dimensionality -- i.e., whether \
cellDim = diffDim -- which we take into account below.\
\>", "Text",
CellChangeTimes->{{3.81288174649824*^9, 3.812881764990966*^9}, {
3.812881800761746*^9, 3.81288187245322*^9}, {3.812907671618938*^9,
3.812907733150786*^9}, {3.812907875227005*^9,
3.812908009796451*^9}},ExpressionUUID->"0eb8c951-aa7e-4e1d-b53f-\
4f43b518a579"],
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{"First", ",", " ",
RowBox[{
"turn", " ", "rPts", " ", "into", " ", "a", " ", "distance", " ", "from",
" ", "the", " ", "wave", " ", "front"}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"rTilde", " ", "=", " ",
RowBox[{"rPts", "-",
RowBox[{
RowBox[{"Last", "[", "rCrit", "]"}], "[",
RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"Then", ",", " ",
RowBox[{
"calculate", " ", "the", " ", "concentration", " ", "according", " ",
"to", " ", "the", " ", "system", " ",
RowBox[{"dimensionality", "."}]}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"concTheory", " ", "=", " ",
RowBox[{"If", "[",
RowBox[{
RowBox[{"cellDim", "\[Equal]", "diffDim"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Cth", " ",
RowBox[{"HeavisideTheta", "[", "rTilde", "]"}],
RowBox[{"Exp", "[",
RowBox[{"-",
FractionBox["rTilde",
RowBox[{"Dc", "/", "vBase"}]]}], "]"}]}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"Cth", "-",
RowBox[{"rTilde", "/", "vBase"}]}], ")"}],
RowBox[{"HeavisideTheta", "[",
RowBox[{"-", "rTilde"}], "]"}]}]}], ",", "\[IndentingNewLine]",
RowBox[{
FractionBox["1",
RowBox[{"Pi", " ", "vBase"}]],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"HeavisideTheta", "[", "rTilde", "]"}],
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox[
RowBox[{"1", "-",
FractionBox["1",
RowBox[{"Sqrt", "[",
RowBox[{"1", "+",
SuperscriptBox["k", "2"]}], "]"}]]}],
SuperscriptBox["k", "2"]],
RowBox[{"Exp", "[",
RowBox[{
RowBox[{"-",
FractionBox[
RowBox[{"vBase", " ", "rTilde"}],
RowBox[{"2", " ", "Dc"}]]}],
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
SuperscriptBox["k", "2"], "+", "1"}], "]"}], "+", "1"}],
")"}]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"k", ",",
RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}]}], "]"}]}],
"+",
RowBox[{
RowBox[{"HeavisideTheta", "[",
RowBox[{"-", "rTilde"}], "]"}],
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox["2",
SuperscriptBox["k", "2"]], "-",
RowBox[{
FractionBox[
RowBox[{"1", "+",
FractionBox["1",
RowBox[{"Sqrt", "[",
RowBox[{"1", "+",
SuperscriptBox["k", "2"]}], "]"}]]}],
SuperscriptBox["k", "2"]],
RowBox[{"Exp", "[",
RowBox[{
FractionBox[
RowBox[{"vBase", " ", "rTilde"}],
RowBox[{"2", " ", "Dc"}]],
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
SuperscriptBox["k", "2"], "+", "1"}], "]"}], "-", "1"}],
")"}]}], "]"}]}]}], ",",
RowBox[{"{",
RowBox[{"k", ",",
RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}]}],
"]"}]}]}], ")"}]}]}], "]"}]}], ";"}]}]}]], "Input",
CellChangeTimes->{{3.812908378909383*^9, 3.812908454754649*^9}, {
3.812908546589994*^9, 3.812908633599968*^9}, 3.812908732165764*^9,
3.812908774048277*^9, {3.812908976886189*^9, 3.8129089922980328`*^9}, {
3.813060610512183*^9, 3.813060643680135*^9}, {3.813060979345771*^9,
3.8130609825406322`*^9}, {3.8130612617776947`*^9, 3.813061266904784*^9}, {
3.813062410077359*^9, 3.813062421704865*^9}, {3.8139308228811913`*^9,
3.8139308232715597`*^9}},
CellLabel->"In[55]:=",ExpressionUUID->"7a29e695-b2e7-4a41-99bc-072959bd17ab"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
RowBox[{"Plot", "!!"}], "!"}], "*)"}], "\[IndentingNewLine]",
RowBox[{"Show", "[",
RowBox[{
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"rPts", ",", "conc"}], "}"}], "]"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<radius\>\"", ",", "\"\<concentration\>\""}], "}"}]}]}],
"]"}], ",", "\[IndentingNewLine]",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"rPts", ",", "concTheory"}], "}"}], "]"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Black"}]}], "]"}], ",",
"\[IndentingNewLine]",
RowBox[{"AspectRatio", "\[Rule]",
RowBox[{"1", "/", "2"}]}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.8128951511017427`*^9, 3.812895162117876*^9}, {
3.812907819014421*^9, 3.8129078230674543`*^9}, {3.812908462013791*^9,
3.812908480948791*^9}, {3.81290879561551*^9, 3.81290886939992*^9}},
CellLabel->"In[57]:=",ExpressionUUID->"815e3bcf-73d1-4f3c-afa8-dd9e15ccf94c"],
Cell[BoxData[
GraphicsBox[{{{},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.011000000000000001`],
AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJwBQQO+/CFib1JlAgAAADMAAAACAAAAYjJVMCqpIz806AwsXTphQLz3y4Lk
+CQ/yhccZmToYEAUvULVnkgmP2NSacXQkWBAboK5J1mYJz94Rexd2DZgQMdH
MHoT6Cg/VaFqQBWvX0AhDafMzTcqP3Q2i2FF6F5AetIdH4iHKz8OfO4RPxle
QNOXlHFC1yw/l+/nnmxCXUAsXQvE/CYuP8aq/xLxY1xAhSKCFrd2Lz+I2UhU
iH5bQPBzfLQ4YzA/cpYtjdmRWkCc1rfdFQsxP7CmLD+XnllASDnzBvOyMT9Y
cxoZxaRYQPWbLjDQWjI/q6EEOmelV0Ci/mlZrQIzPxeeJ7XWn1ZATmGlgoqq
Mz9bkwSqrZRVQPvD4KtnUjQ/XvreMjeEVECoJhzVRPo0P6zsWvbDblNAVIlX
/iGiNT+0gYrhlVRSQADskif/STY/dX2OTek1UUCuTs5Q3PE2P0keSO0EE1BA
WrEJermZNz83nJNpDdhNQAcURaOWQTg/6JwdgoWCS0C0doDMc+k4P7ZyftYQ
JklAYNm79VCROT+4UGvnFcNGQAw89x4uOTo/FgavPWJaREC5njJIC+E6P+ZU
WuYA70FAZgFuceiIOz+CZRAh5xM/QBJkqZrFMDw/UyciWHF5OkC/xuTDotg8
PyaSx0IuPDZAaykg7X+APT9qyz/M3oAyQBiMWxZdKD4/DSTqHPmoLkDF7pY/
OtA+P4cT5/OjXilAclHSaBd4Pz9WlhKZsP4kQA/aBkn6D0A/BWbRq4liIUBm
i6Td6GNAP83cW7XQzxxAvDxCcte3QD9h6sLVXOQXQBLu3wbGC0E/InM21qHT
E0Bon32btF9BP/jlvAT8dhBAv1AbMKOzQT/wF31Lbl0LQBUCucSRB0I/Ni6s
k5HBBkBss1ZZgFtCPx2PHSGa7wJAwmT07W6vQj/sT74Ysoj/PxcWkoJdA0M/
pAFmFspF+j9uxy8XTFdDP8zzMXDE5vU/xXjNqzqrQz9SqRkkZUTyPxsqa0Ap
/0M/nZWdEJx87j9x2wjVF1NEPxLLqLiTc+k/yIymaQanRD+QRs1XlUHlPx0+
RP70+kQ//DWzRCvC4T907+GS405FP8glGm9ert0/NPFmFg==
"]]}, {{}, {}}}, {{}, {{}, {},
{GrayLevel[0], PointSize[0.011000000000000001`], AbsoluteThickness[1.6],
LineBox[CompressedData["
1:eJwt0X0s1HEcB/CLRo2SONeTkPJ4nfMQ55LP5x6EzezI0zHFuQ6TM1eTWS27
TRFLXeWx1UbRISq220WzdR56cNVISCZj0jyUXeah666f33pv333/+Hy/r893
n6+TKCf6rAmFQvEl1uYu8olleiidwUgERsXYtfpOPrN2GKqJuLPEaPuSO3w/
zAV8idhUpaFUrnYVVrnCABH9qgj7Q5kFNnNuICEyGy9Cpx2NA9oAD9L7oErF
gkF7h+JrnlBFeGpaKg5V3snljtDBh/DqL6YgPdmiR+/CIOt7285g0SG5nSrP
i/QXpKfx1/lLHaw0JvwlvG5GMt77rJqgWniT928vJWFYoM5M99yb7C9pTcSV
WgbzY6IP6QVKhVhnyBQ+MfGFd4RnyUjAyNRH8utNvuT7Jhfj8I/mW3N6tB/p
tT+JxccuB4b5G34gJrir2TEYUxJvcKo7RnrCo6eQsqBwNYb7g5bw6ItR2Br5
XvB12R8Mm2CLABOfbS9QVweQ9xvGItHcNqS+gsOCN8T5CGEEduQVDpz/wQIG
0U83Go4pY50rgluBULE58IRQtAxaPcgIZP//Dz6q7/uEWUyxSW82nouSLdLc
ueLjoCDnBbhbrKzpZQaRdUkrG7v7ZjR1o0Fwo036QjN9DLPcHRevFJ4ATuYJ
ff68N+4pS7JLdguGueYevT7uKPYuVQD7YzB45m3Rabe6oSxqMIOWDxBWcpfe
pHRGqy9m4RtWCMX23b80exxRUvpwYi4NIZ1SFD2euA+7griyERWCfkTG59ym
4u6lSbM+Cw5QFQ2ZN5XWmPHgcm3HGQ5M2R4ZtJm3xG7BfmZ9OwekuoKxvBZz
pJqoe26Zc8Egu1qUkGqK59rjhIVJXBCnlGsdOo3wSvx7UdrGhUmqX+3C+jrQ
7BTyZFMeLO+0mp/lrYC034sWEc+DtzlDUJOxDL352mZ2Mw+MpWulqsIF2O+R
he5GHmgEpUpx7HfIHd82TDvFh/xy6/Fd+6bhdVlDplkjH7JrZaysk5NgH8w3
/N7gg2CwMl1UMwYXfk4ppiNDoK6Env105RP8A+OMfYo=
"]]}}, {{}, {}}}},
AspectRatio->NCache[
Rational[1, 2], 0.5],
Axes->{True, True},
AxesLabel->{
FormBox["\"radius\"", TraditionalForm],
FormBox["\"concentration\"", TraditionalForm]},
AxesOrigin->{0.00013957763671897206`, 0},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
ImageSize->{534.12890625, Automatic},
Method->{
"OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
Identity[
Part[#, 1]],
Identity[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
Identity[
Part[#, 1]],
Identity[
Part[#, 2]]}& )}},
PlotRange->{{0.00013957763671897206`, 0.0006502734375}, {
0, 137.8238735439387}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.812895255839679*^9, 3.8128953894540033`*^9, 3.8128954297450867`*^9,
3.8128956757263002`*^9, {3.812908757638755*^9, 3.81290887162022*^9},
3.812909369834551*^9, 3.81292812280097*^9, 3.812929252532555*^9,
3.812930848861638*^9, 3.813060373725398*^9, {3.813060627302177*^9,
3.813060647371336*^9}, {3.813060973429323*^9, 3.8130609861130466`*^9},
3.81306243606868*^9, 3.813068331868132*^9, 3.813070063685548*^9,
3.8139307792921267`*^9, 3.81393082580626*^9, 3.814016124201943*^9,
3.8140169630580893`*^9, 3.814017522143837*^9, 3.8140964152029657`*^9,
3.814186661111508*^9},
CellLabel->"Out[57]=",ExpressionUUID->"ee653c5b-1016-449e-87aa-e016469b49f9"]
}, Open ]],
Cell["\<\
We compute radial gradients of the asymptotic concentration profiles a la \
figure 4 in a similar manner to the above analysis. We may compare this to \
the gradients generated by a static disk (of radius rInv) of cells -- the \
\[OpenCurlyDoubleQuote]simple diffusion\[CloseCurlyDoubleQuote] model -- by \
numerical integration of said source.\
\>", "Text",
CellChangeTimes->{{3.81288174649824*^9, 3.812881764990966*^9}, {
3.812881800761746*^9, 3.81288187245322*^9}, {3.812907671618938*^9,
3.812907733150786*^9}, {3.812907875227005*^9, 3.812908009796451*^9}, {
3.812909610321862*^9, 3.812909622548875*^9}, {3.8129099375762672`*^9,
3.812910042782114*^9}, {3.812927090008402*^9, 3.812927114513781*^9}, {
3.8129277788494883`*^9, 3.812927795736866*^9}, {3.8129278791089354`*^9,
3.8129279040361023`*^9}, {3.812928069163431*^9, 3.8129280950512466`*^9}, {
3.8129287507728786`*^9, 3.812928769381155*^9},
3.814095912912431*^9},ExpressionUUID->"b482c275-92eb-4df9-8b75-\
1a1b1c4b8f56"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"gradRelayTheory", " ", "=", " ",
RowBox[{"If", "[",
RowBox[{
RowBox[{"cellDim", "\[Equal]", "diffDim"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"vBase", " ", "Cth", " ",
RowBox[{"HeavisideTheta", "[", "rTilde", "]"}],
RowBox[{
RowBox[{"Exp", "[",
RowBox[{"-",
FractionBox["rTilde",
RowBox[{"Dc", "/", "vBase"}]]}], "]"}], "/", "Dc"}]}], "-",
RowBox[{
RowBox[{"HeavisideTheta", "[",
RowBox[{"-", "rTilde"}], "]"}], "/", "vBase"}]}], ",",
"\[IndentingNewLine]",
RowBox[{
FractionBox["1",
RowBox[{"2", "Pi", " ", "Dc"}]],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"HeavisideTheta", "[", "rTilde", "]"}],
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox["1",
RowBox[{"Sqrt", "[",
RowBox[{
SuperscriptBox["k", "2"], "+", "1"}], "]"}]],
RowBox[{"Exp", "[",
RowBox[{
RowBox[{"-",
FractionBox[
RowBox[{"vBase", " ", "rTilde"}],
RowBox[{"2", " ", "Dc"}]]}],
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
SuperscriptBox["k", "2"], "+", "1"}], "]"}], "+", "1"}],
")"}]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"k", ",",
RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}]}], "]"}]}],
"+",
RowBox[{
RowBox[{"HeavisideTheta", "[",
RowBox[{"-", "rTilde"}], "]"}],
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{
FractionBox["1",
RowBox[{"Sqrt", "[",
RowBox[{"1", "+",
SuperscriptBox["k", "2"]}], "]"}]],
RowBox[{"Exp", "[",
RowBox[{
FractionBox[
RowBox[{"vBase", " ", "rTilde"}],
RowBox[{"2", " ", "Dc"}]],
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{
SuperscriptBox["k", "2"], "+", "1"}], "]"}], "-", "1"}],
")"}]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"k", ",",
RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}]}],
"]"}]}]}], ")"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"gradSimpleDiff", " ", "=", " ",
RowBox[{"-",
RowBox[{"NIntegrate", "[",
RowBox[{
RowBox[{"gradGF", "[",
RowBox[{"rPts", ",", "R", ",",
RowBox[{"tEval", "-", "tt"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"R", ",", "0", ",", "rInv"}], "}"}], ",",
RowBox[{"{",
RowBox[{"tt", ",", "0", ",", "tEval"}], "}"}], ",",
RowBox[{"Method", "\[Rule]", "nIntOp"}]}], "]"}]}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{"ListLogPlot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"rPts", ",", "gradRelayTheory"}], "}"}], "]"}], ",",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"rPts", ",", "gradSimpleDiff"}], "}"}], "]"}]}], "}"}], ",",
RowBox[{"Joined", "\[Rule]", "True"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<radius\>\"", ",", "\"\<gradient\>\""}], "}"}]}]}],
"]"}], "\[IndentingNewLine]"}], "Input",
CellChangeTimes->{{3.8129100533028183`*^9, 3.8129100671862288`*^9}, {
3.8129113579904747`*^9, 3.812911676530369*^9}, {3.812911713404502*^9,
3.812911750062838*^9}, {3.812927136312624*^9, 3.812927221425005*^9}, {
3.812927263689142*^9, 3.812927290082232*^9}, {3.8129273326531*^9,
3.8129274843556232`*^9}, {3.812927788082411*^9, 3.812927788284029*^9}, {
3.812928110683352*^9, 3.812928145811352*^9}, {3.812928542846753*^9,
3.8129286795064373`*^9}, {3.8129303260026903`*^9, 3.81293036001735*^9},
3.812930457375012*^9, 3.812930619523798*^9, {3.812930768166853*^9,
3.812930783569414*^9}, {3.813060142322863*^9, 3.813060177825492*^9},
3.8130624787081423`*^9, {3.813068080789174*^9, 3.8130680830104856`*^9},
3.8130690702416773`*^9, {3.813069810672179*^9, 3.813069815060051*^9}, {
3.813070104932334*^9, 3.813070105617446*^9}, {3.814016137360138*^9,
3.814016137609561*^9}},
CellLabel->"In[58]:=",ExpressionUUID->"213f45c1-e4a7-4543-8f74-09791c4762e6"],
Cell[BoxData[
GraphicsBox[{{}, {{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.011000000000000001`],
AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwBcQGO/iFib1JlAgAAABYAAAACAAAAv8bkw6LYPD/mLMNEoPopQGspIO1/
gD0/7a1YiCuUKUAYjFsWXSg+P/Mu7su2LSlAxe6WPzrQPj/5r4MPQscoQHJR
0mgXeD8//jAZU81gKEAP2gZJ+g9APwWyrpZY+idAZouk3ehjQD8KM0Ta45Mn
QLw8QnLXt0A/EbTZHW8tJ0AS7t8GxgtBPxc1b2H6xiZAaJ99m7RfQT8etgSl
hWAmQL9QGzCjs0E/Izea6BD6JUAVArnEkQdCPyq4LyyckyVAbLNWWYBbQj8v
OcVvJy0lQMJk9O1ur0I/Nbpas7LGJEAXFpKCXQNDPz078PY9YCRAbscvF0xX
Qz9DvIU6yfkjQMV4zas6q0M/SD0bflSTI0AbKmtAKf9DP06+sMHfLCNAcdsI
1RdTRD9VP0YFa8YiQMiMpmkGp0Q/WsDbSPZfIkAdPkT+9PpEP2JBcYyB+SFA
dO/hkuNORT9nwgbQDJMhQGR0lmU=
"]]},
{RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.011000000000000001`],
AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwBQQO+/CFib1JlAgAAADMAAAACAAAAYjJVMCqpIz9F4vJpz5QqQLz3y4Lk
+CQ/FxquFP12KkAUvULVnkgmP4gTukuJSSpAboK5J1mYJz/bucamLxwqQMdH
MHoT6Cg/R06Gu37vKUAhDafMzTcqP5ZONKlNwylAetIdH4iHKz+VIAz5eJcp
QNOXlHFC1yw/vVOM7+FrKUAsXQvE/CYuP05bt/BtQClAhSKCFrd2Lz8tId3N
BRUpQPBzfLQ4YzA/MBIkGJXpKECc1rfdFQsxP80ohkUKvihASDnzBvOyMT/A
1sQEVJIoQPWbLjDQWjI/M1iraWVmKECi/mlZrQIzP8k6MDgxOihATmGlgoqq
Mz/hD4H8qw0oQPvD4KtnUjQ/AqDvUcvgJ0CoJhzVRPo0P9tAO8+FsydAVIlX
/iGiNT8aohDr0oUnQADskif/STY/ASI51KpXJ0CuTs5Q3PE2P2cVx24GKSdA
WrEJermZNz/KS44q3/kmQAcURaOWQTg/UkoCAi/KJkC0doDMc+k4PwzNumXw
mSZAYNm79VCROT8gD583HmkmQAw89x4uOTo/eZ5UqbM3JkC5njJIC+E6P+yw
4k+sBSZAZgFuceiIOz+nPhsUBNMlQBJkqZrFMDw/oq1DG7efJUC/xuTDotg8
PycO4MvBayVAaykg7X+APT9EcyTPIDclQBiMWxZdKD4/muVS+tABJUDF7pY/
OtA+P4GFrlbPyyRAclHSaBd4Pz/JSZUaGZUkQA/aBkn6D0A/PPGcoatdJEBm
i6Td6GNAP4DfSXuEJSRAvDxCcte3QD9mH1lRoewjQBLu3wbGC0E/Ji9/zv+y
I0Bon32btF9BP57QS/edeCNAv1AbMKOzQT+DIbi/eT0jQBUCucSRB0I/w/hL
WpEBI0Bss1ZZgFtCP7fELr/ixCJAwmT07W6vQj8DsxNqbIciQBcWkoJdA0M/
nv4KpixJIkBuxy8XTFdDPwWXEuIhCiJAxXjNqzqrQz/xaj6gSsohQBsqa0Ap
/0M/TBDJb6WJIUBx2wjVF1NEP/UFH/MwSCFAyIymaQanRD/NsfvY6wUhQB0+
RP70+kQ/gosL39TCIEB07+GS405FP4wqFM7qfiBAgkFYCg==
"]]}}, {{}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{
FormBox["\"radius\"", TraditionalForm],
FormBox["\"gradient\"", TraditionalForm]},
AxesOrigin->{0.00013957763671897206`, 7.853260599956207},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{
Charting`ScaledTicks[{Log, Exp}],
Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
Method->{
"OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
Identity[
Part[#, 1]],
Exp[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
Identity[
Part[#, 1]],
Exp[
Part[#, 2]]}& )}},
PlotRange->{{0.00013957763671897206`, 0.0006502734375}, {7.853260599956207,
13.290644942187223`}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->FrontEndValueCache[{Automatic,
Charting`ScaledTicks[{Log, Exp}]}, {Automatic, {{8.517193191416238,
FormBox["5000", TraditionalForm], {0.01, 0.}}, {9.210340371976184,
FormBox[
TemplateBox[{"\[Times]", "\"\[Times]\"", "1",
TemplateBox[{"10", "4"}, "Superscript", SyntaxForm ->
SuperscriptBox]}, "RowWithSeparators"], TraditionalForm], {0.01,
0.}}, {10.819778284410283`,
FormBox[
TemplateBox[{"\[Times]", "\"\[Times]\"", "5",
TemplateBox[{"10", "4"}, "Superscript", SyntaxForm ->
SuperscriptBox]}, "RowWithSeparators"], TraditionalForm], {0.01,
0.}}, {11.512925464970229`,
FormBox[
TemplateBox[{"\[Times]", "\"\[Times]\"", "1",
TemplateBox[{"10", "5"}, "Superscript", SyntaxForm ->
SuperscriptBox]}, "RowWithSeparators"], TraditionalForm], {0.01,
0.}}, {13.122363377404328`,
FormBox[
TemplateBox[{"\[Times]", "\"\[Times]\"", "5",
TemplateBox[{"10", "5"}, "Superscript", SyntaxForm ->
SuperscriptBox]}, "RowWithSeparators"], TraditionalForm], {0.01,
0.}}, {6.907755278982137,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
7.600902459542082,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.006367567650246,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.294049640102028,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.699514748210191,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.85366542803745,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.987196820661973,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
9.104979856318357,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
9.903487552536127,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
10.308952660644293`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
10.596634733096073`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.002099841204238`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.156250521031495`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.289781913656018`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.407564949312402`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
12.206072645530174`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
12.611537753638338`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
12.89921982609012,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
13.304684934198283`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
13.458835614025542`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
13.592367006650065`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
13.710150042306449`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
13.815510557964274`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
14.22097566607244,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
14.508657738524219`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}}}]]], "Output",
CellChangeTimes->{
3.812929541904025*^9, 3.812930395752589*^9, 3.8129305113118896`*^9,
3.812930659198135*^9, 3.8129309058624477`*^9, 3.8130603801775923`*^9,
3.813062486947851*^9, {3.8130690672852373`*^9, 3.813069077731667*^9},
3.813070073917334*^9, 3.81393083511801*^9, {3.814016127962164*^9,
3.814016141093104*^9}, 3.814016964011634*^9, 3.814017530489252*^9,
3.814096423383594*^9, 3.8141866700328493`*^9},
CellLabel->"Out[60]=",ExpressionUUID->"af3a5f5f-22b2-46a4-9555-21069b82e7dc"]
}, Open ]],
Cell["\<\
Finally, we compute the initiation time for a range of rInv defined by \
[rInvMin, rInvMax] with nrInv total values. We plot the asymptotic forms of \
the initiation times alongside the numerically calculated forms. In this \
section of code, we assume that cellDim = diffDim for purposes of calculating \
the theoretical initiation time.\
\>", "Text",
CellChangeTimes->{{3.81288174649824*^9, 3.812881764990966*^9}, {
3.812881800761746*^9, 3.81288187245322*^9}, {3.812907671618938*^9,
3.812907733150786*^9}, {3.812907875227005*^9, 3.812908009796451*^9}, {
3.812909610321862*^9, 3.812909622548875*^9}, {3.8129099375762672`*^9,
3.812910042782114*^9}, {3.812927090008402*^9, 3.812927114513781*^9}, {
3.8129277788494883`*^9, 3.812927795736866*^9}, {3.8129278791089354`*^9,
3.8129279040361023`*^9}, {3.812928069163431*^9, 3.8129280950512466`*^9}, {
3.8129287507728786`*^9, 3.812928769381155*^9}, {3.8140959210670652`*^9,
3.814095945843891*^9}, {3.814095994644442*^9, 3.814096057041109*^9}, {
3.814096260192081*^9,
3.814096303253516*^9}},ExpressionUUID->"a29a21ba-68df-4007-bc3b-\
2bef71219154"],
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
RowBox[{
"The", " ", "minimum", " ", "initiating", " ", "colony", " ", "sizes",
" ", "for", " ", "1"}], ",", " ", "2", ",", " ",
RowBox[{"and", " ", "3", " ", "dimensions"}]}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"rInvMin", " ", "=", " ",
RowBox[{"{",
RowBox[{
RowBox[{"0.1",
RowBox[{"Dc", "/", "vBase"}]}], ",",
RowBox[{"0.6", " ",
RowBox[{"Dc", "/", "vBase"}]}], ",",
RowBox[{
RowBox[{"Sqrt", "[", "3", "]"}],
RowBox[{"Dc", "/", "vBase"}]}]}], "}"}]}], ";"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"rInvMax", " ", "=", " ",
RowBox[{"10", " ",
RowBox[{"Dc", "/", "vBase"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"nrInv", " ", "=", " ", "25"}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"maximum", " ", "time", " ", "considered", " ", "in", " ", "the", " ",
"numerics"}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"maxInitTime", " ", "=", " ",
RowBox[{"20000", " ",
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"numerical", " ", "time", " ", "resolution"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"initTimeRes", " ", "=", " ",
RowBox[{"0.01", " ",
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}]}]}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
"construct", " ", "an", " ", "array", " ", "with", " ", "all", " ", "the",
" ", "desired", " ", "rInv", " ", "values"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"rInvVals", " ", "=", " ",
RowBox[{
RowBox[{"rInvMin", "[",
RowBox[{"[", "diffDim", "]"}], "]"}], " ",
RowBox[{"Table", "[",
RowBox[{
SuperscriptBox[
RowBox[{"(",
FractionBox["rInvMax",
RowBox[{"rInvMin", "[",
RowBox[{"[", "cellDim", "]"}], "]"}]], ")"}],
FractionBox[
RowBox[{"i", "-", "1"}],
RowBox[{"nrInv", "-", "1"}]]], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "nrInv"}], "}"}]}], "]"}]}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"compute", " ", "the", " ", "initiation", " ", "times"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"inits", " ", "=", " ",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"initTime", "[",
RowBox[{
RowBox[{"rInvVals", "[",
RowBox[{"[", "i", "]"}], "]"}], ",", "maxInitTime", ",",
"initTimeRes"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",",
RowBox[{"Length", "[", "rInvVals", "]"}]}], "}"}]}], "]"}]}],
";"}]}]}]], "Input",
CellChangeTimes->{{3.8140166557952547`*^9, 3.81401667127981*^9}, {
3.8140167167721157`*^9, 3.8140167177636633`*^9}, {3.8140959487537527`*^9,
3.8140959794317923`*^9}, {3.814096047488584*^9, 3.814096109795792*^9}, {
3.814096213428752*^9, 3.814096213771253*^9}, {3.8140963166718597`*^9,
3.814096348075448*^9}, {3.814096413172039*^9, 3.814096413251877*^9}, {
3.81409660536528*^9, 3.814096606583685*^9}, {3.814096686770877*^9,
3.814096696282629*^9}, 3.814097115654194*^9, {3.814111429645339*^9,
3.814111429719967*^9}, 3.814112607442059*^9, {3.8141862137384863`*^9,
3.814186296614746*^9}},
CellLabel->"In[61]:=",ExpressionUUID->"05fe6413-4b85-4e94-b706-56bab029f104"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
"calculate", " ", "the", " ", "theoretical", " ", "initiation", " ",
"times", " ",
RowBox[{"(",
RowBox[{"normalized", " ", "by", " ",
RowBox[{"Dc", "/",
SuperscriptBox["vBase", "2"]}]}], ")"}], " ", "at", " ", "small", " ",
"rInv"}], "*)"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"theoryInitSmall", " ", "=", " ",
RowBox[{
FractionBox[
SuperscriptBox["vBase", "2"], "Dc"],
RowBox[{"{",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"Cth", " ",
RowBox[{
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{"Pi", " ", "Dc"}], "]"}], "/", "2"}], "/",
"rInvVals"}]}], ")"}], "2"], ",",
RowBox[{
RowBox[{"(",
FractionBox[
SuperscriptBox["rInvVals", "2"],
RowBox[{"4", " ", "Dc"}]], ")"}],
RowBox[{"Exp", "[",
SuperscriptBox[
RowBox[{"(",
FractionBox[
RowBox[{"2", "Dc"}],
RowBox[{"vBase", " ", "rInvVals"}]], ")"}], "2"], "]"}]}], ",",
RowBox[{
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"rInvVals",
SuperscriptBox[
RowBox[{"(",
RowBox[{"1", "-",
FractionBox[
RowBox[{"3", " ", "Cth", " ", "Dc"}],
SuperscriptBox["rInvVals", "2"]]}], ")"}],
RowBox[{"-", "1"}]]}], ")"}], "2"], "/", "Pi"}], " ", "Dc"}]}],
"}"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"plot", " ",
RowBox[{"away", "!"}]}], "*)"}], "\[IndentingNewLine]",
RowBox[{"ListLogLogPlot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{
RowBox[{"vBase", " ",
RowBox[{"rInvVals", "/", "Dc"}]}], ",",
RowBox[{
SuperscriptBox["vBase", "2"],
RowBox[{"inits", "/", "Dc"}]}]}], "}"}], "]"}], ",",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{
RowBox[{"vBase", " ",
RowBox[{"rInvVals", "/", "Dc"}]}], ",",
RowBox[{"theoryInitSmall", "[",
RowBox[{"[", "diffDim", "]"}], "]"}]}], "}"}], "]"}], ",",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{
RowBox[{"vBase", " ",
RowBox[{"rInvVals", "/", "Dc"}]}], ",",
RowBox[{
RowBox[{"0.", "rInvVals"}], "+", "2."}]}], "}"}], "]"}]}], "}"}],
",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0.1", ",", "10"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",",
RowBox[{
SuperscriptBox["vBase", "2"],
RowBox[{
RowBox[{"maxInitTime", "/", "Dc"}], "/", "2"}]}]}], "}"}]}],
"}"}]}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<\!\(\*SubscriptBox[\(vr\), \(i\)]\)/D\>\"", ",",
"\"\<\!\(\*SuperscriptBox[\(v\), \(2\)]\)\!\(\*SubscriptBox[\(t\), \
\(init\)]\)/D\>\""}], "}"}]}], ",",
RowBox[{"Joined", "\[Rule]", "True"}]}], "]"}]}]}]], "Input",
CellChangeTimes->{{3.814016948920835*^9, 3.8140169698480997`*^9}, {
3.814017219856649*^9, 3.814017233443769*^9}, {3.814095981125339*^9,
3.814095991745496*^9}, {3.814096362654705*^9, 3.814096490258717*^9}, {
3.814096543424787*^9, 3.814096608970636*^9}, {3.814096644242092*^9,
3.814096645967469*^9}, {3.8140970520527353`*^9, 3.814097086969523*^9}, {
3.814098994599308*^9, 3.814099000939061*^9}, {3.814099056290957*^9,
3.81409907277736*^9}, {3.814186300314896*^9, 3.814186381386406*^9}},
CellLabel->"In[68]:=",ExpressionUUID->"b0e05f92-2b05-4e86-afb5-5a428edf6ffa"],
Cell[BoxData[
GraphicsBox[{{}, {{}, {},
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.011000000000000001`],
AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwBoQFe/iFib1JlAgAAABkAAAACAAAAdVu5fNs84L8WVbW7sWsiQEL1Osy+
MNm/MRDvrYfBHEDewXOjH7DRv5sUvTASWhZA9xxZ9QBfxL9lVZBdiXkRQLbY
Ko8Kd6W/H0o9BBCTC0A2YYdb90azP71QxSJoCwZAXhdS/7mkyD8dKpoz+ewB
QBY/cCj80tM/GulRNi3A/T93cjdRm1PbP8jCIe4KO/k/7FL/PB1q4T+mV7gy
W/D1P6DsYtFsKuU/5fltPWt88z9QhsZlvOroP3eXRiisvPE/AyAq+gur7D9z
CKeXqmrwP9rcRsetNfA/WqNkG4jx7j+zqXiR1RXyP2dm+tD0cu0/jHaqW/31
8z++dRftmmnsP2RD3CUl1vU/AeksZcuc6z8+EA7wTLb3P5z9e04e7uo/F90/
unSW+T/5nbYtpV/qP+6pcYScdvs/TJi/7CPz6T/IdqNOxFb9P0y2lM8thek/
oUPVGOw2/z/JwiHuCjvpPzyIg/GJiwBAmtzAsLgV6T+qbpzWnXsBQHY5Wrk6
8Og/FlW1u7FrAkB2OVq5OvDoP0GWzEw=
"]]},
{RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.011000000000000001`],
AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJwBoQFe/iFib1JlAgAAABkAAAACAAAAUhSB+q5Y4L/HeQlPBWghQEL1Osy+
MNm/ZjxmrTh2GkDewXOjH7DRvzSUtrBeDRRA9xxZ9QBfxL8sRPT9PlsOQLbY
Ko8Kd6W/bcgO3qUJB0A2YYdb90azPyWgNOr4owFAXhdS/7mkyD8A60Pm0Yb7
PxY/cCj80tM/0YwLukYt9j93cjdRm1PbP3DSrQ+puvI/7FL/PB1q4T/T/+vz
NMnwP6DsYtFsKuU/dF82LWsI8D9QhsZlvOroPxlOKlCfOPA/AyAq+gur7D++
kNqlcyfxP9rcRsetNfA/LsflLxGt8j+zqXiR1RXyP15dnnn0qfQ/jHaqW/31
8z9L9gEdMAX3P2RD3CUl1vU/oBYbYgyr+T8+EA7wTLb3P3F+h4Pwi/w/F90/
unSW+T9EHu0zhpv/P+6pcYScdvs/n8LRmwVoAUDIdqNOxFb9Pwh7nLjjEANA
oUPVGOw2/z+ghFzhT8UEQDyIg/GJiwBAzz2V3t+CBkCqbpzWnXsBQH3uDa2q
RwhAFlW1u7FrAkCFRQx9LRIKQCqLxJA=
"]]},
{RGBColor[0.560181, 0.691569, 0.194885], PointSize[0.011000000000000001`],
AbsoluteThickness[1.6], LineBox[CompressedData["
1:eJxTTMoPSmViYGCQBGIQHSTS+GtdxIP97y1//XPSe2bv9NXqzD6Dm3D+vYPF
i+U3XITzv8tEfmWIPwLnb7uh1c9VvhTON0tsj/7uttkexo8TD/q/c8kJOF/M
vkDjz6XLcH55kXng7ODbcP6boP82slkP4fwFb5Iu5mg9hfMD2o6l7nn1As5n
VtD6xb36DZx/647b8bWmH+D8zSsrJl4V/QTn95Stiv779TOcn+J8R1X12lc4
306A74PPtu9wvvhd+10l037C+e9WFrbMKfsN558oW+x3JOwvnL/Q+arEG7P/
cL5NR/PHzm4GBxh/Vd6ca3OrGeF8sdCtuzdmM8H5AL4r4CY=
"]]}}, {{}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{
FormBox[
"\"\\!\\(\\*SubscriptBox[\\(vr\\), \\(i\\)]\\)/D\"", TraditionalForm],
FormBox[
"\"\\!\\(\\*SuperscriptBox[\\(v\\), \\(2\\)]\\)\\!\\(\\*SubscriptBox[\\(t\
\\), \\(init\\)]\\)/D\"", TraditionalForm]},
AxesOrigin->{-2.3025850929940455`, 0.},
DisplayFunction->Identity,
Frame->{{False, False}, {False, False}},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{
Charting`ScaledTicks[{Log, Exp}],
Charting`ScaledFrameTicks[{Log, Exp}]}, {
Charting`ScaledTicks[{Log, Exp}],
Charting`ScaledFrameTicks[{Log, Exp}]}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
Method->{
"OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True,
"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
Exp[
Part[#, 1]],
Exp[
Part[#, 2]]}& ), "CopiedValueFunction" -> ({
Exp[
Part[#, 1]],
Exp[
Part[#, 2]]}& )}},
PlotRange->NCache[{{-2.3025850929940455`,
Log[10]}, {0, 9.210340371976184}}, {{-2.3025850929940455`,
2.302585092994046}, {0, 9.210340371976184}}],
PlotRangeClipping->True,
PlotRangePadding->{{0, 0}, {0, 0}},
Ticks->FrontEndValueCache[{
Charting`ScaledTicks[{Log, Exp}],
Charting`ScaledTicks[{Log, Exp}]}, {{{-2.3025850929940455`,
FormBox[
TagBox[
InterpretationBox[
StyleBox["\"0.1\"", ShowStringCharacters -> False], 0.1, AutoDelete ->
True], NumberForm[#, {
DirectedInfinity[1], 1}]& ], TraditionalForm], {0.01,
0.}}, {-0.6931471805599453,
FormBox[
TagBox[
InterpretationBox[
StyleBox["\"0.5\"", ShowStringCharacters -> False], 0.5, AutoDelete ->
True], NumberForm[#, {
DirectedInfinity[1], 1}]& ], TraditionalForm], {0.01, 0.}}, {0.,
FormBox["1", TraditionalForm], {0.01, 0.}}, {1.6094379124341003`,
FormBox["5", TraditionalForm], {0.01, 0.}}, {2.302585092994046,
FormBox["10", TraditionalForm], {0.01, 0.}}, {-2.995732273553991,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-2.8134107167600364`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-2.659260036932778,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-2.5257286443082556`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-2.4079456086518722`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-1.6094379124341003`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-1.2039728043259361`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.916290731874155,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.5108256237659907,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.35667494393873245`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.2231435513142097,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.10536051565782628`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
0.6931471805599453,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.0986122886681098`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.3862943611198906`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.791759469228055,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.9459101490553132`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
2.0794415416798357`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
2.1972245773362196`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
2.995732273553991,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
3.4011973816621555`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
3.6888794541139363`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
3.912023005428146,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.0943445622221,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.248495242049359,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.382026634673881,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.499809670330265,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.605170185988092,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}}, {{
0.,
FormBox["1", TraditionalForm], {0.01, 0.}}, {2.302585092994046,
FormBox["10", TraditionalForm], {0.01, 0.}}, {4.605170185988092,
FormBox["100", TraditionalForm], {0.01, 0.}}, {6.907755278982137,
FormBox["1000", TraditionalForm], {0.01, 0.}}, {9.210340371976184,
FormBox[
TemplateBox[{"10", "4"}, "Superscript", SyntaxForm -> SuperscriptBox],
TraditionalForm], {0.01, 0.}}, {-2.3025850929940455`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-1.6094379124341003`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-1.2039728043259361`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.916290731874155,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.6931471805599453,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.5108256237659907,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.35667494393873245`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.2231435513142097,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}, {-0.10536051565782628`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
0.6931471805599453,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.0986122886681098`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.3862943611198906`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.6094379124341003`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.791759469228055,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
1.9459101490553132`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
2.0794415416798357`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
2.1972245773362196`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
2.995732273553991,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
3.4011973816621555`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
3.6888794541139363`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
3.912023005428146,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.0943445622221,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.248495242049359,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.382026634673881,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
4.499809670330265,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
5.298317366548036,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
5.703782474656201,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
5.991464547107982,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
6.214608098422191,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
6.396929655216146,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
6.551080335043404,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
6.684611727667927,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
6.802394763324311,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
7.600902459542082,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.006367567650246,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.294049640102028,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.517193191416238,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.699514748210191,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.85366542803745,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
8.987196820661973,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
9.104979856318357,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
9.903487552536127,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
10.308952660644293`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
10.596634733096073`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
10.819778284410283`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.002099841204238`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.156250521031495`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.289781913656018`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.407564949312402`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005, 0.}}, {
11.512925464970229`,
FormBox[
TemplateBox[{0., 0.}, "Spacer2"], TraditionalForm], {0.005,
0.}}}}]]], "Output",
CellChangeTimes->{{3.814016690906822*^9, 3.81401671259899*^9},
3.8140168545793257`*^9, 3.814017057955755*^9, 3.814017197194605*^9,
3.814017233983601*^9, 3.8140358545579033`*^9, 3.8140969881607313`*^9, {
3.814097060337653*^9, 3.8140970872099*^9}, 3.8140982556152773`*^9,
3.814099002665916*^9, 3.814111414005249*^9, 3.814114072984881*^9,
3.8141881403459044`*^9},
CellLabel->"Out[69]=",ExpressionUUID->"3278bf8c-b3f0-4988-93b2-2e2d87ed0644"]
}, Open ]]
}, Open ]]
},
WindowSize->{1440, 855},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"12.1 for Mac OS X x86 (64-bit) (June 19, 2020)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"3e856ba2-27f8-4177-a92f-830098652d9e"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[422, 15, 552, 10, 71, "Subsubsection",ExpressionUUID->"9f649eb5-553d-479c-a38a-504bc2532627"],
Cell[977, 27, 814, 14, 81, "Text",ExpressionUUID->"8110de1f-e0f6-4e12-82cf-66fe630a4d28"],
Cell[1794, 43, 22557, 336, 243, "Input",ExpressionUUID->"aa0f6726-1167-4f75-a1b6-d3e71c46e215",
CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[24354, 381, 978, 15, 81, "Text",ExpressionUUID->"f9c5f87b-a4fe-4512-b47f-8a0fd2619f52"],
Cell[25335, 398, 9159, 281, 634, "Input",ExpressionUUID->"1d5411fc-c466-4794-844c-f005c3e23b2f"],
Cell[34497, 681, 761, 13, 58, "Text",ExpressionUUID->"28ba037a-8215-4606-9435-94131c549906"],
Cell[CellGroupData[{
Cell[35283, 698, 5033, 116, 393, "Input",ExpressionUUID->"d9a87681-af07-467e-9402-8a19d89d974c"],
Cell[40319, 816, 1350, 20, 24, "Print",ExpressionUUID->"bbddf319-3c30-43f2-9362-f3ee0699935d"]
}, Open ]],
Cell[41684, 839, 1153, 18, 81, "Text",ExpressionUUID->"e17f4ff8-50db-4563-944d-7608557cef97"],
Cell[42840, 859, 8091, 174, 620, "Input",ExpressionUUID->"1f0cee4e-a7c2-46d6-bede-6f85096a6f5c",
CellGroupingRules->{GroupTogetherGrouping, 10000.}],
Cell[50934, 1035, 174, 2, 34, "Output",ExpressionUUID->"ad77fcab-354a-4dc8-ae23-ea34b9556d85"],
Cell[51111, 1039, 475, 9, 58, "Text",ExpressionUUID->"b026b330-3d2e-4447-b0e2-fcae81df7333"],
Cell[CellGroupData[{
Cell[51611, 1052, 6562, 158, 534, "Input",ExpressionUUID->"84ac632a-74dc-4374-87e8-ba9f2924b487"],
Cell[CellGroupData[{
Cell[58198, 1214, 479, 8, 24, "Print",ExpressionUUID->"b655968e-f828-45b9-9f09-a3ad4e132d7a"],
Cell[58680, 1224, 481, 8, 24, "Print",ExpressionUUID->"64fb176c-beac-46d7-a4ff-09b6e44d8fc7"],
Cell[59164, 1234, 479, 8, 24, "Print",ExpressionUUID->"86dcb18b-968d-4fd2-a27a-143c9230089c"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[59692, 1248, 2940, 75, 240, "Input",ExpressionUUID->"36cd9677-5a82-481e-ba8a-cd3be4e7b70d"],
Cell[62635, 1325, 2646, 59, 235, "Output",ExpressionUUID->"8ed4f077-da0f-4883-b769-6e6bb29e548f"],
Cell[65284, 1386, 3452, 72, 231, "Output",ExpressionUUID->"dbc7163f-6eb9-4dd9-aa56-adc784cf03c1"]
}, Open ]],
Cell[68751, 1461, 469, 9, 35, "Text",ExpressionUUID->"bd73eded-ac7d-4683-b811-b6086b68897c"],
Cell[69223, 1472, 3069, 73, 304, "Input",ExpressionUUID->"d749fb9b-107e-4226-bd85-8077c148babb"],
Cell[72295, 1547, 495, 9, 35, "Text",ExpressionUUID->"0eb8c951-aa7e-4e1d-b53f-4f43b518a579"],
Cell[72793, 1558, 4349, 110, 309, "Input",ExpressionUUID->"7a29e695-b2e7-4a41-99bc-072959bd17ab"],
Cell[CellGroupData[{
Cell[77167, 1672, 1180, 29, 94, "Input",ExpressionUUID->"815e3bcf-73d1-4f3c-afa8-dd9e15ccf94c"],
Cell[78350, 1703, 4412, 91, 287, "Output",ExpressionUUID->"ee653c5b-1016-449e-87aa-e016469b49f9"]
}, Open ]],
Cell[82777, 1797, 1015, 16, 58, "Text",ExpressionUUID->"b482c275-92eb-4df9-8b75-1a1b1c4b8f56"],
Cell[CellGroupData[{
Cell[83817, 1817, 4797, 120, 352, "Input",ExpressionUUID->"213f45c1-e4a7-4543-8f74-09791c4762e6"],
Cell[88617, 1939, 7871, 179, 232, "Output",ExpressionUUID->"af3a5f5f-22b2-46a4-9555-21069b82e7dc"]
}, Open ]],
Cell[96503, 2121, 1125, 18, 58, "Text",ExpressionUUID->"a29a21ba-68df-4007-bc3b-2bef71219154"],
Cell[97631, 2141, 3666, 95, 355, "Input",ExpressionUUID->"05fe6413-4b85-4e94-b706-56bab029f104"],
Cell[CellGroupData[{
Cell[101322, 2240, 3971, 107, 191, "Input",ExpressionUUID->"b0e05f92-2b05-4e86-afb5-5a428edf6ffa"],
Cell[105296, 2349, 14409, 339, 246, "Output",ExpressionUUID->"3278bf8c-b3f0-4988-93b2-2e2d87ed0644"]
}, Open ]]
}, Open ]]
}
]
*)