(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    165308,       4095]
NotebookOptionsPosition[    149599,       3890]
NotebookOutlinePosition[    150110,       3909]
CellTagsIndexPosition[    150067,       3906]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[TextData[{
 StyleBox["Covid-19 Case Fitting Curves", "Title"],
 "\n",
 StyleBox["Esmerindo Bernardes", "Author"],
 "\n",
 StyleBox["sousa@ifsc.usp.br", "Affiliation"],
 "\n04/2020 - 05/2021 - 07/2022 - 05/2023"
}], "Text",
 CellChangeTimes->{{3.796547237632004*^9, 3.7965472587694883`*^9}, 
   3.7965473940725307`*^9, {3.7965474294289565`*^9, 3.7965474671091633`*^9}, {
   3.7965475246711864`*^9, 3.7965475427308784`*^9}, {3.8311346591403666`*^9, 
   3.831134664636945*^9}, {3.864293943255308*^9, 3.864293948984831*^9}, {
   3.8678300426341963`*^9, 3.867830043402196*^9}, {3.8853032121209507`*^9, 
   3.885303216370919*^9}, {3.8944835980035734`*^9, 3.8944835987001286`*^9}},
 TextAlignment->Center,
 TextJustification->1.,ExpressionUUID->"79c1882b-12d4-4336-94a2-cac77e8e48ae"],

Cell[CellGroupData[{

Cell[TextData[StyleBox["Model function",
 FontColor->RGBColor[0.5, 0, 0.5]]], "Subtitle",
 CellChangeTimes->{{3.867830127019478*^9, 3.8678301312004833`*^9}},
 TextAlignment->Center,ExpressionUUID->"80b7cff6-1a9e-472c-a1ed-a6b22e42e134"],

Cell[TextData[{
 StyleBox["Accumulated cases:",
  FontColor->RGBColor[0, 0, 1]],
 " a linear superposition of hyperbolic tangents \nwhere each hyperbolic \
tangent describes one wave.\n\n",
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{
     SubscriptBox["Z", "N"], "(", "n", ")"}], "=", 
    RowBox[{
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"i", "=", "1"}], "l"], 
      RowBox[{
       SubscriptBox["b", "i"], 
       RowBox[{"tanh", "(", 
        RowBox[{
         RowBox[{
          SubscriptBox["a", "i"], "n"}], "-", 
         SubscriptBox["c", "i"]}], ")"}]}]}], "+", "\[Delta]"}]}], 
   TraditionalForm]],
  FormatType->TraditionalForm,ExpressionUUID->
  "b359aa84-e8c2-4aca-88f3-842524fd89e4"]
}], "Text",
 CellChangeTimes->{{3.867830318549573*^9, 3.867830435697046*^9}, 
   3.8678307541504*^9, {3.8678308425284023`*^9, 3.867830944577354*^9}, {
   3.8678310071053524`*^9, 3.867831070512475*^9}, {3.8678311525075006`*^9, 
   3.8678311597504992`*^9}, {3.8678315790274253`*^9, 3.867831587519425*^9}},
 TextAlignment->Center,
 TextJustification->1,ExpressionUUID->"a6adaaa0-c19e-4e36-9804-31f280d0559d"],

Cell[TextData[{
 StyleBox["Cases per week:",
  FontColor->RGBColor[0, 0, 1]],
 " first derivative of a linear superposition of hyperbolic tangents.\n\n",
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{
     SubscriptBox["\[GothicCapitalZ]", "N"], "(", "n", ")"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"i", "=", "1"}], "l"], 
     FractionBox[
      RowBox[{
       SubscriptBox["a", "i"], 
       SubscriptBox["b", "i"]}], 
      RowBox[{
       SuperscriptBox["cosh", "2"], "(", 
       RowBox[{
        RowBox[{
         SubscriptBox["a", "i"], "n"}], "-", 
        SubscriptBox["c", "i"]}], ")"}]]}]}], TraditionalForm]],
  ExpressionUUID->"652220bb-c954-438e-8dfb-befd4b580dea"]
}], "Text",
 CellChangeTimes->{{3.867830452919488*^9, 3.867830519775468*^9}, 
   3.8678305541014986`*^9, {3.867831320933484*^9, 3.867831323702485*^9}, {
   3.867831356232523*^9, 3.8678314648219013`*^9}, 3.8678315908054276`*^9},
 TextAlignment->Center,ExpressionUUID->"be5c8736-934b-45a1-b64b-475e840fbab2"],

Cell[CellGroupData[{

Cell["Routines", "Section",
 CellChangeTimes->{{3.8709413678476343`*^9, 
  3.870941370279031*^9}},ExpressionUUID->"b45ff5bc-1a50-406c-a509-\
76e4636410a4"],

Cell[CellGroupData[{

Cell["Model functions", "Subsection",
 CellChangeTimes->{{3.870941479130986*^9, 
  3.870941481209756*^9}},ExpressionUUID->"fc505297-3fe8-4104-a3c5-\
30fc749ac8aa"],

Cell[BoxData[
 RowBox[{
  RowBox[{"cshn", "[", 
   RowBox[{"z_", ",", "a_", ",", "b_", ",", "c_"}], "]"}], ":=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", "i", "]"}], "]"}], "*", 
     RowBox[{
      RowBox[{"b", "[", 
       RowBox[{"[", "i", "]"}], "]"}], "/", 
      RowBox[{
       RowBox[{"Cosh", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"a", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "*", "z"}], "-", 
         RowBox[{"c", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "^", "2"}]}]}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"Length", "[", "a", "]"}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8678318895904016`*^9, 3.8678318996774015`*^9}, 
   3.8732921912414303`*^9},
 CellLabel->"In[1]:=",ExpressionUUID->"47d21065-13a5-4b2b-98f6-aeb3fbfd2b51"]
}, Open  ]],

Cell[CellGroupData[{

Cell["RMS", "Subsection",
 CellChangeTimes->{{3.870941553622549*^9, 
  3.8709415549363317`*^9}},ExpressionUUID->"53272786-40cd-4f50-ab15-\
7a386e30b2df"],

Cell[BoxData[
 RowBox[{
  RowBox[{"RMS", "[", 
   RowBox[{"R_", ",", "n_", ",", "m_"}], "]"}], ":=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"FR", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}], "/", 
    RowBox[{"(", 
     RowBox[{"n", "-", "m"}], ")"}]}], "]"}]}]], "Input",
 CellLabel->"In[2]:=",ExpressionUUID->"0bd03155-ebb6-412b-a6cb-6b20ef85d23d"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Joined parameters", "Subsection",
 CellChangeTimes->{{3.9049010362907844`*^9, 
  3.904901050176361*^9}},ExpressionUUID->"dc210d35-2b83-4392-ab6b-\
b8809f2fc4c8"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Tabc", "[", 
   RowBox[{"aux1_", ",", "aux2_", ",", "l1_", ",", "l2_"}], "]"}], ":=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{"1", ";;", "l1"}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["a", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{"l1", "+", "1"}], ";;", 
       RowBox[{"2", "*", "l1"}]}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["c", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"i", "+", "l2"}], ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{
        RowBox[{"2", "*", "l1"}], "+", "1"}], ";;", 
       RowBox[{"3", "*", "l1"}]}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["b", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"i", "+", 
            RowBox[{"2", "*", "l2"}]}], ",", "2"}], "]"}], "]"}]}], "}"}], 
      ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}]}], "]"}]}]], "Input",\
