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

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

(* CreatedBy='Mathematica 10.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     21770,        596]
NotebookOptionsPosition[     16667,        516]
NotebookOutlinePosition[     17100,        533]
CellTagsIndexPosition[     17057,        530]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Program Modularization", "Title",
 CellChangeTimes->{{3.6242127564551*^9, 3.624212773206312*^9}, {
  3.683010423897872*^9, 
  3.683010429203343*^9}},ExpressionUUID->"5116efad-4bef-4544-97c3-\
c253971d72ba"],

Cell[CellGroupData[{

Cell["Package construction by example", "Section",
 CellChangeTimes->{{3.683010459637523*^9, 
  3.6830104672503195`*^9}},ExpressionUUID->"2eaa89f1-c48d-4b81-aeeb-\
ba9aebccc3c7"],

Cell[CellGroupData[{

Cell["\<\
This package will contain the simulation of the straight line parameter fit, 
and a function to extrapolate values using the fitted parameters.
This can be used to check whether the uncertainty in the extrapolated value \
is calculated properly\
\>", "Subsection",
 CellChangeTimes->{{3.6830104782202573`*^9, 3.683010512684145*^9}, {
   3.6830106320313773`*^9, 3.6830106600314064`*^9}, {3.6830107201403217`*^9, 
   3.6830108139623165`*^9}, {3.6830111843948107`*^9, 3.6830112008305235`*^9}, 
   3.8915240157094994`*^9},ExpressionUUID->"ab1ca37d-80d8-4d80-9879-\
b911d9a52972"],

Cell[CellGroupData[{

Cell["\<\
The package must be built in another file. It is recommended to develop all \
lines of code in a notebook and transfer to the package only after being \
completely sure that they work\
\>", "Subsubsection",
 CellChangeTimes->{{3.6830108402881255`*^9, 
  3.683010896126658*^9}},ExpressionUUID->"d0434c09-58ee-495c-aeb5-\
c36194719e08"],

Cell[BoxData[
 RowBox[{
  RowBox[{"a0", "=", 
   RowBox[{"-", "2.5"}]}], ";", " ", 
  RowBox[{"b0", "=", "0.225"}], ";", 
  RowBox[{"\[Sigma]i", "=", "0.3"}], ";", " ", 
  RowBox[{"x", "=", 
   RowBox[{"Range", "[", 
    RowBox[{"4", ",", "24", ",", "4"}], "]"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"straight", " ", "line", " ", "parameters"}], ",", " ", 
     RowBox[{
     "standard", " ", "deviation", " ", "of", " ", "the", " ", "y", " ", 
      "and", " ", "dependent", " ", "values", " ", 
      RowBox[{"list", ":", " ", 
       RowBox[{"they", " ", "will", " ", "be", " ", "parameters"}]}]}]}], " ",
     "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.624212783097642*^9, 3.624212870614581*^9}, {
  3.6830109026972017`*^9, 
  3.6830110305244236`*^9}},ExpressionUUID->"e6100a27-427c-4029-9981-\
8e556e9e473b"],

Cell[BoxData[
 RowBox[{"data", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ",", 
      RowBox[{"RandomVariate", "[", " ", 
       RowBox[{"NormalDistribution", "[", 
        RowBox[{
         RowBox[{"a0", "+", 
          RowBox[{"b0", " ", 
           RowBox[{"x", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}]}], ",", "\[Sigma]i"}], "]"}], 
       " ", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", 
      RowBox[{"Length", "[", "x", "]"}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6242128856113973`*^9, 3.6242129434099207`*^9}, {
  3.6830110461837826`*^9, 
  3.683011050320692*^9}},ExpressionUUID->"4de2b797-2ed0-4948-afe2-\
2afb1a0da9ac"],

Cell[BoxData[
 RowBox[{"ListPlot", "[", "data", "]"}]], "Input",
 CellChangeTimes->{{3.6242129478321247`*^9, 3.6242129534106436`*^9}, {
  3.6830110622016315`*^9, 
  3.6830110630692935`*^9}},ExpressionUUID->"e4fcca9f-9a1e-419a-9147-\
b6802b00546f"],

Cell[BoxData[
 RowBox[{"fit", "=", 
  RowBox[{"LinearModelFit", "[", 
   RowBox[{"data", ",", 
    RowBox[{"{", "z", "}"}], ",", "z", ",", 
    RowBox[{"Weights", "\[Rule]", 
     RowBox[{"1", "/", 
      RowBox[{"ConstantArray", "[", 
       RowBox[{
        RowBox[{"\[Sigma]i", "^", "2"}], ",", 
        RowBox[{"Length", "[", "x", "]"}]}], " ", "]"}]}]}]}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "x", " ", "is", " ", "the", " ", "name", " ", "of", " ", "the", " ", 
     "list", " ", "of", " ", "independent", " ", "values"}], ",", " ", 
    RowBox[{"cannot", " ", "be", " ", "used", " ", "as", " ", "argument"}]}], 
   " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.6242129782874317`*^9, 3.624213076973333*^9}, {
  3.6242131079130683`*^9, 3.624213108756874*^9}, {3.6830110709125433`*^9, 
  3.6830111144421597`*^9}},ExpressionUUID->"c2872320-8c70-41ac-bd56-\
27642bb8c20b"],

Cell[BoxData[
 RowBox[{"\[Chi]2", "=", 
  RowBox[{"Total", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<FitResiduals\>\"", "]"}], "/", "\[Sigma]i"}], 
     ")"}], "^", "2"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6242130864583974`*^9, 
  3.624213158880863*^9}},ExpressionUUID->"151e5462-35b6-41d3-a538-\
edbaceef501f"],

Cell[BoxData[
 RowBox[{"v", "=", 
  RowBox[{"fit", "[", "\"\<CovarianceMatrix\>\"", "]"}]}]], "Input",
 CellChangeTimes->{{3.6242131779398384`*^9, 3.624213185487259*^9}, {
  3.6242132699730735`*^9, 
  3.6242132704731045`*^9}},ExpressionUUID->"6776155a-e327-434e-8780-\
344f0f519bdf"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
This is how we calculate the uncertainty of the extrapolated point at x=40\
\>", "Subsubsection",
 CellChangeTimes->{{3.6830111418488445`*^9, 
  3.6830111638485546`*^9}},ExpressionUUID->"91fc4886-9cba-43d9-aa94-\
87f51932486b"],

Cell[BoxData[
 RowBox[{"\[Sigma]y", "=", 
  RowBox[{"Sqrt", "[", " ", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"1", ",", "40"}], "}"}], ".", "v", ".", 
    RowBox[{"{", 
     RowBox[{"1", ",", "40"}], "}"}]}], " ", "]"}]}]], "Input",
 CellChangeTimes->{{3.6242132095202384`*^9, 
  3.6242132874587126`*^9}},ExpressionUUID->"67d97919-769c-4d20-9dcf-\
70c5b2e6c7e8"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Building the package", "Section",
 CellChangeTimes->{{3.6830112720420914`*^9, 
  3.6830112777182393`*^9}},ExpressionUUID->"59b1e21c-cb1d-4baa-a054-\
32bb76389c34"],

Cell["\<\
This is done in another file: select the drop-down menu File above and select \
New, then Package (.wl).
All comands should be entered in this new file. Save it before starting.\
\>", "Text",
 CellChangeTimes->{{3.6830112854458475`*^9, 3.683011289051996*^9}, {
  3.683011336780436*^9, 
  3.683011433186225*^9}},ExpressionUUID->"ac2cf1c0-e352-45bd-ac3d-\
b95b57f8a1be"],

Cell["\<\
Besides the scoping instructions, there is the ::usage keyword, that can be \
used to document functions in any notebook:\
\>", "Text",
 CellChangeTimes->{{3.683011445081148*^9, 
  3.6830114911826005`*^9}},ExpressionUUID->"f5d2d3b8-cdf8-4fb6-a072-\
989d491e8a80"],

Cell[BoxData[
 RowBox[{
  RowBox[{"func", "[", "x_", "]"}], ":=", 
  RowBox[{"x", "^", "2"}]}]], "Input",
 CellChangeTimes->{{3.6830114944710555`*^9, 
  3.68301150335005*^9}},ExpressionUUID->"7e23e465-d4d5-42eb-905a-\
90b77d4eb958"],

Cell[BoxData[
 RowBox[{
  RowBox[{"func", "::", "usage"}], " ", "=", " ", 
  "\"\<func[x] takes a single argument (either a number or a list) and return \
its square.\>\""}]], "Input",
 CellChangeTimes->{{3.6830115092957315`*^9, 
  3.683011545165801*^9}},ExpressionUUID->"5ef25e4b-1ffe-4160-ac1f-\
939a65197cce"],

Cell["this returns the description entered through usage", "Text",
 CellChangeTimes->{
  3.6830117789212575`*^9},ExpressionUUID->"c6cf4484-f82d-4fd8-8add-\
23ce90804ce3"],

Cell[BoxData[
 RowBox[{"?", "func"}]], "Input",
 CellChangeTimes->{{3.6830117277493415`*^9, 
  3.6830117832421627`*^9}},ExpressionUUID->"206855ee-cc9d-410c-b7fe-\
34d6e7fd6e9c"],

Cell["?? returns also the function code ", "Text",
 CellChangeTimes->{{3.6830117920913696`*^9, 
  3.6830118108987694`*^9}},ExpressionUUID->"81a79f43-a8db-4cb0-9859-\
300f9c4150fd"],

Cell[BoxData[
 RowBox[{"??", "func"}]], "Input",
 CellChangeTimes->{{3.683011814860547*^9, 
  3.683011817009244*^9}},ExpressionUUID->"0f55b598-a73b-4554-8d91-\
60c209b75b1a"]
}, Open  ]],

Cell[CellGroupData[{

Cell["testing the package", "Section",
 CellChangeTimes->{{3.6242142409887433`*^9, 3.6242142474892087`*^9}, {
  3.683011830183536*^9, 
  3.6830118335684566`*^9}},ExpressionUUID->"edb8565f-c68a-4794-90a9-\
ba51172d0da8"],

Cell[BoxData[
 RowBox[{"SetDirectory", "[", " ", 
  RowBox[{"NotebookDirectory", "[", "]"}], " ", "]"}]], "Input",
 CellChangeTimes->{{3.624255594098692*^9, 
  3.624255606130805*^9}},ExpressionUUID->"f84a6db9-fa17-43c0-b2c6-\
c03e52062ef9"],

Cell[BoxData[
 RowBox[{"Needs", "[", "\"\<extrapolation`\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.624214265423806*^9, 
  3.6242142764714746`*^9}},ExpressionUUID->"c7f8860e-ad8c-4a19-b49e-\
28d94efe6820"],

Cell[BoxData[
 RowBox[{"?", "extrapolation`*"}]], "Input",
 CellChangeTimes->{{3.6242142550835075`*^9, 
  3.624214290191208*^9}},ExpressionUUID->"4bf1e9dd-1a4b-436d-bd24-\
cc29455a67de"],

Cell[BoxData[
 RowBox[{"f1", "=", 
  RowBox[{"simul", "[", 
   RowBox[{"a0", ",", "b0", ",", "x", ",", "\[Sigma]i"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6242142965822916`*^9, 
  3.624214313020977*^9}},ExpressionUUID->"7991605f-7f16-444f-9af2-\
bddd6236f124"],

Cell[BoxData[
 RowBox[{"value", "[", 
  RowBox[{"f1", ",", "40"}], "]"}]], "Input",
 CellChangeTimes->{{3.624214318237507*^9, 
  3.624214326425603*^9}},ExpressionUUID->"9d0a1263-6c64-42b0-b4e7-\
686e77c05a73"],

Cell[BoxData[
 RowBox[{
  RowBox[{"u", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"value", "[", " ", 
      RowBox[{
       RowBox[{"simul", "[", 
        RowBox[{"a0", ",", "b0", ",", "x", ",", "\[Sigma]i"}], "]"}], ",", 
       "40"}], " ", "]"}], ",", 
     RowBox[{"{", "1000", "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.624214334832451*^9, 3.6242143872054415`*^9}, 
   3.624255818003399*^9},ExpressionUUID->"02c0f5da-eb78-4600-be36-\
92573e264dcb"],

Cell[BoxData[
 RowBox[{"Mean", "[", "u", "]"}]], "Input",
 CellChangeTimes->{{3.624214396737379*^9, 
  3.6242144039879*^9}},ExpressionUUID->"58c2f45a-b492-47f1-b0f9-a7d75ef5be06"],

Cell[BoxData[
 RowBox[{"StandardDeviation", "[", " ", "u", " ", "]"}]], "Input",
 CellChangeTimes->{{3.624255739454934*^9, 
  3.62425574520535*^9}},ExpressionUUID->"78d17373-701f-418a-ac5e-\
bd19cefc412a"],

Cell[BoxData[
 RowBox[{"\[Sigma]value", "[", 
  RowBox[{"f1", ",", "40"}], "]"}]], "Input",
 CellChangeTimes->{{3.624255718849588*^9, 
  3.6242557305848017`*^9}},ExpressionUUID->"d67d1415-1218-4585-b95e-\
3439cc3cc9e1"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Nest[ ], NestList", "Section",
 CellChangeTimes->{{3.683017228531204*^9, 
  3.6830172337087126`*^9}},ExpressionUUID->"9ca5ec54-9cc0-4873-b655-\
534812f784e1"],

Cell["Nest applies a function repeatedly", "Text",
 CellChangeTimes->{{3.683017257390766*^9, 
  3.683017263479602*^9}},ExpressionUUID->"d498ffda-51ea-4fd6-9882-\
f04a4288fc6c"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Nest", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"#", "^", "2"}], "&"}], ",", "2", ",", "3"}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "last", " ", "parameter", " ", "is", " ", "the", " ", "number", " ", "of", 
    " ", "iterations"}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.6830172662070465`*^9, 3.6830172941124964`*^9}, {
  3.68301745832279*^9, 
  3.6830174811595335`*^9}},ExpressionUUID->"910565be-ff4f-4738-ba30-\
7b3c0ee75ada"],

Cell["\<\
NestList does the same and makes a list with all steps - it helps to \
understand its mechanism\
\>", "Text",
 CellChangeTimes->{{3.6830173046868844`*^9, 
  3.68301735237929*^9}},ExpressionUUID->"4d56c02e-69aa-47a0-8e36-\
e85ec8d06646"],

Cell[BoxData[
 RowBox[{"NestList", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"#", "^", "2"}], "&"}], ",", "2", ",", "3"}], "]"}]], "Input",
 CellChangeTimes->{{3.6830173331876125`*^9, 
  3.6830173338768225`*^9}},ExpressionUUID->"cc6bc007-8a0d-4200-a709-\
1e8943d8f1ba"],

Cell["\<\
Perfect to implementcontinued fractions - remember the algorithm for taking \
the square root? Example for Sqrt[3]\
\>", "Text",
 CellChangeTimes->{{3.6830173590187025`*^9, 3.6830173691936445`*^9}, {
  3.6830174271994786`*^9, 
  3.6830174551094265`*^9}},ExpressionUUID->"9654b588-b976-46c4-a416-\
cf19930bfab3"],

Cell[BoxData[
 RowBox[{
  RowBox[{"NestList", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"#", "+", 
        RowBox[{"3", " ", "/", "#"}]}], ")"}], "/", " ", "2"}], "&"}], ",", 
    "2.", ",", "5"}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "last", " ", "parameter", " ", "is", " ", "the", " ", "number", " ", "of", 
    " ", "iterations"}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.6830173762289267`*^9, 3.6830174236019077`*^9}, {
  3.6830174901044054`*^9, 
  3.6830174937677207`*^9}},ExpressionUUID->"ed3e2321-43b5-4fea-8e21-\
2e1d2449a86a"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Which and Piecewise", "Section",
 CellChangeTimes->{{3.68301750783541*^9, 
  3.683017515407084*^9}},ExpressionUUID->"d2580dc1-1dd1-48a9-a40e-\
be340e60f09d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"selection", "[", 
   RowBox[{"key_", ",", "x_"}], "]"}], ":=", 
  RowBox[{"Which", "[", 
   RowBox[{
    RowBox[{"key", "==", "\"\<k1\>\""}], ",", 
    RowBox[{"x", "^", "2"}], ",", 
    RowBox[{"key", "==", "\"\<k2\>\""}], ",", 
    RowBox[{"Sqrt", "[", "x", "]"}]}], " ", "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"notice", " ", "the"}], "\[Equal]", " ", 
    RowBox[{"for", " ", "comparisons"}]}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.6830175216824017`*^9, 
  3.683017654077058*^9}},ExpressionUUID->"aeb662c5-7536-41e7-8fd2-\
26ba4e80febb"],

Cell[BoxData[
 RowBox[{"selection", "[", 
  RowBox[{"\"\<k1\>\"", ",", "a"}], "]"}]], "Input",
 CellChangeTimes->{{3.683017614213738*^9, 
  3.6830176275806074`*^9}},ExpressionUUID->"27652beb-04d8-4116-8719-\
79f70a33917a"],

