(* 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[ 34265, 816]
NotebookOptionsPosition[ 33004, 768]
NotebookOutlinePosition[ 33341, 783]
CellTagsIndexPosition[ 33298, 780]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", " ", "x"}], "]"}]], "Input",
CellChangeTimes->{{3.6879511125883017`*^9, 3.687951128551201*^9}}],
Cell[BoxData[
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", "x"}], "]"}]], "Output",
CellChangeTimes->{
3.6879511299603863`*^9, {3.687952885427848*^9, 3.687952897510394*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"dBesselJ1", " ", "=", " ",
RowBox[{"D", "[",
RowBox[{
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", " ", "x"}], "]"}], ",", " ", "x"}], "]"}]}]], "Input",
CellChangeTimes->{{3.687951319931106*^9, 3.687951330309391*^9}}],
Cell[BoxData[
RowBox[{
FractionBox["1", "2"], " ",
RowBox[{"(",
RowBox[{
RowBox[{"BesselJ", "[",
RowBox[{"0", ",", "x"}], "]"}], "-",
RowBox[{"BesselJ", "[",
RowBox[{"2", ",", "x"}], "]"}]}], ")"}]}]], "Output",
CellChangeTimes->{
3.6879513307406588`*^9, {3.687952882162366*^9, 3.68795289760646*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", " ", "x"}], "]"}], ",", " ",
RowBox[{"{",
RowBox[{"x", ",", " ", "2"}], "}"}]}], "]"}], " ", "//.", " ",
RowBox[{"x", " ", "\[Rule]", " ", "1"}]}], " ", "//", " ", "N"}], ")"}],
" ", "+", " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", " ", "x"}], "]"}], ",", " ",
RowBox[{"{",
RowBox[{"x", ",", " ", "1"}], "}"}]}], "]"}], " ", "//.", " ",
RowBox[{"x", " ", "\[Rule]", " ", "1"}]}], " ", "//", " ", "N"}],
")"}]}]], "Input",
CellChangeTimes->{{3.687951204980104*^9, 3.6879512349344*^9}}],
Cell[BoxData["0.`"], "Output",
CellChangeTimes->{{3.6879512129131317`*^9, 3.6879512398981333`*^9}, {
3.687952879441258*^9, 3.687952897673266*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"Table", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", " ",
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", " ", "x"}], "]"}], ",", " ", "dBesselJ1"}], "}"}],
",", "\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"x", ",", " ", "1", ",", " ", "2", ",", " ", "0.05"}], "}"}]}],
"\[IndentingNewLine]", "]"}], " ", "//", " ", "TableForm"}]], "Input",
CellChangeTimes->{{3.687951256849247*^9, 3.687951349621264*^9}}],
Cell[BoxData[
TagBox[GridBox[{
{"1.`", "0.44005058574493355`", "0.32514710081303305`"},
{"1.05`", "0.4558967897777311`", "0.30860813759792577`"},
{"1.1`", "0.47090239486629304`", "0.29152893228542653`"},
{"1.15`", "0.48504126615351717`", "0.2739447494580766`"},
{"1.2`", "0.4982890575672156`", "0.25589186295834976`"},
{"1.25`", "0.5106232603198806`", "0.23740747701538092`"},
{"1.3`", "0.5220232474146603`", "0.21852964539638578`"},
{"1.35`", "0.5324703140632868`", "0.19929718876387187`"},
{"1.4`", "0.5419477139308545`", "0.17974961042367837`"},
{"1.45`", "0.5504406911316964`", "0.15992701065240097`"},
{"1.5`", "0.5579365079100995`", "0.1398699997958517`"},
{"1.55`", "0.5644244679492657`", "0.11961961033284683`"},
{"1.6`", "0.5698959352616804`", "0.09921720810083043`"},
{"1.65`", "0.5743443486239402`", "0.07870440288159655`"},
{"1.7000000000000002`", "0.5777652315290233`", "0.058122958546684034`"},
{"1.75`", "0.5801561976389924`", "0.03751470296286935`"},
{"1.8`", "0.5815169517311652`", "0.01692143785857772`"},
{"1.85`", "0.5818492861408044`",
RowBox[{"-", "0.0036151511480226883`"}]},
{"1.9`", "0.5811570727134342`",
RowBox[{"-", "0.024053584159000857`"}]},
{"1.9500000000000002`", "0.5794462502898536`",
RowBox[{"-", "0.044352676045445516`"}]},
{"2.`", "0.5767248077568736`",
RowBox[{"-", "0.06447162473720106`"}]}
},
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
"RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[2.0999999999999996`]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}],
Function[BoxForm`e$,
TableForm[BoxForm`e$]]]], "Output",
CellChangeTimes->{{3.687951293721142*^9, 3.687951349903165*^9},
3.687952897700285*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", " ", "x"}], "]"}], ",", " ", "dBesselJ1"}], "}"}],
",", "\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"x", ",", " ", "1", ",", " ", "2"}], "}"}], ",",
"\[IndentingNewLine]",
RowBox[{"PlotStyle", " ", "\[Rule]", " ",
RowBox[{"{",
RowBox[{"Blue", ",", " ", "Red"}], "}"}]}]}], "\[IndentingNewLine]",
"]"}]], "Input",
CellChangeTimes->{{3.687951491363956*^9, 3.687951518652185*^9}}],
Cell[BoxData[
GraphicsBox[{{}, {},
{RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.],
LineBox[CompressedData["
1:eJwV13k4VO0bB3AS4UURSZssCW8pGkLLLXkRUoQspUL2aJGihCQlS0kilDap
pJTQppB+hDIiFDInhFR2kznPnN89f831uc65nnme77PdR8Uj2H7vNCEhIRFh
ISHBL2t/uKiQ0BCo/lMyXqPVDv5LJrI3CA+BvUalZuC/7XD9ywGDmGlDsPnQ
hLHM8naQ+c/PR0x0CO661NJ2K9uhT8m5RkpyCMY8eOMtBu2QWamfqCQ/BIts
Ol52bGyHaYrD8qs0h2BGADe5YWc7sF/5qHtvHQKvhvRe+wvtcM1GMjb5zhBo
vJl0Xc20w6k0I+MR7hBYxw6Hpnl2AKv6wbqXlsPg3v7FoqOyAxTWtJzcdG0Y
krSXa3ZrdYK0X+KCWUPDsH+LShv7dCesbuVOZliMwKqfqiVC/Z3w74q7Kuey
RyBHbOjUR/NvYLCh99rTXyOg9Xn676nMb3AwtXXXTbNRcDpdblY5/g2qbmwP
eX9pFKKkr8qkqHUBlT6/xblnFJoW/Wdua9UFug2px+6tGQNXneqGGYe64HLF
lRD9hDFwf9IxXnu5C3QmHB/GfB+D3QrF663LuuC8rJOhkdE42GmKWxVRXSDR
Rylonh2HigTpMn8RDjxh1PW/fBuHjMKAjrglHAhpHM2rWjkBXWHZSfJmHPDp
PeIecmYCXPW8zd56ckCdMvGI75yAULaDSXA0B7jpac/bdCZBOllrOTuTA+T8
dYfB2ElwHRWOEirhgJKZ0UbDlknQcvd3/+cjB5hOhaifGlwwGNLNc+3hgOKL
O/w9EVxgK/o2qdEcMLd/x5b9wIWUHLOrsrMoyNW+0ntqyV9o635vqq9OAXvg
8aZdx/6Cfsz245EGFMz2dhY+VPsXdsysm7/GnILTqk+mNqlOgfu+pswWRwpm
Xbli6XpoCpbrPpGQ9qTA6++7bnbFFLRvKJUpDaJgICOvQkqJB/VZNk+GwyiI
G58zHuXHA98PhZdkYyiYivnR2OHPA7HEUAsDdJBc6SPjQB68kHzV5Yp2WOES
OBbEg69qH2tuohf7ZX73CeGBMFcslHWKguftyo22kTyI/FtmYBdLwe9KzYcL
LvEgI0rD5WgcBY4pxn6lb3iw+A1ZHHGOgqJtO5fxy/H/3U8+yEDLK0T92VjJ
g3nT7qsXo5suVx3+WMUDCfe+xt9o++ytJ3vf88AwYfaKXQkU2OZ5Z8o38yAL
lveuSaTAsiylfn8/D64+itzSk0SB8cCArpYcDcyq+JncCxQYvTZP8Z1Nw9uW
ngaZFAoMU2+M3JGn4fbX5xFL0PrrXZ9oKNKwSJWdYY9emVLDUl9AQ9eph1/v
ozWM81Yra9CQVG982vUiBXLxe9cpGNMQPu99aF4qBbK7yrMd1tCgHuBx5yV6
Fmsh/+JaGgZMDlc2oKU7m8rkgAa9O5JlXLS4npnJLDMalqWOv7S4RAG/TXXj
P7bY3lAzh4Me0OqyFPagQV5sTrXYZczzk2mxmycN65fZ9iihhyNuqxV70WAl
7TmwDM1l+/P9fWj4a6GWb48WCx9/0hRIg+M1RjobrfpeUjnvCA0vLA2tl6dT
4OqvP26bSEP/AVb+mgwK3OXTPe4m0dAcolVljd5TNvVR5DwNt4pJlRvaV+7N
vdIUGqq2mkUeQ4c+t9qtmk7DLzebhmfoFMndteM3aADtHYO6VyiovRd/PauU
ht5OdelZmbg+4mDs2zMaOFdTNy5A3/MaM1d7QUNA0bS9muj4Re6Dd1/R8Gr1
F18TtFWK7uqSChqEJcNCgtE1Ya117Doa9p07m1mDLnVKVFb4QEPeov7oJnTe
KtODzh9p6Dlgue0b+syv+4rf2DTkFouVjKEt90R6DH6mIWjbqZOLsij4n6UG
V4yD74tkaAWhi5d8tbKiaAhrZ14fQedOO5+d+J2G1CmfDdHo0y//msr30hAY
uFY0FW2+sj5R5ScNrCHu+VJ0lWKI2tpxGhYsbCkl6KIxzdDICRqy3N0Z0WwK
brE7qismMZ/TfStl0KfOWQRtmsI8R0UdF6HNmHnPnBga8k18r69DV/aW2x6Q
IDC7grI6ii5Njyj4KEmA8W9rjkTnWxnJ6EgRYDkss4xDpxU8qh+QIRAb9oO6
hPYNvWbjJU/AS6FkqBC9U9Mtv0IBn1ddLC5F23+ZI6WiSOBmR73na/Sa9Um1
HUoEfNtdQ+rQMmIRVtuVCWQYt+/rQU8rNbz3dDGBv1vygn6iJ/3GJORVCTjz
G52G0V31ATUN6gSsys06aPSTS66WVtoElJZerZC7SoHzEkPzdfoERtWt/xih
bVpGb2caENi8abb3evSGsw9Fp1YTyDy19q0pWvuXRlWxMfav+7GBNZoUKZit
NCFQ1LhqzBU94s2+mbSBwHd9iVx3dO/cRJFfpgR4tAl4oBuOT6+8+x8Btd2t
S/3QN8xGN6hZE9BVu88+jL48UXA92obA8c3pZUfRCXn+wl2bCYzr9qccQx+W
ot5kbcXxtRQPR6H9y7IX8+wIXPyz6ngM2n2/S5TLNuzvfoXfsWjLpgaY40TA
cafYmXj02tMJ10K2E4hSmfMoAa1raMk0OhMIZuJfJ6HnZb0uS3YjoPLiVmoK
eqbtsUW/dxA4N7LJNRU9XWj1CRt3AiKTDjPS0L88C9ZJ7CEwMj9bMQMdpXXT
eI8HgUeHMkKuoOX/XF79zJNAZOW3Z5novKIEluxeAvkG+39kCfoTHq3r503g
FteKZAvygVCdch/Mi7WfdxXtJRrwr5IfgUDFLs41NPf9Ls0D/gRqa688yhHk
c95hSU0AAfkTOT7X0YudNqmq7MP15TosdgNdNH+9clgQro+ExCSBLTl6C9jB
BAxND/MF/pq7VEnrAIG83Dynm+igwAVzog8SGOxdelngaXqys9sOEWhaMfpa
4LRJ0Vm6hwmYF0iyBdZ+NSV1NpSAU3LQe4HLTv6R4BzB9oWU8wW2t+wWMwoj
YKyzIFTgXuk2kQvhBFrMvTUFDvtUL9R/jMDdI6RC0B+pjApiEkHA7zfnP4Fz
3Eum0k9gHq2zCwXjY6nnTw5FEnDde3G6wNX9OWOW0QQKynebCPLY8fDScM5J
AmsXHPcS5DUUEv+bG0Ng382e/YI8TxlH/twaS6A+NnuvIO+5QiF9eacJiI/c
NBXMR36Vb4/QGdwPc/7OEMyXybmdlPNZXI9q2U8F8+szx6JdPIFA98PGGsH8
T31d07Y7EddjktfKS+ik6ys/lyYRuCBrHnkRXbxsXoPvBQLDb7nNyWirEZn6
NykEKhzedySiO0tE3s9NJXDHcrD+HFps46/K6jQCQQtVfeLQV8SpN4vTCUjt
MpglWM86Hz6/OppBwOBd5vWTaEeXNyWaWQQUq10jItC3gi/eP3Md99NEwdyD
aEP9M3ldNwjYvbabHoyumzp+2/AWgT2rLVoD0KOx3tf6cgmc5Wkb7UWbZhld
tMwnsMv2a6cT+vMenfM5Dwi4iCtJbxPst6VqidwCAqYi99W2oFMeS8XlFRKI
jquTsEBzqr+Fi5cQWOHxz3N9tEe5xs5NpQQmWXLhumjq2T6If4Z5OkcvWS7w
PZ6I1EsCAcYZ69UFTlBMnFlOIEyqKEwW3b11y7U5tdjethez+/D88tqUFr29
jkDcMrd4Ct29ocMzvZ6Anqj3YLvAegGa8xoIvPJ8dZCN7pGPK1zYTIB9+8iZ
Z+je1rJK9W8ENipsaI1Fe7NFc/d2EbC+cVtUcH731ticyeUQ2D55ZJHgfP/x
/IuNZjcBjYfecv7ovqyJ5n/7CWRtOyCyGd2/W6dPb5SAdJGn5Ez0YF+21IYZ
fDC1mf4gBu+bA8U/eVnifLjQ//l4GHoixugnV4IPk31RrGC0kHJzzSMpPriJ
7XBzRcs7SsUtluPDYKJ+3Ar0mvJwYaGFfJBvy9/eLLiPM7ZPlOnx4ZloYIcM
eqbP7Z55LD5Itt3OE0anskabQvX50Kzw2m0M7/urH5Oe6BjyYe7n82Ft6Mei
7/ZfXYf9KS+9eAPddoD1M8KSD6prHE1XojWtZLvWuvNhce21eCOsHzKsCoui
d2H7/3vrrY2WtLY7+243H9Sf9i+fj/5lfUFvqycf5m9hRdBYvzzeLBvr6cuH
H42NymXodXay2vEH+bDNyPfwOrS9s2xIy2k+3J110VNQH1U6F1ouOMOHec6a
OfPQLBe7hXvO8uEG8/Z/4mhF1wvvBs7xodVTobU7jYJ2N1kl/nk+NBoqkUy0
9y7ZMvVMPmQv/+4shj7mLSt+8CEfvmcYuNRiPZcQPhD26BEfHIa5O0vQWUmV
P38X8sHMINP2Jvpl8eGPAUWYz/Go0XA0X+xL2t5nfHC53WSthT5558YSl0o+
NOi0GZzE+jK+b9VGkxY+9GwU46pj/Znu73RChuGD86uvLJ1kCh5sMyA+QgzU
Sn6Wm42uWDvneLkwA6k5em2TWP/+kvkcFjKdARV6tnI52vSJ4+EvEgxYfY1K
tkcP8hwCc+UZ2DloK3IA62eTxG2u67UZ0BFRMcrCevvHw60GQY4M8OJXDd/H
+n3zoxbTGbkMeI9QuU6RFLz3KV4SfIcB9p6keB20pfIl8ZY8fG5W6yCGNkuy
/5B7nwGOLi/r6Qm8fwLrXcwLGUh3N6Hl0Ms0K/fHvmQAZrybqD6O53FOwVWR
TwwEzZ3sVAnH8W9PjA5oYuCOXNGlCfxeEZ8Z6PWpmQHbjbKratHTT2hp32xl
oEy7QCUETbvdemraycDtLZEf3h7F8SleqYvuZ6BuV7/3niN4nnw8+rB/gIF3
96k4Frrv9PYUu0EGFOpikmegu8flnVX+MFAc+tr+QSjO/6fk72/GGEia1Lo7
eRjr5XNB7zQnGHjlG8zUoFtNN989P8mATcLutVnopkLJ4N1TDBxst/A3Qdv7
9dtV8xg4NA3c5NANi6tZKwkDKfXf9bpD8PulNVcxnc9A73OD4afouuTYKYZh
oMDCMC0O/X9qInkK
"]]},
{RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.],
LineBox[CompressedData["
1:eJwVjHk01AkAxx2rhmcd5UjZoocmZtGKNYovUul8YR3ZjWrFjpUpeVvRWjYz
K61IJVk3IYwZSQfVzxFNOTJmNrbXqqSJMXbHMGOcv23/+LzPH5/3PpZHmf7H
NNTU1PZ94n9vPpGgpaYmQ3OkqtX5hQjR1soCL3UZkpUjcVd6RCh5ddLlvIYM
wfwC3ckuEfS2M6KWackgbTo3V8MXYdQs5JmujgxxhXYtFm0i/NHunGFmJEPC
SzuFVqMIGqaTRk5UGUoFXCMiTwTBoyiryAMy9B6eLpEeE6Forw4rs1IGITW5
114uRGoO3U2ukuFlEyX8ZqoQm/kc94e+k+gNTw4ZNxXCeMvAr7uKJqHMyeSn
Fvfjc0aGuYFsElKOjWurQz++HlTN3Ngpx+LylXSrOwLYOdyyvFggx9U02t1I
ugAuXuKixgk5vPm9Pvc6+hB3dTC8zGcKej4PLDR39qGjNDj++bUp8Cdu3rEr
eoHh3DUDIR+mkG+ZHT91uxeb+q4mVm+ZRn5ZUGsMrwfX2/LinX+fRkLwWLzj
o27YKwO5599PY3WaDSe2oQtZhkGudLoCVlWMAym1z6E9OmxMvaDAWkow63Lj
MzSQVs6v3ijgqRNgJq/kI75/qqrDUYmR9rMSbvVTRIlPh8WnKeH92qO1IKcT
VsOeR9OHlFgovCdQY3dAlZvT9Jf9DDw9nDeEsZ9gMavkGylrBpTYB4whZjvM
fOjbXAdmUP4x4CAnug3kkHHyuI0KAhtDdh6jFabNlUtHflbhsYZrNsevBTv8
OwWGvSqo9Pc4WV5/jArbPHGq9SzcL0TR3Q49hEBye1d44iyG0svS7GhNWBkZ
on6qaxZia0efm+vug72+YW7X+jncHj7Gqv/QCIO8PN/QU3MQ0b8UP+lrQMRs
54igbQ4R2wyO13XUQ3Kjqk3XbB4/sKqC9NZw8ZvCRJHMmMdesUehvbIagdlu
jPst8yjfQDd/droCbhLJpo0rFkAzimzO314Kyca3vupHF9BSmh34yC8fXdXp
Jfn3F7CqR1WWGnUF7eLW/Se1F2Fg8CQnxoGNie/r3LWPLMLT5Pg/HP5hvOO/
SaDcW8RbY02b9UMJhHS0QNdr+RLUCv82Yw9mENTdhm+3hi0hZFnbx9zu60Ri
pCEljvupmzIHTSjFRG50UJIeuQRzpuiu+opy4iP3gEtsIInv7C+ZyJhVxD7e
gPfyChKhSW/6HX6pJZ5H3bVmVpKw4G+i7bxcS/iuu0YZqCKRJdzjFFZWS/hc
8u+tqCGxlBqx++LTWmJrTM/BHfUkunu0brzT5xA0avsJ1kMStrOKwqRiDqFb
XFeoKSSRFLCOdba1jkgPzkj5UUSCUjsezhbWERT9mAjhnyQ07oyPZH+oIz5L
2mhbNkii+pzRRLU2l1j4trzRe+jT//Go/KUfl5Ca5nWnjJGQ1OU3rx3mEtEv
znDHJCTkoaDZTHOJUXZwtp+UxKrgilITTR4xojAKsfyXhHvWatcvjHnEa2Hm
+5ZpElwz/QAqnUeEXoztpCpJmKcn1th684hB7323smZIKGg1A7Q9PEJUr8M8
PEeCT66ucTzEI/wZY378eRK5NR7+X0XyiD4L/mbHRRI/pUi7nZg8Yv9ghWnu
Egmqp8ka5zM8ojuTNUeSJIQVnK0uKTziP78phL4=
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{1., 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, 2}, {-0.06447161656765425, 0.581865223202336}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{3.687951519211877*^9, 3.6879528977846212`*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Solve", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"{", "\[IndentingNewLine]",
RowBox[{
RowBox[{"z2next", " ", "\[Equal]", " ",
RowBox[{"z2curr", " ", "+", " ",
RowBox[{"h", " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "z2next"}], " ", "/", " ", "x"}], " ", "-", " ",
RowBox[{
RowBox[{"(",
RowBox[{"1", " ", "-", " ",
RowBox[{"1", "/",
RowBox[{"x", "^", "2"}]}]}], ")"}], " ", "z1next"}]}],
")"}]}]}]}], ",", "\[IndentingNewLine]",
RowBox[{"z1next", " ", "\[Equal]", " ",
RowBox[{"z1curr", " ", "+", " ",
RowBox[{"h", " ", "z2next"}]}]}]}], "\[IndentingNewLine]", "}"}], ",",
"\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"z1next", ",", " ", "z2next"}], "}"}]}], "\[IndentingNewLine]",
"]"}]], "Input",
CellChangeTimes->{{3.68795255574511*^9, 3.6879526211878357`*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{"z1next", "\[Rule]",
FractionBox[
RowBox[{"x", " ",
RowBox[{"(",
RowBox[{
RowBox[{"h", " ", "z1curr"}], "+",
RowBox[{"x", " ", "z1curr"}], "+",
RowBox[{"h", " ", "x", " ", "z2curr"}]}], ")"}]}],
RowBox[{
RowBox[{"-",
SuperscriptBox["h", "2"]}], "+",
RowBox[{"h", " ", "x"}], "+",
SuperscriptBox["x", "2"], "+",
RowBox[{
SuperscriptBox["h", "2"], " ",
SuperscriptBox["x", "2"]}]}]]}], ",",
RowBox[{"z2next", "\[Rule]",
RowBox[{"-",
FractionBox[
RowBox[{
RowBox[{
RowBox[{"-", "h"}], " ", "z1curr"}], "+",
RowBox[{"h", " ",
SuperscriptBox["x", "2"], " ", "z1curr"}], "-",
RowBox[{
SuperscriptBox["x", "2"], " ", "z2curr"}]}],
RowBox[{
RowBox[{"-",
SuperscriptBox["h", "2"]}], "+",
RowBox[{"h", " ", "x"}], "+",
SuperscriptBox["x", "2"], "+",
RowBox[{
SuperscriptBox["h", "2"], " ",
SuperscriptBox["x", "2"]}]}]]}]}]}], "}"}], "}"}]], "Output",
CellChangeTimes->{3.68795262166017*^9, 3.68795289781152*^9}]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"next", "[",
RowBox[{
RowBox[{"{",
RowBox[{"z1curr_", ",", " ", "z2curr_"}], "}"}], ",", " ", "myx_", ",",
" ", "h_"}], "]"}], " ", ":=", " ",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{",
RowBox[{"z1next", ",", " ", "z2next", ",", " ",
RowBox[{"x", " ", "=", " ",
RowBox[{"myx", " ", "+", " ", "h"}]}]}], "}"}], ",",
"\[IndentingNewLine]",
RowBox[{
RowBox[{"z1next", " ", "=", " ",
FractionBox[
RowBox[{"x", " ",
RowBox[{"(",
RowBox[{
RowBox[{"h", " ", "z1curr"}], "+",
RowBox[{"x", " ", "z1curr"}], "+",
RowBox[{"h", " ", "x", " ", "z2curr"}]}], ")"}]}],
RowBox[{
RowBox[{"-",
SuperscriptBox["h", "2"]}], "+",
RowBox[{"h", " ", "x"}], "+",
SuperscriptBox["x", "2"], "+",
RowBox[{
SuperscriptBox["h", "2"], " ",
SuperscriptBox["x", "2"]}]}]]}], ";", "\[IndentingNewLine]",
RowBox[{"z2next", " ", "=", " ",
RowBox[{"-",
FractionBox[
RowBox[{
RowBox[{
RowBox[{"-", "h"}], " ", "z1curr"}], "+",
RowBox[{"h", " ",
SuperscriptBox["x", "2"], " ", "z1curr"}], "-",
RowBox[{
SuperscriptBox["x", "2"], " ", "z2curr"}]}],
RowBox[{
RowBox[{"-",
SuperscriptBox["h", "2"]}], "+",
RowBox[{"h", " ", "x"}], "+",
SuperscriptBox["x", "2"], "+",
RowBox[{
SuperscriptBox["h", "2"], " ",
SuperscriptBox["x", "2"]}]}]]}]}], ";", "\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"z1next", ",", " ", "z2next"}], "}"}]}]}],
"\[IndentingNewLine]", "]"}]}]], "Input",
CellChangeTimes->{{3.687952654994615*^9, 3.6879526669313927`*^9}, {
3.68795270699533*^9, 3.687952774571446*^9}, {3.687952810690868*^9,
3.687952812858673*^9}, {3.68795312977686*^9, 3.687953131712039*^9}}],
Cell[BoxData[{
RowBox[{
RowBox[{"x", " ", "=", " ",
RowBox[{"{", "1", "}"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"z", " ", "=", " ",
RowBox[{"{",
RowBox[{"{",
RowBox[{"0.44005058574493355`", ",", "0.32514710081303305`"}], "}"}],
"}"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.687952819885274*^9, 3.687952867466277*^9}, {
3.687952936314144*^9, 3.687952943521549*^9}, {3.687953115720479*^9,
3.6879531173360987`*^9}}],
Cell[BoxData[
RowBox[{
RowBox[{"h", " ", "=", " ", "0.05"}], ";"}]], "Input",
CellChangeTimes->{{3.68795298243886*^9, 3.687952986832855*^9}}],
Cell[BoxData[
RowBox[{"Do", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"AppendTo", "[",
RowBox[{"z", ",", " ",
RowBox[{"Evaluate", "[",
RowBox[{"next", "[",
RowBox[{
RowBox[{"z", "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}], ",", " ",
RowBox[{"x", "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}], ",", " ", "h"}], "]"}], "]"}]}],
"]"}], ";", "\[IndentingNewLine]",
RowBox[{"AppendTo", "[",
RowBox[{"x", ",", " ",
RowBox[{
RowBox[{"x", "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}], " ", "+", " ", "h"}]}], "]"}]}],
",", "\[IndentingNewLine]", "20"}], "\[IndentingNewLine]", "]"}]], "Input",\
CellChangeTimes->{{3.687952951550411*^9, 3.687952998700535*^9}, {
3.687953031128552*^9, 3.6879530947921257`*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Show", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Plot", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"BesselJ", "[",
RowBox[{"1", ",", " ", "x"}], "]"}], ",", " ", "dBesselJ1"}], "}"}],
",", "\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"x", ",", " ", "1", ",", " ", "2"}], "}"}], ",",
"\[IndentingNewLine]",
RowBox[{"PlotStyle", " ", "\[Rule]", " ",
RowBox[{"{",
RowBox[{"Blue", ",", " ", "Red"}], "}"}]}]}], "\[IndentingNewLine]",
"]"}], ",", "\[IndentingNewLine]",
RowBox[{"ListPlot", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x", "[",
RowBox[{"[", "i", "]"}], "]"}], ",", " ",
RowBox[{"z", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "1"}], "]"}], "]"}]}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"i", ",", " ", "1", ",", " ", "20"}], "}"}]}], "]"}], ",",
"\[IndentingNewLine]",
RowBox[{"PlotStyle", " ", "\[Rule]", " ", "Blue"}]}],
"\[IndentingNewLine]", "]"}], ",", "\[IndentingNewLine]",
RowBox[{"ListPlot", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x", "[",
RowBox[{"[", "i", "]"}], "]"}], ",", " ",
RowBox[{"z", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "2"}], "]"}], "]"}]}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"i", ",", " ", "1", ",", " ", "20"}], "}"}]}], "]"}], ",",
"\[IndentingNewLine]",
RowBox[{"PlotStyle", " ", "\[Rule]", " ", "Red"}]}],
"\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]], "Input",
CellChangeTimes->{{3.687953151463984*^9, 3.687953187447447*^9}, {
3.687953219177075*^9, 3.6879532547988777`*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {},
{RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.],
LineBox[CompressedData["
1:eJwV13k4VO0bB3AS4UURSZssCW8pGkLLLXkRUoQspUL2aJGihCQlS0kilDap
pJTQppB+hDIiFDInhFR2kznPnN89f831uc65nnme77PdR8Uj2H7vNCEhIRFh
ISHBL2t/uKiQ0BCo/lMyXqPVDv5LJrI3CA+BvUalZuC/7XD9ywGDmGlDsPnQ
hLHM8naQ+c/PR0x0CO661NJ2K9uhT8m5RkpyCMY8eOMtBu2QWamfqCQ/BIts
Ol52bGyHaYrD8qs0h2BGADe5YWc7sF/5qHtvHQKvhvRe+wvtcM1GMjb5zhBo
vJl0Xc20w6k0I+MR7hBYxw6Hpnl2AKv6wbqXlsPg3v7FoqOyAxTWtJzcdG0Y
krSXa3ZrdYK0X+KCWUPDsH+LShv7dCesbuVOZliMwKqfqiVC/Z3w74q7Kuey
RyBHbOjUR/NvYLCh99rTXyOg9Xn676nMb3AwtXXXTbNRcDpdblY5/g2qbmwP
eX9pFKKkr8qkqHUBlT6/xblnFJoW/Wdua9UFug2px+6tGQNXneqGGYe64HLF
lRD9hDFwf9IxXnu5C3QmHB/GfB+D3QrF663LuuC8rJOhkdE42GmKWxVRXSDR
Rylonh2HigTpMn8RDjxh1PW/fBuHjMKAjrglHAhpHM2rWjkBXWHZSfJmHPDp
PeIecmYCXPW8zd56ckCdMvGI75yAULaDSXA0B7jpac/bdCZBOllrOTuTA+T8
dYfB2ElwHRWOEirhgJKZ0UbDlknQcvd3/+cjB5hOhaifGlwwGNLNc+3hgOKL
O/w9EVxgK/o2qdEcMLd/x5b9wIWUHLOrsrMoyNW+0ntqyV9o635vqq9OAXvg
8aZdx/6Cfsz245EGFMz2dhY+VPsXdsysm7/GnILTqk+mNqlOgfu+pswWRwpm
Xbli6XpoCpbrPpGQ9qTA6++7bnbFFLRvKJUpDaJgICOvQkqJB/VZNk+GwyiI
G58zHuXHA98PhZdkYyiYivnR2OHPA7HEUAsDdJBc6SPjQB68kHzV5Yp2WOES
OBbEg69qH2tuohf7ZX73CeGBMFcslHWKguftyo22kTyI/FtmYBdLwe9KzYcL
LvEgI0rD5WgcBY4pxn6lb3iw+A1ZHHGOgqJtO5fxy/H/3U8+yEDLK0T92VjJ
g3nT7qsXo5suVx3+WMUDCfe+xt9o++ytJ3vf88AwYfaKXQkU2OZ5Z8o38yAL
lveuSaTAsiylfn8/D64+itzSk0SB8cCArpYcDcyq+JncCxQYvTZP8Z1Nw9uW
ngaZFAoMU2+M3JGn4fbX5xFL0PrrXZ9oKNKwSJWdYY9emVLDUl9AQ9eph1/v
ozWM81Yra9CQVG982vUiBXLxe9cpGNMQPu99aF4qBbK7yrMd1tCgHuBx5yV6
Fmsh/+JaGgZMDlc2oKU7m8rkgAa9O5JlXLS4npnJLDMalqWOv7S4RAG/TXXj
P7bY3lAzh4Me0OqyFPagQV5sTrXYZczzk2mxmycN65fZ9iihhyNuqxV70WAl
7TmwDM1l+/P9fWj4a6GWb48WCx9/0hRIg+M1RjobrfpeUjnvCA0vLA2tl6dT
4OqvP26bSEP/AVb+mgwK3OXTPe4m0dAcolVljd5TNvVR5DwNt4pJlRvaV+7N
vdIUGqq2mkUeQ4c+t9qtmk7DLzebhmfoFMndteM3aADtHYO6VyiovRd/PauU
ht5OdelZmbg+4mDs2zMaOFdTNy5A3/MaM1d7QUNA0bS9muj4Re6Dd1/R8Gr1
F18TtFWK7uqSChqEJcNCgtE1Ya117Doa9p07m1mDLnVKVFb4QEPeov7oJnTe
KtODzh9p6Dlgue0b+syv+4rf2DTkFouVjKEt90R6DH6mIWjbqZOLsij4n6UG
V4yD74tkaAWhi5d8tbKiaAhrZ14fQedOO5+d+J2G1CmfDdHo0y//msr30hAY
uFY0FW2+sj5R5ScNrCHu+VJ0lWKI2tpxGhYsbCkl6KIxzdDICRqy3N0Z0WwK
brE7qismMZ/TfStl0KfOWQRtmsI8R0UdF6HNmHnPnBga8k18r69DV/aW2x6Q
IDC7grI6ii5Njyj4KEmA8W9rjkTnWxnJ6EgRYDkss4xDpxU8qh+QIRAb9oO6
hPYNvWbjJU/AS6FkqBC9U9Mtv0IBn1ddLC5F23+ZI6WiSOBmR73na/Sa9Um1
HUoEfNtdQ+rQMmIRVtuVCWQYt+/rQU8rNbz3dDGBv1vygn6iJ/3GJORVCTjz
G52G0V31ATUN6gSsys06aPSTS66WVtoElJZerZC7SoHzEkPzdfoERtWt/xih
bVpGb2caENi8abb3evSGsw9Fp1YTyDy19q0pWvuXRlWxMfav+7GBNZoUKZit
NCFQ1LhqzBU94s2+mbSBwHd9iVx3dO/cRJFfpgR4tAl4oBuOT6+8+x8Btd2t
S/3QN8xGN6hZE9BVu88+jL48UXA92obA8c3pZUfRCXn+wl2bCYzr9qccQx+W
ot5kbcXxtRQPR6H9y7IX8+wIXPyz6ngM2n2/S5TLNuzvfoXfsWjLpgaY40TA
cafYmXj02tMJ10K2E4hSmfMoAa1raMk0OhMIZuJfJ6HnZb0uS3YjoPLiVmoK
eqbtsUW/dxA4N7LJNRU9XWj1CRt3AiKTDjPS0L88C9ZJ7CEwMj9bMQMdpXXT
eI8HgUeHMkKuoOX/XF79zJNAZOW3Z5novKIEluxeAvkG+39kCfoTHq3r503g
FteKZAvygVCdch/Mi7WfdxXtJRrwr5IfgUDFLs41NPf9Ls0D/gRqa688yhHk
c95hSU0AAfkTOT7X0YudNqmq7MP15TosdgNdNH+9clgQro+ExCSBLTl6C9jB
BAxND/MF/pq7VEnrAIG83Dynm+igwAVzog8SGOxdelngaXqys9sOEWhaMfpa
4LRJ0Vm6hwmYF0iyBdZ+NSV1NpSAU3LQe4HLTv6R4BzB9oWU8wW2t+wWMwoj
YKyzIFTgXuk2kQvhBFrMvTUFDvtUL9R/jMDdI6RC0B+pjApiEkHA7zfnP4Fz
3Eum0k9gHq2zCwXjY6nnTw5FEnDde3G6wNX9OWOW0QQKynebCPLY8fDScM5J
AmsXHPcS5DUUEv+bG0Ng382e/YI8TxlH/twaS6A+NnuvIO+5QiF9eacJiI/c
NBXMR36Vb4/QGdwPc/7OEMyXybmdlPNZXI9q2U8F8+szx6JdPIFA98PGGsH8
T31d07Y7EddjktfKS+ik6ys/lyYRuCBrHnkRXbxsXoPvBQLDb7nNyWirEZn6
NykEKhzedySiO0tE3s9NJXDHcrD+HFps46/K6jQCQQtVfeLQV8SpN4vTCUjt
MpglWM86Hz6/OppBwOBd5vWTaEeXNyWaWQQUq10jItC3gi/eP3Md99NEwdyD
aEP9M3ldNwjYvbabHoyumzp+2/AWgT2rLVoD0KOx3tf6cgmc5Wkb7UWbZhld
tMwnsMv2a6cT+vMenfM5Dwi4iCtJbxPst6VqidwCAqYi99W2oFMeS8XlFRKI
jquTsEBzqr+Fi5cQWOHxz3N9tEe5xs5NpQQmWXLhumjq2T6If4Z5OkcvWS7w
PZ6I1EsCAcYZ69UFTlBMnFlOIEyqKEwW3b11y7U5tdjethez+/D88tqUFr29
jkDcMrd4Ct29ocMzvZ6Anqj3YLvAegGa8xoIvPJ8dZCN7pGPK1zYTIB9+8iZ
Z+je1rJK9W8ENipsaI1Fe7NFc/d2EbC+cVtUcH731ticyeUQ2D55ZJHgfP/x
/IuNZjcBjYfecv7ovqyJ5n/7CWRtOyCyGd2/W6dPb5SAdJGn5Ez0YF+21IYZ
fDC1mf4gBu+bA8U/eVnifLjQ//l4GHoixugnV4IPk31RrGC0kHJzzSMpPriJ
7XBzRcs7SsUtluPDYKJ+3Ar0mvJwYaGFfJBvy9/eLLiPM7ZPlOnx4ZloYIcM
eqbP7Z55LD5Itt3OE0anskabQvX50Kzw2m0M7/urH5Oe6BjyYe7n82Ft6Mei
7/ZfXYf9KS+9eAPddoD1M8KSD6prHE1XojWtZLvWuvNhce21eCOsHzKsCoui
d2H7/3vrrY2WtLY7+243H9Sf9i+fj/5lfUFvqycf5m9hRdBYvzzeLBvr6cuH
H42NymXodXay2vEH+bDNyPfwOrS9s2xIy2k+3J110VNQH1U6F1ouOMOHec6a
OfPQLBe7hXvO8uEG8/Z/4mhF1wvvBs7xodVTobU7jYJ2N1kl/nk+NBoqkUy0
9y7ZMvVMPmQv/+4shj7mLSt+8CEfvmcYuNRiPZcQPhD26BEfHIa5O0vQWUmV
P38X8sHMINP2Jvpl8eGPAUWYz/Go0XA0X+xL2t5nfHC53WSthT5558YSl0o+
NOi0GZzE+jK+b9VGkxY+9GwU46pj/Znu73RChuGD86uvLJ1kCh5sMyA+QgzU
Sn6Wm42uWDvneLkwA6k5em2TWP/+kvkcFjKdARV6tnI52vSJ4+EvEgxYfY1K
tkcP8hwCc+UZ2DloK3IA62eTxG2u67UZ0BFRMcrCevvHw60GQY4M8OJXDd/H
+n3zoxbTGbkMeI9QuU6RFLz3KV4SfIcB9p6keB20pfIl8ZY8fG5W6yCGNkuy
/5B7nwGOLi/r6Qm8fwLrXcwLGUh3N6Hl0Ms0K/fHvmQAZrybqD6O53FOwVWR
TwwEzZ3sVAnH8W9PjA5oYuCOXNGlCfxeEZ8Z6PWpmQHbjbKratHTT2hp32xl
oEy7QCUETbvdemraycDtLZEf3h7F8SleqYvuZ6BuV7/3niN4nnw8+rB/gIF3
96k4Frrv9PYUu0EGFOpikmegu8flnVX+MFAc+tr+QSjO/6fk72/GGEia1Lo7
eRjr5XNB7zQnGHjlG8zUoFtNN989P8mATcLutVnopkLJ4N1TDBxst/A3Qdv7
9dtV8xg4NA3c5NANi6tZKwkDKfXf9bpD8PulNVcxnc9A73OD4afouuTYKYZh
oMDCMC0O/X9qInkK
"]]},
{RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.],
LineBox[CompressedData["
1:eJwVjHk01AkAxx2rhmcd5UjZoocmZtGKNYovUul8YR3ZjWrFjpUpeVvRWjYz
K61IJVk3IYwZSQfVzxFNOTJmNrbXqqSJMXbHMGOcv23/+LzPH5/3PpZHmf7H
NNTU1PZ94n9vPpGgpaYmQ3OkqtX5hQjR1soCL3UZkpUjcVd6RCh5ddLlvIYM
wfwC3ckuEfS2M6KWackgbTo3V8MXYdQs5JmujgxxhXYtFm0i/NHunGFmJEPC
SzuFVqMIGqaTRk5UGUoFXCMiTwTBoyiryAMy9B6eLpEeE6Forw4rs1IGITW5
114uRGoO3U2ukuFlEyX8ZqoQm/kc94e+k+gNTw4ZNxXCeMvAr7uKJqHMyeSn
Fvfjc0aGuYFsElKOjWurQz++HlTN3Ngpx+LylXSrOwLYOdyyvFggx9U02t1I
ugAuXuKixgk5vPm9Pvc6+hB3dTC8zGcKej4PLDR39qGjNDj++bUp8Cdu3rEr
eoHh3DUDIR+mkG+ZHT91uxeb+q4mVm+ZRn5ZUGsMrwfX2/LinX+fRkLwWLzj
o27YKwO5599PY3WaDSe2oQtZhkGudLoCVlWMAym1z6E9OmxMvaDAWkow63Lj
MzSQVs6v3ijgqRNgJq/kI75/qqrDUYmR9rMSbvVTRIlPh8WnKeH92qO1IKcT
VsOeR9OHlFgovCdQY3dAlZvT9Jf9DDw9nDeEsZ9gMavkGylrBpTYB4whZjvM
fOjbXAdmUP4x4CAnug3kkHHyuI0KAhtDdh6jFabNlUtHflbhsYZrNsevBTv8
OwWGvSqo9Pc4WV5/jArbPHGq9SzcL0TR3Q49hEBye1d44iyG0svS7GhNWBkZ
on6qaxZia0efm+vug72+YW7X+jncHj7Gqv/QCIO8PN/QU3MQ0b8UP+lrQMRs
54igbQ4R2wyO13XUQ3Kjqk3XbB4/sKqC9NZw8ZvCRJHMmMdesUehvbIagdlu
jPst8yjfQDd/droCbhLJpo0rFkAzimzO314Kyca3vupHF9BSmh34yC8fXdXp
Jfn3F7CqR1WWGnUF7eLW/Se1F2Fg8CQnxoGNie/r3LWPLMLT5Pg/HP5hvOO/
SaDcW8RbY02b9UMJhHS0QNdr+RLUCv82Yw9mENTdhm+3hi0hZFnbx9zu60Ri
pCEljvupmzIHTSjFRG50UJIeuQRzpuiu+opy4iP3gEtsIInv7C+ZyJhVxD7e
gPfyChKhSW/6HX6pJZ5H3bVmVpKw4G+i7bxcS/iuu0YZqCKRJdzjFFZWS/hc
8u+tqCGxlBqx++LTWmJrTM/BHfUkunu0brzT5xA0avsJ1kMStrOKwqRiDqFb
XFeoKSSRFLCOdba1jkgPzkj5UUSCUjsezhbWERT9mAjhnyQ07oyPZH+oIz5L
2mhbNkii+pzRRLU2l1j4trzRe+jT//Go/KUfl5Ca5nWnjJGQ1OU3rx3mEtEv
znDHJCTkoaDZTHOJUXZwtp+UxKrgilITTR4xojAKsfyXhHvWatcvjHnEa2Hm
+5ZpElwz/QAqnUeEXoztpCpJmKcn1th684hB7323smZIKGg1A7Q9PEJUr8M8
PEeCT66ucTzEI/wZY378eRK5NR7+X0XyiD4L/mbHRRI/pUi7nZg8Yv9ghWnu
Egmqp8ka5zM8ojuTNUeSJIQVnK0uKTziP78phL4=
"]]}}, {{}, {{},
{RGBColor[0, 0, 1], PointSize[0.012833333333333334`], AbsoluteThickness[
1.6], PointBox[CompressedData["
1:eJxTTMoPSmViYGAQAWIQDQEf7HN5pr87qXnH/uwZEPhgH9ZR9ipJ7a79rJkg
8NHe4u1xk3Mi9+zT00Dgk73sG1HplV/u2ZsYg8Bn+/x5didmHL1vzwg274v9
rB+ZZkudH9ifA5v3xX5twKKNO+c8sJ8NNu+rvdbiK8tfvH1gnwE275v9rdQe
WVHzh/amYPO+2+dwm4eqVTy0h7jxh/3kz2Y32Dc+tD8PNu+HfabPgtqVDx/a
zwGb99Oe5+jsvP+cj+wzweb9svc4n6PBpPXI3gxs3m/7tS+latY5PrJnBpv3
x/4J7/5u1sBH9hfA5v2xj35bfpQn4pH9XLB5f+1/qHDeORv2yD4LbN4/+wmM
2nvDfB/Zm4PN+29vdWVi607rR/YAhL+aaw==
"]]}, {}}, {}}, {{}, {{},
{RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
1.6], PointBox[CompressedData["
1:eJxTTMoPSmViYGAQAWIQDQEf7ItTJU6anr9if/YMCHyw91tXuvbz7sv2s2aC
wEf7ow+jzyUuuGSfngYCn+zPvhT/J1l70d7EGAQ+2ze9023fGXTBnhFs3hf7
he9Cth90Pmd/DmzeF3uOttmNaddO288Gm/fV/uuxne/XxJ20zwCb981+/7Wc
k0/uHrM3BZv33X7l3uDUxNAj9hA3/rC/0vPjfcWxg/bnweb9sNdtu57InrrP
fg7YvJ/2amnac7eZ7rTPBJv3y/59yjKlLYxb7M3A5v22t/+3YGbG9LX2zGDz
/tjbuOr2XtZZbH8BbN4f+3kr9PwEJk20nws27699/uXly21eZezPApv3z56p
NUL77Mrp+83B5v23l1hdUd/Ksmw/AGb5npA=
"]]}, {}}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{1., 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, 2}, {-0.06447161656765425, 0.581865223202336}},
PlotRangeClipping->True,
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.05],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.687953146432378*^9, 3.687953187966303*^9, 3.687953219983635*^9, {
3.687953251654253*^9, 3.687953255151575*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"Table", "[", "\[IndentingNewLine]",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x", "[",
RowBox[{"[", "i", "]"}], "]"}], ",", " ",
RowBox[{"z", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "1"}], "]"}], "]"}], ",", " ",
RowBox[{"z", "[",
RowBox[{"[",
RowBox[{"i", ",", " ", "2"}], "]"}], "]"}]}], "}"}], ",",
"\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"i", ",", " ", "1", ",", " ", "20"}], "}"}]}],
"\[IndentingNewLine]", "]"}], " ", "//", " ", "TableForm"}]], "Input",
CellChangeTimes->{{3.687953282490951*^9, 3.6879533023903112`*^9}}],
Cell[BoxData[
TagBox[GridBox[{
{"1", "0.44005058574493355`", "0.32514710081303305`"},
{"1.05`", "0.45546791931895336`", "0.3083466714803976`"},
{"1.1`", "0.47001986650839767`", "0.29103894378888767`"},
{"1.1500000000000001`", "0.4836828960903306`", "0.2732605916386575`"},
{"1.2000000000000002`", "0.4964353518975945`", "0.2550491161452798`"},
{"1.2500000000000002`", "0.5082574904224828`", "0.23644277049776455`"},
{"1.3000000000000003`", "0.5191315144934571`", "0.21748048141948595`"},
{"1.3500000000000003`", "0.5290416028920835`", "0.19820176797252712`"},
{"1.4000000000000004`", "0.5379739358060605`", "0.17864665827953805`"},
{"1.4500000000000004`", "0.5459167160373517`", "0.15885560462582518`"},
{"1.5000000000000004`", "0.5528601859037842`", "0.13886939732865378`"},
{"1.5500000000000005`", "0.5587966397892049`", "0.11872907770841641`"},
{"1.6000000000000005`", "0.563720432312182`", "0.09847585045954385`"},
{"1.6500000000000006`", "0.5676279820968171`", "0.07815099569269866`"},
{"1.7000000000000006`", "0.5705177711418422`", "0.05779578090050252`"},
{"1.7500000000000007`", "0.5723903397960695`", "0.0374513730845458`"},
{"1.8000000000000007`", "0.5732482773595778`", "0.017158751270166854`"},
{"1.8500000000000008`", "0.5730962083408988`",
RowBox[{"-", "0.003041380373583306`"}]},
{"1.9000000000000008`", "0.5719407744109493`",
RowBox[{"-", "0.0231086785989918`"}]},
{"1.9500000000000008`", "0.5697906121046181`",
RowBox[{"-", "0.04300324612662204`"}]}
},
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
"RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[2.0999999999999996`]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}],
Function[BoxForm`e$,
TableForm[BoxForm`e$]]]], "Output",
CellChangeTimes->{{3.687953298055676*^9, 3.68795330274361*^9}}]
}, Open ]]
},
WindowSize->{1054, 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, 156, 3, 32, "Input"],
Cell[739, 27, 177, 4, 32, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[953, 36, 257, 6, 32, "Input"],
Cell[1213, 44, 334, 10, 49, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[1584, 59, 800, 24, 32, "Input"],
Cell[2387, 85, 149, 2, 32, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[2573, 92, 507, 12, 99, "Input"],
Cell[3083, 106, 2002, 41, 383, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[5122, 152, 568, 15, 121, "Input"],
Cell[5693, 169, 7786, 141, 230, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[13516, 315, 979, 25, 165, "Input"],
Cell[14498, 342, 1245, 38, 54, "Output"]
}, Open ]],
Cell[15758, 383, 1968, 54, 198, "Input"],
Cell[17729, 439, 465, 12, 55, "Input"],
Cell[18197, 453, 145, 3, 32, "Input"],
Cell[18345, 458, 894, 25, 121, "Input"],
Cell[CellGroupData[{
Cell[19264, 487, 1941, 50, 341, "Input"],
Cell[21208, 539, 8977, 163, 230, "Output"]
}, Open ]],
Cell[CellGroupData[{
Cell[30222, 707, 661, 18, 99, "Input"],
Cell[30886, 727, 2102, 38, 366, "Output"]
}, Open ]]
}
]
*)
(* End of internal cache information *)
Danny van Dyk