ExpressionUUID->"422317cb-d461-417d-a7f0-3f5a8e767875"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Data", "Section",
 CellChangeTimes->{{3.8709418034969025`*^9, 
  3.870941804124886*^9}},ExpressionUUID->"9c344c06-58c8-4128-8d89-\
219690e1ae32"],

Cell["\<\
Data import (output from Maple).
n (week)
ac (accumulated contamination)
ad (accumulated death)
av (accumulated vaccine)
cw (contamination per week)
dw (death per week)
vw (vaccine per week)\
\>", "Text",
 CellChangeTimes->{{3.864217791423253*^9, 3.8642177986151648`*^9}, {
  3.867832882221464*^9, 3.8678330270923495`*^9}, {3.867833071202381*^9, 
  3.8678331623660336`*^9}, {3.8704433697334566`*^9, 
  3.870443376839944*^9}},ExpressionUUID->"ddbec865-54cf-447b-a02a-\
37d7cdd5e972"],

Cell[BoxData[
 RowBox[{"DataBR", "=", 
  RowBox[{
   RowBox[{"Import", "[", 
    RowBox[{"\"\<D:\\\\sousa\\\\Covid19\\\\GnuPlot\\\\CovidBR.txt\>\"", ",", 
     RowBox[{"{", "\"\<Data\>\"", "}"}]}], "]"}], " ", "//", 
   "Short"}]}]], "Input",
 CellChangeTimes->{{3.875711384956251*^9, 3.875711391116248*^9}, {
   3.893357787545845*^9, 3.893357793121047*^9}, {3.8944841442436514`*^9, 
   3.894484158524002*^9}, 3.8944842082693315`*^9},
 CellLabel->"In[3]:=",ExpressionUUID->"980a1573-88cf-4917-88bc-c1dbb49d1b83"],

Cell[BoxData[
 RowBox[{"nBR", "=", 
  RowBox[{"Length", "[", 
   RowBox[{"DataBR", "[", 
    RowBox[{"[", "1", "]"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8757114099543405`*^9, 3.8757114160843363`*^9}},
 CellLabel->"In[4]:=",ExpressionUUID->"54e70c49-4f26-48b1-abfc-14877e4d3bb9"],

Cell[CellGroupData[{

Cell["Contamination", "Subsection",
 CellChangeTimes->{{3.8709419566949797`*^9, 3.870941959869384*^9}, 
   3.8709420359170327`*^9},ExpressionUUID->"db6dd233-28c1-46ce-99c3-\
4206c96e4105"],

Cell[BoxData[
 RowBox[{"DataAll", "=", "DataBR"}]], "Input",
 CellChangeTimes->{{3.885303593735478*^9, 3.885303601126498*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"7d28f777-8ad8-4b2d-a652-b29e73251dbf"],

Cell[BoxData[
 RowBox[{"n", "=", "nBR"}]], "Input",
 CellChangeTimes->{{3.8770884138340664`*^9, 3.877088416334715*^9}, {
  3.8853036125486765`*^9, 3.885303614955123*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"d8ad62c2-5ddc-41be-b43e-02c255fe0644"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"DataAll", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "i", ",", "5"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.867833607979697*^9, 3.8678336417547903`*^9}, 
   3.867834042482584*^9, 3.867834167808469*^9, {3.869228379395009*^9, 
   3.8692283841724577`*^9}, {3.8692285327297277`*^9, 
   3.8692285335105352`*^9}, {3.8696639875185523`*^9, 3.869663991411394*^9}, {
   3.8707908150861864`*^9, 3.870790820243222*^9}, {3.870891869329871*^9, 
   3.870891887445064*^9}, {3.8708919270279016`*^9, 3.8708919310272264`*^9}, {
   3.877597518035524*^9, 3.8775975217233467`*^9}, {3.8853036238465877`*^9, 
   3.885303627659087*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"3f170062-35df-41b4-ad55-16984b8c95b3"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Datab", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"DataAll", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "i", ",", "5"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870891949547256*^9, 3.870891955747621*^9}, {
  3.8775975418111525`*^9, 3.877597544376419*^9}, {3.885303638987633*^9, 
  3.8853036405189724`*^9}},
 CellLabel->"In[27]:=",ExpressionUUID->"0f1e731c-9ae9-4b51-a5fc-1464839c4b4f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gc", "=", 
   RowBox[{"BarChart", "[", 
    RowBox[{"Data", ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"BarSpacing", "\[Rule]", "Large"}], ",", 
     RowBox[{"ChartLabels", "\[Rule]", 
      RowBox[{"Table", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"ChartStyle", "\[Rule]", "Blue"}], ",", 
     RowBox[{"Frame", "->", 
      RowBox[{"{", 
       RowBox[{"True", ",", "True", ",", "False", ",", "False"}], "}"}]}], 
     ",", 
     RowBox[{"FrameTicks", "->", "False"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "n"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"1.5", "*", 
           SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.8703379978267174`*^9, 3.8703380290143604`*^9}, 
   3.870338677547329*^9, 3.8707912196612453`*^9, {3.8708776637366166`*^9, 
   3.8708777320344996`*^9}, {3.870877808093382*^9, 3.8708778242955365`*^9}, 
   3.870878218443729*^9, {3.87089197651674*^9, 3.870891987903799*^9}, 
   3.8708920591276574`*^9, 3.8708921042647977`*^9, 3.8708921494281025`*^9, {
   3.870893360390463*^9, 3.8708933823700247`*^9}, {3.8770874161586857`*^9, 
   3.8770874163488903`*^9}, 3.8770884440135345`*^9, 3.877597567357426*^9, 
   3.8853036629569817`*^9, {3.8853037007712545`*^9, 3.8853037754966154`*^9}, {
   3.8853046244511027`*^9, 3.8853046355775223`*^9}, {3.893357868051364*^9, 
   3.8933578813730707`*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"93adb5a8-31bc-47c2-9007-b0ca9312aeba"],

Cell[BoxData["gc"], "Input",
 CellChangeTimes->{{3.894484265540986*^9, 3.894484265626023*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"bafb72b9-4a4a-4f67-9032-3cc670c401e7"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Deaths", "Subsection",
 CellChangeTimes->{{3.870942006800601*^9, 
  3.8709420080954022`*^9}},ExpressionUUID->"251fb970-e571-4ea5-af54-\
71a1ac49be81"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"DataAll", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "i", ",", "6"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.867833607979697*^9, 3.8678336417547903`*^9}, 
   3.867834042482584*^9, 3.867834167808469*^9, {3.869228379395009*^9, 
   3.8692283841724577`*^9}, {3.8692285327297277`*^9, 
   3.8692285335105352`*^9}, {3.8696639875185523`*^9, 3.869663991411394*^9}, {
   3.8707908150861864`*^9, 3.870790820243222*^9}, {3.870891869329871*^9, 
   3.870891887445064*^9}, {3.8708919270279016`*^9, 3.8708919310272264`*^9}, {
   3.87094213776725*^9, 3.8709421385938425`*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"575b742e-f6ff-42b1-9c1f-dd0ac86b7797"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Datab", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"DataAll", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "i", ",", "6"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870891949547256*^9, 3.870891955747621*^9}, {
  3.870942154736216*^9, 3.870942155364955*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"74bb772b-acfd-4111-8f3b-66b66a492dce"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gd", "=", 
   RowBox[{"BarChart", "[", 
    RowBox[{"Data", ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"BarSpacing", "\[Rule]", "Large"}], ",", 
     RowBox[{"ChartLabels", "\[Rule]", 
      RowBox[{"Table", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ",", 
     RowBox[{"ChartStyle", "\[Rule]", "Blue"}], ",", 
     RowBox[{"Frame", "->", 
      RowBox[{"{", 
       RowBox[{"True", ",", "True", ",", "False", ",", "False"}], "}"}]}], 
     ",", 
     RowBox[{"FrameTicks", "->", "False"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "140"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"24", "*", 
           SuperscriptBox["10", "3"]}]}], "}"}]}], "}"}]}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.8703379978267174`*^9, 3.8703380290143604`*^9}, 
   3.870338677547329*^9, 3.8707912196612453`*^9, {3.8708776637366166`*^9, 
   3.8708777320344996`*^9}, {3.870877808093382*^9, 3.8708778242955365`*^9}, 
   3.870878218443729*^9, {3.87089197651674*^9, 3.870891987903799*^9}, 
   3.8708920591276574`*^9, 3.8708921042647977`*^9, 3.8708921494281025`*^9, {
   3.870893360390463*^9, 3.8708933823700247`*^9}, {3.8709421691568766`*^9, 
   3.87094219654668*^9}, 3.870942257629178*^9, 3.871053120411022*^9, {
   3.8710551323925123`*^9, 3.8710551571969447`*^9}},
 CellLabel->"In[13]:=",ExpressionUUID->"73f1e0f4-fc18-48e7-80f0-9136d9878e48"],

Cell[BoxData["gd"], "Input",
 CellChangeTimes->{{3.904896674136001*^9, 3.9048966742610025`*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"30dcb392-f692-4697-be0c-86025840ddf8"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Contaminations", "Section",
 CellChangeTimes->{{3.870941623488143*^9, 3.870941624781828*^9}, {
  3.8710566572013187`*^9, 
  3.871056661255905*^9}},ExpressionUUID->"c72d7ca5-a1cb-4eec-a1db-\
1087298cfd1f"],

Cell[CellGroupData[{

Cell["First part", "Subsection",
 CellChangeTimes->{{3.8709441005548306`*^9, 
  3.870944103121689*^9}},ExpressionUUID->"d4d8a613-d3ca-4a7e-b8dd-\
99e76939338b"],

Cell[BoxData[
 RowBox[{"m1", "=", " ", "96", " ", 
  RowBox[{"(*", "BR", "*)"}]}]], "Input",
 CellChangeTimes->{{3.885303853688653*^9, 3.885303876799182*^9}, {
  3.8853039412389402`*^9, 3.8853039554738045`*^9}, {3.885304648156492*^9, 
  3.8853046494377527`*^9}},
 CellLabel->"In[80]:=",ExpressionUUID->"a3644006-c726-4b08-a19d-da718b514fa2"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data1", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"DataAll", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "i", ",", "5"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "m1"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.867833607979697*^9, 3.8678336417547903`*^9}, 
   3.867834042482584*^9, 3.867834167808469*^9, {3.869228379395009*^9, 
   3.8692283841724577`*^9}, {3.8692285327297277`*^9, 
   3.8692285335105352`*^9}, {3.8696639875185523`*^9, 3.869663991411394*^9}, {
   3.8707908150861864`*^9, 3.870790820243222*^9}, {3.870891869329871*^9, 
   3.870891887445064*^9}, {3.8708919270279016`*^9, 3.8708919310272264`*^9}, {
   3.87094213776725*^9, 3.8709421385938425`*^9}, {3.8709444109197607`*^9, 
   3.8709444181136146`*^9}, {3.8770888680926695`*^9, 3.877088868188154*^9}, {
   3.8853039688966703`*^9, 3.8853039690060587`*^9}},
 CellLabel->"In[81]:=",ExpressionUUID->"6738a859-1b22-42cf-b082-a48fda9110c8"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data1b", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"DataAll", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "i", ",", "5"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "m1"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870891949547256*^9, 3.870891955747621*^9}, {
  3.870942154736216*^9, 3.870942155364955*^9}, {3.870944422170905*^9, 
  3.870944429072955*^9}, {3.877088871350814*^9, 3.877088871466403*^9}, {
  3.8853039785222006`*^9, 3.885303978647207*^9}},
 CellLabel->"In[82]:=",ExpressionUUID->"e2d9512f-534c-404e-8f63-766e589a7778"],

Cell[BoxData[
 RowBox[{"l1", "=", "10"}]], "Input",
 CellChangeTimes->{{3.870941700057722*^9, 3.870941703626641*^9}, {
   3.870943576800144*^9, 3.870943577707052*^9}, {3.8709442099463186`*^9, 
   3.870944220344717*^9}, {3.8853047181917753`*^9, 3.8853047186448994`*^9}, {
   3.885305799588354*^9, 3.8853058001039734`*^9}, {3.885305854624651*^9, 
   3.885305855921527*^9}, {3.8853060259949713`*^9, 3.8853060266668434`*^9}, {
   3.885307432116104*^9, 3.8853074332411003`*^9}, {3.8853083274055295`*^9, 
   3.8853083303742857`*^9}, {3.885308405159259*^9, 3.8853084058000755`*^9}, {
   3.885308511759242*^9, 3.8853085125717373`*^9}, {3.8853086247009425`*^9, 
   3.885308624810318*^9}, {3.885308672032172*^9, 3.8853086730790586`*^9}, {
   3.885308729004796*^9, 3.885308729082919*^9}, {3.8853089268103523`*^9, 
   3.885308927685171*^9}, {3.885309847597536*^9, 3.8853098479098406`*^9}, {
   3.8853102226079473`*^9, 3.8853102227173147`*^9}, {3.8933597295444546`*^9, 
   3.893359730105548*^9}, {3.8933601503149147`*^9, 3.893360151534171*^9}, {
   3.893360202131283*^9, 3.8933602026804967`*^9}, {3.8933602651385827`*^9, 
   3.893360265853343*^9}, {3.8933603229097652`*^9, 3.893360323027627*^9}, {
   3.8933621421327477`*^9, 3.893362142680525*^9}, 3.904897675661442*^9, {
   3.904897739908194*^9, 3.9048977404549136`*^9}},
 CellLabel->
  "In[101]:=",ExpressionUUID->"919a2a64-fcd6-4a02-bb5f-c4a64de8d9f1"],

Cell[BoxData[
 RowBox[{
  RowBox[{"mf1", "=", 
   RowBox[{"cshn", "[", 
    RowBox[{"z", ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867832129979463*^9, 3.867832158572497*^9}, 
   3.8708918194752913`*^9, {3.8709442342350492`*^9, 3.8709442476434917`*^9}},
 CellLabel->
  "In[102]:=",ExpressionUUID->"e31eb85c-ab8e-4ea3-9b33-b8a95000e4ad"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars1", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867834570829482*^9, 3.8678345775014815`*^9}, {
  3.870942389746482*^9, 3.8709424187653246`*^9}, {3.8709442605776134`*^9, 
  3.870944290983293*^9}},
 CellLabel->
  "In[103]:=",ExpressionUUID->"6ab1d067-a4b2-4f2e-b517-208c21e9bf9d"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cons1", "=", 
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.01", "<", 
         SubscriptBox["a", "i"], "<", "0.96"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"1", "<", 
         SubscriptBox["c", "i"], "<", "80"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox["10", "1"], "<", 
         SubscriptBox["b", "i"], "<", 
         SuperscriptBox["10", "7"]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          SubscriptBox["c", "i"], "/", 
          SubscriptBox["a", "i"]}], "<", "m1"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], ";"}],
   " "}]], "Input",
 CellChangeTimes->{{3.885304247672944*^9, 3.8853042565635853`*^9}, {
   3.8853048040867057`*^9, 3.885304804680252*^9}, {3.8853048555117745`*^9, 
   3.885304855683449*^9}, {3.8853049250797567`*^9, 3.8853049591151047`*^9}, {
   3.8853050568741612`*^9, 3.8853050570460224`*^9}, {3.8853051106279087`*^9, 
   3.885305114737282*^9}, {3.8853051826457787`*^9, 3.885305182802038*^9}, {
   3.8853052256177626`*^9, 3.8853052259615107`*^9}, {3.885305275559865*^9, 
   3.885305276700673*^9}, {3.88530533367214*^9, 3.885305334172145*^9}, {
   3.885305393300015*^9, 3.885305393971732*^9}, {3.885305436287636*^9, 
   3.885305436944083*^9}, {3.8853055231201305`*^9, 3.88530552393264*^9}, {
   3.8853056116411304`*^9, 3.8853056117661314`*^9}, {3.8853056549250956`*^9, 
   3.8853056557065487`*^9}, {3.8853057007558928`*^9, 3.885305701365469*^9}, {
   3.8853059187539897`*^9, 3.8853059193008614`*^9}, {3.8853059656476965`*^9, 
   3.885305966303751*^9}, {3.885306080029665*^9, 3.8853060817640305`*^9}, {
   3.8853061308621893`*^9, 3.8853061753032146`*^9}, {3.885306219009184*^9, 
   3.8853062213529425`*^9}, {3.885306362576228*^9, 3.88530636415434*^9}, {
   3.88530770743559*^9, 3.8853077085920277`*^9}, {3.8853077393285723`*^9, 
   3.8853077399223213`*^9}, {3.885307795330293*^9, 3.885307795971113*^9}, 
   3.885307847645425*^9, {3.8853079078509197`*^9, 3.8853079084759207`*^9}, {
   3.885307971884968*^9, 3.8853079767755923`*^9}, {3.8853080104656305`*^9, 
   3.8853080119812536`*^9}, {3.8853080908430767`*^9, 
   3.8853080918431015`*^9}, {3.885308137814732*^9, 3.8853081628955746`*^9}, {
   3.885308194132041*^9, 3.8853081959133015`*^9}, {3.8853082425723615`*^9, 
   3.8853082431350527`*^9}, {3.885308284966072*^9, 3.885308285622325*^9}, {
   3.8853085298534937`*^9, 3.8853085299628735`*^9}, {3.8853088056951523`*^9, 
   3.8853088061951623`*^9}, {3.8853088724323473`*^9, 3.8853088728852577`*^9}, 
   3.8853089959898567`*^9, {3.8853097202176213`*^9, 3.885309720858058*^9}, {
   3.8853097586577654`*^9, 3.8853097626421356`*^9}, {3.885309922804079*^9, 
   3.885309928819701*^9}, {3.885309970477724*^9, 3.885309971165225*^9}, {
   3.8933589572245283`*^9, 3.893358958659233*^9}, {3.8933591358220215`*^9, 
   3.8933591359182024`*^9}, {3.8933591800516405`*^9, 3.893359180176862*^9}, {
   3.893359254006399*^9, 3.893359254914632*^9}, {3.893359330267268*^9, 
   3.893359330742404*^9}, {3.8933594355298257`*^9, 3.8933594366323586`*^9}, {
   3.893359499054765*^9, 3.893359499627181*^9}, {3.8933595697846518`*^9, 
   3.893359576007849*^9}, {3.8933596397204404`*^9, 3.893359639811404*^9}, {
   3.893359696387888*^9, 3.8933596975696697`*^9}, {3.8933598029658904`*^9, 
   3.8933598368125577`*^9}, {3.893359943201434*^9, 3.893359944314442*^9}, {
   3.893359991164522*^9, 3.8933599919295316`*^9}, {3.8933600397556343`*^9, 
   3.893360040801405*^9}, {3.893360083908104*^9, 3.893360084568057*^9}, {
   3.893360393979486*^9, 3.8933603945074177`*^9}, {3.893360452735924*^9, 
   3.8933604529123516`*^9}, {3.8933604956407957`*^9, 
   3.8933604964254675`*^9}, {3.893360548497713*^9, 3.8933605499286695`*^9}, {
   3.8933621827043734`*^9, 3.8933622016356516`*^9}, {3.9048977814003024`*^9, 
   3.904897781963*^9}, {3.904897846600715*^9, 3.9048978466161814`*^9}, {
   3.9048978821233764`*^9, 3.9048978956419225`*^9}, {3.904897934805147*^9, 
   3.9048979355869474`*^9}, {3.904897965921009*^9, 3.90489796648332*^9}},
 CellLabel->
  "In[127]:=",ExpressionUUID->"cd1e4f9f-3a61-44bb-ac03-b564a3cecf7f"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"t0", ",", "lsf1"}], "}"}], "=", 
   RowBox[{"AbsoluteTiming", "[", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"Data1", ",", 
      RowBox[{"{", 
       RowBox[{"mf1", ",", "Cons1"}], "}"}], ",", "Pars1", ",", 
      RowBox[{"{", "z", "}"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "40000"}]}], "]"}], "]"}]}], ";", 
  " ", "t0"}]], "Input",
 CellChangeTimes->{{3.8696649023379245`*^9, 3.8696649531392283`*^9}, {
   3.8697258207695074`*^9, 3.869725821480483*^9}, {3.8697265516411*^9, 
   3.8697265520491714`*^9}, {3.869732294992208*^9, 3.86973229683278*^9}, {
   3.869733165401655*^9, 3.8697331661405754`*^9}, {3.8697373580290446`*^9, 
   3.869737358697692*^9}, {3.869744460132951*^9, 3.869744470596495*^9}, {
   3.8697523910705624`*^9, 3.869752393040555*^9}, {3.869755922118797*^9, 
   3.8697559238834257`*^9}, {3.8697561207431426`*^9, 3.869756121493045*^9}, {
   3.869756160011915*^9, 3.869756168565979*^9}, {3.86975627964835*^9, 
   3.869756281294915*^9}, 3.869838664700287*^9, {3.870265616754628*^9, 
   3.870265616882077*^9}, {3.87026726926554*^9, 3.870267269435379*^9}, {
   3.870892909198313*^9, 3.870892910305993*^9}, 3.8708936939969606`*^9, 
   3.870894335050465*^9, {3.8709443474842157`*^9, 3.870944356400553*^9}, 
   3.870944445627025*^9, {3.8853050280109634`*^9, 3.8853050281517916`*^9}, {
   3.893359313922022*^9, 3.8933593147862034`*^9}, {3.8944845308631086`*^9, 
   3.8944845322624216`*^9}, {3.9048975997234635`*^9, 3.904897599801575*^9}},
 CellLabel->
  "In[128]:=",ExpressionUUID->"3ed0b713-0c5b-4ea2-b5a3-8672433c3d83"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"FR", "=", 
   RowBox[{"lsf1", "[", "\"\<FitResiduals\>\"", "]"}]}], ";"}], "\n", 
 RowBox[{"lsf1", "[", 
  RowBox[{"\"\<EstimatedVariance\>\"", ",", 
   RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"RMS", "[", 
       RowBox[{"#", ",", "m1", ",", 
        RowBox[{"3", "*", "l1"}]}], "]"}], "&"}], ")"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.867834878616255*^9, 3.867834878772255*^9}, 
   3.8692243637591076`*^9, {3.869224444478524*^9, 3.8692244539570875`*^9}, {
   3.869229367657488*^9, 3.869229370946748*^9}, {3.870944594624638*^9, 
   3.8709446032108283`*^9}, {3.8770890473344407`*^9, 
   3.8770890602566085`*^9}, {3.8853040677297254`*^9, 3.8853040677922115`*^9}},
 CellLabel->
  "In[129]:=",ExpressionUUID->"9105c868-3b2a-4597-8225-05d441f766fc"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"lsf1", "[", "i", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"i", ",", "1", ",", "100"}], "}"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"PlotRange", "->", "Full"}], ",", 
   RowBox[{"Epilog", "->", 
    RowBox[{"Point", "[", "Data1b", "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8708070875569153`*^9, 3.8708070877519207`*^9}, {
  3.8708927694670935`*^9, 3.8708927701268067`*^9}, {3.870944630431551*^9, 
  3.8709446429311743`*^9}, {3.8853047007061033`*^9, 3.8853047014404097`*^9}},
 CellLabel->
  "In[131]:=",ExpressionUUID->"c4d5c03b-8dd4-42f8-9de4-afb63c3067bf"],