Cell[BoxData[
 RowBox[{"selection", "[", 
  RowBox[{"\"\<k2\>\"", ",", "a"}], "]"}]], "Input",
 CellChangeTimes->{{3.6830176676596775`*^9, 
  3.683017667793172*^9}},ExpressionUUID->"31b49bac-6699-43fc-9f69-\
efc948963ad9"],

Cell["Like Cases in C. For functions, use Piecewise", "Text",
 CellChangeTimes->{{3.6830176798764915`*^9, 
  3.6830176917224517`*^9}},ExpressionUUID->"ffc3af29-461e-440c-bdcc-\
9ba04d24870a"],

Cell[BoxData[
 RowBox[{
  RowBox[{"f", "[", "x_", "]"}], ":=", 
  RowBox[{"Piecewise", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "x"}], ",", " ", 
        RowBox[{"x", "<", "0"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"x", "^", "2"}], ",", 
        RowBox[{"0", "<", "x", "<", "1"}]}], "}"}]}], "}"}], ",", "1"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.6830176967297063`*^9, 
  3.6830177651030564`*^9}},ExpressionUUID->"109f0d7f-88e4-40e3-beaf-\
74f8788e200b"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"f", "[", "x", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "2"}], ",", "3"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.683017768097169*^9, 
  3.6830177817959127`*^9}},ExpressionUUID->"7ff268af-1579-421c-a8b4-\
76105893ead7"],

