Newer
Older
Lecture_repo / Lectures_my / NumMet / 2016 / Lecture10 / fdm.nb
@Danny van Dyk Danny van Dyk on 14 Nov 2016 18 KB Add missing files for lecture 10
(* Content-type: application/vnd.wolfram.mathematica *)

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

(* CreatedBy='Mathematica 10.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     18407,        488]
NotebookOptionsPosition[     17562,        454]
NotebookOutlinePosition[     17898,        469]
CellTagsIndexPosition[     17855,        466]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"h", " ", "=", " ", "0.1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"A", " ", "=", " ", 
   RowBox[{
    RowBox[{"1", "/", 
     RowBox[{"h", "^", "2"}]}], " ", 
    RowBox[{"{", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "2"}], ",", " ", 
        RowBox[{"+", "1"}], ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", 
        " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0"}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"+", "1"}], ",", " ", 
        RowBox[{"-", "2"}], ",", " ", 
        RowBox[{"+", "1"}], ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", 
        " ", "0", ",", " ", "0", ",", " ", "0"}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"0", ",", " ", 
        RowBox[{"+", "1"}], ",", " ", 
        RowBox[{"-", "2"}], ",", " ", 
        RowBox[{"+", "1"}], ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", 
        " ", "0", ",", " ", "0"}], "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"0", ",", " ", "0", ",", " ", 
        RowBox[{"+", "1"}], ",", " ", 
        RowBox[{"-", "2"}], ",", " ", 
        RowBox[{"+", "1"}], ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", 
        " ", "0"}], "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"0", ",", " ", "0", ",", " ", "0", ",", 
        RowBox[{"+", "1"}], ",", " ", 
        RowBox[{"-", "2"}], ",", " ", 
        RowBox[{"+", "1"}], ",", " ", "0", ",", " ", "0", ",", " ", "0"}], 
       "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"0", ",", " ", "0", ",", " ", "0", ",", "0", ",", " ", 
        RowBox[{"+", "1"}], ",", " ", 
        RowBox[{"-", "2"}], ",", " ", 
        RowBox[{"+", "1"}], ",", " ", "0", ",", " ", "0"}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
       "0", ",", " ", "0", ",", " ", "0", ",", "0", ",", " ", "0", ",", " ", 
        RowBox[{"+", "1"}], ",", " ", 
        RowBox[{"-", "2"}], ",", " ", 
        RowBox[{"+", "1"}], ",", " ", "0"}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
       "0", ",", " ", "0", ",", " ", "0", ",", "0", ",", " ", "0", ",", " ", 
        "0", ",", 
        RowBox[{"+", "1"}], ",", "  ", 
        RowBox[{"-", "2"}], ",", 
        RowBox[{"+", "1"}]}], "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
       "0", ",", " ", "0", ",", " ", "0", ",", "0", ",", " ", "0", ",", " ", 
        "0", ",", "0", ",", 
        RowBox[{"+", "1"}], ",", 
        RowBox[{"-", "2"}]}], "}"}]}], "\[IndentingNewLine]", "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"b", " ", "=", " ", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"4", " ", "+", 
       RowBox[{"1", " ", "/", " ", 
        RowBox[{"h", "^", "2"}]}]}], ",", " ", "4", ",", " ", "4", ",", " ", 
      "4", ",", " ", "4", ",", " ", "4", ",", " ", "4", ",", " ", "4", ",", 
      " ", 
      RowBox[{"4", " ", "-", " ", 
       RowBox[{"1", "/", 
        RowBox[{"h", "^", "2"}]}]}]}], "}"}]}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"y", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
    "y1", ",", " ", "y2", ",", " ", "y3", ",", " ", "y4", ",", " ", "y5", ",",
      " ", "y6", ",", " ", "y7", ",", " ", "y8", ",", " ", "y9"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "A", "]"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "b", "]"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "y", "]"}]}], "Input",
 CellChangeTimes->{{3.6881134166143703`*^9, 3.6881135429801407`*^9}, {
  3.688113614373023*^9, 3.688113701931562*^9}, {3.6881137320757027`*^9, 
  3.688113732475355*^9}, {3.688113846083762*^9, 3.688113902986618*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`", "0.`", 
      "0.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`", "0.`", 
      "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`", "0.`", 
      "0.`", "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`", "0.`", 
      "0.`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`", "0.`", 
      "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`", "0.`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}], "99.99999999999999`"},
     {"0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "0.`", "99.99999999999999`", 
      RowBox[{"-", "199.99999999999997`"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.688113521414591*^9, 3.688113621437024*^9, {3.688113674139687*^9, 
   3.6881137023408117`*^9}, 3.688113732904488*^9, {3.688113873518072*^9, 
   3.688113903424461*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {"103.99999999999999`"},
      {"4"},
      {"4"},
      {"4"},
      {"4"},
      {"4"},
      {"4"},
      {"4"},
      {
       RowBox[{"-", "95.99999999999999`"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.688113521414591*^9, 3.688113621437024*^9, {3.688113674139687*^9, 
   3.6881137023408117`*^9}, 3.688113732904488*^9, {3.688113873518072*^9, 
   3.688113903426567*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {"y1"},
      {"y2"},
      {"y3"},
      {"y4"},
      {"y5"},
      {"y6"},
      {"y7"},
      {"y8"},
      {"y9"}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
        "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.688113521414591*^9, 3.688113621437024*^9, {3.688113674139687*^9, 
   3.6881137023408117`*^9}, 3.688113732904488*^9, {3.688113873518072*^9, 
   3.688113903427741*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"solution", " ", "=", " ", 
  RowBox[{"Solve", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"A", ".", "y"}], " ", "\[Equal]", " ", "b"}], ",", 
    "\[IndentingNewLine]", "y"}], "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.688113680416813*^9, 3.688113704899498*^9}, {
   3.688113751098772*^9, 3.688113757674786*^9}, 3.68811379211481*^9, {
   3.688113965625009*^9, 3.688113967024735*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"y1", "\[Rule]", 
     RowBox[{"-", "0.9800000000000001`"}]}], ",", 
    RowBox[{"y2", "\[Rule]", 
     RowBox[{"-", "0.9200000000000002`"}]}], ",", 
    RowBox[{"y3", "\[Rule]", 
     RowBox[{"-", "0.8200000000000004`"}]}], ",", 
    RowBox[{"y4", "\[Rule]", 
     RowBox[{"-", "0.6800000000000004`"}]}], ",", 
    RowBox[{"y5", "\[Rule]", 
     RowBox[{"-", "0.5000000000000003`"}]}], ",", 
    RowBox[{"y6", "\[Rule]", 
     RowBox[{"-", "0.28000000000000025`"}]}], ",", 
    RowBox[{"y7", "\[Rule]", 
     RowBox[{"-", "0.020000000000000177`"}]}], ",", 
    RowBox[{"y8", "\[Rule]", "0.28`"}], ",", 
    RowBox[{"y9", "\[Rule]", "0.6199999999999999`"}]}], "}"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.688113758061264*^9, 3.688113792694109*^9, {3.688113881155346*^9, 
   3.688113906365438*^9}, 3.6881139677847967`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"y", " ", "//.", " ", 
  RowBox[{"%", "[", 
   RowBox[{"[", "1", "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6881137940686827`*^9, 3.688113799594346*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.9800000000000001`"}], ",", 
   RowBox[{"-", "0.9200000000000002`"}], ",", 
   RowBox[{"-", "0.8200000000000004`"}], ",", 
   RowBox[{"-", "0.6800000000000004`"}], ",", 
   RowBox[{"-", "0.5000000000000003`"}], ",", 
   RowBox[{"-", "0.28000000000000025`"}], ",", 
   RowBox[{"-", "0.020000000000000177`"}], ",", "0.28`", ",", 
   "0.6199999999999999`"}], "}"}]], "Output",
 CellChangeTimes->{{3.688113883193458*^9, 3.6881139083347187`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"points", " ", "=", " ", 
  RowBox[{"Table", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"0", " ", "+", " ", 
       RowBox[{"h", " ", "*", " ", "i"}]}], ",", " ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"y", " ", "//.", " ", 
         RowBox[{"solution", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "[", 
       RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"i", ",", " ", "1", ",", " ", "9"}], "}"}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.6881139220814857`*^9, 3.6881140228081217`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0.1`", ",", 
     RowBox[{"-", "0.9800000000000001`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.2`", ",", 
     RowBox[{"-", "0.9200000000000002`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.30000000000000004`", ",", 
     RowBox[{"-", "0.8200000000000004`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.4`", ",", 
     RowBox[{"-", "0.6800000000000004`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.5`", ",", 
     RowBox[{"-", "0.5000000000000003`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.6000000000000001`", ",", 
     RowBox[{"-", "0.28000000000000025`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.7000000000000001`", ",", 
     RowBox[{"-", "0.020000000000000177`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.8`", ",", "0.28`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.9`", ",", "0.6199999999999999`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.688113947360079*^9, 3.688113984829784*^9, 
  3.688114023090886*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Plot", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"2", 
       RowBox[{"x", "^", "2"}]}], " ", "-", " ", "1"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"x", ",", " ", "0", ",", " ", "1"}], "}"}]}], 
    "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{
   "ListPlot", "[", "\[IndentingNewLine]", "points", "\[IndentingNewLine]", 
    "]"}]}], "\[IndentingNewLine]", "]"}]], "Input",
 CellChangeTimes->{{3.688113987402186*^9, 3.6881140297520638`*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
     1.], LineBox[CompressedData["
1:eJwV1gk0lGsYB3DRIhO35aamuCKmVJZESenRviFblluJQombVEpuWSNliUwq
DdGEQkYokry2myUyBhlhzIzljmlm3rpZKtF9+875znd+5zvv8/6/8zzvOZ/m
MT87T3k5OblUcv967vMUtVYOH9zy9eevC6N/8x2P9VGNgZ32I+WXVSkd8UXU
HaBx5JvPFPHLimDfizRrGDIbO/uDWH/j75mpVAegmX5hfSOWWzCc9FPDBXgO
n7TGiFslr8Pdaa6gkih9+5nY9nhDqXuyO7yQiHMkxGbsdQW3qB7QeUKEhoiX
m6dl1zC84Iz8kBKfeE6O4oMRDW8IK++P4RKPqJ5L1mH6gFKyYBubuDeiN86R
dhooN/sM6oiPf2VW1zicgcjsXrvXxJ36LlsOJ/vDxb7u3ELiyvs1JnHU81Dz
lPszhdiYE/hM51AAqO/o/B5H/ERRX6+CcQE+fetYHUp8K+CuNta4BFM1bQbH
iGfmWaVfcw+CjVyOvB1xkFBebRnzb5BQOJStv/Y/8NcCW1owlLeykRoxN1Lr
5vCJEOAeZ5+aTUzj+cyocggF34Xs/aNTGKX8rhnmkhwGk+9bWI3EEWExgTeo
V6G6jf3Th7hCzeuFY9lVSHdv/dOO+GuJxYjWoUgYU+F83EDsKxv1e8WIAkp9
G0uO2OGQ2wmJxnV4X8RF0ZMYxY+ZZZZWXocA8w+NJ4nrExf2X3W/AXFD3ZQ9
xOb1ja7qzBjQudNnMJ2YZmLiaE2Lh3F/0c+zPzByY/9GX1IXD87j4jJL4vs+
4tahEzdhLUOaq0Os2NqShRwSwGrHF4OOCYwuz4rUck5OBB0VucpVxB7nMfU6
lQ7u0RpyBd8wWqdwK6jwDB2mv5FrDCSWv2Xc3V1HhyJ2YJYF8QPWJYb+hdtg
ZKV3p/krRh/EChrtnGQoEWY79YxjdMB9sbZG7D0QhvQ7NI9i9Mensqt7+u/B
pvi1e64RS4KPDPqbpUCat6q1BfENxsOsWlEKRF97lMgawai2c43uqZ0MEB7l
jId/IfNnvVX/+VQalISft5j8RL6ntz+ed/AB1LobHHlM3OkThWc9fQC9ymGx
dsRltZ7OT/XSYW2zSOsRxshZmEOV60kHN1GD/mYZmRd1E0am6UPgJXCq4SOp
R9+b8enzI3DOmeymD2A0YSqYyFiSCdvGAsS6xLLeQEf77Zlw9vBn3Yp+jNpW
PKa8oGdC9sXPmwaEGKW9mnnh7w1Z4DGsmkcTkHntr9k380o26JSNlV/tJf0x
2jKyVDEH6uYqZWR2YGTf2WHdbJgDMZmDLcuId13+60mwSw7ccag1uN+O0eo3
DFdBTg6ETYsJj2vDaMxloj7LKheeJJufP9mKUWxYaerapDywyN8QPd6EUQnb
aPfuP1iQcCvkbFANmadhyZ41u1hg6Z/B51Vj5CWfvW/eaRbk97OvbSNWNV5q
3f2aBScjnMNmVWH06m1YTol+AeTmK2lGV2C01ODUs/3vCuBs8ff5Hi9Jf0fN
0DnlQphxZXZVVD7JK7/YwkuzEO6we204TzEqVBmtdDYphJYSRFUnvrGCVW1+
pBBMtuY5FuSS/rksfzPrKXkPK680P8boXjml+b5lEeyS5F4WPiT9iej5UBNT
DBvkhTsvJmM0l09bcyKjGIRFN67dvY1Rw2b/K0olxaD9kNpXSif1xmYssxUW
w2ljccvYLZLX28CTZ/ocFhk6lp64iRHPOhx/HXwOlREuj9SjyXwvWTld36IE
TjOUh7sCMfJ7dk7v7shLWJe0WofuRPq7WdDWoFgGDYGyvDhHcl7rrYMm1MoA
4t6ejDqI0ck+3TrXnWXQtK39coA96acy300nuQw8+Rqq+w5gZONtSS/c8Aru
DT5Qbt6FkaGmzkRTUDkoLt4+18wEo/8S3tdPm4ZAXl3Vo+A3cr6r3qCWzdUg
vDnp6FEkQ96X9iq7SGrBjpIq2LBPhqpmG/3M+6sOXja4Spd2S1EWs3ZqiU4j
6Cf6Lfd3kCIL1dcP7fObYNWilCmvBgnSq+spfbOqBToZPi8ZxhLUbmEveyJj
w86H62etjfyIfJV8OB2GHNjdqH7hWLsYveU6vVaJbwO/zStlPjPFyH6Xs2Bg
oB28XxjaVugNo/XpMU7p2u/h9rvu5nJHEeJPJbWfDOkENgTL+x39Fw3ubWGq
1nBh5QGa02GfIVRwr+9Sq+YHyFLYLTvnMIgyoxWu6wZ1w41C/I+C8QAaPf6n
2OtVD4wPC4wF1H5E2RFl92wRD9rDt09WdgnQ5bmDcwSH+6Apy3NtSAYfzalY
LddpxYdk+kRq/j88ZH7hwGTNdz4E1j5O2zTWjRiPYpnLGQIwXeC+JGxzF3I9
aNw8e78QHESVc+ir36PnHZ5JFWNCCOpPUJFQ2tDGpgL5oeR+sDGZWqGuyEbL
dnqNJu0YgIE6xfx0qwZkYmySWDw0AIayLtWv9dWo77+4BmbkILANyJHwKkWS
dbaRVKMh4FN8l7a4MBHto554WtsQhMb1tYQOJEKtdgDnTNi/wK/wO9SztRhu
tncct9MWQfrMvBDa8yrA9pWJ8rUi4Jq0spkHG6Dw4MLtG/2H4YxNlKNomA3M
5ohU4XwxnO8qoEu/t8N43eG/51eLwfSDh1ZPSBfwE9Ho3NMfwUbTqblnDg/G
WjdRw1UkUOmuEdpuJIAzdK3KS+USWJ9tZuvI6IdFzVpzB1ylYB0rC9YTD4J9
R94PtSkpvAi8IpGYiuC/oxm+PTkyMNu0xaC7TgzIjRWpYonBsk6JrWkqhbIW
3/08awz7Ap+v0NpO1m1ZNS/fFgPN2Onycmsp5KllMqycMHg1hajpekghhZtS
FHsMw8oQT1OLBCkE2EQJlC5hSJGG81JFUlhtccR8VhYGU9mMAdUUGdBYSxQ6
H2MI/Z74YX2mDLT+4NZn5WLImrfwnVOBDBb/sHPY9QxDkJtcHuONDGaW7vGJ
LMfgEeqz3eiLDAQGxvcU2kiehiR+LMndk/bZta0Dw/yIhKxykourzNJmcjHw
dl/1lpJc7z7qFmzjYfjx2WbIluQqz9aoCxsm9euDX9FIrlLV3lhbCYYtXeu8
Xcm+RZEpdpqY1BnpW3CH7JtzfCGvcgRDmqW6G4Xsm8nhMBPGMdRH50/bTeqm
b03wdvuOobFlfXoEqcsosDIwnMSQp/liUxVZd0eDMkr+08AjeE37r+f/qB5u
WA==
      "]]}}, {{}, {{}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.012833333333333334`], AbsoluteThickness[1.6], 
      PointBox[{{0.1, -0.9800000000000001}, {0.2, -0.9200000000000002}, {
       0.30000000000000004`, -0.8200000000000004}, {
       0.4, -0.6800000000000004}, {0.5, -0.5000000000000003}, {
       0.6000000000000001, -0.28000000000000025`}, {
       0.7000000000000001, -0.020000000000000177`}, {0.8, 0.28}, {0.9, 
       0.6199999999999999}}]}, {}}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None},
  PlotRange->{{0, 1}, {-0.9999999999999991, 0.9999999183673478}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.688114030148334*^9}]
}, Open  ]]
},
WindowSize->{862, 1179},
WindowMargins->{{0, Automatic}, {Automatic, 19}},
FrontEndVersion->"10.3 for Linux x86 (64-bit) (October 9, 2015)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 3874, 94, 385, "Input"],
Cell[4457, 118, 1996, 43, 179, "Output"],
Cell[6456, 163, 967, 31, 179, "Output"],
Cell[7426, 196, 918, 30, 179, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8381, 231, 451, 9, 99, "Input"],
Cell[8835, 242, 897, 23, 55, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9769, 270, 181, 4, 32, "Input"],
Cell[9953, 276, 499, 11, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10489, 292, 667, 17, 99, "Input"],
Cell[11159, 311, 1044, 29, 55, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12240, 345, 609, 15, 209, "Input"],
Cell[12852, 362, 4694, 89, 224, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)