Cell[BoxData[
 RowBox[{"pars1", "=", 
  RowBox[{"lsf1", "[", "\"\<BestFitParameters\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.8709447037162905`*^9, 3.870944706183713*^9}},
 CellLabel->
  "In[132]:=",ExpressionUUID->"57a46e14-df0e-4679-80dc-0c271f6540c5"],

Cell[BoxData[
 RowBox[{"Pars1", "=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["a", "i"], ",", 
        RowBox[{
         SubscriptBox["a", "i"], "/.", "pars1"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["c", "i"], ",", 
        RowBox[{
         SubscriptBox["c", "i"], "/.", "pars1"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["b", "i"], ",", 
        RowBox[{
         SubscriptBox["b", "i"], "/.", "pars1"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8703383358095846`*^9, 3.8703383359032965`*^9}, 
   3.8703385368297215`*^9, 3.8708069708168297`*^9, {3.8709447711598516`*^9, 
   3.870944777123497*^9}, {3.8770995071274195`*^9, 3.8770995162602825`*^9}},
 CellLabel->
  "In[133]:=",ExpressionUUID->"1f2b68f9-5d7a-4ed7-ac66-7cf93ec5aafc"],

Cell[BoxData[
 RowBox[{"Pars1BR", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "1"], ",", "0.3442544314568885`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "2"], ",", "0.06265025791695183`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "3"], ",", "0.11751913671825302`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "4"], ",", "0.3695226437922263`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "5"], ",", "0.19378419843330846`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "6"], ",", "0.29894114249384784`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "7"], ",", "0.18752557321622793`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "8"], ",", "0.21221509339360087`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "9"], ",", "0.5971713010553253`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "10"], ",", "0.9600000099999672`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "1"], ",", "14.38120931929726`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "2"], ",", "4.351620378460686`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "3"], ",", "3.7478183168053985`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "4"], ",", "17.54671424061865`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "5"], ",", "10.903826746993905`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "6"], ",", "19.457759706853334`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "7"], ",", "4.446031475614506`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "8"], ",", "3.516623652679805`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "9"], ",", "41.47475464006512`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "10"], ",", "71.06632302071381`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "1"], ",", "633664.3243725311`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "2"], ",", "4.0044840715635186`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "3"], ",", "1.135564175304713`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "4"], ",", "686639.3094733317`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "5"], ",", "1.9440307014502825`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "6"], ",", "549142.436488662`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "7"], ",", "1.178598948495137`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "8"], ",", "680689.1254867393`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "9"], ",", "370152.93649694487`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "10"], ",", "79019.55848215359`"}], "}"}]}], 
   "}"}]}]], "Input",
 CellChangeTimes->{{3.8933587887212152`*^9, 3.8933588162926664`*^9}, 
   3.8933621305306177`*^9, 3.8933622464930716`*^9, 3.89448461630927*^9, 
   3.904897657203971*^9, 3.904898012945341*^9},
 CellLabel->
  "In[134]:=",ExpressionUUID->"8e8b57ef-6f28-48ad-be27-d4173bae4035"],

Cell["Go Maple!", "Text",
 CellChangeTimes->{{3.8697246415680885`*^9, 
  3.869724644647354*^9}},ExpressionUUID->"11bab72b-c843-498d-9164-\
6321269a8360"],

Cell[BoxData[
 RowBox[{
  RowBox[{"ParsMaple", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "i"], "/.", "pars1"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["c", "i"], "/.", "pars1"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["b", "i"], "/.", "pars1"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867848620243645*^9, 3.8678486971683884`*^9}, {
   3.8679151951023393`*^9, 3.8679151971583385`*^9}, 3.869226041175602*^9, 
   3.87761538957613*^9, {3.8944847109876146`*^9, 3.8944847651704254`*^9}},
 CellLabel->
  "In[176]:=",ExpressionUUID->"1fa3207b-9d89-47ac-bfa1-bdd2a6414e32"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Export", "[", 
   RowBox[{
   "\"\<D:\\\\sousa\\\\Covid19\\\\GnuPlot\\\\ParsBRp1.txt\>\"", ",", 
    "ParsMaple", ",", "\"\<CSV\>\""}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.867915697525957*^9, 3.8679157280658646`*^9}, 
   3.8679157628765945`*^9, {3.8679164509181013`*^9, 3.867916451654106*^9}, 
   3.867916808995467*^9, 3.867917142053749*^9, {3.867917632604063*^9, 
   3.8679176413030043`*^9}, {3.867929514447388*^9, 3.867929516560387*^9}, {
   3.867930359137655*^9, 3.867930360405655*^9}, {3.869042378475341*^9, 
   3.8690423865541596`*^9}, {3.8690498646245327`*^9, 3.869049868272987*^9}, 
   3.869226047759435*^9, {3.869664545710162*^9, 3.869664546845023*^9}, {
   3.8697257951606436`*^9, 3.8697257954476876`*^9}, {3.869751738310028*^9, 
   3.869751744380788*^9}, 3.8697558833468914`*^9, {3.869836556856716*^9, 
   3.869836556925766*^9}, {3.8703384927006083`*^9, 3.8703384969338017`*^9}, {
   3.870339143253948*^9, 3.8703391440818763`*^9}, {3.87036072023104*^9, 
   3.870360723804949*^9}, 3.870365323681324*^9, {3.870448556460325*^9, 
   3.870448556651752*^9}, {3.870879015485437*^9, 3.8708790192284355`*^9}, {
   3.870893538346167*^9, 3.870893539412013*^9}, 3.87104557858768*^9, {
   3.8771008479686832`*^9, 3.877100850810518*^9}, {3.894484730195363*^9, 
   3.89448473581075*^9}, 3.894484802723254*^9},
 CellLabel->
  "In[177]:=",ExpressionUUID->"bec225a2-9ef8-4ad0-8be2-c60327f8742f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf1", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"lsf1", "[", "i", "]"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "n"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"1.5", "*", 
           SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", "Red", "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{
        "{", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \(a\)]\)(n)\>\"", 
         "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Contamination/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870893042085024*^9, 3.870893071424758*^9}, {
  3.8708932351373615`*^9, 3.870893263626501*^9}, {3.870893298766883*^9, 
  3.8708932988554955`*^9}, {3.8709636554052777`*^9, 3.8709637057698674`*^9}, {
  3.870963784521944*^9, 3.870963790720209*^9}, {3.8709638331309447`*^9, 
  3.8709638335118675`*^9}, {3.8770995365775566`*^9, 3.877099540165415*^9}, {
  3.8853043973507404`*^9, 3.8853044067417297`*^9}, {3.885304530573325*^9, 
  3.885304533760824*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"e9485626-019c-4a45-b231-29b0cbfff40f"],

Cell[BoxData[
 RowBox[{"gfc1", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf1", ",", "gc"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8708930927308626`*^9, 3.8708930987028093`*^9}, {
  3.8709637151085815`*^9, 3.870963719673771*^9}, {3.870963876897517*^9, 
  3.8709638841272526`*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"7f856a76-b2fa-4725-a20a-a52ee624559e"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Second part", "Subsection",
 CellChangeTimes->{{3.8709441005548306`*^9, 3.870944103121689*^9}, {
  3.8709639383725348`*^9, 
  3.87096394221332*^9}},ExpressionUUID->"898d8142-4872-4e3e-b06a-\
f54f94a176e7"],

Cell[BoxData[
 RowBox[{"n2", "=", "96"}]], "Input",
 CellChangeTimes->{{3.8775980142897835`*^9, 3.877598017741531*^9}, {
  3.9048969295306463`*^9, 3.9048969314843936`*^9}},
 CellLabel->"In[34]:=",ExpressionUUID->"cc359a44-f488-4c9d-8a11-5974229044a0"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data2", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Data", "[", 
      RowBox[{"[", "i", "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "n2", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.867833607979697*^9, 3.8678336417547903`*^9}, 
   3.867834042482584*^9, 3.867834167808469*^9, {3.869228379395009*^9, 
   3.8692283841724577`*^9}, {3.8692285327297277`*^9, 
   3.8692285335105352`*^9}, {3.8696639875185523`*^9, 3.869663991411394*^9}, {
   3.8707908150861864`*^9, 3.870790820243222*^9}, {3.870891869329871*^9, 
   3.870891887445064*^9}, {3.8708919270279016`*^9, 3.8708919310272264`*^9}, {
   3.87094213776725*^9, 3.8709421385938425`*^9}, {3.8709444109197607`*^9, 
   3.8709444181136146`*^9}, {3.8709639486815453`*^9, 3.870963961045286*^9}, {
   3.8770885127216387`*^9, 3.87708851281291*^9}, {3.8775979842738657`*^9, 
   3.877597985965374*^9}, {3.8775980310197105`*^9, 3.877598050100524*^9}, {
   3.90489691349618*^9, 3.9048969190130215`*^9}},
 CellLabel->"In[35]:=",ExpressionUUID->"591fa2db-2ae7-4594-b898-ef97114f8d2a"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data2b", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"Data", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "n2", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870891949547256*^9, 3.870891955747621*^9}, {
  3.870942154736216*^9, 3.870942155364955*^9}, {3.870944422170905*^9, 
  3.870944429072955*^9}, {3.870963971662084*^9, 3.870963980655541*^9}, {
  3.8770885165320873`*^9, 3.877088516656339*^9}, {3.8775980687937145`*^9, 
  3.8775980895430183`*^9}, {3.877598144366598*^9, 3.8775981627457943`*^9}, {
  3.904896922622941*^9, 3.9048969243264837`*^9}},
 CellLabel->"In[36]:=",ExpressionUUID->"2e275371-0b11-4809-8d63-0fc8848217d7"],

Cell[BoxData[
 RowBox[{"l1", "=", "10"}]], "Input",
 CellChangeTimes->{{3.877087804782305*^9, 3.8770878106827483`*^9}, {
  3.904898029875766*^9, 3.9048980303447547`*^9}},
 CellLabel->
  "In[147]:=",ExpressionUUID->"96b5e056-1cbc-4e35-b45f-902d0be5a8cf"],

Cell[BoxData[
 RowBox[{"l2", "=", "8"}]], "Input",
 CellChangeTimes->{{3.870941700057722*^9, 3.870941703626641*^9}, {
  3.870943576800144*^9, 3.870943577707052*^9}, {3.8709442099463186`*^9, 
  3.870944220344717*^9}, {3.87096398743102*^9, 3.870963987521147*^9}, {
  3.877598277243244*^9, 3.877598278296707*^9}, {3.8775999474783506`*^9, 
  3.8775999485783377`*^9}, {3.8776043932778254`*^9, 3.877604394020425*^9}, {
  3.90489706315802*^9, 3.904897063830085*^9}},
 CellLabel->
  "In[148]:=",ExpressionUUID->"c8410cea-0d57-4af9-8faa-456d7cfb36f1"],

Cell[BoxData[
 RowBox[{
  RowBox[{"mf2", "=", 
   RowBox[{"cshn", "[", 
    RowBox[{"z", ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867832129979463*^9, 3.867832158572497*^9}, 
   3.8708918194752913`*^9, {3.8709442342350492`*^9, 3.8709442476434917`*^9}, {
   3.8709639958539886`*^9, 3.870964006063445*^9}, {3.8709640539212065`*^9, 
   3.8709641051602545`*^9}},
 CellLabel->
  "In[149]:=",ExpressionUUID->"7b9b9d83-1d7b-4250-a272-5215425484e7"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars2", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867834570829482*^9, 3.8678345775014815`*^9}, {
  3.870942389746482*^9, 3.8709424187653246`*^9}, {3.8709442605776134`*^9, 
  3.870944290983293*^9}, {3.870964012793183*^9, 3.8709640234502068`*^9}, {
  3.8709641225539055`*^9, 3.870964129481372*^9}, {3.870964170592956*^9, 
  3.8709641774432645`*^9}},
 CellLabel->
  "In[150]:=",ExpressionUUID->"4d5ff2a2-02e8-4315-9690-f16b5cc9b4e7"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Cons2", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"0.01", "<", 
        SubscriptBox["a", "i"], "<", "0.9"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"1", "<", 
        SubscriptBox["c", "i"], "<", "80"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SuperscriptBox["10", "1"], "<", 
        SubscriptBox["b", "i"], "<", 
        SuperscriptBox["10", "7"]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.870943197443465*^9, 3.870943201625883*^9}, {
  3.870944283544839*^9, 3.8709443201375895`*^9}, {3.8709640311805105`*^9, 
  3.8709640348758936`*^9}, {3.870964133634597*^9, 3.8709641425285974`*^9}, {
  3.870964184136936*^9, 3.87096419317742*^9}},
 CellLabel->"In[74]:=",ExpressionUUID->"a358bd6e-2fd3-441d-aa07-e16a19918c5a"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cons2", "=", 
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.01", "<", 
         SubscriptBox["a", "i"], "<", "0.6"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.1", "<", 
         SubscriptBox["c", "i"], "<", "20"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox["10", "1"], "<", 
         SubscriptBox["b", "i"], "<", 
         SuperscriptBox["10", "7"]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          SubscriptBox["c", "i"], "/", 
          SubscriptBox["a", "i"]}], "<", "n"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], ";"}], 
  " "}]], "Input",
 CellChangeTimes->{{3.8702680472882824`*^9, 3.870268062880148*^9}, {
   3.870274527961333*^9, 3.8702745289124813`*^9}, {3.870274712729985*^9, 
   3.8702747138006673`*^9}, {3.8703380760634556`*^9, 3.870338085777099*^9}, {
   3.870338134356515*^9, 3.870338134781331*^9}, {3.870338715435112*^9, 
   3.870338740776287*^9}, {3.870359410836882*^9, 3.8703594114353127`*^9}, {
   3.870359494778078*^9, 3.870359503140952*^9}, 3.8703645089762597`*^9, {
   3.8704435176581383`*^9, 3.8704435249828515`*^9}, 3.8704465109072294`*^9, {
   3.8704471172830877`*^9, 3.870447117490363*^9}, {3.870448573777196*^9, 
   3.8704485746457186`*^9}, {3.8707868280045576`*^9, 3.8707869329824204`*^9}, 
   3.870892974501359*^9, {3.87089364835896*^9, 3.8708936623592963`*^9}, {
   3.870894444232833*^9, 3.8708944531899633`*^9}, {3.8708945052751646`*^9, 
   3.870894509928249*^9}, 3.8709424914766407`*^9, {3.8709443307752733`*^9, 
   3.870944339866149*^9}, {3.8709641474240932`*^9, 3.8709641588991756`*^9}, {
   3.8709641994977407`*^9, 3.8709642292079687`*^9}, {3.8770874902344704`*^9, 
   3.8770874967737436`*^9}, {3.8770878904304543`*^9, 3.877087890547224*^9}, {
   3.877087940184823*^9, 3.877087948274808*^9}, {3.877088035860699*^9, 
   3.877088036610381*^9}, {3.877088087761205*^9, 3.877088087856166*^9}, {
   3.8770886360573907`*^9, 3.8770886361904974`*^9}, 3.877088700692484*^9, {
   3.8775984010866833`*^9, 3.877598401251259*^9}, {3.8775988951982093`*^9, 
   3.8775988969679337`*^9}, {3.8775993225083437`*^9, 
   3.8775993226093397`*^9}, {3.8775993833177752`*^9, 3.877599384018194*^9}, {
   3.877599445377966*^9, 3.877599445448095*^9}, {3.877599518818184*^9, 
   3.8775995189683037`*^9}, {3.8775995637978663`*^9, 3.877599563951803*^9}, {
   3.877599806898203*^9, 3.8775998458581743`*^9}, {3.8775999986891384`*^9, 
   3.877600000848173*^9}, {3.877600043609413*^9, 3.877600044098344*^9}, {
   3.877600076618606*^9, 3.877600076898345*^9}, {3.8776043460862846`*^9, 
   3.8776043469278917`*^9}, {3.9048970820224433`*^9, 3.904897082132004*^9}, {
   3.9048972075790234`*^9, 3.9048972165341406`*^9}, {3.9048972677643414`*^9, 
   3.90489726784247*^9}, {3.904897326495208*^9, 3.9048973359503546`*^9}, {
   3.9048973872419024`*^9, 3.9048973918989*^9}, {3.9048974252655125`*^9, 
   3.9048974260779753`*^9}},
 CellLabel->
  "In[151]:=",ExpressionUUID->"67ec7e92-a3a0-4bca-8809-a0da020b024c"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"t0", ",", "lsf2"}], "}"}], "=", 
   RowBox[{"AbsoluteTiming", "[", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"Data2", ",", 
      RowBox[{"{", 
       RowBox[{"mf2", ",", "Cons2"}], "}"}], ",", "Pars2", ",", 
      RowBox[{"{", "z", "}"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "40000"}]}], "]"}], "]"}]}], ";", 
  " ", "t0"}]], "Input",
 CellChangeTimes->{{3.8696649023379245`*^9, 3.8696649531392283`*^9}, {
   3.8697258207695074`*^9, 3.869725821480483*^9}, {3.8697265516411*^9, 
   3.8697265520491714`*^9}, {3.869732294992208*^9, 3.86973229683278*^9}, {
   3.869733165401655*^9, 3.8697331661405754`*^9}, {3.8697373580290446`*^9, 
   3.869737358697692*^9}, {3.869744460132951*^9, 3.869744470596495*^9}, {
   3.8697523910705624`*^9, 3.869752393040555*^9}, {3.869755922118797*^9, 
   3.8697559238834257`*^9}, {3.8697561207431426`*^9, 3.869756121493045*^9}, {
   3.869756160011915*^9, 3.869756168565979*^9}, {3.86975627964835*^9, 
   3.869756281294915*^9}, 3.869838664700287*^9, {3.870265616754628*^9, 
   3.870265616882077*^9}, {3.87026726926554*^9, 3.870267269435379*^9}, {
   3.870892909198313*^9, 3.870892910305993*^9}, 3.8708936939969606`*^9, 
   3.870894335050465*^9, {3.8709443474842157`*^9, 3.870944356400553*^9}, 
   3.870944445627025*^9, {3.8709642380156755`*^9, 3.870964251855894*^9}, {
   3.90489744861427*^9, 3.904897448692395*^9}},
 CellLabel->
  "In[152]:=",ExpressionUUID->"fac56909-b1c8-42cb-a95a-1f5224613c2e"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"FR", "=", 
   RowBox[{"lsf2", "[", "\"\<FitResiduals\>\"", "]"}]}], ";"}], "\n", 
 RowBox[{"lsf2", "[", 
  RowBox[{"\"\<EstimatedVariance\>\"", ",", 
   RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"RMS", "[", 
       RowBox[{"#", ",", 
        RowBox[{"n", "-", "n2", "+", "1"}], ",", 
        RowBox[{"3", "*", "l2"}]}], "]"}], "&"}], ")"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.867834878616255*^9, 3.867834878772255*^9}, 
   3.8692243637591076`*^9, {3.869224444478524*^9, 3.8692244539570875`*^9}, {
   3.869229367657488*^9, 3.869229370946748*^9}, {3.870944594624638*^9, 
   3.8709446032108283`*^9}, {3.8709642777743073`*^9, 3.870964292671986*^9}, 
   3.8770875898825474`*^9, {3.877087660813163*^9, 3.87708766603794*^9}, {
   3.8770877489645805`*^9, 3.8770877548889027`*^9}, {3.9048969857132273`*^9, 
   3.904896987573018*^9}},
 CellLabel->
  "In[153]:=",ExpressionUUID->"26ebda72-b06a-45e0-8ca4-96517c355222"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"lsf2", "[", 
    RowBox[{"i", "-", "n2"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"i", ",", "90", ",", "n"}], "}"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"PlotRange", "->", "Full"}], ",", 
   RowBox[{"Epilog", "->", 
    RowBox[{"Point", "[", "Data2b", "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8708070875569153`*^9, 3.8708070877519207`*^9}, {
  3.8708927694670935`*^9, 3.8708927701268067`*^9}, {3.870944630431551*^9, 
  3.8709446429311743`*^9}, {3.870964310082224*^9, 3.8709643691971035`*^9}, {
  3.8770879812225513`*^9, 3.8770879820541553`*^9}, {3.9048970040610304`*^9, 
  3.9048970348328724`*^9}},
 CellLabel->
  "In[155]:=",ExpressionUUID->"94d53b18-115d-43c0-a73d-0bd67be707cc"],

Cell[BoxData[
 RowBox[{"pars2", "=", 
  RowBox[{"lsf2", "[", "\"\<BestFitParameters\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.8709447037162905`*^9, 3.870944706183713*^9}, {
  3.870964389688529*^9, 3.8709643930406733`*^9}},
 CellLabel->
  "In[156]:=",ExpressionUUID->"8efde151-d94d-42ca-bf34-b0f3a89d6f74"],

Cell[BoxData[
 RowBox[{"Pars2", "=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["a", "i"], ",", 
        RowBox[{
         SubscriptBox["a", "i"], "/.", "pars2"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"l1", "+", "1"}], ",", 
        RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["c", "i"], ",", 
        RowBox[{
         SubscriptBox["c", "i"], "/.", "pars2"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"l1", "+", "1"}], ",", 
        RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["b", "i"], ",", 
        RowBox[{
         SubscriptBox["b", "i"], "/.", "pars2"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"l1", "+", "1"}], ",", 
        RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8703383358095846`*^9, 3.8703383359032965`*^9}, 
   3.8703385368297215`*^9, 3.8708069708168297`*^9, {3.8709447711598516`*^9, 
   3.870944777123497*^9}, {3.870964415028035*^9, 3.870964484931066*^9}},
 CellLabel->
  "In[157]:=",ExpressionUUID->"efeb5087-f51d-4c95-bc6d-86748962196b"],

Cell["Go Maple!", "Text",
 CellChangeTimes->{{3.8697246415680885`*^9, 
  3.869724644647354*^9}},ExpressionUUID->"7c5c632c-b2ea-4e74-97fc-\
7435fb030ba7"],

Cell[BoxData[
 RowBox[{
  RowBox[{"ParsMaple2", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["a", 
         RowBox[{"i", "+", "l1"}]], "/.", "pars2"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["c", 
         RowBox[{"i", "+", "l1"}]], "/.", "pars2"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["b", 
         RowBox[{"i", "+", "l1"}]], "/.", "pars2"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867848620243645*^9, 3.8678486971683884`*^9}, {
   3.8679151951023393`*^9, 3.8679151971583385`*^9}, 3.869226041175602*^9, 
   3.87761538957613*^9, {3.8944847109876146`*^9, 3.8944847651704254`*^9}, {
   3.9049012180033426`*^9, 3.904901265245146*^9}, {3.9049013145776997`*^9, 
   3.904901332301627*^9}},
 CellLabel->
  "In[180]:=",ExpressionUUID->"76ece0fb-aa74-4bee-b252-701dfbcb7121"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Export", "[", 
   RowBox[{
   "\"\<D:\\\\sousa\\\\Covid19\\\\GnuPlot\\\\ParsBRp2.txt\>\"", ",", 
    "ParsMaple2", ",", "\"\<CSV\>\""}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.867915697525957*^9, 3.8679157280658646`*^9}, 
   3.8679157628765945`*^9, {3.8679164509181013`*^9, 3.867916451654106*^9}, 
   3.867916808995467*^9, 3.867917142053749*^9, {3.867917632604063*^9, 
   3.8679176413030043`*^9}, {3.867929514447388*^9, 3.867929516560387*^9}, {
   3.867930359137655*^9, 3.867930360405655*^9}, {3.869042378475341*^9, 
   3.8690423865541596`*^9}, {3.8690498646245327`*^9, 3.869049868272987*^9}, 
   3.869226047759435*^9, {3.869664545710162*^9, 3.869664546845023*^9}, {
   3.8697257951606436`*^9, 3.8697257954476876`*^9}, {3.869751738310028*^9, 
   3.869751744380788*^9}, 3.8697558833468914`*^9, {3.869836556856716*^9, 
   3.869836556925766*^9}, {3.8703384927006083`*^9, 3.8703384969338017`*^9}, {
   3.870339143253948*^9, 3.8703391440818763`*^9}, {3.87036072023104*^9, 
   3.870360723804949*^9}, 3.870365323681324*^9, {3.870448556460325*^9, 
   3.870448556651752*^9}, {3.870879015485437*^9, 3.8708790192284355`*^9}, {
   3.870893538346167*^9, 3.870893539412013*^9}, 3.87104557858768*^9, {
   3.8771008479686832`*^9, 3.877100850810518*^9}, {3.894484730195363*^9, 
   3.89448473581075*^9}, 3.894484802723254*^9, {3.9049012673941545`*^9, 
   3.9049012766681585`*^9}},
 CellLabel->
  "In[181]:=",ExpressionUUID->"fdfd2541-f325-4eaa-935c-4e3380c865cd"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf2", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"lsf2", "[", 
       RowBox[{"i", "-", "74"}], "]"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "69", ",", "150"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "150"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"6", "*", 
           SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", "Red", "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{
        "{", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \(b\)]\)(n)\>\"", 
         "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Contamination/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870893042085024*^9, 3.870893071424758*^9}, {
  3.8708932351373615`*^9, 3.870893263626501*^9}, {3.870893298766883*^9, 
  3.8708932988554955`*^9}, {3.8709636554052777`*^9, 3.8709637057698674`*^9}, {
  3.870963784521944*^9, 3.870963790720209*^9}, {3.8709638331309447`*^9, 
  3.8709638335118675`*^9}, {3.8709645215678253`*^9, 3.870964549699584*^9}, {
  3.877088815433672*^9, 3.877088821650298*^9}},
 CellLabel->
  "In[128]:=",ExpressionUUID->"93a2d356-f1fe-4709-9c12-be1d27061769"],

Cell[BoxData[
 RowBox[{"gfc2", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf2", ",", "gc"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8708930927308626`*^9, 3.8708930987028093`*^9}, {
  3.8709637151085815`*^9, 3.870963719673771*^9}, {3.870963876897517*^9, 
  3.8709638841272526`*^9}, {3.870964560150072*^9, 3.870964564545*^9}},
 CellLabel->
  "In[129]:=",ExpressionUUID->"7091e121-da3f-4501-a591-a4e36a5e2a4b"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf12", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"lsf1", "[", "i", "]"}], ",", 
       RowBox[{"lsf2", "[", 
        RowBox[{"i", "-", "74"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "140"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "140"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"6", "*", 
           SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Red", ",", "Brown"}], "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
         "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \(a\)]\)(n)\>\"", 
          ",", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \
\(b\)]\)(n)\>\""}], "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Contamination/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870964711371406*^9, 3.870964733600583*^9}, {
  3.870964783638962*^9, 3.870964880301608*^9}, {3.8709649138296537`*^9, 
  3.87096493651463*^9}, {3.870965147951193*^9, 3.8709651903825254`*^9}},
 CellLabel->
  "In[131]:=",ExpressionUUID->"cb8d7e11-0f44-4528-b617-63db8ef0482e"],

Cell[BoxData[
 RowBox[{"gfc12", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf12", ",", "gc"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.87096488840285*^9, 3.8709648918741126`*^9}},
 CellLabel->
  "In[132]:=",ExpressionUUID->"78b307fa-af54-40f7-b693-914721405c0d"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Joint parameters", "Subsection",
 CellChangeTimes->{{3.871043208282177*^9, 
  3.8710432117624273`*^9}},ExpressionUUID->"9bca91aa-a2f1-4d56-b884-\
426f694047c5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Tabc", "[", 
   RowBox[{"aux1_", ",", "aux2_", ",", "l1_", ",", "l2_"}], "]"}], ":=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{"1", ";;", "l1"}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["a", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{"l1", "+", "1"}], ";;", 
       RowBox[{"2", "*", "l1"}]}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["c", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"i", "+", "l2"}], ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{
        RowBox[{"2", "*", "l1"}], "+", "1"}], ";;", 
       RowBox[{"3", "*", "l1"}]}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["b", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"i", "+", 
            RowBox[{"2", "*", "l2"}]}], ",", "2"}], "]"}], "]"}]}], "}"}], 
      ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.87104457892337*^9, 3.87104471153677*^9}, 
   3.871044747799478*^9, {3.8710448500107675`*^9, 3.8710448741324635`*^9}, {
   3.9048991727469015`*^9, 3.904899184765229*^9}},
 CellLabel->
  "In[171]:=",ExpressionUUID->"fc152712-ee36-400a-a424-9767fbd1c621"],

Cell[BoxData[
 RowBox[{"Pars", "=", 
  RowBox[{"Tabc", "[", 
   RowBox[{"Pars1", ",", "Pars2", ",", "l1", ",", "l2"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.9048988835789995`*^9, 3.904898946091647*^9}, {
  3.9048990392197084`*^9, 3.9048990392822037`*^9}, {3.9048992058163624`*^9, 
  3.9048992070823984`*^9}},
 CellLabel->
  "In[173]:=",ExpressionUUID->"976ae3e0-57b6-4f17-bc69-1484359b4d6c"],

Cell[BoxData[
 RowBox[{"ParsM", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"Pars", "[", 
     RowBox[{"[", 
      RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"3", "*", 
       RowBox[{"(", 
        RowBox[{"l1", "+", "l2"}], ")"}]}]}], "}"}]}], "]"}]}]], "Input",
 CellLabel->
  "In[174]:=",ExpressionUUID->"28e2fc1b-f9f1-42ed-a30a-75f579303514"],

Cell[BoxData[
 RowBox[{"ParsM", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"Pars1BR", "[", 
     RowBox[{"[", 
      RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"3", "*", 
       RowBox[{"(", 
        RowBox[{"l1", "+", "0"}], ")"}]}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.893362493814966*^9, 3.893362503429675*^9}},
 CellLabel->
  "In[220]:=",ExpressionUUID->"3f6348e1-0fef-41ca-ad70-bcefe4ea95bd"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Export", "[", 
   RowBox[{
   "\"\<D:\\\\sousa\\\\Covid19\\\\GnuPlot\\\\ParsBR.txt\>\"", ",", "ParsM", 
    ",", "\"\<CSV\>\""}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.8770997847964945`*^9, 3.877099793386468*^9}, {
   3.893362388166809*^9, 3.893362391239275*^9}, 3.8933624263378887`*^9, {
   3.8933625183862476`*^9, 3.893362531747977*^9}, 3.904899277940061*^9},
 CellLabel->
  "In[175]:=",ExpressionUUID->"864be332-a77f-4354-8f69-a39f487ee6ba"],

Cell[BoxData[
 RowBox[{"ParsBRc", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "1"], ",", "0.3442544314568885`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "2"], ",", "0.06265025791695183`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "3"], ",", "0.11751913671825302`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "4"], ",", "0.3695226437922263`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "5"], ",", "0.19378419843330846`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "6"], ",", "0.29894114249384784`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "7"], ",", "0.18752557321622793`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "8"], ",", "0.21221509339360087`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "9"], ",", "0.5971713010553253`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "10"], ",", "0.9600000099999672`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "11"], ",", "0.6000000100054333`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "12"], ",", "0.22157316504842214`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "13"], ",", "0.19779096864618154`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "14"], ",", "0.6000000100018034`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "15"], ",", "0.3250359893137727`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "16"], ",", "0.05092284557133669`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "17"], ",", "0.6000000099999716`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "18"], ",", "0.6000000099999978`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "1"], ",", "14.38120931929726`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "2"], ",", "4.351620378460686`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "3"], ",", "3.7478183168053985`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "4"], ",", "17.54671424061865`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "5"], ",", "10.903826746993905`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "6"], ",", "19.457759706853334`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "7"], ",", "4.446031475614506`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "8"], ",", "3.516623652679805`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "9"], ",", "41.47475464006512`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "10"], ",", "71.06632302071381`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "11"], ",", "3.3242044737711267`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "12"], ",", "2.7965316615865374`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "13"], ",", "5.633124600538097`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "14"], ",", "4.027765699621836`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "15"], ",", "16.940564466764428`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "16"], ",", "3.022304386922389`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "17"], ",", "2.518949546658343`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "18"], ",", "5.436739265157799`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "1"], ",", "633664.3243725311`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "2"], ",", "4.0044840715635186`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "3"], ",", "1.135564175304713`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "4"], ",", "686639.3094733317`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "5"], ",", "1.9440307014502825`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "6"], ",", "549142.436488662`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "7"], ",", "1.178598948495137`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "8"], ",", "680689.1254867393`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "9"], ",", "370152.93649694487`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "10"], ",", "79019.55848215359`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "11"], ",", "920994.5623687214`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "12"], ",", "1.0226601270489899`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "13"], ",", "1.9668264057100106`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "14"], ",", "1.351941244276504`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "15"], ",", "692029.5451003822`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "16"], ",", "1.1248643751510964`*^6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "17"], ",", "135827.8438653382`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "18"], ",", "678825.3872858201`"}], "}"}]}], 
   "}"}]}]], "Input",
 CellChangeTimes->{{3.9048993332656107`*^9, 
  3.904899365053023*^9}},ExpressionUUID->"42406004-0047-4213-ba5e-\
c5607d601621"],

Cell[BoxData[
 RowBox[{
  RowBox[{"ParsUSc", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "1"], ",", "0.22307323041929009`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "2"], ",", "0.46626237661958264`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "3"], ",", "0.2858617973432077`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "4"], ",", "0.5889348660805614`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "5"], ",", "0.269154465636806`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "6"], ",", "0.06696604631700814`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "7"], ",", "0.48462232049929177`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "8"], ",", "0.5429243341906766`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "9"], ",", "0.5877996581269637`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "10"], ",", "0.15870457585807324`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "11"], ",", "0.290831710931846`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "12"], ",", "0.21220602021402282`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "13"], ",", "0.028658433031011843`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "14"], ",", "0.6065802948734511`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "15"], ",", "0.6616591166626528`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "16"], ",", "0.1255880312581142`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["a", "17"], ",", "0.3655277215373987`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "1"], ",", "5.942158564428258`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "2"], ",", "20.380190456077532`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "3"], ",", "18.574452856249945`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "4"], ",", "30.15739461022833`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "5"], ",", "3.668747115167499`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "6"], ",", "3.310702150300672`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "7"], ",", "22.988822305378054`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "8"], ",", "29.385332300047775`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "9"], ",", "0.3655277215373987`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "10"], ",", "19.96380498296387`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "11"], ",", "9.082681348053972`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "12"], ",", "14.090047144010972`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "13"], ",", "2.457470772507497`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "14"], ",", "2.1105418187161815`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "15"], ",", "18.994858171374325`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "16"], ",", "19.422408277503614`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["c", "17"], ",", "3.298842077560371`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "1"], ",", "1.6296647872806403`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "2"], ",", "1.2552303394426382`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "3"], ",", "886397.2717791495`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "4"], ",", "1.775234431633412`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "5"], ",", "760606.4057041183`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "6"], ",", "8.244728731029832`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "7"], ",", "1.883600316780816`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "8"], ",", "737239.1395470635`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "9"], ",", "3.298842077560371`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "10"], ",", "2.8775094990698076`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "11"], ",", "639997.8081759373`"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "12"], ",", "4.347049924845669`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "13"], ",", "1.5626173188986466`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "14"], ",", "4.371205197098493`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "15"], ",", "8.040848471546578`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "16"], ",", "5.221108990394659`*^6"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SubscriptBox["b", "17"], ",", "5.80930989835473`*^6"}], "}"}]}], 
    "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.877099829935419*^9, 3.8770998365362988`*^9}, {
  3.8776047876883574`*^9, 3.877604789572139*^9}},
 CellLabel->
  "In[138]:=",ExpressionUUID->"d0792500-13e2-4b2d-b73d-d3eb0e0c7e5f"]
}, Open  ]],

Cell[CellGroupData[{

Cell["All", "Subsection",
 CellChangeTimes->{{3.8709441947786517`*^9, 
  3.8709441956603746`*^9}},ExpressionUUID->"8cf1d87c-89f4-4242-a4a4-\
d1ac79da9065"],

Cell[BoxData[
 RowBox[{"l", "=", "16"}]], "Input",
 CellChangeTimes->{{3.870941700057722*^9, 3.870941703626641*^9}, {
  3.870943576800144*^9, 3.870943577707052*^9}, {3.877606072903012*^9, 
  3.8776060749367027`*^9}, {3.8855624544166975`*^9, 3.8855624545104866`*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"fbd3b7f1-0622-47fe-819b-e3174ee7a4af"],