Cell[BoxData[
 RowBox[{"D", "[", 
  RowBox[{
   RowBox[{"f", "[", "z", "]"}], " ", ",", "z"}], "]"}]], "Input",
 CellChangeTimes->{{3.6830178526314187`*^9, 
  3.683017852864189*^9}},ExpressionUUID->"2f3cfa8e-deae-4824-8b67-\
0acdee71453e"],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{"f", "[", "z", "]"}], " ", ",", "z"}], "]"}], "/.", 
    RowBox[{"z", "\[Rule]", "x"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.683017793616143*^9, 3.683017825020549*^9}, {
  3.6830178571793385`*^9, 
  3.6830178575213637`*^9}},ExpressionUUID->"26b1c59a-0a43-49ab-b6bf-\
954bee8aa647"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{"f", "[", "x", "]"}], " ", ",", "x"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", 
      RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "This", " ", "does", " ", "not", " ", "work", " ", "and", " ", "takes", 
    " ", "a", " ", "while", " ", "do", " ", "understand", " ", "why"}], " ", 
   "*)"}]}]], "Input",
 CellChangeTimes->{{3.6830178660321045`*^9, 
  3.6830178931358004`*^9}},ExpressionUUID->"73121a2a-6ac5-4989-acc2-\
53bcc71954cb"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1141.2, 576.6},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
Magnification:>1.3 Inherited,
FrontEndVersion->"13.1 for Microsoft Windows (64-bit) (August 22, 2022)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"af8b7212-e3c6-4f93-892a-d62ac80155a1"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 212, 4, 126, "Title",ExpressionUUID->"5116efad-4bef-4544-97c3-c253971d72ba"],
Cell[CellGroupData[{
Cell[817, 30, 178, 3, 86, "Section",ExpressionUUID->"2eaa89f1-c48d-4b81-aeeb-ba9aebccc3c7"],
Cell[CellGroupData[{
Cell[1020, 37, 585, 10, 176, "Subsection",ExpressionUUID->"ab1ca37d-80d8-4d80-9879-b911d9a52972"],
Cell[CellGroupData[{
Cell[1630, 51, 344, 7, 91, "Subsubsection",ExpressionUUID->"d0434c09-58ee-495c-aeb5-c36194719e08"],
Cell[1977, 60, 833, 21, 87, "Input",ExpressionUUID->"e6100a27-427c-4029-9981-8e556e9e473b"],
Cell[2813, 83, 770, 22, 36, "Input",ExpressionUUID->"4de2b797-2ed0-4948-afe2-2afb1a0da9ac"],
Cell[3586, 107, 247, 5, 36, "Input",ExpressionUUID->"e4fcca9f-9a1e-419a-9147-b6802b00546f"],
Cell[3836, 114, 903, 21, 62, "Input",ExpressionUUID->"c2872320-8c70-41ac-bd56-27642bb8c20b"],
Cell[4742, 137, 357, 10, 36, "Input",ExpressionUUID->"151e5462-35b6-41d3-a538-edbaceef501f"],
Cell[5102, 149, 283, 6, 36, "Input",ExpressionUUID->"6776155a-e327-434e-8780-344f0f519bdf"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5422, 160, 237, 5, 58, "Subsubsection",ExpressionUUID->"91fc4886-9cba-43d9-aa94-87f51932486b"],
Cell[5662, 167, 366, 10, 36, "Input",ExpressionUUID->"67d97919-769c-4d20-9dcf-70c5b2e6c7e8"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[6089, 184, 169, 3, 86, "Section",ExpressionUUID->"59b1e21c-cb1d-4baa-a054-32bb76389c34"],
Cell[6261, 189, 378, 8, 75, "Text",ExpressionUUID->"ac2cf1c0-e352-45bd-ac3d-b95b57f8a1be"],
Cell[6642, 199, 273, 6, 45, "Text",ExpressionUUID->"f5d2d3b8-cdf8-4fb6-a072-989d491e8a80"],
Cell[6918, 207, 232, 6, 36, "Input",ExpressionUUID->"7e23e465-d4d5-42eb-905a-90b77d4eb958"],
Cell[7153, 215, 312, 7, 36, "Input",ExpressionUUID->"5ef25e4b-1ffe-4160-ac1f-939a65197cce"],
Cell[7468, 224, 170, 3, 45, "Text",ExpressionUUID->"c6cf4484-f82d-4fd8-8add-23ce90804ce3"],
Cell[7641, 229, 177, 4, 36, "Input",ExpressionUUID->"206855ee-cc9d-410c-b7fe-34d6e7fd6e9c"],
Cell[7821, 235, 180, 3, 45, "Text",ExpressionUUID->"81a79f43-a8db-4cb0-9859-300f9c4150fd"],
Cell[8004, 240, 174, 4, 36, "Input",ExpressionUUID->"0f55b598-a73b-4554-8d91-60c209b75b1a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8215, 249, 219, 4, 86, "Section",ExpressionUUID->"edb8565f-c68a-4794-90a9-ba51172d0da8"],
Cell[8437, 255, 240, 5, 36, "Input",ExpressionUUID->"f84a6db9-fa17-43c0-b2c6-c03e52062ef9"],
Cell[8680, 262, 207, 4, 36, "Input",ExpressionUUID->"c7f8860e-ad8c-4a19-b49e-28d94efe6820"],
Cell[8890, 268, 186, 4, 36, "Input",ExpressionUUID->"4bf1e9dd-1a4b-436d-bd24-cc29455a67de"],
Cell[9079, 274, 264, 6, 36, "Input",ExpressionUUID->"7991605f-7f16-444f-9af2-bddd6236f124"],
Cell[9346, 282, 209, 5, 36, "Input",ExpressionUUID->"9d0a1263-6c64-42b0-b4e7-686e77c05a73"],
Cell[9558, 289, 482, 13, 36, "Input",ExpressionUUID->"02c0f5da-eb78-4600-be36-92573e264dcb"],
Cell[10043, 304, 179, 3, 36, "Input",ExpressionUUID->"58c2f45a-b492-47f1-b0f9-a7d75ef5be06"],
Cell[10225, 309, 205, 4, 36, "Input",ExpressionUUID->"78d17373-701f-418a-ac5e-bd19cefc412a"],
Cell[10433, 315, 219, 5, 36, "Input",ExpressionUUID->"d67d1415-1218-4585-b95e-3439cc3cc9e1"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10689, 325, 164, 3, 86, "Section",ExpressionUUID->"9ca5ec54-9cc0-4873-b655-534812f784e1"],
Cell[10856, 330, 176, 3, 45, "Text",ExpressionUUID->"d498ffda-51ea-4fd6-9882-f04a4288fc6c"],
Cell[11035, 335, 487, 13, 36, "Input",ExpressionUUID->"910565be-ff4f-4738-ba30-7b3c0ee75ada"],
Cell[11525, 350, 246, 6, 45, "Text",ExpressionUUID->"4d56c02e-69aa-47a0-8e36-e85ec8d06646"],
Cell[11774, 358, 268, 7, 36, "Input",ExpressionUUID->"cc6bc007-8a0d-4200-a709-1e8943d8f1ba"],
Cell[12045, 367, 321, 7, 45, "Text",ExpressionUUID->"9654b588-b976-46c4-a416-cf19930bfab3"],
Cell[12369, 376, 594, 17, 36, "Input",ExpressionUUID->"ed3e2321-43b5-4fea-8e21-2e1d2449a86a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13000, 398, 163, 3, 86, "Section",ExpressionUUID->"d2580dc1-1dd1-48a9-a40e-be340e60f09d"],
Cell[13166, 403, 605, 16, 36, "Input",ExpressionUUID->"aeb662c5-7536-41e7-8fd2-26ba4e80febb"],
Cell[13774, 421, 222, 5, 36, "Input",ExpressionUUID->"27652beb-04d8-4116-8719-79f70a33917a"],
Cell[13999, 428, 222, 5, 36, "Input",ExpressionUUID->"31b49bac-6699-43fc-9f69-efc948963ad9"],
Cell[14224, 435, 191, 3, 45, "Text",ExpressionUUID->"ffc3af29-461e-440c-bdcc-9ba04d24870a"],
Cell[14418, 440, 560, 18, 36, "Input",ExpressionUUID->"109f0d7f-88e4-40e3-beaf-74f8788e200b"],
Cell[14981, 460, 317, 9, 36, "Input",ExpressionUUID->"7ff268af-1579-421c-a8b4-76105893ead7"],
Cell[15301, 471, 239, 6, 36, "Input",ExpressionUUID->"2f3cfa8e-deae-4824-8b67-0acdee71453e"],
Cell[15543, 479, 488, 14, 36, "Input",ExpressionUUID->"26b1c59a-0a43-49ab-b6bf-954bee8aa647"],
Cell[16034, 495, 605, 17, 36, "Input",ExpressionUUID->"73121a2a-6ac5-4989-acc2-53bcc71954cb"]
}, Open  ]]
}, Open  ]]
}
]
*)