diff --git a/Lectures_my/NumMet/2016/Lecture10/fdm.nb b/Lectures_my/NumMet/2016/Lecture10/fdm.nb new file mode 100644 index 0000000..f9dbafc --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture10/fdm.nb @@ -0,0 +1,497 @@ +(* 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 *) diff --git a/Lectures_my/NumMet/2016/Lecture10/fixed-point.nb b/Lectures_my/NumMet/2016/Lecture10/fixed-point.nb new file mode 100644 index 0000000..771a9d4 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture10/fixed-point.nb @@ -0,0 +1,474 @@ +(* 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[ 15417, 465] +NotebookOptionsPosition[ 13779, 402] +NotebookOutlinePosition[ 14115, 417] +CellTagsIndexPosition[ 14072, 414] +WindowFrame->Normal*) + +(* Beginning of Notebook Content *) +Notebook[{ +Cell[BoxData[{ + RowBox[{ + RowBox[{ + RowBox[{"f", "[", + RowBox[{"x_", ",", " ", "y_"}], "]"}], " ", ":=", " ", + RowBox[{ + RowBox[{"-", "2.1"}], "y"}]}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"h", " ", "=", " ", "0.5"}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"x0", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"y0", " ", "=", " ", "1"}], ";"}], "\[IndentingNewLine]", + RowBox[{ + RowBox[{"yNextValues", " ", "=", " ", + RowBox[{"{", "}"}]}], ";"}]}], "Input", + CellChangeTimes->{{3.688109047059428*^9, 3.6881091615507298`*^9}, { + 3.688109196246677*^9, 3.688109226023087*^9}, {3.6881093786779633`*^9, + 3.688109383397826*^9}, {3.688110265128317*^9, 3.6881102693437967`*^9}, { + 3.688110342415762*^9, 3.6881103425034227`*^9}, {3.688110560365991*^9, + 3.688110561005869*^9}}], + +Cell["explicit Euler step with h 0.5", "Text", + CellChangeTimes->{{3.688109042375373*^9, 3.688109081215054*^9}, + 3.688109115526836*^9}], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{"x0", ",", " ", "y0"}], "]"}]}]}]}], "]"}]], "Input", + CellChangeTimes->{{3.688109023442931*^9, 3.688109034352015*^9}, + 3.6881090887992973`*^9, {3.688109133567192*^9, 3.6881092158464327`*^9}}], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{"-", "0.050000000000000044`"}], "}"}]], "Output", + CellChangeTimes->{{3.6881092163858*^9, 3.6881092272861643`*^9}, + 3.688109387188862*^9, 3.688110272361723*^9, {3.688110344468812*^9, + 3.688110350385796*^9}, 3.688110565279694*^9}] +}, Open ]], + +Cell["implicit Euler step", "Text", + CellChangeTimes->{{3.688109251600007*^9, 3.6881092747580843`*^9}}], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input", + CellChangeTimes->{{3.688109281639778*^9, 3.688109302422203*^9}}], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`"}], "}"}]], "Output",\ + + CellChangeTimes->{3.688109302822234*^9, 3.688109390123254*^9, + 3.688110272427383*^9, 3.688110350460581*^9, 3.6881105653414927`*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`", ",", + RowBox[{"-", "0.10512500000000014`"}]}], "}"}]], "Output", + CellChangeTimes->{3.688109310311212*^9, 3.688109390988352*^9, + 3.6881102724889593`*^9, 3.68811035053411*^9, 3.688110565396331*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`", ",", + RowBox[{"-", "0.10512500000000014`"}], ",", "1.11038125`"}], + "}"}]], "Output", + CellChangeTimes->{3.6881093159400263`*^9, 3.688109392546154*^9, + 3.688110272540285*^9, 3.6881103506078377`*^9, 3.688110565400127*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`", ",", + RowBox[{"-", "0.10512500000000014`"}], ",", "1.11038125`", ",", + RowBox[{"-", "0.16590031250000026`"}]}], "}"}]], "Output", + CellChangeTimes->{3.688109473510303*^9, 3.6881102726069508`*^9, + 3.688110350612899*^9, 3.6881105654640217`*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`", ",", + RowBox[{"-", "0.10512500000000014`"}], ",", "1.11038125`", ",", + RowBox[{"-", "0.16590031250000026`"}], ",", "1.1741953281250002`"}], + "}"}]], "Output", + CellChangeTimes->{3.6881095000395107`*^9, 3.688110272668214*^9, + 3.688110350675782*^9, 3.68811056553791*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`", ",", + RowBox[{"-", "0.10512500000000014`"}], ",", "1.11038125`", ",", + RowBox[{"-", "0.16590031250000026`"}], ",", "1.1741953281250002`", ",", + RowBox[{"-", "0.2329050945312503`"}]}], "}"}]], "Output", + CellChangeTimes->{3.688109503620624*^9, 3.688110272735523*^9, + 3.688110350738181*^9, 3.68811056561063*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`", ",", + RowBox[{"-", "0.10512500000000014`"}], ",", "1.11038125`", ",", + RowBox[{"-", "0.16590031250000026`"}], ",", "1.1741953281250002`", ",", + RowBox[{"-", "0.2329050945312503`"}], ",", "1.2445503492578127`"}], + "}"}]], "Output", + CellChangeTimes->{3.68810950655245*^9, 3.688110272805271*^9, + 3.688110350802301*^9, 3.6881105656837397`*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"AppendTo", "[", + RowBox[{"yNextValues", ",", " ", + RowBox[{"y0", " ", "+", " ", + RowBox[{"h", " ", "*", " ", + RowBox[{"f", "[", + RowBox[{ + RowBox[{"x0", " ", "+", " ", "h"}], ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", + RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}]}]}], "]"}]], "Input"], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"-", "0.050000000000000044`"}], ",", "1.0525`", ",", + RowBox[{"-", "0.10512500000000014`"}], ",", "1.11038125`", ",", + RowBox[{"-", "0.16590031250000026`"}], ",", "1.1741953281250002`", ",", + RowBox[{"-", "0.2329050945312503`"}], ",", "1.2445503492578127`", ",", + RowBox[{"-", "0.3067778667207035`"}]}], "}"}]], "Output", + CellChangeTimes->{3.688109510119643*^9, 3.688110272870085*^9, + 3.688110350875483*^9, 3.688110565688951*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"xValues", " ", "=", " ", + RowBox[{"{", " ", + RowBox[{ + "0.5", ",", " ", "0.5", ",", " ", "0.5", ",", " ", "0.5", ",", " ", "0.5", + ",", " ", "0.5", ",", " ", "0.5", ",", " ", "0.5", ",", " ", "0.5"}], + "}"}]}]], "Input", + CellChangeTimes->{{3.6881094113321342`*^9, 3.688109427557271*^9}, { + 3.688109529300817*^9, 3.688109540588752*^9}}], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + "0.5`", ",", "0.5`", ",", "0.5`", ",", "0.5`", ",", "0.5`", ",", "0.5`", + ",", "0.5`", ",", "0.5`", ",", "0.5`"}], "}"}]], "Output", + CellChangeTimes->{ + 3.688109429506762*^9, {3.688109526925215*^9, 3.688109541022224*^9}, + 3.6881102729397783`*^9, 3.688110350945903*^9, 3.68811056575429*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"pValues", " ", "=", " ", + RowBox[{"Table", "[", + RowBox[{ + RowBox[{"{", + RowBox[{"i", ",", " ", + RowBox[{"yNextValues", "[", + RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", " ", + RowBox[{"{", + RowBox[{"i", ",", " ", "1", ",", " ", "9"}], "}"}]}], "]"}]}]], "Input", + CellChangeTimes->{{3.688109435316019*^9, 3.688109477381198*^9}, { + 3.688109551381194*^9, 3.6881095636526012`*^9}}], + +Cell[BoxData[ + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{"1", ",", + RowBox[{"-", "0.050000000000000044`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"2", ",", "1.0525`"}], "}"}], ",", + RowBox[{"{", + RowBox[{"3", ",", + RowBox[{"-", "0.10512500000000014`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"4", ",", "1.11038125`"}], "}"}], ",", + RowBox[{"{", + RowBox[{"5", ",", + RowBox[{"-", "0.16590031250000026`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"6", ",", "1.1741953281250002`"}], "}"}], ",", + RowBox[{"{", + RowBox[{"7", ",", + RowBox[{"-", "0.2329050945312503`"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"8", ",", "1.2445503492578127`"}], "}"}], ",", + RowBox[{"{", + RowBox[{"9", ",", + RowBox[{"-", "0.3067778667207035`"}]}], "}"}]}], "}"}]], "Output", + CellChangeTimes->{{3.688109451265951*^9, 3.68810947771984*^9}, { + 3.688109560290724*^9, 3.688109563868795*^9}, 3.688110272995986*^9, + 3.6881103510253687`*^9, 3.6881105658340073`*^9}] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Show", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Plot", "[", "\[IndentingNewLine]", + RowBox[{ + RowBox[{"Exp", "[", + RowBox[{"-", "0.5"}], "]"}], ",", "\[IndentingNewLine]", + RowBox[{"{", + RowBox[{"i", ",", " ", "1", ",", " ", "9"}], "}"}]}], + "\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]", + RowBox[{ + "ListPlot", "[", "\[IndentingNewLine]", "pValues", "\[IndentingNewLine]", + "]"}]}], "\[IndentingNewLine]", "]"}]], "Input", + CellChangeTimes->{{3.6881094064823713`*^9, 3.688109464029291*^9}, { + 3.6881095680930147`*^9, 3.6881096166203403`*^9}}], + +Cell[BoxData[ + GraphicsBox[{{{}, {}, + {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ + 1.], LineBox[CompressedData[" +1:eJxTTMoPSmViYGAwAWIQvaTxsjYDwwd7rmn5fzZlPLaX2NGxfIbjRzi/552t +smHzJzifMTVIo8bmK5yvu+XMxYnHv8P5axzz3j50/QXnb7+xU8zn1x84v9Ij ++0lDw384/0nC4/vpcowOMP7WPUvadMqY4PyHLDPrr21jhvMvBybV9wqxwvkK +3XEXNcPZ4PzpyXkxyVPZ4fwrDsHvVr7jgPNzuLIvXTXggvNP3wjfy9fHDecH +u0U8fPKEB843W9AdvkCFD85/8G/ylYx6fjj/85Lngd2JAnD+zTPv97wQEITz +Q440HhHZi+D/MdQt3RwrBOdzKDmvvMEtDOc7dhS1FaxD8Pn3aTNc9xWB8x3K +/P8e/oXgz1/Ss1h5jiicnxhqcpbTWwzO33E1dfK+bwi+zZkNTM+micP5yq5p +Xye7SMD5FiamE7c8Q/Affeo9ubhVEs5/bxzYKmkkBedrvtZ9xXgZwT+uUnqp +oFEazp905WpykIoMnP8p+MBEpiMI/tZQUWfLQlk4f9nZ5rmPhOTg/F/HY6qF +DiH4jyfu/yqQJw/n/30xl8eRXQHOt/QSfGATh+D3pAlyFK1H8FdmhdXx/Ufw +f6wPMMsLVYTz4zZcd2JfhuCLLlg3j/kygn+mv/XX//8IPgDh0JgR + "]]}}, {{}, {{}, + {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ + 0.012833333333333334`], AbsoluteThickness[1.6], + PointBox[{{1., -0.050000000000000044`}, {2., 1.0525}, { + 3., -0.10512500000000014`}, {4., 1.11038125}, { + 5., -0.16590031250000026`}, {6., 1.1741953281250002`}, { + 7., -0.2329050945312503}, {8., 1.2445503492578127`}, { + 9., -0.3067778667207035}}]}, {}}, {}}}, + AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], + Axes->{True, True}, + AxesLabel->{None, None}, + AxesOrigin->{2., 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->{{1, 9}, {0., 1.2130613194252668`}}, + PlotRangeClipping->True, + PlotRangePadding->{{ + Scaled[0.02], + Scaled[0.02]}, { + Scaled[0.05], + Scaled[0.05]}}, + Ticks->{Automatic, Automatic}]], "Output", + CellChangeTimes->{{3.68810946439671*^9, 3.6881094785983686`*^9}, { + 3.6881095647315187`*^9, 3.6881095801719913`*^9}, 3.688109617483754*^9, + 3.6881102730755672`*^9, 3.68811035114158*^9, 3.688110565938201*^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[558, 20, 843, 20, 121, "Input"], +Cell[1404, 42, 138, 2, 31, "Text"], +Cell[CellGroupData[{ +Cell[1567, 48, 376, 8, 32, "Input"], +Cell[1946, 58, 274, 5, 32, "Output"] +}, Open ]], +Cell[2235, 66, 103, 1, 31, "Text"], +Cell[CellGroupData[{ +Cell[2363, 71, 434, 11, 32, "Input"], +Cell[2800, 84, 255, 6, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[3092, 95, 368, 10, 32, "Input"], +Cell[3463, 107, 300, 6, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[3800, 118, 368, 10, 32, "Input"], +Cell[4171, 130, 326, 7, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[4534, 142, 368, 10, 32, "Input"], +Cell[4905, 154, 349, 7, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[5291, 166, 368, 10, 32, "Input"], +Cell[5662, 178, 377, 8, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[6076, 191, 368, 10, 32, "Input"], +Cell[6447, 203, 419, 8, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[6903, 216, 368, 10, 32, "Input"], +Cell[7274, 228, 452, 9, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[7763, 242, 368, 10, 32, "Input"], +Cell[8134, 254, 495, 9, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[8666, 268, 375, 8, 32, "Input"], +Cell[9044, 278, 339, 7, 32, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[9420, 290, 440, 11, 32, "Input"], +Cell[9863, 303, 1024, 28, 55, "Output"] +}, Open ]], +Cell[CellGroupData[{ +Cell[10924, 336, 628, 14, 209, "Input"], +Cell[11555, 352, 2208, 47, 241, "Output"] +}, Open ]] +} +] +*) + +(* End of internal cache information *) diff --git a/Lectures_my/NumMet/2016/Lecture10/images/BG_lower.png b/Lectures_my/NumMet/2016/Lecture10/images/BG_lower.png new file mode 100755 index 0000000..073034a --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture10/images/BG_lower.png Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture10/images/Shooting_method_error.pdf b/Lectures_my/NumMet/2016/Lecture10/images/Shooting_method_error.pdf new file mode 100644 index 0000000..c965752 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture10/images/Shooting_method_error.pdf Binary files differ diff --git a/Lectures_my/NumMet/2016/Lecture10/images/Shooting_method_trajectories.pdf b/Lectures_my/NumMet/2016/Lecture10/images/Shooting_method_trajectories.pdf new file mode 100644 index 0000000..dec4329 --- /dev/null +++ b/Lectures_my/NumMet/2016/Lecture10/images/Shooting_method_trajectories.pdf Binary files differ