Cell[BoxData[
 RowBox[{
  RowBox[{"mf", "=", 
   RowBox[{"cshn", "[", 
    RowBox[{"z", ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867832129979463*^9, 3.867832158572497*^9}, 
   3.8708918194752913`*^9},
 CellLabel->"In[11]:=",ExpressionUUID->"4931fa6b-24a0-4de9-8358-9929a13a1039"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867834570829482*^9, 3.8678345775014815`*^9}, {
  3.870942389746482*^9, 3.8709424187653246`*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"1018d7f9-8778-4c69-8363-10123b67c0e6"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Cons", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"0.01", "<", 
        SubscriptBox["a", "i"], "<", "0.9"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"1", "<", 
        SubscriptBox["c", "i"], "<", "80"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SuperscriptBox["10", "1"], "<", 
        SubscriptBox["b", "i"], "<", 
        SuperscriptBox["10", "7"]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.870943197443465*^9, 3.870943201625883*^9}},
 CellLabel->"In[18]:=",ExpressionUUID->"acb2c888-3248-4660-90d6-24e25e67ccd9"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cons", "=", 
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.01", "<", 
         SubscriptBox["a", "i"], "<", "2.5"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"1", "<", 
         SubscriptBox["c", "i"], "<", "80"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox["10", "1"], "<", 
         SubscriptBox["b", "i"], "<", 
         SuperscriptBox["10", "7"]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          SubscriptBox["c", "i"], "/", 
          SubscriptBox["a", "i"]}], "<", 
         RowBox[{"n", "+", "2"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], ";"}], 
  " "}]], "Input",
 CellChangeTimes->{{3.8702680472882824`*^9, 3.870268062880148*^9}, {
   3.870274527961333*^9, 3.8702745289124813`*^9}, {3.870274712729985*^9, 
   3.8702747138006673`*^9}, {3.8703380760634556`*^9, 3.870338085777099*^9}, {
   3.870338134356515*^9, 3.870338134781331*^9}, {3.870338715435112*^9, 
   3.870338740776287*^9}, {3.870359410836882*^9, 3.8703594114353127`*^9}, {
   3.870359494778078*^9, 3.870359503140952*^9}, 3.8703645089762597`*^9, {
   3.8704435176581383`*^9, 3.8704435249828515`*^9}, 3.8704465109072294`*^9, {
   3.8704471172830877`*^9, 3.870447117490363*^9}, {3.870448573777196*^9, 
   3.8704485746457186`*^9}, {3.8707868280045576`*^9, 3.8707869329824204`*^9}, 
   3.870892974501359*^9, {3.87089364835896*^9, 3.8708936623592963`*^9}, {
   3.870894444232833*^9, 3.8708944531899633`*^9}, {3.8708945052751646`*^9, 
   3.870894509928249*^9}, 3.8709424914766407`*^9, 3.8710454439888363`*^9, {
   3.877100748287959*^9, 3.8771007483712378`*^9}, {3.877606147034653*^9, 
   3.8776061544831514`*^9}, {3.8776062168979917`*^9, 3.877606217031622*^9}, {
   3.8776062821348906`*^9, 3.8776062931383066`*^9}, {3.8776063744335375`*^9, 
   3.8776063745873184`*^9}, {3.8776064238679185`*^9, 
   3.8776064239865093`*^9}, {3.8776064782443686`*^9, 
   3.8776065050666475`*^9}, {3.8776065400088463`*^9, 3.877606540769662*^9}, {
   3.8776120040949388`*^9, 3.877612008213312*^9}, {3.8776120669677286`*^9, 
   3.8776120680432634`*^9}},
 CellLabel->"In[13]:=",ExpressionUUID->"add0452d-0f32-4abe-9d7c-2a8acf596353"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cons", "=", 
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.01", "<", 
         SubscriptBox["a", "i"], "<", "4.6"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"1", "<", 
         SubscriptBox["c", "i"], "<", "280"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox["10", "1"], "<", 
         SubscriptBox["b", "i"], "<", 
         SuperscriptBox["10", "7"]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          SubscriptBox["c", "i"], "/", 
          SubscriptBox["a", "i"]}], "<", 
         RowBox[{"n", "+", "2"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], ";"}], 
  " "}]], "Input",
 CellChangeTimes->{{3.885562747591427*^9, 3.8855627537177215`*^9}, {
  3.8855628740212135`*^9, 3.885562875021226*^9}, {3.885563093999267*^9, 
  3.885563094108673*^9}, {3.885563179349765*^9, 3.8855631833501177`*^9}, {
  3.8855632999828453`*^9, 3.8855633004985065`*^9}, {3.8855634051452284`*^9, 
  3.8855634064264812`*^9}, {3.8855674821771092`*^9, 3.8855674827862844`*^9}, {
  3.8855675430256763`*^9, 3.8855675437288013`*^9}, {3.885567598983508*^9, 
  3.8855675994836683`*^9}, {3.8855677636516743`*^9, 3.885567764182722*^9}, {
  3.8855694766221066`*^9, 3.885569478153169*^9}},
 CellLabel->"In[71]:=",ExpressionUUID->"f476cebe-9218-4d02-acf6-89cd076a668b"],

Cell[BoxData[
 RowBox[{"Cons", "//", "Short"}]], "Input",
 CellChangeTimes->{{3.870942368113262*^9, 3.87094241536086*^9}, {
  3.8709653784232564`*^9, 3.8709653794744134`*^9}},
 CellLabel->"In[99]:=",ExpressionUUID->"4a9d94fb-fd4b-4a8a-8546-fb8d42ab4015"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"t0", ",", "lsf"}], "}"}], "=", 
   RowBox[{"AbsoluteTiming", "[", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"Data", ",", 
      RowBox[{"{", 
       RowBox[{"mf", ",", "Cons"}], "}"}], ",", "Pars", ",", 
      RowBox[{"{", "z", "}"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "40000"}]}], "]"}], "]"}]}], ";", 
  " ", "t0"}]], "Input",
 CellChangeTimes->{{3.8696649023379245`*^9, 3.8696649531392283`*^9}, {
   3.8697258207695074`*^9, 3.869725821480483*^9}, {3.8697265516411*^9, 
   3.8697265520491714`*^9}, {3.869732294992208*^9, 3.86973229683278*^9}, {
   3.869733165401655*^9, 3.8697331661405754`*^9}, {3.8697373580290446`*^9, 
   3.869737358697692*^9}, {3.869744460132951*^9, 3.869744470596495*^9}, {
   3.8697523910705624`*^9, 3.869752393040555*^9}, {3.869755922118797*^9, 
   3.8697559238834257`*^9}, {3.8697561207431426`*^9, 3.869756121493045*^9}, {
   3.869756160011915*^9, 3.869756168565979*^9}, {3.86975627964835*^9, 
   3.869756281294915*^9}, 3.869838664700287*^9, {3.870265616754628*^9, 
   3.870265616882077*^9}, {3.87026726926554*^9, 3.870267269435379*^9}, {
   3.870892909198313*^9, 3.870892910305993*^9}, 3.8708936939969606`*^9, 
   3.870894335050465*^9, {3.8710451286298876`*^9, 3.871045128809559*^9}},
 CellLabel->"In[72]:=",ExpressionUUID->"d87c187a-9670-4a1c-a83a-fe19bb88633c"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"FR", "=", 
   RowBox[{"lsf", "[", "\"\<FitResiduals\>\"", "]"}]}], ";"}], "\n", 
 RowBox[{"lsf", "[", 
  RowBox[{"\"\<EstimatedVariance\>\"", ",", 
   RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"RMS", "[", 
       RowBox[{"#", ",", "n", ",", 
        RowBox[{"3", "*", "l"}]}], "]"}], "&"}], ")"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.867834878616255*^9, 3.867834878772255*^9}, 
   3.8692243637591076`*^9, {3.869224444478524*^9, 3.8692244539570875`*^9}, {
   3.869229367657488*^9, 3.869229370946748*^9}},
 CellLabel->"In[73]:=",ExpressionUUID->"9c5934b0-179f-45d6-97c9-3fc57e8e5efb"],

Cell[BoxData[
 RowBox[{"lsf", "[", "z", "]"}]], "Input",
 CellChangeTimes->{{3.877614596006403*^9, 3.877614601881892*^9}, {
  3.877614653475318*^9, 3.877614658158794*^9}, {3.877614710801098*^9, 
  3.8776147185943856`*^9}, {3.8776147638477893`*^9, 3.877614774376788*^9}},
 CellLabel->"In[65]:=",ExpressionUUID->"86ba14f5-0ec5-4759-bfcb-1e5b167a4368"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"lsf", "[", "i", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"i", ",", "1", ",", "160"}], "}"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"PlotRange", "->", "Full"}], ",", 
   RowBox[{"Epilog", "->", 
    RowBox[{"Point", "[", "Datab", "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8708070875569153`*^9, 3.8708070877519207`*^9}, {
  3.8708927694670935`*^9, 3.8708927701268067`*^9}, {3.877100788988387*^9, 
  3.8771007899602*^9}, {3.885562630509494*^9, 3.885562631212621*^9}},
 CellLabel->"In[66]:=",ExpressionUUID->"22a9ea35-674f-41a0-91e3-d0ab3d537c4f"],

Cell[BoxData[
 RowBox[{"pars", "=", 
  RowBox[{"lsf", "[", "\"\<BestFitParameters\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.8709447037162905`*^9, 3.870944706183713*^9}, {
  3.870966985723645*^9, 3.8709669903137703`*^9}},
 CellLabel->
  "In[216]:=",ExpressionUUID->"597b583c-03fa-42e3-b48e-cb6d8b30e33d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         SubscriptBox["a", "i"], ",", 
         RowBox[{
          SubscriptBox["a", "i"], "/.", "pars"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         SubscriptBox["c", "i"], ",", 
         RowBox[{
          SubscriptBox["c", "i"], "/.", "pars"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         SubscriptBox["b", "i"], ",", 
         RowBox[{
          SubscriptBox["b", "i"], "/.", "pars"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.8703383358095846`*^9, 3.8703383359032965`*^9}, 
   3.8703385368297215`*^9, 3.8708069708168297`*^9, {3.8709447711598516`*^9, 
   3.870944777123497*^9}, {3.870967026265746*^9, 3.8709670433658247`*^9}, 
   3.87104556011141*^9, 3.877100930672902*^9},
 CellLabel->
  "In[208]:=",ExpressionUUID->"dc637d05-948b-4e9c-902b-2879c62737a3"],

Cell["Go Maple!", "Text",
 CellChangeTimes->{{3.8697246415680885`*^9, 
  3.869724644647354*^9}},ExpressionUUID->"f1614df4-cb52-445d-bbc7-\
8abd06148c4a"],

Cell[BoxData[
 RowBox[{"ParsMaple", "=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       SubscriptBox["a", "i"], "/.", "pars"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       SubscriptBox["c", "i"], "/.", "pars"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       SubscriptBox["b", "i"], "/.", "pars"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.867848620243645*^9, 3.8678486971683884`*^9}, {
   3.8679151951023393`*^9, 3.8679151971583385`*^9}, 3.869226041175602*^9, 
   3.87761538957613*^9},
 CellLabel->
  "In[217]:=",ExpressionUUID->"097d0611-bab5-45cf-953e-507ade9e8796"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Export", "[", 
   RowBox[{
   "\"\<D:\\\\sousa\\\\pesquisa\\\\Covid19\\\\GnuPlot\\\\ParsUSc.txt\>\"", 
    ",", "ParsMaple", ",", "\"\<CSV\>\""}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.867915697525957*^9, 3.8679157280658646`*^9}, 
   3.8679157628765945`*^9, {3.8679164509181013`*^9, 3.867916451654106*^9}, 
   3.867916808995467*^9, 3.867917142053749*^9, {3.867917632604063*^9, 
   3.8679176413030043`*^9}, {3.867929514447388*^9, 3.867929516560387*^9}, {
   3.867930359137655*^9, 3.867930360405655*^9}, {3.869042378475341*^9, 
   3.8690423865541596`*^9}, {3.8690498646245327`*^9, 3.869049868272987*^9}, 
   3.869226047759435*^9, {3.869664545710162*^9, 3.869664546845023*^9}, {
   3.8697257951606436`*^9, 3.8697257954476876`*^9}, {3.869751738310028*^9, 
   3.869751744380788*^9}, 3.8697558833468914`*^9, {3.869836556856716*^9, 
   3.869836556925766*^9}, {3.8703384927006083`*^9, 3.8703384969338017`*^9}, {
   3.870339143253948*^9, 3.8703391440818763`*^9}, {3.87036072023104*^9, 
   3.870360723804949*^9}, 3.870365323681324*^9, {3.870448556460325*^9, 
   3.870448556651752*^9}, {3.870879015485437*^9, 3.8708790192284355`*^9}, {
   3.870893538346167*^9, 3.870893539412013*^9}, 3.87104557858768*^9, {
   3.8771008479686832`*^9, 3.877100850810518*^9}},
 CellLabel->
  "In[210]:=",ExpressionUUID->"55c9bd54-93f6-4d7c-aa86-732147512f0c"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"lsf", "[", "i", "]"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "150"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "150"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"6", "*", 
           SuperscriptBox["10", "6"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", "Red", "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{
        "{", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \
\(147\)]\)(n)\>\"", "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Contamination/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870893042085024*^9, 3.870893071424758*^9}, {
  3.8708932351373615`*^9, 3.870893263626501*^9}, {3.870893298766883*^9, 
  3.8708932988554955`*^9}, {3.8709636554052777`*^9, 3.8709637057698674`*^9}, {
  3.870963784521944*^9, 3.870963790720209*^9}, {3.8709638331309447`*^9, 
  3.8709638335118675`*^9}, {3.8709670491579113`*^9, 3.870967064159026*^9}, {
  3.877612148546461*^9, 3.8776121645978775`*^9}},
 CellLabel->
  "In[204]:=",ExpressionUUID->"78076fd0-a60f-4a22-8db5-5bfac7cb10b4"],

Cell[BoxData[
 RowBox[{"gfc", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf", ",", "gc"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8708930927308626`*^9, 3.8708930987028093`*^9}, {
  3.8709637151085815`*^9, 3.870963719673771*^9}, {3.870963876897517*^9, 
  3.8709638841272526`*^9}, {3.8709670733660917`*^9, 3.8709670763221116`*^9}, {
  3.871045252774228*^9, 3.8710452607917747`*^9}, {3.8776122076375575`*^9, 
  3.8776122156493917`*^9}},
 CellLabel->
  "In[205]:=",ExpressionUUID->"7bbb9ee2-3ea8-479c-b137-646676e9d40b"]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{

Cell["Deaths", "Section",
 CellChangeTimes->{{3.870941623488143*^9, 
  3.870941624781828*^9}},ExpressionUUID->"d335295a-360c-4ea8-8148-\
12a2a7be3e4e"],

Cell[CellGroupData[{

Cell["All", "Subsection",
 CellChangeTimes->{{3.8709441947786517`*^9, 
  3.8709441956603746`*^9}},ExpressionUUID->"cd9101aa-571f-406e-b312-\
a006b1c92f75"],

Cell[BoxData[
 RowBox[{"l", "=", "16"}]], "Input",
 CellChangeTimes->{{3.870941700057722*^9, 3.870941703626641*^9}, {
  3.870943576800144*^9, 3.870943577707052*^9}},
 CellLabel->
  "In[133]:=",ExpressionUUID->"bed05695-006c-4eb6-878c-65a5ce4266ec"],

Cell[BoxData[
 RowBox[{
  RowBox[{"mf", "=", 
   RowBox[{"cshn", "[", 
    RowBox[{"z", ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867832129979463*^9, 3.867832158572497*^9}, 
   3.8708918194752913`*^9},
 CellLabel->
  "In[134]:=",ExpressionUUID->"1879b5b4-edfd-4554-9894-8a66c4454fb6"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867834570829482*^9, 3.8678345775014815`*^9}, {
  3.870942389746482*^9, 3.8709424187653246`*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"63a1013e-93a8-4fb9-8ce3-cb23f701b3b0"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Cons", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"0.01", "<", 
        SubscriptBox["a", "i"], "<", "0.9"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"1", "<", 
        SubscriptBox["c", "i"], "<", "80"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SuperscriptBox["10", "1"], "<", 
        SubscriptBox["b", "i"], "<", 
        SuperscriptBox["10", "7"]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.870943197443465*^9, 3.870943201625883*^9}},
 CellLabel->"In[18]:=",ExpressionUUID->"2e13a3de-eed2-4653-a637-d329f7322a8d"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cons", "=", 
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.01", "<", 
         SubscriptBox["a", "i"], "<", "0.96"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"1", "<", 
         SubscriptBox["c", "i"], "<", "80"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox["10", "1"], "<", 
         SubscriptBox["b", "i"], "<", 
         SuperscriptBox["10", "7"]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          SubscriptBox["c", "i"], "/", 
          SubscriptBox["a", "i"]}], "<", "n"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], ";"}], 
  " "}]], "Input",
 CellChangeTimes->{{3.8702680472882824`*^9, 3.870268062880148*^9}, {
   3.870274527961333*^9, 3.8702745289124813`*^9}, {3.870274712729985*^9, 
   3.8702747138006673`*^9}, {3.8703380760634556`*^9, 3.870338085777099*^9}, {
   3.870338134356515*^9, 3.870338134781331*^9}, {3.870338715435112*^9, 
   3.870338740776287*^9}, {3.870359410836882*^9, 3.8703594114353127`*^9}, {
   3.870359494778078*^9, 3.870359503140952*^9}, 3.8703645089762597`*^9, {
   3.8704435176581383`*^9, 3.8704435249828515`*^9}, 3.8704465109072294`*^9, {
   3.8704471172830877`*^9, 3.870447117490363*^9}, {3.870448573777196*^9, 
   3.8704485746457186`*^9}, {3.8707868280045576`*^9, 3.8707869329824204`*^9}, 
   3.870892974501359*^9, {3.87089364835896*^9, 3.8708936623592963`*^9}, {
   3.870894444232833*^9, 3.8708944531899633`*^9}, {3.8708945052751646`*^9, 
   3.870894509928249*^9}, 3.8709424914766407`*^9, 3.8710454439888363`*^9},
 CellLabel->
  "In[135]:=",ExpressionUUID->"b4d5036b-b7eb-45ad-aa98-48ef609fe014"],

Cell[BoxData[
 RowBox[{"Cons", "//", "Short"}]], "Input",
 CellChangeTimes->{{3.870942368113262*^9, 3.87094241536086*^9}, {
  3.8709653784232564`*^9, 3.8709653794744134`*^9}},
 CellLabel->"In[99]:=",ExpressionUUID->"512de362-116b-49b3-a5a3-9376312c59ce"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"t0", ",", "lsf"}], "}"}], "=", 
   RowBox[{"AbsoluteTiming", "[", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"Data", ",", 
      RowBox[{"{", 
       RowBox[{"mf", ",", "Cons"}], "}"}], ",", "Pars", ",", 
      RowBox[{"{", "z", "}"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "40000"}]}], "]"}], "]"}]}], ";", 
  " ", "t0"}]], "Input",
 CellChangeTimes->{{3.8696649023379245`*^9, 3.8696649531392283`*^9}, {
   3.8697258207695074`*^9, 3.869725821480483*^9}, {3.8697265516411*^9, 
   3.8697265520491714`*^9}, {3.869732294992208*^9, 3.86973229683278*^9}, {
   3.869733165401655*^9, 3.8697331661405754`*^9}, {3.8697373580290446`*^9, 
   3.869737358697692*^9}, {3.869744460132951*^9, 3.869744470596495*^9}, {
   3.8697523910705624`*^9, 3.869752393040555*^9}, {3.869755922118797*^9, 
   3.8697559238834257`*^9}, {3.8697561207431426`*^9, 3.869756121493045*^9}, {
   3.869756160011915*^9, 3.869756168565979*^9}, {3.86975627964835*^9, 
   3.869756281294915*^9}, 3.869838664700287*^9, {3.870265616754628*^9, 
   3.870265616882077*^9}, {3.87026726926554*^9, 3.870267269435379*^9}, {
   3.870892909198313*^9, 3.870892910305993*^9}, 3.8708936939969606`*^9, 
   3.870894335050465*^9, {3.8710451286298876`*^9, 3.871045128809559*^9}},
 CellLabel->
  "In[136]:=",ExpressionUUID->"e2c58dbd-4bd7-45a8-bcda-e9bcc73e5b5e"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"FR", "=", 
   RowBox[{"lsf", "[", "\"\<FitResiduals\>\"", "]"}]}], ";"}], "\n", 
 RowBox[{"lsf", "[", 
  RowBox[{"\"\<EstimatedVariance\>\"", ",", 
   RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"RMS", "[", 
       RowBox[{"#", ",", "n", ",", 
        RowBox[{"3", "*", "l"}]}], "]"}], "&"}], ")"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.867834878616255*^9, 3.867834878772255*^9}, 
   3.8692243637591076`*^9, {3.869224444478524*^9, 3.8692244539570875`*^9}, {
   3.869229367657488*^9, 3.869229370946748*^9}},
 CellLabel->
  "In[137]:=",ExpressionUUID->"c867c94f-59e0-4ce0-aed0-5a3a04d1156c"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"lsf", "[", "i", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"i", ",", "1", ",", 
     RowBox[{"n", "+", "10"}]}], "}"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"PlotRange", "->", "Full"}], ",", 
   RowBox[{"Epilog", "->", 
    RowBox[{"Point", "[", "Datab", "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8708070875569153`*^9, 3.8708070877519207`*^9}, {
  3.8708927694670935`*^9, 3.8708927701268067`*^9}},
 CellLabel->
  "In[139]:=",ExpressionUUID->"de9ca565-4254-45cb-9c47-ad85f95c5d08"],

Cell[BoxData[
 RowBox[{"pars", "=", 
  RowBox[{"lsf", "[", "\"\<BestFitParameters\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.8709447037162905`*^9, 3.870944706183713*^9}, {
  3.870966985723645*^9, 3.8709669903137703`*^9}},
 CellLabel->
  "In[140]:=",ExpressionUUID->"82be80d3-9d7c-4d4d-b14a-843b3fbc4133"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         SubscriptBox["a", "i"], ",", 
         RowBox[{
          SubscriptBox["a", "i"], "/.", "pars"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         SubscriptBox["c", "i"], ",", 
         RowBox[{
          SubscriptBox["c", "i"], "/.", "pars"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         SubscriptBox["b", "i"], ",", 
         RowBox[{
          SubscriptBox["b", "i"], "/.", "pars"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.8703383358095846`*^9, 3.8703383359032965`*^9}, 
   3.8703385368297215`*^9, 3.8708069708168297`*^9, {3.8709447711598516`*^9, 
   3.870944777123497*^9}, {3.870967026265746*^9, 3.8709670433658247`*^9}, 
   3.87104556011141*^9},
 CellLabel->
  "In[141]:=",ExpressionUUID->"e7acdc00-60c1-43ac-a816-4a2138f6c147"],

Cell["Go Maple!", "Text",
 CellChangeTimes->{{3.8697246415680885`*^9, 
  3.869724644647354*^9}},ExpressionUUID->"9fccfbc0-db4c-4aae-aaae-\
9360d4e2cefd"],

Cell[BoxData[
 RowBox[{
  RowBox[{"ParsMaple", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["a", "i"], "/.", "pars"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["c", "i"], "/.", "pars"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SubscriptBox["b", "i"], "/.", "pars"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867848620243645*^9, 3.8678486971683884`*^9}, {
   3.8679151951023393`*^9, 3.8679151971583385`*^9}, 3.869226041175602*^9},
 CellLabel->
  "In[142]:=",ExpressionUUID->"5ad82bff-0fbd-4410-b058-72e8687a6875"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Export", "[", 
   RowBox[{
   "\"\<D:\\\\sousa\\\\pesquisa\\\\Covid19\\\\GnuPlot\\\\parsUSall3.txt\>\"", 
    ",", "ParsMaple", ",", "\"\<CSV\>\""}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.867915697525957*^9, 3.8679157280658646`*^9}, 
   3.8679157628765945`*^9, {3.8679164509181013`*^9, 3.867916451654106*^9}, 
   3.867916808995467*^9, 3.867917142053749*^9, {3.867917632604063*^9, 
   3.8679176413030043`*^9}, {3.867929514447388*^9, 3.867929516560387*^9}, {
   3.867930359137655*^9, 3.867930360405655*^9}, {3.869042378475341*^9, 
   3.8690423865541596`*^9}, {3.8690498646245327`*^9, 3.869049868272987*^9}, 
   3.869226047759435*^9, {3.869664545710162*^9, 3.869664546845023*^9}, {
   3.8697257951606436`*^9, 3.8697257954476876`*^9}, {3.869751738310028*^9, 
   3.869751744380788*^9}, 3.8697558833468914`*^9, {3.869836556856716*^9, 
   3.869836556925766*^9}, {3.8703384927006083`*^9, 3.8703384969338017`*^9}, {
   3.870339143253948*^9, 3.8703391440818763`*^9}, {3.87036072023104*^9, 
   3.870360723804949*^9}, 3.870365323681324*^9, {3.870448556460325*^9, 
   3.870448556651752*^9}, {3.870879015485437*^9, 3.8708790192284355`*^9}, {
   3.870893538346167*^9, 3.870893539412013*^9}, 3.87104557858768*^9, {
   3.8732958250200577`*^9, 3.8732958252666883`*^9}},
 CellLabel->
  "In[143]:=",ExpressionUUID->"65396194-ba47-41a1-bfd2-30afc764a0dd"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"lsf", "[", "i", "]"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "140"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "140"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"24", "*", 
           SuperscriptBox["10", "3"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", "Red", "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{
        "{", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \
\(140\)]\)(n)\>\"", "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Deaths/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870893042085024*^9, 3.870893071424758*^9}, {
  3.8708932351373615`*^9, 3.870893263626501*^9}, {3.870893298766883*^9, 
  3.8708932988554955`*^9}, {3.8709636554052777`*^9, 3.8709637057698674`*^9}, {
  3.870963784521944*^9, 3.870963790720209*^9}, {3.8709638331309447`*^9, 
  3.8709638335118675`*^9}, {3.8709670491579113`*^9, 3.870967064159026*^9}, {
  3.8710567021854973`*^9, 3.871056709206452*^9}, {3.8732958527142677`*^9, 
  3.8732958663303823`*^9}},
 CellLabel->
  "In[144]:=",ExpressionUUID->"ba9550fc-392b-41cd-9870-c14997dbafed"],

Cell[BoxData[
 RowBox[{"gfd", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf", ",", "gd"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8708930927308626`*^9, 3.8708930987028093`*^9}, {
  3.8709637151085815`*^9, 3.870963719673771*^9}, {3.870963876897517*^9, 
  3.8709638841272526`*^9}, {3.8709670733660917`*^9, 3.8709670763221116`*^9}, {
  3.871045252774228*^9, 3.8710452607917747`*^9}},
 CellLabel->
  "In[145]:=",ExpressionUUID->"4b9c3dd4-19c8-4b6c-9f60-581ef19a4b93"]
}, Closed]],

Cell[CellGroupData[{

Cell["First part", "Subsection",
 CellChangeTimes->{{3.8709441005548306`*^9, 
  3.870944103121689*^9}},ExpressionUUID->"66afd28b-f580-40f2-ab1d-\
a1e947869f46"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data1", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"DataAll", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "i", ",", "6"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "75"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.867833607979697*^9, 3.8678336417547903`*^9}, 
   3.867834042482584*^9, 3.867834167808469*^9, {3.869228379395009*^9, 
   3.8692283841724577`*^9}, {3.8692285327297277`*^9, 
   3.8692285335105352`*^9}, {3.8696639875185523`*^9, 3.869663991411394*^9}, {
   3.8707908150861864`*^9, 3.870790820243222*^9}, {3.870891869329871*^9, 
   3.870891887445064*^9}, {3.8708919270279016`*^9, 3.8708919310272264`*^9}, {
   3.87094213776725*^9, 3.8709421385938425`*^9}, {3.8709444109197607`*^9, 
   3.8709444181136146`*^9}},
 CellLabel->"In[5]:=",ExpressionUUID->"91f8ecca-2e58-44b4-b65c-f9330de3ae66"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data1b", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"DataAll", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "i", ",", "6"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "75"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870891949547256*^9, 3.870891955747621*^9}, {
  3.870942154736216*^9, 3.870942155364955*^9}, {3.870944422170905*^9, 
  3.870944429072955*^9}},
 CellLabel->"In[6]:=",ExpressionUUID->"0034a340-1fa2-49cf-94a0-5ee247f6966a"],

Cell[BoxData[
 RowBox[{"l1", "=", "8"}]], "Input",
 CellChangeTimes->{{3.870941700057722*^9, 3.870941703626641*^9}, {
  3.870943576800144*^9, 3.870943577707052*^9}, {3.8709442099463186`*^9, 
  3.870944220344717*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"10542a1f-a53e-43f1-b326-90c5553b8c2a"],

Cell[BoxData[
 RowBox[{
  RowBox[{"mf1", "=", 
   RowBox[{"cshn", "[", 
    RowBox[{"z", ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867832129979463*^9, 3.867832158572497*^9}, 
   3.8708918194752913`*^9, {3.8709442342350492`*^9, 3.8709442476434917`*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"a8fbbeb5-886e-45cd-936b-c45a66e7aa9e"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars1", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867834570829482*^9, 3.8678345775014815`*^9}, {
  3.870942389746482*^9, 3.8709424187653246`*^9}, {3.8709442605776134`*^9, 
  3.870944290983293*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"0280e91f-a6d9-406f-ab00-3b1749ed26b5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Cons1", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"0.01", "<", 
        SubscriptBox["a", "i"], "<", "0.9"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"1", "<", 
        SubscriptBox["c", "i"], "<", "80"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SuperscriptBox["10", "1"], "<", 
        SubscriptBox["b", "i"], "<", 
        SuperscriptBox["10", "7"]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.870943197443465*^9, 3.870943201625883*^9}, {
  3.870944283544839*^9, 3.8709443201375895`*^9}},
 CellLabel->"In[40]:=",ExpressionUUID->"1bc98aba-dac9-4ea4-8490-749971a7b616"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cons1", "=", 
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.01", "<", 
         SubscriptBox["a", "i"], "<", "0.96"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"1", "<", 
         SubscriptBox["c", "i"], "<", "90"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox["10", "1"], "<", 
         SubscriptBox["b", "i"], "<", 
         SuperscriptBox["10", "6"]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{
          SubscriptBox["c", "i"], "/", 
          SubscriptBox["a", "i"]}], "<", "n"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}], ";"}],
   " "}]], "Input",
 CellChangeTimes->{{3.8702680472882824`*^9, 3.870268062880148*^9}, {
   3.870274527961333*^9, 3.8702745289124813`*^9}, {3.870274712729985*^9, 
   3.8702747138006673`*^9}, {3.8703380760634556`*^9, 3.870338085777099*^9}, {
   3.870338134356515*^9, 3.870338134781331*^9}, {3.870338715435112*^9, 
   3.870338740776287*^9}, {3.870359410836882*^9, 3.8703594114353127`*^9}, {
   3.870359494778078*^9, 3.870359503140952*^9}, 3.8703645089762597`*^9, {
   3.8704435176581383`*^9, 3.8704435249828515`*^9}, 3.8704465109072294`*^9, {
   3.8704471172830877`*^9, 3.870447117490363*^9}, {3.870448573777196*^9, 
   3.8704485746457186`*^9}, {3.8707868280045576`*^9, 3.8707869329824204`*^9}, 
   3.870892974501359*^9, {3.87089364835896*^9, 3.8708936623592963`*^9}, {
   3.870894444232833*^9, 3.8708944531899633`*^9}, {3.8708945052751646`*^9, 
   3.870894509928249*^9}, 3.8709424914766407`*^9, {3.8709443307752733`*^9, 
   3.870944339866149*^9}, {3.871056838055213*^9, 3.8710568522492886`*^9}, {
   3.871056909373443*^9, 3.8710569098387823`*^9}, {3.8710576380394773`*^9, 
   3.871057642402381*^9}, {3.871057806961446*^9, 3.8710578087198954`*^9}, {
   3.871057864428071*^9, 3.8710578653910975`*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"cc34ef7c-6dbe-40fa-b47a-f9f441dcb0d3"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"t0", ",", "lsf1"}], "}"}], "=", 
   RowBox[{"AbsoluteTiming", "[", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"Data1", ",", 
      RowBox[{"{", 
       RowBox[{"mf1", ",", "Cons1"}], "}"}], ",", "Pars1", ",", 
      RowBox[{"{", "z", "}"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "40000"}]}], "]"}], "]"}]}], ";", 
  " ", "t0"}]], "Input",
 CellChangeTimes->{{3.8696649023379245`*^9, 3.8696649531392283`*^9}, {
   3.8697258207695074`*^9, 3.869725821480483*^9}, {3.8697265516411*^9, 
   3.8697265520491714`*^9}, {3.869732294992208*^9, 3.86973229683278*^9}, {
   3.869733165401655*^9, 3.8697331661405754`*^9}, {3.8697373580290446`*^9, 
   3.869737358697692*^9}, {3.869744460132951*^9, 3.869744470596495*^9}, {
   3.8697523910705624`*^9, 3.869752393040555*^9}, {3.869755922118797*^9, 
   3.8697559238834257`*^9}, {3.8697561207431426`*^9, 3.869756121493045*^9}, {
   3.869756160011915*^9, 3.869756168565979*^9}, {3.86975627964835*^9, 
   3.869756281294915*^9}, 3.869838664700287*^9, {3.870265616754628*^9, 
   3.870265616882077*^9}, {3.87026726926554*^9, 3.870267269435379*^9}, {
   3.870892909198313*^9, 3.870892910305993*^9}, 3.8708936939969606`*^9, 
   3.870894335050465*^9, {3.8709443474842157`*^9, 3.870944356400553*^9}, 
   3.870944445627025*^9, {3.8710567765875664`*^9, 3.8710567767085204`*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"ec5a62f9-cb07-4b16-b983-142ac3963670"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"FR", "=", 
   RowBox[{"lsf1", "[", "\"\<FitResiduals\>\"", "]"}]}], ";"}], "\n", 
 RowBox[{"lsf1", "[", 
  RowBox[{"\"\<EstimatedVariance\>\"", ",", 
   RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"RMS", "[", 
       RowBox[{"#", ",", "75", ",", 
        RowBox[{"3", "*", "l1"}]}], "]"}], "&"}], ")"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.867834878616255*^9, 3.867834878772255*^9}, 
   3.8692243637591076`*^9, {3.869224444478524*^9, 3.8692244539570875`*^9}, {
   3.869229367657488*^9, 3.869229370946748*^9}, {3.870944594624638*^9, 
   3.8709446032108283`*^9}, {3.8732927012769575`*^9, 3.8732927091266403`*^9}, 
   3.873292894230918*^9},
 CellLabel->"In[17]:=",ExpressionUUID->"9ea867a0-6970-4605-a870-2b09a8a9bc87"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"lsf1", "[", "i", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"i", ",", "1", ",", "80"}], "}"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"PlotRange", "->", "Full"}], ",", 
   RowBox[{"Epilog", "->", 
    RowBox[{"Point", "[", "Data1b", "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8708070875569153`*^9, 3.8708070877519207`*^9}, {
  3.8708927694670935`*^9, 3.8708927701268067`*^9}, {3.870944630431551*^9, 
  3.8709446429311743`*^9}},
 CellLabel->"In[19]:=",ExpressionUUID->"e2c5f4bc-43a9-4131-9092-e9138138ce44"],

Cell[BoxData[
 RowBox[{"pars1", "=", 
  RowBox[{"lsf1", "[", "\"\<BestFitParameters\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.8709447037162905`*^9, 3.870944706183713*^9}},
 CellLabel->"In[16]:=",ExpressionUUID->"3ac0b1c9-30ff-4706-8b2f-0daa32fc1b7e"],

Cell[BoxData[
 RowBox[{"Pars1", "=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["a", "i"], ",", 
        RowBox[{
         SubscriptBox["a", "i"], "/.", "pars1"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["c", "i"], ",", 
        RowBox[{
         SubscriptBox["c", "i"], "/.", "pars1"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["b", "i"], ",", 
        RowBox[{
         SubscriptBox["b", "i"], "/.", "pars1"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l1"}], "}"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8703383358095846`*^9, 3.8703383359032965`*^9}, 
   3.8703385368297215`*^9, 3.8708069708168297`*^9, {3.8709447711598516`*^9, 
   3.870944777123497*^9}, {3.8710579208737564`*^9, 3.871057930010552*^9}},
 CellLabel->
  "In[122]:=",ExpressionUUID->"db259c83-f156-4b2b-aca3-de1d01a49e95"],

Cell[BoxData[
 RowBox[{"Pars1", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "1"], ",", "0.1381608070121996`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "2"], ",", "0.14638418501069994`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "3"], ",", "0.5279302641157294`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "4"], ",", "0.9157789286701654`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "5"], ",", "0.2173116437906986`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "6"], ",", "0.29520749024480714`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "7"], ",", "0.5834548838763242`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["a", "8"], ",", "0.08556712721338508`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "1"], ",", "4.016173455884751`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "2"], ",", "6.657301100425359`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "3"], ",", "24.93065167827863`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "4"], ",", "47.00255809930365`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "5"], ",", "11.681824205982881`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "6"], ",", "4.65227584335838`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "7"], ",", "7.29846129780184`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["c", "8"], ",", "5.436268385557547`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "1"], ",", "48876.959694935125`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "2"], ",", "49993.15567563796`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "3"], ",", "12649.417137695442`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "4"], ",", "4748.186990122024`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "5"], ",", "77120.98352739475`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "6"], ",", "36936.084869910555`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "7"], ",", "18751.91815379085`"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      SubscriptBox["b", "8"], ",", "58988.11396158941`"}], "}"}]}], 
   "}"}]}]], "Input",
 CellChangeTimes->{{3.8732929644874673`*^9, 3.8732929705500402`*^9}},
 CellLabel->
  "In[121]:=",ExpressionUUID->"8cd3b1de-68b1-4652-a162-6b090fc37590"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf1", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"lsf1", "[", "i", "]"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "140"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "140"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"24", "*", 
           SuperscriptBox["10", "3"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", "Red", "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{
        "{", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \(a\)]\)(n)\>\"", 
         "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Deaths/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870893042085024*^9, 3.870893071424758*^9}, {
  3.8708932351373615`*^9, 3.870893263626501*^9}, {3.870893298766883*^9, 
  3.8708932988554955`*^9}, {3.8709636554052777`*^9, 3.8709637057698674`*^9}, {
  3.870963784521944*^9, 3.870963790720209*^9}, {3.8709638331309447`*^9, 
  3.8709638335118675`*^9}, {3.8710577516420135`*^9, 3.871057767453901*^9}},
 CellLabel->"In[21]:=",ExpressionUUID->"9970887a-0428-45a5-8181-dd5e56a8e31a"],

Cell[BoxData[
 RowBox[{"gfd1", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf1", ",", "gd"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8708930927308626`*^9, 3.8708930987028093`*^9}, {
  3.8709637151085815`*^9, 3.870963719673771*^9}, {3.870963876897517*^9, 
  3.8709638841272526`*^9}, {3.8710577808461533`*^9, 3.8710577874358983`*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"a17d4ce1-fe07-4bad-a2e9-3ea827c395ef"]
}, Closed]],

Cell[CellGroupData[{

Cell["Second part", "Subsection",
 CellChangeTimes->{{3.8709441005548306`*^9, 3.870944103121689*^9}, {
  3.8709639383725348`*^9, 
  3.87096394221332*^9}},ExpressionUUID->"96d408fb-f1c8-434a-90e1-\
fdea94cb786e"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data2", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"DataAll", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "i", ",", "6"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "75", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.867833607979697*^9, 3.8678336417547903`*^9}, 
   3.867834042482584*^9, 3.867834167808469*^9, {3.869228379395009*^9, 
   3.8692283841724577`*^9}, {3.8692285327297277`*^9, 
   3.8692285335105352`*^9}, {3.8696639875185523`*^9, 3.869663991411394*^9}, {
   3.8707908150861864`*^9, 3.870790820243222*^9}, {3.870891869329871*^9, 
   3.870891887445064*^9}, {3.8708919270279016`*^9, 3.8708919310272264`*^9}, {
   3.87094213776725*^9, 3.8709421385938425`*^9}, {3.8709444109197607`*^9, 
   3.8709444181136146`*^9}, {3.8709639486815453`*^9, 3.870963961045286*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"4b4c2ec0-15a2-45f9-88a7-f00cad117989"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Data2b", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"DataAll", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "i", ",", "6"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "75", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870891949547256*^9, 3.870891955747621*^9}, {
  3.870942154736216*^9, 3.870942155364955*^9}, {3.870944422170905*^9, 
  3.870944429072955*^9}, {3.870963971662084*^9, 3.870963980655541*^9}},
 CellLabel->"In[27]:=",ExpressionUUID->"fd4630b0-589d-48fc-8e5f-201a22f34466"],

Cell[BoxData[
 RowBox[{"l2", "=", "8"}]], "Input",
 CellChangeTimes->{{3.870941700057722*^9, 3.870941703626641*^9}, {
  3.870943576800144*^9, 3.870943577707052*^9}, {3.8709442099463186`*^9, 
  3.870944220344717*^9}, {3.87096398743102*^9, 3.870963987521147*^9}},
 CellLabel->"In[28]:=",ExpressionUUID->"8e6d0dfd-6454-47a3-9b72-2eece97625b7"],

Cell[BoxData[
 RowBox[{
  RowBox[{"mf2", "=", 
   RowBox[{"cshn", "[", 
    RowBox[{"z", ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867832129979463*^9, 3.867832158572497*^9}, 
   3.8708918194752913`*^9, {3.8709442342350492`*^9, 3.8709442476434917`*^9}, {
   3.8709639958539886`*^9, 3.870964006063445*^9}, {3.8709640539212065`*^9, 
   3.8709641051602545`*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"42480e1c-df2c-4c77-a303-7306d3526e34"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Pars2", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["a", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["c", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       SubscriptBox["b", "i"], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.867834570829482*^9, 3.8678345775014815`*^9}, {
  3.870942389746482*^9, 3.8709424187653246`*^9}, {3.8709442605776134`*^9, 
  3.870944290983293*^9}, {3.870964012793183*^9, 3.8709640234502068`*^9}, {
  3.8709641225539055`*^9, 3.870964129481372*^9}, {3.870964170592956*^9, 
  3.8709641774432645`*^9}},
 CellLabel->"In[43]:=",ExpressionUUID->"c8d0907d-8e73-4d52-b6a2-e861bf04cc74"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Cons2", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"0.01", "<", 
        SubscriptBox["a", "i"], "<", "0.9"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"1", "<", 
        SubscriptBox["c", "i"], "<", "80"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        SuperscriptBox["10", "1"], "<", 
        SubscriptBox["b", "i"], "<", 
        SuperscriptBox["10", "7"]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"l1", "+", "1"}], ",", 
         RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.870943197443465*^9, 3.870943201625883*^9}, {
  3.870944283544839*^9, 3.8709443201375895`*^9}, {3.8709640311805105`*^9, 
  3.8709640348758936`*^9}, {3.870964133634597*^9, 3.8709641425285974`*^9}, {
  3.870964184136936*^9, 3.87096419317742*^9}},
 CellLabel->"In[74]:=",ExpressionUUID->"3e77d8a9-d9a3-4ff7-bc6a-f522b8824ea1"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cons2", "=", 
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.01", "<", 
         SubscriptBox["a", "i"], "<", "0.98"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"0.1", "<", 
         SubscriptBox["c", "i"], "<", "90"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         SuperscriptBox["10", "1"], "<", 
         SubscriptBox["b", "i"], "<", 
         SuperscriptBox["10", "6"]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"5", "<", 
         RowBox[{
          SubscriptBox["c", "i"], "/", 
          SubscriptBox["a", "i"]}], "<", "n"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", 
          RowBox[{"l1", "+", "1"}], ",", 
          RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}], ";"}], 
  " "}]], "Input",
 CellChangeTimes->{{3.8702680472882824`*^9, 3.870268062880148*^9}, {
   3.870274527961333*^9, 3.8702745289124813`*^9}, {3.870274712729985*^9, 
   3.8702747138006673`*^9}, {3.8703380760634556`*^9, 3.870338085777099*^9}, {
   3.870338134356515*^9, 3.870338134781331*^9}, {3.870338715435112*^9, 
   3.870338740776287*^9}, {3.870359410836882*^9, 3.8703594114353127`*^9}, {
   3.870359494778078*^9, 3.870359503140952*^9}, 3.8703645089762597`*^9, {
   3.8704435176581383`*^9, 3.8704435249828515`*^9}, 3.8704465109072294`*^9, {
   3.8704471172830877`*^9, 3.870447117490363*^9}, {3.870448573777196*^9, 
   3.8704485746457186`*^9}, {3.8707868280045576`*^9, 3.8707869329824204`*^9}, 
   3.870892974501359*^9, {3.87089364835896*^9, 3.8708936623592963`*^9}, {
   3.870894444232833*^9, 3.8708944531899633`*^9}, {3.8708945052751646`*^9, 
   3.870894509928249*^9}, 3.8709424914766407`*^9, {3.8709443307752733`*^9, 
   3.870944339866149*^9}, {3.8709641474240932`*^9, 3.8709641588991756`*^9}, {
   3.8709641994977407`*^9, 3.8709642292079687`*^9}, {3.871057997425127*^9, 
   3.871058010808834*^9}, {3.871058136627864*^9, 3.8710581374740925`*^9}, {
   3.8710581983406124`*^9, 3.871058207632083*^9}, {3.871058250276015*^9, 
   3.871058281720549*^9}, {3.8710583333025413`*^9, 3.8710583397348056`*^9}, {
   3.8732934766766386`*^9, 3.8732934847178607`*^9}, {3.8732935860422187`*^9, 
   3.8732935869682226`*^9}, {3.8732936571064262`*^9, 3.873293662371318*^9}, {
   3.873293740177701*^9, 3.8732937558818274`*^9}, {3.8732938163835344`*^9, 
   3.873293832654558*^9}, 3.8732938918912926`*^9, {3.8732939687956233`*^9, 
   3.8732939749546776`*^9}, {3.8732940343213444`*^9, 
   3.8732940521766996`*^9}, {3.8732940823778353`*^9, 3.87329413179287*^9}, {
   3.8732942159247503`*^9, 3.8732942203154693`*^9}, {3.8732942887336044`*^9, 
   3.8732942893794236`*^9}},
 CellLabel->
  "In[109]:=",ExpressionUUID->"9e7cfe2d-4a5c-4b8d-9589-ad0a92c450fb"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"t0", ",", "lsf2"}], "}"}], "=", 
   RowBox[{"AbsoluteTiming", "[", 
    RowBox[{"NonlinearModelFit", "[", 
     RowBox[{"Data2", ",", 
      RowBox[{"{", 
       RowBox[{"mf2", ",", "Cons2"}], "}"}], ",", "Pars2", ",", 
      RowBox[{"{", "z", "}"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "30000"}]}], "]"}], "]"}]}], ";", 
  " ", "t0"}]], "Input",
 CellChangeTimes->{{3.8696649023379245`*^9, 3.8696649531392283`*^9}, {
   3.8697258207695074`*^9, 3.869725821480483*^9}, {3.8697265516411*^9, 
   3.8697265520491714`*^9}, {3.869732294992208*^9, 3.86973229683278*^9}, {
   3.869733165401655*^9, 3.8697331661405754`*^9}, {3.8697373580290446`*^9, 
   3.869737358697692*^9}, {3.869744460132951*^9, 3.869744470596495*^9}, {
   3.8697523910705624`*^9, 3.869752393040555*^9}, {3.869755922118797*^9, 
   3.8697559238834257`*^9}, {3.8697561207431426`*^9, 3.869756121493045*^9}, {
   3.869756160011915*^9, 3.869756168565979*^9}, {3.86975627964835*^9, 
   3.869756281294915*^9}, 3.869838664700287*^9, {3.870265616754628*^9, 
   3.870265616882077*^9}, {3.87026726926554*^9, 3.870267269435379*^9}, {
   3.870892909198313*^9, 3.870892910305993*^9}, 3.8708936939969606`*^9, 
   3.870894335050465*^9, {3.8709443474842157`*^9, 3.870944356400553*^9}, 
   3.870944445627025*^9, {3.8709642380156755`*^9, 3.870964251855894*^9}},
 CellLabel->
  "In[110]:=",ExpressionUUID->"7a0df496-b5a8-4638-8bff-f764498bb1f9"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"FR", "=", 
   RowBox[{"lsf2", "[", "\"\<FitResiduals\>\"", "]"}]}], ";"}], "\n", 
 RowBox[{"lsf2", "[", 
  RowBox[{"\"\<EstimatedVariance\>\"", ",", 
   RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"RMS", "[", 
       RowBox[{"#", ",", 
        RowBox[{"n", "-", "74"}], ",", 
        RowBox[{"3", "*", "l2"}]}], "]"}], "&"}], ")"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.867834878616255*^9, 3.867834878772255*^9}, 
   3.8692243637591076`*^9, {3.869224444478524*^9, 3.8692244539570875`*^9}, {
   3.869229367657488*^9, 3.869229370946748*^9}, {3.870944594624638*^9, 
   3.8709446032108283`*^9}, {3.8709642777743073`*^9, 3.870964292671986*^9}, {
   3.8732931134856486`*^9, 3.873293133160251*^9}},
 CellLabel->
  "In[111]:=",ExpressionUUID->"7d833c5d-f61f-4775-9b56-d89bc598521b"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"lsf2", "[", 
    RowBox[{"i", "-", "74"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"i", ",", "69", ",", "140"}], "}"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", "Thick"}], "}"}]}], ",", 
   RowBox[{"PlotRange", "->", "Full"}], ",", 
   RowBox[{"Epilog", "->", 
    RowBox[{"Point", "[", "Data2b", "]"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8708070875569153`*^9, 3.8708070877519207`*^9}, {
  3.8708927694670935`*^9, 3.8708927701268067`*^9}, {3.870944630431551*^9, 
  3.8709446429311743`*^9}, {3.870964310082224*^9, 3.8709643691971035`*^9}},
 CellLabel->
  "In[113]:=",ExpressionUUID->"998e1c25-cff4-4a1f-99e3-4afdc0f5ee1d"],

Cell[BoxData[
 RowBox[{"pars2", "=", 
  RowBox[{"lsf2", "[", "\"\<BestFitParameters\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.8709447037162905`*^9, 3.870944706183713*^9}, {
  3.870964389688529*^9, 3.8709643930406733`*^9}},
 CellLabel->
  "In[114]:=",ExpressionUUID->"7fe546e9-3de1-493c-8231-218e620333f3"],

Cell[BoxData[
 RowBox[{"Pars2", "=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["a", "i"], ",", 
        RowBox[{
         SubscriptBox["a", "i"], "/.", "pars2"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"l1", "+", "1"}], ",", 
        RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["c", "i"], ",", 
        RowBox[{
         SubscriptBox["c", "i"], "/.", "pars2"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"l1", "+", "1"}], ",", 
        RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["b", "i"], ",", 
        RowBox[{
         SubscriptBox["b", "i"], "/.", "pars2"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"l1", "+", "1"}], ",", 
        RowBox[{"l1", "+", "l2"}]}], "}"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8703383358095846`*^9, 3.8703383359032965`*^9}, 
   3.8703385368297215`*^9, 3.8708069708168297`*^9, {3.8709447711598516`*^9, 
   3.870944777123497*^9}, {3.870964415028035*^9, 3.870964484931066*^9}},
 CellLabel->
  "In[115]:=",ExpressionUUID->"64cf2943-8185-47e1-b396-99cff3998851"],

Cell[BoxData[
 RowBox[{"Pars2", "="}]], "Input",
 CellChangeTimes->{{3.8732946889226494`*^9, 
  3.8732946933004026`*^9}},ExpressionUUID->"35b9f060-be60-4061-89db-\
dc30a01ab22e"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf2", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"lsf2", "[", 
       RowBox[{"i", "-", "74"}], "]"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "69", ",", "140"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "140"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"24", "*", 
           SuperscriptBox["10", "3"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", "Red", "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{
        "{", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \(b\)]\)(n)\>\"", 
         "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Deaths/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870893042085024*^9, 3.870893071424758*^9}, {
  3.8708932351373615`*^9, 3.870893263626501*^9}, {3.870893298766883*^9, 
  3.8708932988554955`*^9}, {3.8709636554052777`*^9, 3.8709637057698674`*^9}, {
  3.870963784521944*^9, 3.870963790720209*^9}, {3.8709638331309447`*^9, 
  3.8709638335118675`*^9}, {3.8709645215678253`*^9, 3.870964549699584*^9}, {
  3.8732933413224607`*^9, 3.873293352454122*^9}},
 CellLabel->
  "In[116]:=",ExpressionUUID->"548568e7-4e57-4b8e-8a5d-f3fcb236bb65"],

Cell[BoxData[
 RowBox[{"gfd2", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf2", ",", "gd"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8708930927308626`*^9, 3.8708930987028093`*^9}, {
  3.8709637151085815`*^9, 3.870963719673771*^9}, {3.870963876897517*^9, 
  3.8709638841272526`*^9}, {3.870964560150072*^9, 3.870964564545*^9}, {
  3.873293220343669*^9, 3.8732932314936385`*^9}},
 CellLabel->
  "In[117]:=",ExpressionUUID->"f62348a3-aa7a-4568-acf5-a6f840dc5307"],

Cell[BoxData[
 RowBox[{
  RowBox[{"gf12", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"lsf1", "[", "i", "]"}], ",", 
       RowBox[{"lsf2", "[", 
        RowBox[{"i", "-", "74"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "140"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "140"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0", ",", 
          RowBox[{"24", "*", 
           SuperscriptBox["10", "3"]}]}], "}"}]}], "}"}]}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Red", ",", "Brown"}], "}"}]}], ",", 
     RowBox[{"GridLines", "->", "Automatic"}], ",", 
     RowBox[{"Frame", "->", "True"}], ",", 
     RowBox[{"ImageSize", "->", "800"}], ",", 
     RowBox[{"ImagePadding", "->", "60"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"Placed", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
         "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \(a\)]\)(n)\>\"", 
          ",", "\"\<\!\(\*SubscriptBox[\(\[ScriptCapitalZ]\), \
\(b\)]\)(n)\>\""}], "}"}], ",", "Below"}], "]"}]}], ",", 
     RowBox[{"FrameLabel", "->", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"\"\<Deaths/week\>\"", ",", "None"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"\"\<n (week)\>\"", ",", "None"}], "}"}]}], "}"}]}]}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.870964711371406*^9, 3.870964733600583*^9}, {
  3.870964783638962*^9, 3.870964880301608*^9}, {3.8709649138296537`*^9, 
  3.87096493651463*^9}, {3.870965147951193*^9, 3.8709651903825254`*^9}, {
  3.873293389452095*^9, 3.873293396580164*^9}},
 CellLabel->
  "In[118]:=",ExpressionUUID->"b11329de-3c58-4555-b3d4-ff126824b21e"],

Cell[BoxData[
 RowBox[{"gfd12", "=", 
  RowBox[{"Overlay", "[", 
   RowBox[{"{", 
    RowBox[{"gf12", ",", "gd"}], "}"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.87096488840285*^9, 3.8709648918741126`*^9}, {
  3.873293407505667*^9, 3.8732934158881917`*^9}},
 CellLabel->
  "In[119]:=",ExpressionUUID->"77836a2c-fac0-4f84-bb7a-56b14d7bcab4"]
}, Closed]],

Cell[CellGroupData[{

Cell["Joint parameters", "Subsection",
 CellChangeTimes->{{3.871043208282177*^9, 
  3.8710432117624273`*^9}},ExpressionUUID->"f63965c0-0312-4951-a10f-\
2e08913e4d42"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Tabc", "[", 
   RowBox[{"aux1_", ",", "aux2_", ",", "l1_", ",", "l2_"}], "]"}], ":=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{"1", ";;", "l1"}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["a", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{"l1", "+", "1"}], ";;", 
       RowBox[{"2", "*", "l1"}]}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["c", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"i", "+", "l1"}], ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}], ",", 
    RowBox[{"aux1", "[", 
     RowBox[{"[", 
      RowBox[{
       RowBox[{
        RowBox[{"2", "*", "l1"}], "+", "1"}], ";;", 
       RowBox[{"3", "*", "l1"}]}], "]"}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        SubscriptBox["b", 
         RowBox[{"i", "+", "l1"}]], ",", 
        RowBox[{"aux2", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"i", "+", 
            RowBox[{"2", "*", "l1"}]}], ",", "2"}], "]"}], "]"}]}], "}"}], 
      ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "l2"}], "}"}]}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.87104457892337*^9, 3.87104471153677*^9}, 
   3.871044747799478*^9, {3.8710448500107675`*^9, 3.8710448741324635`*^9}},
 CellLabel->
  "In[123]:=",ExpressionUUID->"e5adb368-73bb-44cc-b62a-e93d1cca2c6e"],

Cell[BoxData[
 RowBox[{"Pars", "=", 
  RowBox[{"Tabc", "[", 
   RowBox[{"Pars1", ",", "Pars2", ",", "l1", ",", "l2"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8710447678753657`*^9, 3.8710447920978208`*^9}, {
  3.873294763084342*^9, 3.8732947711704564`*^9}},
 CellLabel->
  "In[130]:=",ExpressionUUID->"9d92ad06-e5d8-41f7-9a52-7bf087179751"],

Cell[BoxData[
 RowBox[{"ParsM", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"Pars", "[", 
     RowBox[{"[", 
      RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"3", "*", 
       RowBox[{"(", 
        RowBox[{"l1", "+", "l2"}], ")"}]}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.8732949860616693`*^9, 3.873295007979524*^9}, 
   3.8732950843915386`*^9, {3.8732951590278554`*^9, 3.873295170659943*^9}, {
   3.8732952059978976`*^9, 3.8732952709208856`*^9}, {3.873295319556636*^9, 
   3.8732953247722964`*^9}},
 CellLabel->
  "In[131]:=",ExpressionUUID->"454c593b-a8c0-4ba8-8bd0-14bdf10290a3"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Export", "[", 
   RowBox[{
   "\"\<D:\\\\sousa\\\\pesquisa\\\\Covid19\\\\GnuPlot\\\\teste0.txt\>\"", ",",
     "ParsM", ",", "\"\<CSV\>\""}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.8732953716610055`*^9, 3.8732953847095757`*^9}},
 CellLabel->
  "In[132]:=",ExpressionUUID->"b2cbd333-a4e4-4eef-9685-6fe25eb047b3"]
}, Closed]]
}, Closed]]
}, Open  ]]
},
WindowSize->{951, 691.1999999999999},
WindowMargins->{{-1009.8000000000001`, Automatic}, {
  Automatic, 49.199999999999996`}},
TaggingRules-><|"TryRealOnly" -> False|>,
Magnification:>1.,
FrontEndVersion->"13.3 for Microsoft Windows (64-bit) (July 24, 2023)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"c8bdd9ec-4051-4b6a-93a9-0270cef2231b"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 782, 15, 146, "Text",ExpressionUUID->"79c1882b-12d4-4336-94a2-cac77e8e48ae"],
Cell[CellGroupData[{
Cell[1365, 39, 236, 3, 53, "Subtitle",ExpressionUUID->"80b7cff6-1a9e-472c-a1ed-a6b22e42e134"],
Cell[1604, 44, 1137, 30, 99, "Text",ExpressionUUID->"a6adaaa0-c19e-4e36-9804-31f280d0559d"],
Cell[2744, 76, 1023, 27, 85, "Text",ExpressionUUID->"be5c8736-934b-45a1-b64b-475e840fbab2"],
Cell[CellGroupData[{
Cell[3792, 107, 155, 3, 67, "Section",ExpressionUUID->"b45ff5bc-1a50-406c-a509-76e4636410a4"],
Cell[CellGroupData[{
Cell[3972, 114, 163, 3, 54, "Subsection",ExpressionUUID->"fc505297-3fe8-4104-a3c5-30fc749ac8aa"],
Cell[4138, 119, 882, 25, 28, "Input",ExpressionUUID->"47d21065-13a5-4b2b-98f6-aeb3fbfd2b51"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5057, 149, 153, 3, 54, "Subsection",ExpressionUUID->"53272786-40cd-4f50-ab15-7a386e30b2df"],
Cell[5213, 154, 514, 15, 28, "Input",ExpressionUUID->"0bd03155-ebb6-412b-a6cb-6b20ef85d23d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5764, 174, 167, 3, 54, "Subsection",ExpressionUUID->"dc210d35-2b83-4392-ab6b-b8809f2fc4c8"],
Cell[5934, 179, 1816, 57, 48, "Input",ExpressionUUID->"422317cb-d461-417d-a7f0-3f5a8e767875"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[7799, 242, 151, 3, 53, "Section",ExpressionUUID->"9c344c06-58c8-4128-8d89-219690e1ae32"],
Cell[7953, 247, 492, 14, 196, "Text",ExpressionUUID->"ddbec865-54cf-447b-a02a-37d7cdd5e972"],
Cell[8448, 263, 512, 10, 28, "Input",ExpressionUUID->"980a1573-88cf-4917-88bc-c1dbb49d1b83"],
Cell[8963, 275, 291, 6, 28, "Input",ExpressionUUID->"54e70c49-4f26-48b1-abfc-14877e4d3bb9"],
Cell[CellGroupData[{
Cell[9279, 285, 188, 3, 54, "Subsection",ExpressionUUID->"db6dd233-28c1-46ce-99c3-4206c96e4105"],
Cell[9470, 290, 204, 3, 28, "Input",ExpressionUUID->"7d28f777-8ad8-4b2d-a652-b29e73251dbf"],
Cell[9677, 295, 249, 4, 28, "Input",ExpressionUUID->"d8ad62c2-5ddc-41be-b43e-02c255fe0644"],
Cell[9929, 301, 898, 18, 28, "Input",ExpressionUUID->"3f170062-35df-41b4-ad55-16984b8c95b3"],
Cell[10830, 321, 600, 15, 84, "Input",ExpressionUUID->"0f1e731c-9ae9-4b51-a5fc-1464839c4b4f"],
Cell[11433, 338, 1774, 39, 53, "Input",ExpressionUUID->"93adb5a8-31bc-47c2-9007-b0ca9312aeba"],
Cell[13210, 379, 173, 2, 28, "Input",ExpressionUUID->"bafb72b9-4a4a-4f67-9032-3cc670c401e7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13420, 386, 156, 3, 54, "Subsection",ExpressionUUID->"251fb970-e571-4ea5-af54-71a1ac49be81"],
Cell[13579, 391, 845, 17, 28, "Input",ExpressionUUID->"575b742e-f6ff-42b1-9c1f-dd0ac86b7797"],
Cell[14427, 410, 547, 14, 28, "Input",ExpressionUUID->"74bb772b-acfd-4111-8f3b-66b66a492dce"],
Cell[14977, 426, 1641, 37, 53, "Input",ExpressionUUID->"73f1e0f4-fc18-48e7-80f0-9136d9878e48"],
Cell[16621, 465, 175, 2, 28, "Input",ExpressionUUID->"30dcb392-f692-4697-be0c-86025840ddf8"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[16845, 473, 210, 4, 53, "Section",ExpressionUUID->"c72d7ca5-a1cb-4eec-a1db-1087298cfd1f"],
Cell[CellGroupData[{
Cell[17080, 481, 160, 3, 54, "Subsection",ExpressionUUID->"d4d8a613-d3ca-4a7e-b8dd-99e76939338b"],
Cell[17243, 486, 341, 6, 28, "Input",ExpressionUUID->"a3644006-c726-4b08-a19d-da718b514fa2"],
Cell[17587, 494, 1003, 19, 28, "Input",ExpressionUUID->"6738a859-1b22-42cf-b082-a48fda9110c8"],
Cell[18593, 515, 695, 16, 28, "Input",ExpressionUUID->"e2d9512f-534c-404e-8f63-766e589a7778"],
Cell[19291, 533, 1392, 20, 28, "Input",ExpressionUUID->"919a2a64-fcd6-4a02-bb5f-c4a64de8d9f1"],
Cell[20686, 555, 847, 24, 28, "Input",ExpressionUUID->"e31eb85c-ab8e-4ea3-9b33-b8a95000e4ad"],
Cell[21536, 581, 865, 25, 28, "Input",ExpressionUUID->"6ab1d067-a4b2-4f2e-b517-208c21e9bf9d"],
Cell[22404, 608, 4554, 81, 53, "Input",ExpressionUUID->"cd1e4f9f-3a61-44bb-ac03-b564a3cecf7f"],
Cell[26961, 691, 1619, 29, 28, "Input",ExpressionUUID->"3ed0b713-0c5b-4ea2-b5a3-8672433c3d83"],
Cell[28583, 722, 843, 18, 48, "Input",ExpressionUUID->"9105c868-3b2a-4597-8225-05d441f766fc"],
Cell[29429, 742, 710, 16, 28, "Input",ExpressionUUID->"c4d5c03b-8dd4-42f8-9de4-afb63c3067bf"],
Cell[30142, 760, 259, 5, 28, "Input",ExpressionUUID->"57a46e14-df0e-4679-80dc-0c271f6540c5"],
Cell[30404, 767, 1241, 35, 48, "Input",ExpressionUUID->"1f2b68f9-5d7a-4ed7-ac66-7cf93ec5aafc"],
Cell[31648, 804, 3516, 99, 162, "Input",ExpressionUUID->"8e8b57ef-6f28-48ad-be27-d4173bae4035"],
Cell[35167, 905, 153, 3, 35, "Text",ExpressionUUID->"11bab72b-c843-498d-9164-6321269a8360"],
Cell[35323, 910, 1020, 28, 28, "Input",ExpressionUUID->"1fa3207b-9d89-47ac-bfa1-bdd2a6414e32"],
Cell[36346, 940, 1432, 23, 28, "Input",ExpressionUUID->"bec225a2-9ef8-4ad0-8be2-c60327f8742f"],
Cell[37781, 965, 1860, 45, 72, "Input",ExpressionUUID->"e9485626-019c-4a45-b231-29b0cbfff40f"],
Cell[39644, 1012, 389, 8, 28, "Input",ExpressionUUID->"7f856a76-b2fa-4725-a20a-a52ee624559e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40070, 1025, 211, 4, 54, "Subsection",ExpressionUUID->"898d8142-4872-4e3e-b06a-f54f94a176e7"],
Cell[40284, 1031, 251, 4, 28, "Input",ExpressionUUID->"cc359a44-f488-4c9d-8a11-5974229044a0"],
Cell[40538, 1037, 1110, 20, 28, "Input",ExpressionUUID->"591fa2db-2ae7-4594-b898-ef97114f8d2a"],
Cell[41651, 1059, 804, 17, 28, "Input",ExpressionUUID->"2e275371-0b11-4809-8d63-0fc8848217d7"],
Cell[42458, 1078, 253, 5, 28, "Input",ExpressionUUID->"96b5e056-1cbc-4e35-b45f-902d0be5a8cf"],
Cell[42714, 1085, 542, 9, 28, "Input",ExpressionUUID->"c8410cea-0d57-4af9-8faa-456d7cfb36f1"],
Cell[43259, 1096, 1139, 32, 28, "Input",ExpressionUUID->"7b9b9d83-1d7b-4250-a272-5215425484e7"],
Cell[44401, 1130, 1201, 33, 28, "Input",ExpressionUUID->"4d5ff2a2-02e8-4315-9690-f16b5cc9b4e7"],
Cell[45605, 1165, 1328, 36, 53, "Input",ExpressionUUID->"a358bd6e-2fd3-441d-aa07-e16a19918c5a"],
Cell[46936, 1203, 3714, 75, 53, "Input",ExpressionUUID->"67ec7e92-a3a0-4bca-8809-a0da020b024c"],
Cell[50653, 1280, 1512, 28, 28, "Input",ExpressionUUID->"fac56909-b1c8-42cb-a95a-1f5224613c2e"],
Cell[52168, 1310, 1005, 21, 48, "Input",ExpressionUUID->"26ebda72-b06a-45e0-8ca4-96517c355222"],
Cell[53176, 1333, 839, 19, 28, "Input",ExpressionUUID->"94d53b18-115d-43c0-a73d-0bd67be707cc"],
Cell[54018, 1354, 310, 6, 28, "Input",ExpressionUUID->"8efde151-d94d-42ca-bf34-b0f3a89d6f74"],
Cell[54331, 1362, 1417, 41, 48, "Input",ExpressionUUID->"efeb5087-f51d-4c95-bc6d-86748962196b"],
Cell[55751, 1405, 153, 3, 35, "Text",ExpressionUUID->"7c5c632c-b2ea-4e74-97fc-7435fb030ba7"],
Cell[55907, 1410, 1218, 33, 28, "Input",ExpressionUUID->"76ece0fb-aa74-4bee-b252-701dfbcb7121"],
Cell[57128, 1445, 1487, 24, 28, "Input",ExpressionUUID->"fdfd2541-f325-4eaa-935c-4e3380c865cd"],
Cell[58618, 1471, 1843, 46, 72, "Input",ExpressionUUID->"93a2d356-f1fe-4709-9c12-be1d27061769"],
Cell[60464, 1519, 436, 9, 28, "Input",ExpressionUUID->"7091e121-da3f-4501-a591-a4e36a5e2a4b"],
Cell[60903, 1530, 1830, 48, 72, "Input",ExpressionUUID->"cb8d7e11-0f44-4528-b617-63db8ef0482e"],
Cell[62736, 1580, 290, 7, 28, "Input",ExpressionUUID->"78b307fa-af54-40f7-b693-914721405c0d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[63063, 1592, 166, 3, 54, "Subsection",ExpressionUUID->"9bca91aa-a2f1-4d56-b884-426f694047c5"],
Cell[63232, 1597, 2034, 61, 48, "Input",ExpressionUUID->"fc152712-ee36-400a-a424-9767fbd1c621"],
Cell[65269, 1660, 394, 8, 28, "Input",ExpressionUUID->"976ae3e0-57b6-4f17-bc69-1484359b4d6c"],
Cell[65666, 1670, 425, 13, 28, "Input",ExpressionUUID->"28e2fc1b-f9f1-42ed-a30a-75f579303514"],
Cell[66094, 1685, 493, 14, 28, "Input",ExpressionUUID->"3f6348e1-0fef-41ca-ad70-bcefe4ea95bd"],
Cell[66590, 1701, 487, 10, 28, "Input",ExpressionUUID->"864be332-a77f-4354-8f69-a39f487ee6ba"],
Cell[67080, 1713, 5919, 169, 276, "Input",ExpressionUUID->"42406004-0047-4213-ba5e-c5607d601621"],
Cell[73002, 1884, 5858, 162, 276, "Input",ExpressionUUID->"d0792500-13e2-4b2d-b73d-d3eb0e0c7e5f"]
}, Open  ]],
Cell[CellGroupData[{
Cell[78897, 2051, 155, 3, 54, "Subsection",ExpressionUUID->"8cf1d87c-89f4-4242-a4a4-d1ac79da9065"],
Cell[79055, 2056, 345, 5, 28, "Input",ExpressionUUID->"fbd3b7f1-0622-47fe-819b-e3174ee7a4af"],
Cell[79403, 2063, 789, 23, 28, "Input",ExpressionUUID->"4931fa6b-24a0-4de9-8358-9929a13a1039"],
Cell[80195, 2088, 806, 23, 28, "Input",ExpressionUUID->"1018d7f9-8778-4c69-8363-10123b67c0e6"],
Cell[81004, 2113, 938, 27, 31, "Input",ExpressionUUID->"acb2c888-3248-4660-90d6-24e25e67ccd9"],
Cell[81945, 2142, 2626, 56, 31, "Input",ExpressionUUID->"add0452d-0f32-4abe-9d7c-2a8acf596353"],
Cell[84574, 2200, 1781, 44, 31, "Input",ExpressionUUID->"f476cebe-9218-4d02-acf6-89cd076a668b"],
Cell[86358, 2246, 254, 4, 28, "Input",ExpressionUUID->"4a9d94fb-fd4b-4a8a-8546-fb8d42ab4015"],
Cell[86615, 2252, 1380, 25, 28, "Input",ExpressionUUID->"d87c187a-9670-4a1c-a83a-fe19bb88633c"],
Cell[87998, 2279, 679, 15, 48, "Input",ExpressionUUID->"9c5934b0-179f-45d6-97c9-3fc57e8e5efb"],
Cell[88680, 2296, 349, 5, 28, "Input",ExpressionUUID->"86ba14f5-0ec5-4759-bfcb-1e5b167a4368"],
Cell[89032, 2303, 696, 15, 28, "Input",ExpressionUUID->"22a9ea35-674f-41a0-91e3-d0ab3d537c4f"],
Cell[89731, 2320, 308, 6, 28, "Input",ExpressionUUID->"597b583c-03fa-42e3-b48e-cb6d8b30e33d"],
Cell[90042, 2328, 1329, 38, 28, "Input",ExpressionUUID->"dc637d05-948b-4e9c-902b-2879c62737a3"],
Cell[91374, 2368, 153, 3, 35, "Text",ExpressionUUID->"f1614df4-cb52-445d-bbc7-8abd06148c4a"],
Cell[91530, 2373, 923, 26, 28, "Input",ExpressionUUID->"097d0611-bab5-45cf-953e-507ade9e8796"],
Cell[92456, 2401, 1372, 22, 28, "Input",ExpressionUUID->"55c9bd54-93f6-4d7c-aa86-732147512f0c"],
Cell[93831, 2425, 1807, 45, 53, "Input",ExpressionUUID->"78076fd0-a60f-4a22-8db5-5bfac7cb10b4"],
Cell[95641, 2472, 545, 11, 28, "Input",ExpressionUUID->"7bbb9ee2-3ea8-479c-b137-646676e9d40b"]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell[96235, 2489, 151, 3, 53, "Section",ExpressionUUID->"d335295a-360c-4ea8-8148-12a2a7be3e4e"],
Cell[CellGroupData[{
Cell[96411, 2496, 155, 3, 54, "Subsection",ExpressionUUID->"cd9101aa-571f-406e-b312-a006b1c92f75"],
Cell[96569, 2501, 248, 5, 27, "Input",ExpressionUUID->"bed05695-006c-4eb6-878c-65a5ce4266ec"],
Cell[96820, 2508, 793, 24, 27, "Input",ExpressionUUID->"1879b5b4-edfd-4554-9894-8a66c4454fb6"],
Cell[97616, 2534, 806, 23, 27, "Input",ExpressionUUID->"63a1013e-93a8-4fb9-8ce3-cb23f701b3b0"],
Cell[98425, 2559, 938, 27, 32, "Input",ExpressionUUID->"2e13a3de-eed2-4653-a637-d329f7322a8d"],
Cell[99366, 2588, 2079, 48, 32, "Input",ExpressionUUID->"b4d5036b-b7eb-45ad-aa98-48ef609fe014"],
Cell[101448, 2638, 254, 4, 27, "Input",ExpressionUUID->"512de362-116b-49b3-a5a3-9376312c59ce"],
Cell[101705, 2644, 1384, 26, 27, "Input",ExpressionUUID->"e2c58dbd-4bd7-45a8-bcda-e9bcc73e5b5e"],
Cell[103092, 2672, 683, 16, 48, "Input",ExpressionUUID->"c867c94f-59e0-4ce0-aed0-5a3a04d1156c"],
Cell[103778, 2690, 632, 16, 27, "Input",ExpressionUUID->"de9ca565-4254-45cb-9c47-ad85f95c5d08"],
Cell[104413, 2708, 308, 6, 27, "Input",ExpressionUUID->"82be80d3-9d7c-4d4d-b14a-843b3fbc4133"],
Cell[104724, 2716, 1307, 38, 27, "Input",ExpressionUUID->"e7acdc00-60c1-43ac-a816-4a2138f6c147"],
Cell[106034, 2756, 153, 3, 33, "Text",ExpressionUUID->"9fccfbc0-db4c-4aae-aaae-9360d4e2cefd"],
Cell[106190, 2761, 939, 27, 27, "Input",ExpressionUUID->"5ad82bff-0fbd-4410-b058-72e8687a6875"],
Cell[107132, 2790, 1377, 22, 27, "Input",ExpressionUUID->"65396194-ba47-41a1-bfd2-30afc764a0dd"],
Cell[108512, 2814, 1854, 46, 72, "Input",ExpressionUUID->"ba9550fc-392b-41cd-9870-c14997dbafed"],
Cell[110369, 2862, 492, 10, 27, "Input",ExpressionUUID->"4b9c3dd4-19c8-4b6c-9f60-581ef19a4b93"]
}, Closed]],
Cell[CellGroupData[{
Cell[110898, 2877, 160, 3, 38, "Subsection",ExpressionUUID->"66afd28b-f580-40f2-ab1d-a1e947869f46"],
Cell[111061, 2882, 900, 18, 27, "Input",ExpressionUUID->"91f8ecca-2e58-44b4-b65c-f9330de3ae66"],
Cell[111964, 2902, 597, 15, 27, "Input",ExpressionUUID->"0034a340-1fa2-49cf-94a0-5ee247f6966a"],
Cell[112564, 2919, 294, 5, 27, "Input",ExpressionUUID->"10542a1f-a53e-43f1-b326-90c5553b8c2a"],
Cell[112861, 2926, 842, 23, 27, "Input",ExpressionUUID->"a8fbbeb5-886e-45cd-936b-c45a66e7aa9e"],
Cell[113706, 2951, 860, 24, 27, "Input",ExpressionUUID->"0280e91f-a6d9-406f-ab00-3b1749ed26b5"],
Cell[114569, 2977, 993, 28, 32, "Input",ExpressionUUID->"1bc98aba-dac9-4ea4-8490-749971a7b616"],
Cell[115565, 3007, 2360, 51, 53, "Input",ExpressionUUID->"cc34ef7c-6dbe-40fa-b47a-f9f441dcb0d3"],
Cell[117928, 3060, 1461, 26, 27, "Input",ExpressionUUID->"ec5a62f9-cb07-4b16-b983-142ac3963670"],
Cell[119392, 3088, 811, 17, 48, "Input",ExpressionUUID->"9ea867a0-6970-4605-a870-2b09a8a9bc87"],
Cell[120206, 3107, 655, 15, 27, "Input",ExpressionUUID->"e2c5f4bc-43a9-4131-9092-e9138138ce44"],
Cell[120864, 3124, 255, 4, 27, "Input",ExpressionUUID->"3ac0b1c9-30ff-4706-8b2f-0daa32fc1b7e"],
Cell[121122, 3130, 1239, 35, 27, "Input",ExpressionUUID->"db259c83-f156-4b2b-aca3-de1d01a49e95"],
Cell[122364, 3167, 2746, 79, 123, "Input",ExpressionUUID->"8cd3b1de-68b1-4652-a162-6b090fc37590"],
Cell[125113, 3248, 1754, 43, 72, "Input",ExpressionUUID->"9970887a-0428-45a5-8181-dd5e56a8e31a"],
Cell[126870, 3293, 439, 8, 27, "Input",ExpressionUUID->"a17d4ce1-fe07-4bad-a2e9-3ea827c395ef"]
}, Closed]],
Cell[CellGroupData[{
Cell[127346, 3306, 211, 4, 38, "Subsection",ExpressionUUID->"96d408fb-f1c8-434a-90e1-fdea94cb786e"],
Cell[127560, 3312, 949, 18, 28, "Input",ExpressionUUID->"4b4c2ec0-15a2-45f9-88a7-f00cad117989"],
Cell[128512, 3332, 644, 15, 28, "Input",ExpressionUUID->"fd4630b0-589d-48fc-8e5f-201a22f34466"],
Cell[129159, 3349, 340, 5, 28, "Input",ExpressionUUID->"8e6d0dfd-6454-47a3-9b72-2eece97625b7"],
Cell[129502, 3356, 1135, 31, 28, "Input",ExpressionUUID->"42480e1c-df2c-4c77-a303-7306d3526e34"],
Cell[130640, 3389, 1197, 32, 28, "Input",ExpressionUUID->"c8d0907d-8e73-4d52-b6a2-e861bf04cc74"],
Cell[131840, 3423, 1328, 36, 32, "Input",ExpressionUUID->"3e77d8a9-d9a3-4ff7-bc6a-f522b8824ea1"],
Cell[133171, 3461, 3281, 69, 53, "Input",ExpressionUUID->"9e7cfe2d-4a5c-4b8d-9589-ad0a92c450fb"],
Cell[136455, 3532, 1463, 27, 28, "Input",ExpressionUUID->"7a0df496-b5a8-4638-8bff-f764498bb1f9"],
Cell[137921, 3561, 868, 19, 48, "Input",ExpressionUUID->"7d833c5d-f61f-4775-9b56-d89bc598521b"],
Cell[138792, 3582, 735, 17, 27, "Input",ExpressionUUID->"998e1c25-cff4-4a1f-99e3-4afdc0f5ee1d"],
Cell[139530, 3601, 310, 6, 27, "Input",ExpressionUUID->"7fe546e9-3de1-493c-8231-218e620333f3"],
Cell[139843, 3609, 1417, 41, 48, "Input",ExpressionUUID->"64cf2943-8185-47e1-b396-99cff3998851"],
Cell[141263, 3652, 178, 4, 27, "Input",ExpressionUUID->"35b9f060-be60-4061-89db-dc30a01ab22e"],
Cell[141444, 3658, 1839, 46, 72, "Input",ExpressionUUID->"548568e7-4e57-4b8e-8a5d-f3fcb236bb65"],
Cell[143286, 3706, 487, 10, 27, "Input",ExpressionUUID->"f62348a3-aa7a-4568-acf5-a6f840dc5307"],
Cell[143776, 3718, 1873, 49, 72, "Input",ExpressionUUID->"b11329de-3c58-4555-b3d4-ff126824b21e"],
Cell[145652, 3769, 341, 8, 27, "Input",ExpressionUUID->"77836a2c-fac0-4f84-bb7a-56b14d7bcab4"]
}, Closed]],
Cell[CellGroupData[{
Cell[146030, 3782, 166, 3, 38, "Subsection",ExpressionUUID->"f63965c0-0312-4951-a10f-2e08913e4d42"],
Cell[146199, 3787, 1982, 60, 48, "Input",ExpressionUUID->"e5adb368-73bb-44cc-b62a-e93d1cca2c6e"],
Cell[148184, 3849, 341, 7, 28, "Input",ExpressionUUID->"9d92ad06-e5d8-41f7-9a52-7bf087179751"],
Cell[148528, 3858, 675, 17, 28, "Input",ExpressionUUID->"454c593b-a8c0-4ba8-8bd0-14bdf10290a3"],
Cell[149206, 3877, 353, 8, 28, "Input",ExpressionUUID->"b2cbd333-a4e4-4eef-9685-6fe25eb047b3"]
}, Closed]]
}, Closed]]
}, Open  ]]
}
]
*)