(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.0' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 209600, 5723]*) (*NotebookOutlinePosition[ 210336, 5748]*) (* CellTagsIndexPosition[ 210292, 5744]*) (*WindowFrame->Normal*) Notebook[{ Cell["Math 221 - Numerical PDE's, Lab 2", "Subtitle"], Cell["\<\ We've seen how finite difference formulas may be derived by various \ means in Lab. 1. Today we'll look at applying these formulas to generate ODE \ solvers and test the solvers.\ \>", "Text"], Cell[BoxData[ \(\(Off[General::spell1];\)\)], "Input", InitializationCell->True], Cell[CellGroupData[{ Cell["0. Define operators", "Section", InitializationCell->True], Cell[CellGroupData[{ Cell["Finite difference operators", "Subsection", InitializationCell->True], Cell[BoxData[{ \(\(dforward[q_[t_]] := q[t + k] - q[t];\)\), "\[IndentingNewLine]", \(\(dbackward[q_[t_]] := q[t] - q[t - k];\)\), "\[IndentingNewLine]", \(\(dcentered[q_[t_]] := q[t + k/2] - q[t - k/2];\)\), "\[IndentingNewLine]", \(dforward[q_[t_], n_] := If[n \[Equal] 1, \ dforward[q[t]], \ dforward[q[t + k] - q[t], n - 1]]\), "\[IndentingNewLine]", \(\(dforward[a_\ q_[t_], n_]\ = \ a\ dforward[q[t], n];\)\ \), "\[IndentingNewLine]", \(\(dforward[a_\ + \ b_, n_] = dforward[a, n] + dforward[b, n];\)\), "\[IndentingNewLine]", \(\(dbackward[q_[t_], n_] := If[n \[Equal] 1, \ dbackward[q[t]], \ dbackward[q[t] - q[t - k], n - 1]];\)\), "\[IndentingNewLine]", \(\(dbackward[a_\ q_[t_], n_] = a\ dbackward[q[t], n];\)\), "\[IndentingNewLine]", \(\(dbackward[a_\ + \ b_, n_] = dbackward[a, n] + dbackward[b, n];\)\), "\[IndentingNewLine]", \(\(dcentered[q_[t_], n_] := If[n \[Equal] 1, \ dcentered[q[t]], \ dcentered[q[t + k/2] - q[t - k/2], n - 1]];\)\), "\[IndentingNewLine]", \(\(dcentered[a_\ q_[t_], n_] = a\ dcentered[q[t], n];\)\), "\[IndentingNewLine]", \(\(dcentered[a_\ + \ b_, n_] = dcentered[a, n] + dcentered[b, n];\)\)}], "Input", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell["Tests", "Subsection", InitializationCell->True], Cell[CellGroupData[{ Cell[BoxData[{ \(dforward[q[t], 2]\), "\[IndentingNewLine]", \(dbackward[q[t], 2]\), "\[IndentingNewLine]", \(dcentered[q[t], 2]\)}], "Input", InitializationCell->True], Cell[BoxData[ \(q[t] - 2\ q[k + t] + q[2\ k + t]\)], "Output"], Cell[BoxData[ \(q[t] + q[\(-2\)\ k + t] - 2\ q[\(-k\) + t]\)], "Output"], Cell[BoxData[ \(\(-2\)\ q[t] + q[\(-k\) + t] + q[k + t]\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Operator series approximations of first derivative", "Subsection", InitializationCell->True], Cell[BoxData[{ \(s1forward[q_[t_], n_] := Simplify[\ Sum[\ Coefficient[Series[Log[1 + x], {x, 0, n}], x^j]\ dforward[q[t], j]\ , {j, 1, n}]/k]\), "\[IndentingNewLine]", \(s1backward[q_[t_], n_] := Simplify[\ Sum[\ Coefficient[Series[\(-Log[1 - x]\), {x, 0, n}], x^j]\ dbackward[ q[t], j]\ , {j, 1, n}]/k]\), "\[IndentingNewLine]", \(s1centered[q_[t_], n_] := Simplify[\ Sum[\ Coefficient[Series[2\ ArcSinh[x/2], {x, 0, n}], x^j]\ dcentered[ q[t], j]\ , {j, 1, n}]/k]\)}], "Input", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell["Tests", "Subsection", InitializationCell->True], Cell[CellGroupData[{ Cell[BoxData[ \(f4 = s1forward[q[t], 4]\)], "Input", InitializationCell->True], Cell[BoxData[ \(\(-\(\(25\ q[t] - 48\ q[k + t] + 36\ q[2\ k + t] - 16\ q[3\ k + t] + 3\ q[4\ k + t]\)\/\(12\ k\)\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Series[f4, {k, 0, 4}]\)], "Input", InitializationCell->True], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{ SuperscriptBox["q", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{\(1\/5\), " ", RowBox[{ SuperscriptBox["q", TagBox[\((5)\), Derivative], MultilineFunction->None], "[", "t", "]"}], " ", \(k\^4\)}], "+", InterpretationBox[\(O[k]\^5\), SeriesData[ k, 0, {}, 0, 5, 1], Editable->False]}], SeriesData[ k, 0, { Derivative[ 1][ q][ t], 0, 0, 0, Times[ Rational[ -1, 5], Derivative[ 5][ q][ t]]}, 0, 5, 1], Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(b4 = s1backward[q[t], 4]\)], "Input", InitializationCell->True], Cell[BoxData[ \(\(25\ q[t] + 3\ q[\(-4\)\ k + t] - 16\ q[\(-3\)\ k + t] + 36\ q[\(-2\)\ \ k + t] - 48\ q[\(-k\) + t]\)\/\(12\ k\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Series[b4, {k, 0, 4}]\)], "Input", InitializationCell->True], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{ SuperscriptBox["q", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{\(1\/5\), " ", RowBox[{ SuperscriptBox["q", TagBox[\((5)\), Derivative], MultilineFunction->None], "[", "t", "]"}], " ", \(k\^4\)}], "+", InterpretationBox[\(O[k]\^5\), SeriesData[ k, 0, {}, 0, 5, 1], Editable->False]}], SeriesData[ k, 0, { Derivative[ 1][ q][ t], 0, 0, 0, Times[ Rational[ -1, 5], Derivative[ 5][ q][ t]]}, 0, 5, 1], Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(c3 = s1centered[q[t], 3]\)], "Input", InitializationCell->True], Cell[BoxData[ \(\(q[\(-\(\(3\ k\)\/2\)\) + t] - 27\ q[\(-\(k\/2\)\) + t] + 27\ q[k\/2 + \ t] - q[\(3\ k\)\/2 + t]\)\/\(24\ k\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Series[c3, {k, 0, 4}]\)], "Input", InitializationCell->True], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{ SuperscriptBox["q", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{\(3\/640\), " ", RowBox[{ SuperscriptBox["q", TagBox[\((5)\), Derivative], MultilineFunction->None], "[", "t", "]"}], " ", \(k\^4\)}], "+", InterpretationBox[\(O[k]\^5\), SeriesData[ k, 0, {}, 0, 5, 1], Editable->False]}], SeriesData[ k, 0, { Derivative[ 1][ q][ t], 0, 0, 0, Times[ Rational[ -3, 640], Derivative[ 5][ q][ t]]}, 0, 5, 1], Editable->False]], "Output"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ 1. Numerically solve an ODE by finite difference approximation of \ the derivative\ \>", "Section"], Cell[CellGroupData[{ Cell["Choose an interesting problem", "Subsection"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[k, \[Lambda], a]\), "\[IndentingNewLine]", \(\(fmodel[t_, q_, \[Lambda]_, a_]\ = \ \(-\[Lambda]\)\ \ q[t]\ + \ a\ Sin[t];\)\), "\[IndentingNewLine]", \(\(q0 = 1;\)\), "\[IndentingNewLine]", \(\(T = 8 \[Pi];\)\), "\[IndentingNewLine]", \(qmodel[t_, \[Lambda]_, a_]\ = \ q[t]\ /. \ First[\ DSolve[{\(q'\)[t]\ \ \[Equal] \ fmodel[t, q, \[Lambda], a], q[0] \[Equal] q0}, q[t], t]\ ]\)}], "Input"], Cell[BoxData[ \(8\ \[Pi]\)], "Output"], Cell[BoxData[ \(\(\[ExponentialE]\^\(\(-t\)\ \[Lambda]\)\ \((1 + a + \[Lambda]\^2 - a\ \ \[ExponentialE]\^\(t\ \[Lambda]\)\ Cos[t] + a\ \[ExponentialE]\^\(t\ \ \[Lambda]\)\ \[Lambda]\ Sin[t])\)\)\/\(1 + \[Lambda]\^2\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[ Evaluate[ Table[qmodel[t, \[Lambda], 1], {\[Lambda], \(-0.05\), 0.05, 0.05}]], {t, 0, T}, PlotRange \[Rule] All, PlotStyle \[Rule] {Thickness[0.007]}];\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.037894 0.0484518 0.0786127 [ [.21328 .03595 -3 -9 ] [.21328 .03595 3 0 ] [.40275 .03595 -6 -9 ] [.40275 .03595 6 0 ] [.59222 .03595 -6 -9 ] [.59222 .03595 6 0 ] [.78169 .03595 -6 -9 ] [.78169 .03595 6 0 ] [.97116 .03595 -6 -9 ] [.97116 .03595 6 0 ] [.01131 .12706 -6 -4.5 ] [.01131 .12706 0 4.5 ] [.01131 .20568 -6 -4.5 ] [.01131 .20568 0 4.5 ] [.01131 .28429 -6 -4.5 ] [.01131 .28429 0 4.5 ] [.01131 .3629 -6 -4.5 ] [.01131 .3629 0 4.5 ] [.01131 .44152 -6 -4.5 ] [.01131 .44152 0 4.5 ] [.01131 .52013 -6 -4.5 ] [.01131 .52013 0 4.5 ] [.01131 .59874 -6 -4.5 ] [.01131 .59874 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21328 .04845 m .21328 .0547 L s [(5)] .21328 .03595 0 1 Mshowa .40275 .04845 m .40275 .0547 L s [(10)] .40275 .03595 0 1 Mshowa .59222 .04845 m .59222 .0547 L s [(15)] .59222 .03595 0 1 Mshowa .78169 .04845 m .78169 .0547 L s [(20)] .78169 .03595 0 1 Mshowa .97116 .04845 m .97116 .0547 L s [(25)] .97116 .03595 0 1 Mshowa .125 Mabswid .0617 .04845 m .0617 .0522 L s .0996 .04845 m .0996 .0522 L s .13749 .04845 m .13749 .0522 L s .17539 .04845 m .17539 .0522 L s .25117 .04845 m .25117 .0522 L s .28907 .04845 m .28907 .0522 L s .32696 .04845 m .32696 .0522 L s .36486 .04845 m .36486 .0522 L s .44064 .04845 m .44064 .0522 L s .47854 .04845 m .47854 .0522 L s .51643 .04845 m .51643 .0522 L s .55433 .04845 m .55433 .0522 L s .63011 .04845 m .63011 .0522 L s .66801 .04845 m .66801 .0522 L s .7059 .04845 m .7059 .0522 L s .7438 .04845 m .7438 .0522 L s .81958 .04845 m .81958 .0522 L s .85748 .04845 m .85748 .0522 L s .89537 .04845 m .89537 .0522 L s .93327 .04845 m .93327 .0522 L s .25 Mabswid 0 .04845 m 1 .04845 L s .02381 .12706 m .03006 .12706 L s [(1)] .01131 .12706 1 0 Mshowa .02381 .20568 m .03006 .20568 L s [(2)] .01131 .20568 1 0 Mshowa .02381 .28429 m .03006 .28429 L s [(3)] .01131 .28429 1 0 Mshowa .02381 .3629 m .03006 .3629 L s [(4)] .01131 .3629 1 0 Mshowa .02381 .44152 m .03006 .44152 L s [(5)] .01131 .44152 1 0 Mshowa .02381 .52013 m .03006 .52013 L s [(6)] .01131 .52013 1 0 Mshowa .02381 .59874 m .03006 .59874 L s [(7)] .01131 .59874 1 0 Mshowa .125 Mabswid .02381 .06417 m .02756 .06417 L s .02381 .0799 m .02756 .0799 L s .02381 .09562 m .02756 .09562 L s .02381 .11134 m .02756 .11134 L s .02381 .14279 m .02756 .14279 L s .02381 .15851 m .02756 .15851 L s .02381 .17423 m .02756 .17423 L s .02381 .18995 m .02756 .18995 L s .02381 .2214 m .02756 .2214 L s .02381 .23712 m .02756 .23712 L s .02381 .25284 m .02756 .25284 L s .02381 .26857 m .02756 .26857 L s .02381 .30001 m .02756 .30001 L s .02381 .31573 m .02756 .31573 L s .02381 .33146 m .02756 .33146 L s .02381 .34718 m .02756 .34718 L s .02381 .37863 m .02756 .37863 L s .02381 .39435 m .02756 .39435 L s .02381 .41007 m .02756 .41007 L s .02381 .42579 m .02756 .42579 L s .02381 .45724 m .02756 .45724 L s .02381 .47296 m .02756 .47296 L s .02381 .48868 m .02756 .48868 L s .02381 .50441 m .02756 .50441 L s .02381 .53585 m .02756 .53585 L s .02381 .55157 m .02756 .55157 L s .02381 .5673 m .02756 .5673 L s .02381 .58302 m .02756 .58302 L s .02381 .03273 m .02756 .03273 L s .02381 .01701 m .02756 .01701 L s .02381 .00128 m .02756 .00128 L s .02381 .61446 m .02756 .61446 L s .25 Mabswid .02381 0 m .02381 .61803 L s .007 w .02381 .12706 m .02499 .12722 L .02605 .12744 L .02729 .12776 L .02846 .12814 L .03053 .129 L .03279 .1302 L .0379 .13394 L .04262 .1386 L .05224 .15145 L .06244 .16928 L .0842 .21638 L .10458 .26152 L .11392 .27925 L .12281 .29311 L .1277 .29926 L .13305 .30463 L .1354 .30653 L .1379 .30823 L .14005 .30942 L .14241 .31045 L .14376 .3109 L .14502 .31123 L .14635 .31148 L .14709 .31158 L .14777 .31165 L .14898 .31171 L .1503 .31168 L .15155 .31156 L .1527 .31138 L .15402 .31109 L .15526 .31073 L .1581 .30962 L .16085 .30814 L .16383 .30612 L .16895 .30166 L .17435 .29574 L .184 .28248 L .22211 .21685 L .23219 .20242 L .23715 .19665 L .24172 .19223 L .24698 .18834 L .24976 .18684 L .25117 .18622 L .25269 .18568 L .25397 .18531 L .25515 .18506 L .25632 .18487 L .25744 .18477 L Mistroke .25874 .18473 L .25994 .18478 L .26126 .18492 L .262 .18504 L .26267 .18518 L .265 .18585 L .26632 .18636 L .26756 .18692 L .27272 .19018 L .27718 .19415 L .28195 .19953 L .29286 .21586 L .30279 .23486 L .32352 .28207 L .34268 .32601 L .35229 .34502 L .36131 .35974 L .36621 .36623 L .3707 .37115 L .37569 .37539 L .37831 .37709 L .38106 .37847 L .3823 .37897 L .38365 .3794 L .38491 .37973 L .38608 .37995 L .38714 .38008 L .38829 .38016 L .3895 .38017 L .39063 .38011 L .39196 .37994 L .39321 .37971 L .39434 .37942 L .39556 .37904 L .3983 .37791 L .40085 .37652 L .40508 .37353 L .40975 .36929 L .41923 .35807 L .44017 .325 L .4594 .29334 L .46914 .28008 L .47822 .27077 L .48279 .26744 L .48536 .26601 L .48769 .26499 L .48893 .26456 L .49024 .2642 L .49093 .26404 L .49168 .2639 L Mistroke .49237 .2638 L .49302 .26373 L .49428 .26365 L .49544 .26365 L .49608 .26369 L .49678 .26375 L .49805 .26393 L .4987 .26406 L .49938 .26421 L .50061 .26456 L .50337 .26565 L .50576 .26694 L .50838 .2687 L .51322 .27294 L .51774 .27802 L .52787 .29304 L .53872 .31395 L .55854 .36048 L .5801 .41276 L .59014 .43395 L .60081 .45237 L .60992 .46395 L .6151 .46865 L .61762 .47043 L .61996 .47177 L .62224 .47278 L .6234 .47319 L .62466 .47356 L .62539 .47373 L .62605 .47386 L .62736 .47405 L .62855 .47414 L .6298 .47416 L .63103 .4741 L .63235 .47394 L .6336 .47371 L .63474 .47343 L .63752 .47248 L .64016 .47122 L .64476 .46826 L .64902 .46468 L .65867 .45405 L .67752 .42658 L .69787 .3958 L .70332 .38883 L .70906 .38251 L .71438 .37777 L .71698 .37589 L .71934 .37445 L .72157 .37333 L Mistroke .72399 .37239 L .72532 .372 L .72655 .37172 L .72773 .37153 L .72897 .3714 L .73022 .37136 L .73092 .37137 L .73156 .3714 L .73281 .37153 L .73399 .37173 L .73536 .37206 L .73664 .37246 L .73953 .37369 L .74217 .37523 L .7446 .37698 L .75009 .38213 L .75493 .388 L .76004 .3955 L .77864 .4325 L .8166 .52965 L .8361 .57165 L .84172 .58094 L .84704 .58835 L .85178 .59374 L .8543 .59612 L .85701 .59831 L .85961 .60003 L .86199 .60129 L .86427 .60221 L .86541 .60256 L .86666 .60288 L .86735 .60301 L .86801 .60312 L .86923 .60326 L .87055 .60332 L .87129 .60331 L .87197 .60329 L .87312 .60319 L .87437 .603 L .87569 .60272 L .87692 .60239 L .87929 .60153 L .88155 .60048 L .88579 .59789 L .89036 .5943 L .89533 .58952 L .91668 .56241 L .93602 .5368 L .94572 .52687 L .95118 .52274 L Mistroke .95377 .5212 L .9562 .52004 L .95832 .51925 L .95942 .51893 L .9606 .51865 L .9619 .51842 L .96313 .51828 L .96433 .51822 L .96544 .51823 L .9666 .51831 L .96769 .51845 L .96888 .51868 L .97017 .51901 L .97279 .51997 L .97519 .52117 L .97619 .52177 L Mfstroke .02381 .12706 m .02499 .1271 L .02605 .1272 L .02729 .1274 L .02846 .12765 L .03053 .1283 L .03279 .12926 L .03527 .13063 L .0379 .13244 L .04262 .13655 L .05205 .14791 L .06244 .1645 L .08274 .20445 L .09394 .22739 L .10458 .24748 L .11464 .26348 L .11995 .27036 L .12562 .27629 L .12831 .27857 L .13085 .28038 L .13333 .28182 L .1356 .28285 L .13777 .28358 L .13886 .28385 L .14005 .28407 L .14125 .28422 L .14254 .28429 L .14376 .28427 L .14487 .28418 L .14605 .28401 L .14715 .28379 L .14837 .28346 L .14967 .28302 L .15234 .28184 L .15478 .28043 L .15925 .27705 L .16395 .27242 L .17387 .2594 L .18471 .24102 L .20596 .19825 L .22527 .16102 L .23491 .14618 L .23999 .13985 L .24536 .13444 L .25006 .13088 L .25273 .12936 L .25401 .12877 L .25522 .12828 L .25632 .12792 L .25751 .12759 L Mistroke .2586 .12736 L .25961 .12721 L .26081 .1271 L .26189 .12706 L .26314 .12711 L .26432 .12722 L .26549 .12742 L .26674 .1277 L .26899 .12843 L .27167 .12966 L .27411 .13111 L .27942 .13531 L .28444 .14056 L .2933 .15253 L .30292 .16879 L .32234 .20758 L .33266 .22863 L .34376 .24936 L .35435 .2657 L .35898 .27145 L .36405 .2766 L .36637 .27854 L .36886 .28032 L .371 .28159 L .37332 .2827 L .37469 .28322 L .37595 .28361 L .37716 .2839 L .37844 .28412 L .37958 .28424 L .38065 .28429 L .38182 .28427 L .38307 .28417 L .38429 .28399 L .38498 .28385 L .38563 .28369 L .38835 .2828 L .3907 .2817 L .39323 .2802 L .399 .27554 L .40429 .26985 L .41408 .25611 L .4244 .23803 L .44435 .19766 L .46248 .16255 L .47187 .14774 L .47682 .14132 L .48207 .1357 L .48473 .13336 L .48767 .13119 L Mistroke .49037 .12959 L .49283 .12847 L .49405 .12803 L .49537 .12765 L .49609 .12748 L .49676 .12735 L .49745 .12724 L .49809 .12716 L .49934 .12708 L .50047 .12707 L .50169 .12714 L .50301 .12731 L .50425 .12756 L .50538 .12786 L .50793 .12878 L .51029 .12994 L .51252 .13132 L .51752 .13532 L .52296 .14106 L .53391 .15649 L .54403 .17443 L .5643 .21556 L .57422 .23541 L .58332 .25186 L .592 .2651 L .5966 .27089 L .6015 .27601 L .60409 .27825 L .60644 .27998 L .60884 .28145 L .61104 .28254 L .61352 .28346 L .61483 .2838 L .6155 .28394 L .61621 .28407 L .61741 .28422 L .61853 .28428 L .61922 .28429 L .61987 .28427 L .62109 .28418 L .62234 .28399 L .62304 .28385 L .62369 .2837 L .62614 .28292 L .62836 .28193 L .63071 .2806 L .6356 .27691 L .64095 .27152 L .65048 .25877 L .65935 .24386 L Mistroke .69833 .16651 L .70854 .14978 L .71392 .14252 L .71963 .13621 L .72461 .13199 L .72744 .13015 L .7301 .12881 L .73135 .12831 L .73253 .12791 L .73361 .12762 L .73476 .12737 L .736 .12718 L .73715 .12709 L .73786 .12707 L .73852 .12707 L .73977 .12714 L .74093 .12728 L .74217 .12752 L .74335 .12782 L .74442 .12816 L .74704 .12924 L .74952 .13061 L .75446 .13428 L .75984 .13965 L .76943 .15244 L .77831 .16732 L .79855 .20761 L .81756 .24516 L .8264 .25981 L .83572 .27206 L .84059 .27691 L .84317 .27901 L .84593 .28087 L .84847 .28224 L .85078 .28319 L .8519 .28354 L .8531 .28384 L .85425 .28406 L .8553 .2842 L .85655 .28428 L .85774 .28428 L .8588 .28421 L .85997 .28407 L .86123 .28383 L .86257 .28348 L .86499 .28261 L .86744 .28141 L .86974 .27998 L .87405 .27659 L .87922 .27132 L Mistroke .88397 .2654 L .89469 .24876 L .91514 .20884 L .93422 .17056 L .94365 .15432 L .95257 .14185 L .9575 .13643 L .9601 .13405 L .96286 .13188 L .96517 .13036 L .96773 .12902 L .96885 .12854 L .97005 .12809 L .97109 .12778 L .97223 .12749 L .97323 .1273 L .97416 .12718 L .97514 .12709 L .97619 .12706 L Mfstroke .02381 .12706 m .02499 .12698 L .02605 .12697 L .02729 .12704 L .02846 .12717 L .02954 .12737 L .03053 .1276 L .03279 .12833 L .03506 .12933 L .03756 .13074 L .04262 .13454 L .04755 .13941 L .05285 .14581 L .06244 .15995 L .08405 .19888 L .09469 .21836 L .10458 .23463 L .11412 .24747 L .11945 .25306 L .12429 .25702 L .12706 .25877 L .12961 .26004 L .13106 .26061 L .13241 .26104 L .13382 .26139 L .13457 .26153 L .13537 .26165 L .13607 .26172 L .13671 .26177 L .13742 .26179 L .13817 .26179 L .13943 .26171 L .14015 .26163 L .14083 .26153 L .14208 .26127 L .14322 .26097 L .1458 .26001 L .14829 .25875 L .15054 .25732 L .15561 .25314 L .1603 .24812 L .16473 .24241 L .18543 .20571 L .20375 .16555 L .21416 .1429 L .22389 .12357 L .23428 .10631 L .24377 .09451 L .24626 .09214 L .24897 .08991 L Mistroke .25157 .08815 L .25392 .08686 L .25665 .08574 L .25816 .08529 L .25884 .08513 L .25957 .08499 L .26088 .08481 L .26209 .08472 L .26339 .08472 L .26413 .08476 L .26481 .08482 L .26608 .08501 L .26672 .08513 L .26742 .0853 L .26984 .08606 L .27119 .08663 L .27265 .08734 L .27529 .08892 L .28047 .09302 L .28521 .09788 L .29426 .10978 L .30377 .12526 L .32378 .16293 L .34193 .19541 L .35205 .20968 L .35668 .21488 L .36158 .21932 L .36408 .22114 L .36682 .22278 L .36932 .22393 L .37044 .22433 L .37165 .2247 L .37238 .22489 L .37306 .22503 L .3737 .22515 L .3744 .22525 L .37561 .22536 L .37691 .22539 L .37764 .22536 L .37831 .22532 L .37904 .22524 L .37982 .22512 L .38112 .22485 L .38255 .22445 L .38508 .22347 L .38782 .22201 L .39281 .21836 L .39766 .21358 L .40213 .20817 L .42308 .17241 L Mistroke .44341 .12875 L .46245 .09024 L .4718 .07506 L .48069 .06402 L .48559 .05956 L .48818 .05771 L .49093 .05612 L .49225 .05551 L .49348 .05502 L .49469 .05462 L .49579 .05432 L .49685 .0541 L .49796 .05394 L .49908 .05384 L .5003 .05381 L .50161 .05387 L .50282 .05401 L .50345 .05411 L .50414 .05425 L .50555 .05461 L .50811 .05554 L .51052 .05673 L .515 .05974 L .51983 .0641 L .53074 .07774 L .54075 .094 L .56129 .13328 L .5713 .15241 L .58039 .16817 L .58924 .18108 L .5939 .18665 L .59883 .1915 L .60134 .19352 L .6041 .19538 L .60664 .19674 L .60896 .19769 L .61018 .19808 L .61134 .19837 L .61257 .1986 L .61326 .1987 L .6139 .19876 L .61507 .19882 L .61614 .19881 L .61737 .19872 L .61852 .19856 L .61969 .19833 L .62094 .19799 L .6232 .19718 L .62588 .19585 L .62832 .19431 L Mistroke .63363 .18989 L .63865 .18443 L .64752 .17204 L .65714 .15525 L .67658 .11502 L .698 .07109 L .70778 .05469 L .71317 .04723 L .7183 .04134 L .72305 .03703 L .72556 .03522 L .72824 .03365 L .73071 .03255 L .73296 .03183 L .73404 .03158 L .73521 .03139 L .73632 .03127 L .73733 .03122 L .7386 .03124 L .73979 .03134 L .74086 .0315 L .74203 .03174 L .7433 .03208 L .74464 .03254 L .74706 .03361 L .75138 .03628 L .75615 .04032 L .7668 .05309 L .77681 .06901 L .81548 .14281 L .82586 .15901 L .83552 .17037 L .84039 .17447 L .84311 .17624 L .84559 .17753 L .84691 .17808 L .84831 .17856 L .84904 .17877 L .84984 .17897 L .85058 .17912 L .85126 .17923 L .85194 .17932 L .85259 .17938 L .85328 .17942 L .85402 .17943 L .85467 .17942 L .85537 .17938 L .85661 .17925 L .85786 .17902 L .85923 .17868 L Mistroke .86065 .17822 L .86198 .17769 L .86692 .1749 L .86973 .17273 L .8723 .17041 L .87806 .16404 L .88842 .14893 L .89823 .13118 L .9366 .05277 L .94651 .03659 L .95581 .02501 L .96048 .02073 L .96288 .01896 L .9655 .01737 L .96822 .01611 L .9697 .01559 L .97108 .01521 L .97238 .01496 L .97355 .0148 L .9742 .01475 L .97491 .01472 L .97619 .01473 L Mfstroke 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{562.063, 347.375}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg7oo1000@Goo1000;goo000=Ool0 0`00Oomoo`1dOol4000GOol4001LOol4000/Ool4001;Ool4000iOol30011Ool3000`Ool000eoo`03 001oogoo07Aoo`<001Uoo`@005]oo`<002ioo`<004]oo`<003Yoo`@003moo`@0031oo`003Goo00<0 07ooOol0Lgoo10006Woo0`00FWoo1000;Woo1000BGoo1000>goo0`00?goo0`00goo10007oo0`00=7oo000=Ool00`00Oomoo`1^Ool4000T Ool3001AOol4000fOol30013Ool30013Ool3000gOol4000dOol000eoo`03001oogoo06ioo`<002Eo o`@0051oo`<003Moo`@0045oo`@004=oo`@003Ioo`<003Eoo`003Goo00<007ooOol0KGoo10009Woo 0`00Cgoo1000>7oo0`00@Goo0`00AGoo0`00=Goo1000=Goo000=Ool00`00Oomoo`1/Ool4000WOol4 001>Ool3000iOol4000oOol40015Ool4000dOol3000fOol00003Ool0000000<000Moo`03001oogoo 06aoo`<002Uoo`@004aoo`@003Yoo`<003moo`<004Moo`<003=oo`@003Ioo`000goo00<007ooOol0 1Woo1P00JGoo10004Goo1@0057oo0`00C7oo0`006Woo1`006Woo1000?Goo100087oo1@008Woo1000 goo0`0077oo2@0000=oo`0000001P007Goo 1000;Woo100087oo2@003goo0002Ool20009Ool01000Oomoogoo1`00HGoo10002goo1P002Goo1@00 3goo0`00B7oo0`005Goo1P002Goo1@005Woo1000>Goo10006goo1@0027oo1`007Goo0`00;Woo0`00 87oo1@0057oo000=Ool00`00Oomoo`03Ool6001POol3000Ool40016Ool4000D Ool6000;Ool5000FOol3000iOol3000JOol6000=Ool4000LOol4000/Ool4000OOol5000EOol000eo o`03001oogoo00Aoo`H005moo`<000]oo`@000moo`@000ioo`<004Ioo`<001Aoo`D000moo`@001Eo o`@003Moo`@001Yoo`@0011oo`@001aoo`<002aoo`<001moo`@001Moo`003Goo00<007ooOol01Goo 1P00GGoo10002Woo10004Goo10003Goo1000AGoo0`004goo10004Woo10005Goo0`00=goo0`006Woo 10004Woo10006goo1000:Woo10007goo0`0067oo000=Ool00`00Oomoo`06Ool6001LOol3000:Ool4 000COol4000=Ool30014Ool4000BOol4000DOol4000DOol4000eOol4000IOol4000DOol4000KOol3 000ZOol3000OOol4000HOol000eoo`03001oogoo00Moo`H005Yoo`@000Uoo`@001Eoo`@000aoo`@0 04=oo`<001=oo`<001Ioo`@001Aoo`<003Eoo`<001Uoo`@001Ioo`@001Yoo`@002Qoo`@001ioo`@0 01Uoo`003Goo0`0027oo1P00FGoo0`002Goo10005goo100037oo0`00@Woo10004Woo0`0067oo1000 4goo1000Goo100047oo10008goo10004goo 0`00:Woo0`005goo100097oo0`006goo0`0077oo0`007Woo100087oo000=Ool00`00Oomoo`0?Ool6 001=Ool40007Ool3000UOol4000:Ool4000hOol3000@Ool4000UOol3000COol4000XOol4000GOol3 000UOol4000JOol4000JOol4000NOol3000QOol000eoo`03001oogoo00moo`L004aoo`<000Moo`@0 02Ioo`@000Yoo`<003Moo`@0011oo`<002Ioo`@001=oo`<002Qoo`<001Moo`@002Ioo`<001]oo`@0 01Qoo`@001ioo`@0025oo`003Goo00<007ooOol047oo1P00Bgoo10001goo0`00:7oo0`002Woo1000 =Woo0`0047oo10009goo10004Woo10009Woo10005Woo10009goo10006goo0`0067oo0`007Woo1000 8Woo000=Ool00`00Oomoo`0@Ool7001:Ool30007Ool4000XOol4000:Ool3000fOol3000@Ool3000Y Ool3000COol4000UOol3000GOol3000YOol4000JOol4000FOol3000OOol3000SOol000eoo`03001o ogoo015oo`H004Uoo`@000Moo`<002Yoo`@000Uoo`@003Aoo`@000moo`@002Uoo`@001=oo`<002Ao o`@001Ioo`@002Yoo`<001]oo`@001Aoo`@001ioo`@002=oo`003Goo00<007ooOol04Woo1P00B7oo 0`001goo1000:goo0`002Woo0`00=7oo0`0047oo0`00:goo0`004goo10008goo0`005goo0`00:goo 10006goo10004Woo10007Woo100097oo000=Ool00`00Oomoo`0BOol70016Ool40007Ool3000/Ool4 0009Ool4000bOol4000?Ool4000[Ool4000COol4000QOol4000FOol4000/Ool3000LOol4000?Ool5 000OOol3000UOol000eoo`03001oogoo01=oo`H004Ioo`<000Moo`@002eoo`<000Yoo`<0039oo`<0 011oo`<002eoo`@001=oo`<0025oo`<001Moo`<002eoo`@001]oo`D000eoo`D001moo`@002Eoo`00 3Goo0`004goo1`00A7oo10001goo0`00;Woo10002Goo1000<7oo10003goo1000;Woo0`004goo1000 7goo0`005goo1000;Woo0`007Goo1@001goo00<007ooOol01@0087oo0`009Woo000=Ool00`00Oomo o`0DOol60014Ool30007Ool4000_Ool3000:Ool3000_Ool4000@Ool3000_Ool4000COol3000NOol4 000GOol3000_Ool4000MOol70003Ool7000QOol4000VOol000eoo`03001oogoo01Aoo`L004=oo`<0 00Ioo`@0031oo`@000Uoo`@002ioo`<0011oo`@0031oo`<001=oo`@001aoo`@001Moo`@0031oo`<0 01ioo`l0029oo`<002Moo`003Goo00<007ooOol05Goo1P00@Woo10001Woo0`00goo10001Woo0`003goo1P00 2goo10004goo10002Woo10009Woo0`004Goo0`00=goo10005Goo10003goo1@006Woo0`00=goo1000 B7oo0`00:Woo000=Ool00`00Oomoo`0HOol7000kOol30006Ool4000>Ool5000>Ool4000COol3000; Ool3000UOol4000@Ool4000hOol3000FOol4000=Ool5000JOol4000hOol30017Ool4000ZOol000eo o`03001oogoo01Qoo`P003Uoo`@000Ioo`<000ioo`D0011oo`@0019oo`@000Yoo`@002=oo`@0015o o`<003Uoo`<001Moo`H000Uoo`D001]oo`<003Uoo`@004Ioo`<002]oo`003Goo00<007ooOol06Goo 1`00>Goo0`001Woo10003Woo0`004goo10004Woo0`002goo0`008goo0`004Goo1000>Goo10005goo 4P006goo1000>Woo0`00AGoo1000:goo000=Ool00`00Oomoo`0IOol8000hOol30006Ool3000>Ool3 000EOol4000AOol4000:Ool4000QOol4000AOol3000kOol3000HOol?000MOol3000kOol40014Ool3 000/Ool000eoo`03001oogoo01Yoo`L003Moo`@000Eoo`@000eoo`@001Ioo`@0015oo`<000]oo`@0 021oo`<0015oo`@003]oo`@001Yoo`/001eoo`@003aoo`<004=oo`@002aoo`0000=oo`0000000P00 27oo00<007ooOol06Woo2000=Woo0`001Woo0`003Goo100067oo100047oo10002goo0`007goo1000 4Goo0`00?Goo0`007Goo0`008Woo0`00?Goo1000@Woo0`00;Goo00000goo001oo`0:Ool00`00Oomo o`0KOol7000eOol40005Ool4000=Ool3000JOol4000@Ool3000Woo1000A7oo0`00>goo1000<7oo000=Ool00`00Oomoo`0NOol8000]Ool40006Ool3 000;Ool4000TOol3000?Ool4000=Ool4000BOol4000COol30015Ool3000jOol30015Ool4000jOol3 000aOol000eoo`03001oogoo01ioo`T002aoo`<000Ioo`@000]oo`<002Eoo`@000moo`<000ioo`D0 00moo`@001=oo`@004Eoo`@003Qoo`@004Ioo`<003Uoo`@0035oo`003Goo00<007ooOol07goo2000 :goo10001Woo0`002goo10009Woo0`003goo10003Woo1@003Goo100057oo0`00Agoo0`00>7oo0`00 Agoo1000>7oo0`00Ool4000BOol:000HOol4001:Ool4000cOol3001;Ool4 000dOol3000dOol000eoo`03001oogoo025oo`T002Eoo`@000Ioo`@000Yoo`@002aoo`<000moo`<0 03Aoo`<004aoo`<0039oo`@004aoo`<003=oo`@003Aoo`003Goo00<007ooOol08Goo2P0097oo0`00 1goo0`002goo0`00;Goo10003Woo1000Ool4000`Ool4001>Ool4000_Ool3001?Ool4000_Ool4000fOol000eoo`03 001oogoo029oo`H00004Ool00000000QOol40006Ool4000:Ool3000aOol3000?Ool3000`Ool3001@ Ool3000^Ool4001@Ool3000_Ool3000gOol000eoo`<002=oo`X001moo`@000Moo`<000Yoo`@0035o o`@000ioo`@002ioo`@0051oo`@002eoo`<0055oo`@002eoo`@003Moo`003Goo00<007ooOol08goo 2`007Woo0`001goo10002Woo0`007oo000=Ool00`00Oomoo`0TOol600000goo00000002000LOol40007Ool3000:Ool4000cOol3 000?Ool4000/Ool4001BOol4000[Ool3001COol4000[Ool4000hOol000eoo`03001oogoo02Aoo`L0 0003Ool00000008001]oo`<000Moo`@000Yoo`<003Aoo`@000moo`<002]oo`@005Aoo`<002Yoo`@0 05Aoo`<002]oo`<003Uoo`003Goo00<007ooOol09Goo1P000Woo10006Goo10001goo0`002Woo1000 =Goo0`003goo1000:Woo0`00EGoo1000:Goo0`00EGoo1000:Goo1000>Goo000=Ool00`00Oomoo`0U Ool70002Ool3000HOol40007Ool40009Ool4000fOol4000?Ool3000YOol4001FOol4000WOol4001F Ool3000YOol3000jOol000eoo`03001oogoo02Eoo`L000=oo`<001Ioo`@000Qoo`<000Yoo`<003Qo o`<000moo`@002Qoo`<005Qoo`<002Ioo`@005Moo`@002Moo`@003Yoo`003Goo00<007ooOol09Woo 1`000Woo100057oo100027oo10002Goo1000>7oo10003goo0`009goo1000F7oo10009Goo0`00FGoo 10009Woo0`00>goo000=Ool00`00Oomoo`0VOol80002Ool4000BOol40008Ool4000:Ool3000jOol3 000?Ool4000UOol4001JOol3000TOol4001JOol3000UOol4000kOol000eoo`<002Moo`<00004Ool0 00000003Ool4000@Ool40009Ool3000:Ool4000jOol3000@Ool4000TOol3001KOol4000SOol3001K Ool4000SOol4000lOol000eoo`03001oogoo02Moo`P000=oo`D000aoo`D000Uoo`@000Uoo`@003]o o`@0011oo`<002=oo`@005aoo`@0025oo`@005aoo`<002=oo`<003eoo`003Goo00<007ooOol0:7oo 0`0000=oo`0000000P000goo1@002Woo1@002Woo0`002Woo0`00?Goo0`0047oo10008Woo0`00GWoo 0`0087oo1000GGoo10008Goo1000?Goo000=Ool00`00Oomoo`0XOol30002Ool30004Ool70004Ool7 000:Ool40009Ool4000mOol4000@Ool4000POol3001OOol4000OOol3001OOol4000OOol4000nOol0 00eoo`03001oogoo02Qoo`@00003Ool00000008000Eoo`h000aoo`<000Yoo`<003moo`<0015oo`<0 01moo`@0061oo`<001ioo`@0061oo`<001moo`<003moo`003Goo00<007ooOol0:Goo0`000Woo1000 1Goo300037oo10002Goo1000?goo100047oo10007Goo1000:Goo2@00;goo100077oo1000HGoo1000 7Goo1000?goo000=Ool00`00Oomoo`0YOol40002Ool30008Ool6000>Ool4000:Ool30011Ool3000A Ool4000KOol4000XOol<000_Ool4000JOol4001SOol4000KOol40010Ool000eoo`03001oogoo02Yo o`<0009oo`@001Yoo`@000Yoo`@0045oo`<0019oo`<001Yoo`@002Mooa0002ioo`<001Yoo`<006Eo o`<001]oo`<0045oo`003Goo00<007ooOol0:Woo10000Woo10006Goo0`002Woo1000@Woo10004Goo 100067oo10009goo1P001goo1@00;Goo100067oo1000IGoo10006Goo1000@Goo000=Ool3000[Ool4 0002Ool3000HOol4000:Ool30014Ool3000BOol4000GOol3000WOol5000:Ool5000]Ool4000FOol4 001WOol4000GOol40012Ool000eoo`03001oogoo02aoo`<000=oo`<001Ioo`@000Yoo`@004Aoo`@0 019oo`@001Eoo`<002Moo`@000ioo`@002eoo`<001Eoo`@006Uoo`@001Eoo`@004=oo`003Goo00<0 07ooOol0;7oo10000Woo100057oo10002goo0`00AWoo0`004goo10004goo10009Woo100047oo1000 ;7oo10004goo1000Jgoo10004goo1@00@goo000=Ool00`00Oomoo`0]Ool30003Ool5000@Ool6000: Ool40016Ool4000COol4000AOol4000VOol4000BOol4000/Ool4000AOol4001/Ool5000AOol40015 Ool000eoo`03001oogoo02eoo`@000=oo`D000ioo`D000]oo`@004Qoo`<001Aoo`D000ioo`@002Io o`@001Aoo`@002aoo`D000ioo`@006moo`@000moo`@004Ioo`003Goo00<007ooOol0;Woo10000goo 1P002goo1@0037oo0`00BGoo0`005Goo1@002goo1@009Woo10005Woo1000;7oo1@002goo1@00LGoo 100037oo1@00Agoo0002Ool20009Ool00`00Oomoo`0_Ool30005Ool50008Ool5000=Ool40019Ool4 000EOol50009Ool5000WOol3000HOol4000/Ool70007Ool5001cOol60008Ool50018Ool00005Ool0 07ooOol00008Ool00`00Oomoo`0_Ool40005Ool@000>Ool3001;Ool3000GOol@000WOol3000JOol4 000]Ool@001eOolA0019Ool000=oo`03001oogoo00Moo`D002ioo`@000Ioo`d000ioo`@004]oo`@0 01Moo`d002Qoo`@001]oo`<002ioo`d007Qoo`h004]oo`0017oo00<007ooOol01Woo00<007ooOol0 Ool4001JOol4 000oOol4000ZOol4003oOol5Ool000eoo`03001oogoo03Uoo`D000]oo`D005aoo`<003moo`<002ao o`<00?moo`Eoo`003Goo00<007ooOol0>goo1P0017oo2@00G7oo1000?Goo1000;7oo1000ogoo17oo 000=Ool00`00Oomoo`0lOol@001OOol3000lOol4000^Ool3003oOol4Ool000eoo`03001oogoo03eo o`d0065oo`@003]oo`<002moo`@00?moo`=oo`003Goo00<007ooOol0@7oo1P00IWoo0`00>Woo1000 <7oo0`00ogoo0goo000=Ool3002/Ool4000iOol3000aOol4003oOol2Ool000eoo`03001oogoo0:eo o`<003Qoo`@0039oo`<00?moo`9oo`003Goo00<007ooOol0[Goo0`00=goo1000Woo0`00oWoo000=Ool3002`Ool4000_Ool4000jOol3 003nOol000eoo`03001oogoo0;5oo`<002moo`<003]oo`@00?eoo`003Goo00<007ooOol0/Goo1000 ;Goo1000?7oo0`00oGoo000=Ool00`00Oomoo`2bOol3000/Ool4000mOol4003lOol000eoo`03001o ogoo0;9oo`@002]oo`<003moo`<00?aoo`003Goo00<007ooOol0/goo0`00:Woo1000?goo1000ngoo 000=Ool00`00Oomoo`2cOol4000XOol40011Ool3003kOol000eoo`03001oogoo0;Aoo`<002Qoo`<0 049oo`<00?]oo`003Goo00<007ooOol0]7oo10009Woo1000@Woo1000nWoo000=Ool3002eOol3000U Ool40014Ool3003jOol000eoo`03001oogoo0;Eoo`@002Aoo`<004Eoo`@00?Uoo`003Goo00<007oo Ool0]Woo0`008goo1000AWoo0`00nGoo000=Ool00`00Oomoo`2fOol4000ROol30017Ool4003hOol0 00eoo`03001oogoo0;Moo`@0021oo`@004Qoo`<00?Qoo`003Goo00<007ooOol0^7oo0`007goo1000 BGoo0`00n7oo000=Ool00`00Oomoo`2hOol4000MOol4001:Ool4003gOol000=oo`<000Moo`03001o ogoo0;Uoo`<001aoo`@004aoo`<00?Moo`0017oo00<007ooOol01Woo00<007ooOol0^Goo10006Woo 1000CGoo1000mWoo00000goo000000030007Ool5002hOol4000HOol4001?Ool3003fOol00005Ool0 07ooOol00008Ool00`00Oomoo`2kOol4000FOol4001@Ool3003fOol0009oo`03001oo`0000Qoo`03 001oogoo0;aoo`@001Aoo`@0055oo`@00?Eoo`000goo0P0027oo00<007ooOol0_Goo10004Woo1000 Dgoo0`00Egoo00<007oo00001@00UWoo000=Ool00`00Oomoo`2nOol4000@Ool4001DOol4001EOol; 002DOol000eoo`03001oogoo0;moo`@000ioo`@005Ioo`<005Aoo`l0095oo`003Goo00<007ooOol0 `7oo1@002Woo1@00Egoo1000DWoo1P0017oo1`00T7oo000=Ool00`00Oomoo`31Ool50008Ool5001I Ool3001@Ool50009Ool6002?Ool000eoo`<00<5ooa4005Yoo`<004moo`D000eoo`@008ioo`003Goo 00<007ooOol0a7oo3000G7oo1000CGoo1@003goo1000SGoo000=Ool00`00Oomoo`35Ool:001NOol3 001Goo0`00:Woo1000PGoo000=Ool00`00 Oomoo`3oOoleOol4000gOol4000[Ool30021Ool000eoo`03001oogoo0?moocIoo`<003Moo`<002ao o`<0085oo`003Goo00<007ooOol0ogoo=Woo1000=Goo1000;7oo1000P7oo000=Ool00`00Oomoo`3o OolgOol3000dOol4000^Ool30020Ool000eoo`03001oogoo0?moocMoo`@003=oo`<002moo`@007mo o`003Goo00<007ooOol0ogoo>7oo0`00Goo1000;Woo1000goo0`00:goo1000=goo0`00O7oo000=Ool00`00Oomoo`3oOolkOol4000YOol4000hOol4001kOol0 00eoo`03001oogoo0?moocaoo`<002Uoo`<003Yoo`<007]oo`000Woo0P002Goo00<007ooOol0ogoo ?7oo10009goo1000>Woo0`00Ngoo00001Goo001oogoo000027oo00<007ooOol0ogoo?Goo0`009Woo 1000>goo1000NWoo0004Ool00`00Oomoo`06Ool5003oOolkOol4000TOol4000mOol3001jOol0009o o`8000Uoo`03001oogoo0?moocioo`@002=oo`<003ioo`@007Uoo`000Woo00<007ooOol027oo00<0 07ooOol0ogoo?goo0`008Woo1000?goo0`00NGoo0002Ool30008Ool00`00Oomoo`3oOoloOol4000P Ool40010Ool3001iOol000eoo`03001oogoo0?mood1oo`<0021oo`<0045oo`@007Qoo`003Goo00<0 07ooOol0ogoo@7oo10007Woo1000@Woo0`00N7oo000=Ool00`00Oomoo`3oOom1Ool3000MOol40013 Ool3001hOol000eoo`03001oogoo0?mood5oo`@001aoo`<004Aoo`@007Moo`003Goo0`00ogoo@Woo 10006Woo0`00AWoo0`00Mgoo000=Ool00`00Oomoo`3oOom3Ool4000HOol40016Ool4001fOol000eo o`03001oogoo0?moodAoo`@001Ioo`@004Qoo`<007Ioo`003Goo00<007ooOol0ogooAGoo0`005Goo 1000BGoo0`00MWoo000=Ool00`00Oomoo`3oOom6Ool3000COol4001:Ool4001eOol000eoo`03001o ogoo0?moodIoo`@0015oo`@004aoo`<007Eoo`003Goo00<007ooOol0ogooAgoo10003Woo1@00CGoo 1000M7oo000=Ool00`00Oomoo`3oOom8Ool4000;Ool6001?Ool3001dOol000eoo`03001oogoo0?mo odUoo`H000Eoo`P0051oo`<007Aoo`003Goo0`00ogooBWoo4000DWoo1000Lgoo000=Ool00`00Oomo o`3oOom;Ool=001EOol3001cOol000eoo`03001oogoo0?moodioo`L005Qoo`@0079oo`003Goo00<0 07ooOol0ogoo[Woo0`00LWoo000=Ool00`00Oomoo`3oOon^Ool3001bOol000eoo`03001oogoo0?mo ojioo`@0075oo`003Goo00<007ooOol0ogoo[goo0`00LGoo000=Ool00`00Oomoo`3oOon_Ool3001a Ool000eoo`03001oogoo0?moojmoo`@0071oo`003Goo0`00ogoo/7oo0`00L7oo000=Ool00`00Oomo o`3oOon`Ool4001_Ool000eoo`03001oogoo0?mook5oo`<006moo`003Goo00<007ooOol0ogoo/Goo 0`00Kgoo000=Ool00`00Oomoo`3oOonaOol4001^Ool000eoo`03001oogoo0?mook9oo`<006ioo`00 3Goo00<007ooOol0ogoo/Woo1000KGoo000=Ool00`00Oomoo`3oOoncOol3001]Ool000eoo`03001o ogoo0?mook=oo`<006eoo`003Goo0`00ogoo/goo1000K7oo000=Ool00`00Oomoo`3oOondOol3001/ Ool000eoo`03001oogoo0?mookAoo`@006]oo`003Goo00<007ooOol0ogoo]Goo0`00Jgoo000=Ool0 0`00Oomoo`3oOoneOol3001[Ool000eoo`03001oogoo0?mookEoo`@006Yoo`003Goo00<007ooOol0 ogoo]Woo0`00JWoo000017oo000000002Goo00<007ooOol0ogoo]Woo0`00D7oo2@004Goo00001Goo 001oogoo000027oo00<007ooOol0ogoo]Woo1000C7oo3P003goo00001Goo001oogoo000027oo1@00 ogoo]Goo0`00Bgoo40003Woo000017oo000000002Goo00<007ooOol0ogoo]goo1000BGoo1P001goo 1@003Goo0002Ool00`00Oomoo`08Ool00`00Oomoo`3oOonhOol30017Ool6000;Ool3000=Ool0009o o`<000Qoo`03001oogoo0?mookQoo`<004Ioo`D000ioo`03001oogoo00aoo`003Goo00<007ooOol0 ogoo^7oo1000A7oo1@007Woo000=Ool00`00Oomoo`3oOoniOol30013Ool4000POol000eoo`03001o ogoo0?mookUoo`@0045oo`@0025oo`003Goo00<007ooOol0ogoo^Woo0`00@7oo10008Woo000=Ool3 003oOonjOol4000nOol4000SOol000eoo`03001oogoo0?mook]oo`<003ioo`<002Aoo`003Goo00<0 07ooOol0ogoo^goo1000?7oo100097oo000=Ool00`00Oomoo`3oOonlOol3000kOol4000UOol000eo o`03001oogoo0?mookaoo`<003Yoo`@002Ioo`003Goo00<007ooOol0ogoo_7oo1000>Goo0`009goo 000=Ool00`00Oomoo`3oOonmOol3000hOol4000WOol000eoo`03001oogoo0?mookeoo`@003Ioo`@0 02Qoo`003Goo00<007ooOol0ogoo_Woo0`00=Woo0`00:Goo000=Ool3003oOonnOol4000dOol4000Y Ool000eoo`03001oogoo0?mookmoo`<003=oo`@002Yoo`003Goo00<007ooOol0ogoo_goo10007oo000=Ool00`00Oomoo`3oOoo;Ool3000HOol4000i Ool000eoo`03001oogoo0?mool]oo`@001Ioo`@003Yoo`003Goo00<007ooOol0ogooc7oo10004goo 1@00>goo000=Ool00`00Oomoo`3oOoo=Ool4000AOol5000lOol000=oo`03001oogoo00Moo`03001o ogoo0?mooleoo`H000eoo`H003eoo`000goo00<007ooOol01goo00<007ooOol0ogoocgoo1@002goo 1@00?goo0003Ool00`00Oomoo`07Ool5003oOoo>Ool50008Ool60010Ool000Aoo`03001oogoo00Io o`03001oogoo0?moom9oo`l0049oo`0000Eoo`00Oomoo`0000Qoo`03001oogoo0?moom=oo`d004=o o`0000=oo`0000000P0027oo00<007ooOol0ogooe7oo2P00AGoo000=Ool00`00Oomoo`3oOoooOolT Ool000eoo`03001oogoo0?mooomoobAoo`003Goo00<007ooOol0ogooogoo97oo000=Ool00`00Oomo o`3oOoooOolTOol000eoo`<00?mooomoobAoo`003Goo00<007ooOol0ogooogoo97oo000=Ool00`00 Oomoo`3oOoooOolTOol00001\ \>"], ImageRangeCache->{{{0, 561.063}, {346.375, 0}} -> {-0.642867, -0.616375, \ 0.0470866, 0.0226973}}] }, Open ]], Cell[BoxData[ \(qmax = \(-1\); \ qmin = 8;\)], "Input"] }, Closed]], Cell[CellGroupData[{ Cell["First-order forward difference approximation", "Subsection"], Cell[CellGroupData[{ Cell["Find the update formula", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[k, a, \[Lambda]]; \ order = 1;\), "\[IndentingNewLine]", \(s1forward[q[t], order]\)}], "Input"], Cell[BoxData[ \(\(\(-q[t]\) + q[k + t]\)\/k\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(fmodel[t, q, \[Lambda], a]\)], "Input"], Cell[BoxData[ \(\(-\[Lambda]\)\ q[t] + a\ Sin[t]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(fq\)\(\ \)\(=\)\(\ \)\(First[\ Solve[s1forward[q[t], order]\ \[Equal] \ fmodel[t, q, \[Lambda], a], q[order\ k + t]]\ ]\)\(\ \)\)\)], "Input"], Cell[BoxData[ \({q[k + t] \[Rule] q[t] - k\ \[Lambda]\ q[t] + a\ k\ Sin[t]}\)], "Output"] }, Open ]], Cell["Define the update rule for the sequence of values Qf2", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[Q, n]\), "\[IndentingNewLine]", \(UpdateRule\ = \ Simplify[\ \(\(fq\ \ /. \ {t \[Rule] n\ k}\)\ /. \ q[t_] \[Rule] Q[t/k]\)\ /. \ n \[Rule] n - order\ ]\)}], "Input"], Cell[BoxData[ \({Q[n] \[Rule] \((1 - k\ \[Lambda])\)\ Q[\(-1\) + n] + a\ k\ Sin[k\ \((\(-1\) + n)\)]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[Qupdate];\)\), "\[IndentingNewLine]", \(Qupdate[n_]\ = \ Q[n]\ /. \ UpdateRule\)}], "Input"], Cell[BoxData[ \(\((1 - k\ \[Lambda])\)\ Q[\(-1\) + n] + a\ k\ Sin[k\ \((\(-1\) + n)\)]\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Carry out the computation for specific numerical parameters", \ "Subsubsection"], Cell[BoxData[{ \(\(Clear[Q, k, \[Lambda], a];\)\), "\[IndentingNewLine]", \(k = 1. ; \ \[Lambda] = 0.05; \ a = 1. ;\), "\[IndentingNewLine]", \(Q[n_]\ := \ \(Q[n]\ = \ \ Qupdate[n]\)\)}], "Input"], Cell["Set initial conditions", "Text"], Cell[BoxData[ \(For[n = 0, n < order, \(n++\), Q[n] = qmodel[n\ k, \[Lambda], a]]\)], "Input"], Cell["Carry out computations and compare with exact solution", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Qvalues = Table[{k\ n, Q[n]}, {n, 0, T/k}];\)\), "\[IndentingNewLine]", \(\(Qplot\ = \ ListPlot[Qvalues, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(qplot\ = \ Plot[qmodel[t, \[Lambda], a], {t, 0, T}, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(Show[Qplot, qplot, DisplayFunction \[Rule] $DisplayFunction, PlotRange \[Rule] {{0, T}, {qmin, qmax}}];\)\)}], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -1.14383e-017 0.0397887 0.0686704 0.0686704 [ [.19894 .05617 -3 -9 ] [.19894 .05617 3 0 ] [.39789 .05617 -6 -9 ] [.39789 .05617 6 0 ] [.59683 .05617 -6 -9 ] [.59683 .05617 6 0 ] [.79577 .05617 -6 -9 ] [.79577 .05617 6 0 ] [.99472 .05617 -6 -9 ] [.99472 .05617 6 0 ] [-0.0125 .20601 -6 -4.5 ] [-0.0125 .20601 0 4.5 ] [-0.0125 .34335 -6 -4.5 ] [-0.0125 .34335 0 4.5 ] [-0.0125 .48069 -6 -4.5 ] [-0.0125 .48069 0 4.5 ] [-0.0125 .61803 -6 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .19894 .06867 m .19894 .07492 L s [(5)] .19894 .05617 0 1 Mshowa .39789 .06867 m .39789 .07492 L s [(10)] .39789 .05617 0 1 Mshowa .59683 .06867 m .59683 .07492 L s [(15)] .59683 .05617 0 1 Mshowa .79577 .06867 m .79577 .07492 L s [(20)] .79577 .05617 0 1 Mshowa .99472 .06867 m .99472 .07492 L s [(25)] .99472 .05617 0 1 Mshowa .125 Mabswid .03979 .06867 m .03979 .07242 L s .07958 .06867 m .07958 .07242 L s .11937 .06867 m .11937 .07242 L s .15915 .06867 m .15915 .07242 L s .23873 .06867 m .23873 .07242 L s .27852 .06867 m .27852 .07242 L s .31831 .06867 m .31831 .07242 L s .3581 .06867 m .3581 .07242 L s .43768 .06867 m .43768 .07242 L s .47746 .06867 m .47746 .07242 L s .51725 .06867 m .51725 .07242 L s .55704 .06867 m .55704 .07242 L s .63662 .06867 m .63662 .07242 L s .67641 .06867 m .67641 .07242 L s .7162 .06867 m .7162 .07242 L s .75599 .06867 m .75599 .07242 L s .83556 .06867 m .83556 .07242 L s .87535 .06867 m .87535 .07242 L s .91514 .06867 m .91514 .07242 L s .95493 .06867 m .95493 .07242 L s .25 Mabswid 0 .06867 m 1 .06867 L s 0 .20601 m .00625 .20601 L s [(2)] -0.0125 .20601 1 0 Mshowa 0 .34335 m .00625 .34335 L s [(4)] -0.0125 .34335 1 0 Mshowa 0 .48069 m .00625 .48069 L s [(6)] -0.0125 .48069 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(8)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .10301 m .00375 .10301 L s 0 .13734 m .00375 .13734 L s 0 .17168 m .00375 .17168 L s 0 .24035 m .00375 .24035 L s 0 .27468 m .00375 .27468 L s 0 .30902 m .00375 .30902 L s 0 .37769 m .00375 .37769 L s 0 .41202 m .00375 .41202 L s 0 .44636 m .00375 .44636 L s 0 .51503 m .00375 .51503 L s 0 .54936 m .00375 .54936 L s 0 .5837 m .00375 .5837 L s 0 .03434 m .00375 .03434 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .008 w 0 .13734 Mdot .03979 .13391 Mdot .07958 .18843 Mdot .11937 .24488 Mdot .15915 .24576 Mdot .19894 .18494 Mdot .23873 .11328 Mdot .27852 .09186 Mdot .31831 .13581 Mdot .3581 .2004 Mdot .39789 .22211 Mdot .43768 .17708 Mdot .47746 .10299 Mdot .51725 .06443 Mdot .55704 .09349 Mdot .59683 .16028 Mdot .63662 .20035 Mdot .67641 .174 Mdot .7162 .10271 Mdot .75599 .04944 Mdot .79577 .06069 Mdot .83556 .12378 Mdot .87535 .17848 Mdot .91514 .17238 Mdot .95493 .10909 Mdot .99472 .04488 Mdot .5 Mabswid 0 .13734 m .00123 .13727 L .00236 .13726 L .00365 .13732 L .00488 .13744 L .00601 .1376 L .00706 .13781 L .00942 .13844 L .01181 .13932 L .01444 .14055 L .01975 .14387 L .02493 .14813 L .0305 .15371 L .04057 .16606 L .06325 .20007 L .07443 .21709 L .08481 .2313 L .09483 .24251 L .10043 .24741 L .1055 .25086 L .10841 .25239 L .11109 .2535 L .11261 .25399 L .11403 .25437 L .11551 .25468 L .1163 .2548 L .11714 .2549 L .11787 .25497 L .11855 .25501 L .11929 .25503 L .12008 .25503 L .12141 .25496 L .12216 .25489 L .12287 .2548 L .12418 .25458 L .12538 .25431 L .12809 .25347 L .1307 .25237 L .13307 .25112 L .13839 .24747 L .14331 .24308 L .14796 .2381 L .1697 .20604 L .18893 .17096 L .19986 .15117 L .21009 .13429 L .22099 .11921 L .23096 .1089 L .23357 .10683 L .23642 .10489 L Mistroke .23915 .10335 L .24161 .10222 L .24449 .10124 L .24606 .10085 L .24679 .10071 L .24755 .10059 L .24893 .10043 L .2502 .10035 L .25156 .10035 L .25233 .10039 L .25305 .10044 L .25438 .1006 L .25506 .10071 L .2558 .10086 L .25833 .10152 L .25975 .10202 L .26128 .10264 L .26406 .10402 L .26949 .1076 L .27447 .11185 L .28397 .12224 L .29396 .13576 L .31497 .16867 L .33403 .19704 L .34465 .20951 L .34952 .21405 L .35466 .21793 L .35729 .21952 L .36017 .22095 L .36279 .22195 L .36396 .22231 L .36523 .22263 L .366 .22279 L .36671 .22292 L .36738 .22302 L .36812 .22311 L .3694 .2232 L .37075 .22323 L .37152 .22321 L .37222 .22317 L .37299 .2231 L .37381 .223 L .37518 .22276 L .37668 .22241 L .37933 .22155 L .38221 .22028 L .38745 .21709 L .39254 .21292 L .39724 .20819 L .41923 .17695 L Mistroke .44058 .13881 L .46057 .10517 L .47039 .09191 L .47972 .08227 L .48487 .07838 L .48759 .07676 L .49048 .07537 L .49186 .07483 L .49316 .07441 L .49442 .07406 L .49558 .0738 L .49669 .07361 L .49786 .07346 L .49903 .07338 L .50032 .07335 L .50169 .0734 L .50296 .07352 L .50362 .07361 L .50435 .07373 L .50583 .07405 L .50852 .07486 L .51104 .0759 L .51575 .07853 L .52082 .08234 L .53228 .09426 L .54278 .10846 L .56435 .14277 L .57486 .15948 L .58441 .17325 L .5937 .18452 L .59859 .18939 L .60377 .19362 L .60641 .19539 L .60931 .19701 L .61197 .19821 L .61441 .19904 L .61569 .19937 L .61691 .19963 L .61819 .19983 L .61892 .19991 L .61959 .19997 L .62083 .20002 L .62194 .20001 L .62323 .19993 L .62445 .1998 L .62567 .19959 L .62699 .1993 L .62935 .19859 L .63217 .19743 L .63473 .19608 L Mistroke .64031 .19222 L .64558 .18746 L .65489 .17663 L .665 .16196 L .68541 .12682 L .7079 .08844 L .71817 .07412 L .72383 .0676 L .72922 .06246 L .73421 .05869 L .73684 .05711 L .73965 .05574 L .74225 .05478 L .7446 .05415 L .74574 .05393 L .74697 .05376 L .74813 .05366 L .7492 .05362 L .75053 .05364 L .75178 .05372 L .75291 .05386 L .75413 .05407 L .75546 .05437 L .75688 .05477 L .75942 .0557 L .76395 .05804 L .76896 .06156 L .78014 .07272 L .79065 .08663 L .83126 .1511 L .84215 .16525 L .85229 .17517 L .85741 .17875 L .86027 .1803 L .86287 .18142 L .86426 .1819 L .86573 .18232 L .86649 .18251 L .86734 .18268 L .86811 .18281 L .86882 .18291 L .86954 .18299 L .87022 .18304 L .87094 .18308 L .87173 .18309 L .87241 .18308 L .87314 .18304 L .87444 .18292 L .87575 .18273 L .87719 .18243 L Mistroke .87868 .18203 L .88008 .18156 L .88526 .17912 L .88822 .17723 L .89091 .17521 L .89696 .16964 L .90784 .15644 L .91815 .14093 L .95843 .07244 L .96884 .05831 L .9786 .04819 L .9835 .04446 L .98602 .04291 L .98878 .04152 L .99163 .04042 L .99319 .03996 L .99463 .03964 L .996 .03941 L .99723 .03927 L .99791 .03923 L .99866 .0392 L 1 .03921 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{641, 396.188}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol2001mOol00`00Oomoo`0WOol00`00 Oomoo`1XOol00`00Oomoo`0MOol000eoo`03001oogoo0?moobeoo`03001oogoo011oo`03001oogoo 07Uoo`03001oogoo02Uoo`03001oogoo06Moo`03001oogoo01eoo`003Goo00<007ooOol0ogoo:goo 0P0057oo0P00N7oo00<007ooOol0:Woo00<007ooOol0IWoo00<007ooOol07Woo000=Ool00`00Oomo o`3oOolYOol2000HOol00`00Oomoo`1eOol00`00Oomoo`0[Ool00`00Oomoo`1UOol00`00Oomoo`0N Ool000eoo`03001oogoo0?moobQoo`03001oogoo01Uoo`03001oogoo07=oo`03001oogoo02eoo`03 001oogoo06=oo`03001oogoo01moo`003Goo00<007ooOol0ogoo9goo00<007ooOol06goo00<007oo Ool0LGoo00<007ooOol0;Woo00<007ooOol0HWoo00<007ooOol087oo000=Ool00`00Oomoo`3oOolV Ool00`00Oomoo`0MOol00`00Oomoo`1`Ool00`00Oomoo`0_Ool00`00Oomoo`1QOol00`00Oomoo`0P Ool000eoo`03001oogoo0?moobIoo`03001oogoo01eoo`03001oogoo06moo`03001oogoo035oo`03 001oogoo05moo`03001oogoo025oo`003Goo00<007ooOol0ZWoo0P00N7oo00<007ooOol07goo00<0 07ooOol0KGoo00<007ooOol0Ool2001LOol00`00Oomoo`0cOol00`00Oomoo`1M Ool00`00Oomoo`0ROol000eoo`03001oogoo0:Uoo`@007Eoo`03001oogoo02=oo`03001oogoo00ao o`@005Yoo`03001oogoo03Eoo`03001oogoo05]oo`03001oogoo02=oo`003Goo00<007ooOol0ZWoo 0P00MGoo00<007ooOol09Goo00<007ooOol02goo1000FWoo00<007ooOol0=Goo00<007ooOol0Fgoo 00<007ooOol08goo000=Ool00`00Oomoo`3oOolROol00`00Oomoo`0VOol00`00Oomoo`0;Ool2001J Ool00`00Oomoo`0gOol00`00Oomoo`1IOol00`00Oomoo`0TOol000eoo`03001oogoo0?moob5oo`03 001oogoo02Qoo`03001oogoo06Eoo`03001oogoo03Uoo`03001oogoo05Qoo`03001oogoo02Aoo`00 3Goo00<007ooOol0ogoo87oo00<007ooOol0:Goo00<007ooOol0IGoo00<007ooOol0>Goo00<007oo Ool0Egoo00<007ooOol09Goo000=Ool00`00Oomoo`2FOol9001oOol00`00Oomoo`05Ool2000TOol0 0`00Oomoo`1SOol00`00Oomoo`06Ool2000cOol00`00Oomoo`1FOol00`00Oomoo`0UOol000eoo`@0 095oo`@000Uoo`8007eoo`03001oogoo00Aoo`@002Aoo`03001oogoo069oo`03001oogoo00Eoo`@0 039oo`03001oogoo05Eoo`03001oogoo02Ioo`003Goo00<007ooOol0TGoo00<007ooOol03Goo0P00 NWoo00<007ooOol01Goo10009Goo00<007ooOol0H7oo00<007ooOol01Woo1000Ool0 0`00Oomoo`07Ool4000NOol000eoo`03001oogoo08]oo`03001oogoo00=oo`8001Eoo`03001oogoo 06moo`03001oogoo03Aoo`03001oogoo05Yoo`03001oogoo04=oo`03001oogoo04eoo`03001oogoo 00Uoo`8001moo`003Goo00<007ooOol0RWoo00<007ooOol00goo10005Goo00<007ooOol0KGoo00<0 07ooOol0=Woo00<007ooOol0FGoo00<007ooOol0A7oo00<007ooOol0C7oo00<007ooOol0:Woo000= Ool00`00Oomoo`29Ool00`00Oomoo`04Ool4000FOol00`00Oomoo`1/Ool00`00Oomoo`0fOol00`00 Oomoo`1HOol00`00Oomoo`16Ool00`00Oomoo`1:Ool00`00Oomoo`0[Ool000eoo`03001oogoo08Qo o`03001oogoo00Ioo`8001Qoo`03001oogoo06Yoo`03001oogoo03Qoo`03001oogoo05Ioo`03001o ogoo04Moo`03001oogoo04Yoo`03001oogoo02]oo`003Goo00<007ooOol0Qgoo00<007ooOol08Woo 00<007ooOol0J7oo00<007ooOol0>Woo00<007ooOol0EGoo00<007ooOol0B7oo00<007ooOol0B7oo 00<007ooOol0;7oo000=Ool00`00Oomoo`26Ool00`00Oomoo`0TOol00`00Oomoo`1WOol00`00Oomo o`0jOol00`00Oomoo`1DOol00`00Oomoo`1:Ool00`00Oomoo`17Ool00`00Oomoo`0/Ool000eoo`03 001oogoo08Ioo`03001oogoo02Eoo`03001oogoo06Eoo`03001oogoo03aoo`03001oogoo05=oo`03 001oogoo04Yoo`03001oogoo00Yoo`8003Yoo`03001oogoo02eoo`003Goo00<007ooOol0QGoo00<0 07ooOol09goo00<007ooOol0I7oo00<007ooOol0?Goo00<007ooOol0DGoo00<007ooOol0C7oo00<0 07ooOol027oo1000>7oo00<007ooOol0;Woo000=Ool00`00Oomoo`24Ool00`00Oomoo`0YOol00`00 Oomoo`1ROol00`00Oomoo`0nOol00`00Oomoo`1AOol00`00Oomoo`1=Ool00`00Oomoo`07Ool4000h Ool00`00Oomoo`0^Ool000eoo`03001oogoo08=oo`03001oogoo02Yoo`03001oogoo069oo`03001o ogoo03moo`03001oogoo04moo`03001oogoo04ioo`03001oogoo00Qoo`8003Qoo`03001oogoo02mo o`003Goo00<007ooOol0Pgoo00<007ooOol0:goo00<007ooOol0H7oo00<007ooOol0@7oo00<007oo Ool0CWoo00<007ooOol0D7oo00<007ooOol0@Goo00<007ooOol0;goo000=Ool00`00Oomoo`22Ool0 0`00Oomoo`0]Ool00`00Oomoo`1NOol00`00Oomoo`12Ool00`00Oomoo`1=Ool00`00Oomoo`1AOol0 0`00Oomoo`0oOol00`00Oomoo`0`Ool000eoo`03001oogoo01Eoo`8006Yoo`03001oogoo02moo`03 001oogoo05eoo`03001oogoo04=oo`03001oogoo04]oo`03001oogoo059oo`03001oogoo03ioo`03 001oogoo035oo`003Goo00<007ooOol057oo1000JGoo00<007ooOol0;goo00<007ooOol0G7oo00<0 07ooOol0A7oo00<007ooOol0Bgoo00<007ooOol0Dgoo00<007ooOol0?Goo00<007ooOol0Goo00<007ooOol07oo00<007ooOol0E7oo00<007ooOol0 Bgoo00<007ooOol0@goo00<007ooOol0Fgoo00<007ooOol0=Goo00<007ooOol0=Goo000=Ool00`00 Oomoo`09Ool2001aOol00`00Oomoo`0iOol00`00Oomoo`1COol00`00Oomoo`1=Ool00`00Oomoo`11 Ool00`00Oomoo`1LOol00`00Oomoo`0eOol00`00Oomoo`0eOol000eoo`03001oogoo00]oo`8006mo o`03001oogoo03Yoo`03001oogoo059oo`03001oogoo04eoo`03001oogoo045oo`03001oogoo05eo o`03001oogoo03=oo`03001oogoo03Ioo`003Goo00<007ooOol03Goo00<007ooOol0Jgoo00<007oo Ool0?7oo00<007ooOol0D7oo00<007ooOol0Cgoo00<007ooOol0?goo00<007ooOol0Ggoo00<007oo Ool0Ool00`00Oomoo`1ZOol00`00Oomoo`0lOol0 0`00Oomoo`1@Ool00`00Oomoo`1@Ool00`00Oomoo`0nOol00`00Oomoo`1OOol00`00Oomoo`0`Ool0 0`00Oomoo`0hOol000eoo`03001oogoo00ioo`03001oogoo06Uoo`03001oogoo03ioo`03001oogoo 04ioo`03001oogoo055oo`03001oogoo03eoo`03001oogoo065oo`03001oogoo02moo`03001oogoo 03Qoo`003Goo00<007ooOol03goo00<007ooOol0Igoo00<007ooOol0?goo00<007ooOol0CGoo00<0 07ooOol0Dgoo00<007ooOol0?7oo00<007ooOol0HWoo00<007ooOol0;Goo00<007ooOol0>Goo000= Ool00`00Oomoo`0@Ool00`00Oomoo`1VOol00`00Oomoo`10Ool00`00Oomoo`17oo 00<007ooOol0IWoo00<007ooOol0:Goo00<007ooOol0>goo000=Ool00`00Oomoo`0COol00`00Oomo o`1QOol00`00Oomoo`14Ool00`00Oomoo`18Ool00`00Oomoo`1HOol00`00Oomoo`09Ool2000[Ool0 0`00Oomoo`1WOol00`00Oomoo`0XOol00`00Oomoo`0lOol000eoo`03001oogoo01=oo`03001oogoo 065oo`03001oogoo04Eoo`03001oogoo04Moo`03001oogoo05Uoo`03001oogoo00Moo`@002Yoo`03 001oogoo06Qoo`03001oogoo02Ioo`03001oogoo03eoo`003Goo00<007ooOol057oo00<007ooOol0 Ggoo00<007ooOol0AWoo00<007ooOol0AWoo00<007ooOol0FWoo00<007ooOol01goo1000:Goo00<0 07ooOol0JWoo00<007ooOol09Goo00<007ooOol0?Goo000=Ool00`00Oomoo`0EOol00`00Oomoo`1N Ool00`00Oomoo`17Ool00`00Oomoo`15Ool00`00Oomoo`1KOol00`00Oomoo`07Ool2000YOol00`00 Oomoo`1[Ool00`00Oomoo`0TOol00`00Oomoo`0nOol000eoo`03001oogoo01Ioo`03001oogoo05ao o`03001oogoo04Uoo`03001oogoo04=oo`03001oogoo05eoo`03001oogoo035oo`03001oogoo06ao o`03001oogoo029oo`03001oogoo03moo`003Goo00<007ooOol05goo00<007ooOol0Fgoo00<007oo Ool0BGoo00<007ooOol0@Woo00<007ooOol0GWoo00<007ooOol0<7oo00<007ooOol0KWoo00<007oo Ool087oo00<007ooOol0@7oo000=Ool00`00Oomoo`0HOol00`00Oomoo`1IOol00`00Oomoo`1;Ool0 0`00Oomoo`11Ool00`00Oomoo`1OOol00`00Oomoo`0^Ool00`00Oomoo`1`Ool00`00Oomoo`0OOol0 0`00Oomoo`10Ool000eoo`03001oogoo01Qoo`03001oogoo05Uoo`03001oogoo04aoo`03001oogoo 03moo`03001oogoo065oo`03001oogoo02eoo`03001oogoo075oo`8001ioo`03001oogoo00Moo`80 03Qoo`003Goo100067oo00<007ooOol0Egoo00<007ooOol0CGoo00<007ooOol0?goo00<007ooOol0 HGoo00<007ooOol0;7oo00<007ooOol01goo0P00Jgoo00<007ooOol06Woo00<007ooOol01goo1000 =goo000=Ool00`00Oomoo`0JOol00`00Oomoo`1EOol00`00Oomoo`1?Ool00`00Oomoo`0mOol00`00 Oomoo`1SOol00`00Oomoo`0ZOol00`00Oomoo`07Ool4001[Ool00`00Oomoo`0GOol2000:Ool4000g Ool000eoo`03001oogoo01Yoo`03001oogoo05Eoo`03001oogoo051oo`03001oogoo03aoo`03001o ogoo00Ioo`8005aoo`03001oogoo02Uoo`03001oogoo00Moo`@006aoo`03001oogoo01Eoo`03001o ogoo00]oo`8003Qoo`003Goo00<007ooOol06goo00<007ooOol0Dgoo00<007ooOol0DGoo00<007oo Ool0>goo00<007ooOol01Woo1000Fgoo00<007ooOol0:7oo00<007ooOol02Goo0P00KWoo0P002Woo 0P001goo0P00B7oo000=Ool00`00Oomoo`0LOol00`00Oomoo`1BOol00`00Oomoo`1BOol00`00Oomo o`0iOol00`00Oomoo`07Ool4001LOol00`00Oomoo`0VOol00`00Oomoo`1lOol00`00Oomoo`06Ool4 0005Ool00`00Oomoo`18Ool000eoo`03001oogoo01aoo`03001oogoo055oo`03001oogoo05Aoo`03 001oogoo03Qoo`03001oogoo00Qoo`8005ioo`03001oogoo02Aoo`03001oogoo07ioo`8000Ioo`@0 00=oo`8004]oo`003Goo00<007ooOol07Goo00<007ooOol0D7oo00<007ooOol0E7oo00<007ooOol0 =goo00<007ooOol0JWoo00<007ooOol08Woo00<007ooOol0PGoo0P001Goo0P0000Aoo`00000004eo o`003Goo00<007ooOol07Woo00<007ooOol0CWoo00<007ooOol01Woo0P00CWoo00<007ooOol0=Goo 00<007ooOol0Jgoo00<007ooOol08Woo00<007ooOol0Pgoo2000D7oo000=Ool00`00Oomoo`0NOol0 0`00Oomoo`1=Ool00`00Oomoo`06Ool4001>Ool00`00Oomoo`0cOol00`00Oomoo`1]Ool00`00Oomo o`0POol00`00Oomoo`3LOol000eoo`03001oogoo01moo`03001oogoo00aoo`8003ioo`03001oogoo 00Ioo`@004ioo`03001oogoo03=oo`03001oogoo06ioo`03001oogoo01ioo`03001oogoo0=eoo`00 3Goo00<007ooOol087oo00<007ooOol02Woo1000?7oo00<007ooOol027oo0P00D7oo00<007ooOol0 1oo`003Goo00<007ooOol08Woo 00<007ooOol0Agoo00<007ooOol0GGoo00<007ooOol0;Goo00<007ooOol0MWoo00<007ooOol05Goo 00<007ooOol0hGoo000=Ool00`00Oomoo`0ROol00`00Oomoo`16Ool00`00Oomoo`1OOol00`00Oomo o`0[Ool00`00Oomoo`1hOol00`00Oomoo`0BOol2003TOol000eoo`03001oogoo02=oo`03001oogoo 04Aoo`03001oogoo061oo`03001oogoo02]oo`03001oogoo07Uoo`03001oogoo00moo`800>Ioo`00 3Goo00<007ooOol097oo00<007ooOol0@goo00<007ooOol0HGoo00<007ooOol0:Goo00<007ooOol0 Ngoo0P003Goo0P00j7oo000=Ool00`00Oomoo`0TOol00`00Oomoo`12Ool00`00Oomoo`1SOol00`00 Oomoo`0:Ool2000KOol00`00Oomoo`1nOol50004Ool40002Ool2003VOol000eoo`03001oogoo02Eo o`03001oogoo045oo`03001oogoo06=oo`03001oogoo00Uoo`@001Uoo`03001oogoo08Aoo`D000Ao o`@00>Eoo`0000=oo`0000000P0027oo00<007ooOol09Woo00<007ooOol0?goo00<007ooOol0IGoo 00<007ooOol027oo10006Goo00<007ooOol0SGoo1000iGoo00000goo001oo`0:Ool00`00Oomoo`0V Ool00`00Oomoo`0oOol00`00Oomoo`1VOol00`00Oomoo`08Ool2000IOol00`00Oomoo`2?Ool2003V Ool0009oo`03001oogoo00Qoo`D002Eoo`03001oogoo03eoo`03001oogoo06Qoo`03001oogoo025o o`03001oogoo0?moogUoo`000goo00<007ooOol01goo00<007ooOol0:7oo00<007ooOol0>goo00<0 07ooOol0JGoo00<007ooOol08Goo00<007ooOol0ogooNGoo00001Goo001oogoo000027oo00<007oo Ool0:7oo00<007ooOol0>goo00<007ooOol0JWoo00<007ooOol07goo00<007ooOol0ogooNWoo0002 Ool20009Ool00`00Oomoo`0YOol00`00Oomoo`0iOol00`00Oomoo`1/Ool00`00Oomoo`0MOol00`00 Oomoo`3oOomkOol000eoo`03001oogoo02Yoo`03001oogoo03Moo`03001oogoo06ioo`8001aoo`03 001oogoo0?moogaoo`003Goo00<007ooOol0:goo00<007ooOol0=Woo00<007ooOol0L7oo00<007oo Ool067oo00<007ooOol0ogooOGoo000=Ool00`00Oomoo`0[Ool00`00Oomoo`0eOol00`00Oomoo`1b Ool00`00Oomoo`0EOol2003oOon0Ool000eoo`03001oogoo02aoo`03001oogoo03=oo`03001oogoo 07Aoo`03001oogoo01=oo`03001oogoo0?mooh1oo`003Goo00<007ooOol0;Goo00<007ooOol07oo00<007ooOol07Goo00<007ooOol0ogooogoo67oo 000=Ool00`00Oomoo`0iOol00`00Oomoo`0;Ool2000>Ool00`00Oomoo`06Ool2003oOoooOolAOol0 00eoo`03001oogoo03Yoo`8000Yoo`@000]oo`8000Qoo`@00?mooomooa1oo`003Goo00<007ooOol0 ?7oo00<007ooOol01goo10002Woo00<007ooOol027oo1000ogooogoo47oo000=Ool00`00Oomoo`0m Ool20008Ool2000:Ool00`00Oomoo`0:Ool2003oOoooOolAOol000eoo`03001oogoo03moo`03001o ogoo00moo`800?mooomoob1oo`003Goo00<007ooOol0@7oo0P003Goo0P00ogooogoo8Woo000=Ool0 0`00Oomoo`12Ool30008Ool2003oOoooOolTOol000eoo`03001oogoo04Aoo`T00?mooomoobIoo`00 3Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mo oomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03 001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomc Ool000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomo o`3oOoooOomcOol000eoo`@00?mooomoog9oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00 Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo 000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0 ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo 00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomo og=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001o ogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol0 00eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3o OoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo1000ogooogooLWoo000=Ool00`00Oomo o`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000= Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogoo ogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<0 07ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=o o`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo 0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eo o`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo0003Ool30007Ool00`00Oomo o`3oOoooOomcOol000Aoo`03001oogoo00Ioo`03001oogoo0?mooomoog=oo`0000=oo`0000000`00 1goo1@00ogooogooLGoo00001Goo001oogoo000027oo00<007ooOol0ogooogooLgoo0002Ool00`00 Ool00008Ool00`00Oomoo`3oOoooOomcOol000=oo`8000Qoo`03001oogoo0?mooomoog=oo`003Goo 00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomo og=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001o ogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol0 00eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3o OoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool0 0`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogoo Lgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`@00?mooomoog9oo`003Goo00<007ooOol0ogoo ogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<0 07ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=o o`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo 0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eo o`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOooo OomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00 Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo 000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo1000ogooogoo LWoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007oo Ool0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`00 3Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mo oomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03 001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomc Ool000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomo o`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000= Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo1000ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0 ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo 00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomo og=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001o ogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol0 00eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3o OoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool0 0`00Oomoo`3oOoooOomcOol00004Ool000000009Ool00`00Oomoo`3oOoooOomcOol00005Ool007oo Ool00008Ool00`00Oomoo`3oOoooOomcOol00005Ool007ooOol00008Ool5003oOoooOomaOol00004 Ool000000009Ool00`00Oomoo`3oOoooOomcOol0009oo`03001oogoo00Qoo`03001oogoo0?mooomo og=oo`000Woo0`0027oo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eo o`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOooo OomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00 Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo 000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0 ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo 00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomo og=oo`003Goo1000ogooogooLWoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mo oomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03 001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomc Ool000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomo o`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000= Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogoo ogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<0 07ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=o o`003Goo00<007ooOol0ogooogooLgoo000=Ool4003oOoooOombOol000eoo`03001oogoo0?mooomo og=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001o ogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol0 00eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3o OoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool0 0`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogoo Lgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007oo Ool0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`00 3Goo00<007ooOol0ogooogooLgoo000=Ool4003oOoooOombOol000eoo`03001oogoo0?mooomoog=o o`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo 0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol000eo o`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00Oomoo`3oOooo OomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo000=Ool00`00 Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0ogooogooLgoo 000=Ool00`00Oomoo`3oOoooOomcOol000eoo`03001oogoo0?mooomoog=oo`003Goo00<007ooOol0 ogooogooLgoo000=Ool00`00Oomoo`3oOoooOomcOol0009oo`8000Uoo`03001oogoo0?mooomoog=o o`0000Eoo`00Oomoo`0000Qoo`03001oogoo0?mooomoog=oo`0000Eoo`00Oomoo`0000Qoo`D00?mo oomoog5oo`000Woo0P00ogooogooOgoo00001Goo001oogoo0000ogooogooOWoo0002Ool2003oOooo OomoOol00?mooomooh=oo`00ogooogooPgoo003oOoooOon3Ool00?mooomooh=oo`00\ \>"], ImageRangeCache->{{{0, 640}, {395.188, 0}} -> {-0.556144, -1.06249, \ 0.0403097, 0.0233561}}] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["First-order backward difference approximation", "Subsection"], Cell[CellGroupData[{ Cell["Find the update formula", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[k, a, \[Lambda]]; \ order = 1;\), "\[IndentingNewLine]", \(s1backward[q[t], order]\)}], "Input"], Cell[BoxData[ \(\(q[t] - q[\(-k\) + t]\)\/k\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(fmodel[t, q, \[Lambda], a]\)], "Input"], Cell[BoxData[ \(\(-\[Lambda]\)\ q[t] + a\ Sin[t]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(fq\)\(\ \)\(=\)\(\ \)\(First[\ Solve[s1backward[q[t], order]\ \[Equal] \ fmodel[t, q, \[Lambda], a], q[t]]\ ]\)\(\ \)\)\)], "Input"], Cell[BoxData[ \({q[t] \[Rule] \(q[\(-k\) + t] + a\ k\ Sin[t]\)\/\(1 + k\ \ \[Lambda]\)}\)], "Output"] }, Open ]], Cell["Define the update rule for the sequence of values Qf2", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[Q, n]\), "\[IndentingNewLine]", \(UpdateRule\ = \ Simplify[\ \(fq\ \ /. \ {t \[Rule] n\ k}\)\ /. \ q[t_] \[Rule] Q[t/k]\ \ ]\)}], "Input"], Cell[BoxData[ \({Q[n] \[Rule] \(Q[\(-1\) + n] + a\ k\ Sin[k\ n]\)\/\(1 + k\ \ \[Lambda]\)}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[Qupdate];\)\), "\[IndentingNewLine]", \(Qupdate[n_]\ = \ Q[n]\ /. \ UpdateRule\)}], "Input"], Cell[BoxData[ \(\(Q[\(-1\) + n] + a\ k\ Sin[k\ n]\)\/\(1 + k\ \[Lambda]\)\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Carry out the computation for specific numerical parameters", \ "Subsubsection"], Cell[BoxData[{ \(\(Clear[Q, k, \[Lambda], a];\)\), "\[IndentingNewLine]", \(k = 1. ; \ \[Lambda] = 0.05; \ a = 1. ;\), "\[IndentingNewLine]", \(Q[n_]\ := \ \(Q[n]\ = \ \ Qupdate[n]\)\)}], "Input"], Cell["Set initial conditions", "Text"], Cell[BoxData[ \(For[n = 0, n < order, \(n++\), Q[n] = qmodel[n\ k, \[Lambda], a]]\)], "Input"], Cell["Carry out computations and compare with exact solution", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Qvalues = Table[{k\ n, Q[n]}, {n, 0, T/k}];\)\), "\[IndentingNewLine]", \(\(Qplot\ = \ ListPlot[Qvalues, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(qplot\ = \ Plot[qmodel[t, \[Lambda], a], {t, 0, T}, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(Show[Qplot, qplot, DisplayFunction \[Rule] $DisplayFunction, PlotRange \[Rule] {{0, T}, {qmin, qmax}}];\)\)}], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -1.14383e-017 0.0397887 0.0686704 0.0686704 [ [.19894 .05617 -3 -9 ] [.19894 .05617 3 0 ] [.39789 .05617 -6 -9 ] [.39789 .05617 6 0 ] [.59683 .05617 -6 -9 ] [.59683 .05617 6 0 ] [.79577 .05617 -6 -9 ] [.79577 .05617 6 0 ] [.99472 .05617 -6 -9 ] [.99472 .05617 6 0 ] [-0.0125 .20601 -6 -4.5 ] [-0.0125 .20601 0 4.5 ] [-0.0125 .34335 -6 -4.5 ] [-0.0125 .34335 0 4.5 ] [-0.0125 .48069 -6 -4.5 ] [-0.0125 .48069 0 4.5 ] [-0.0125 .61803 -6 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .19894 .06867 m .19894 .07492 L s [(5)] .19894 .05617 0 1 Mshowa .39789 .06867 m .39789 .07492 L s [(10)] .39789 .05617 0 1 Mshowa .59683 .06867 m .59683 .07492 L s [(15)] .59683 .05617 0 1 Mshowa .79577 .06867 m .79577 .07492 L s [(20)] .79577 .05617 0 1 Mshowa .99472 .06867 m .99472 .07492 L s [(25)] .99472 .05617 0 1 Mshowa .125 Mabswid .03979 .06867 m .03979 .07242 L s .07958 .06867 m .07958 .07242 L s .11937 .06867 m .11937 .07242 L s .15915 .06867 m .15915 .07242 L s .23873 .06867 m .23873 .07242 L s .27852 .06867 m .27852 .07242 L s .31831 .06867 m .31831 .07242 L s .3581 .06867 m .3581 .07242 L s .43768 .06867 m .43768 .07242 L s .47746 .06867 m .47746 .07242 L s .51725 .06867 m .51725 .07242 L s .55704 .06867 m .55704 .07242 L s .63662 .06867 m .63662 .07242 L s .67641 .06867 m .67641 .07242 L s .7162 .06867 m .7162 .07242 L s .75599 .06867 m .75599 .07242 L s .83556 .06867 m .83556 .07242 L s .87535 .06867 m .87535 .07242 L s .91514 .06867 m .91514 .07242 L s .95493 .06867 m .95493 .07242 L s .25 Mabswid 0 .06867 m 1 .06867 L s 0 .20601 m .00625 .20601 L s [(2)] -0.0125 .20601 1 0 Mshowa 0 .34335 m .00625 .34335 L s [(4)] -0.0125 .34335 1 0 Mshowa 0 .48069 m .00625 .48069 L s [(6)] -0.0125 .48069 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(8)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .10301 m .00375 .10301 L s 0 .13734 m .00375 .13734 L s 0 .17168 m .00375 .17168 L s 0 .24035 m .00375 .24035 L s 0 .27468 m .00375 .27468 L s 0 .30902 m .00375 .30902 L s 0 .37769 m .00375 .37769 L s 0 .41202 m .00375 .41202 L s 0 .44636 m .00375 .44636 L s 0 .51503 m .00375 .51503 L s 0 .54936 m .00375 .54936 L s 0 .5837 m .00375 .5837 L s 0 .03434 m .00375 .03434 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .008 w 0 .13734 Mdot .03979 .1891 Mdot .07958 .24284 Mdot .11937 .24377 Mdot .15915 .18594 Mdot .19894 .11764 Mdot .23873 .09704 Mdot .27852 .13865 Mdot .31831 .20002 Mdot .3581 .22072 Mdot .39789 .1779 Mdot .43768 .1073 Mdot .47746 .07037 Mdot .51725 .09777 Mdot .55704 .16117 Mdot .59683 .19929 Mdot .63662 .17424 Mdot .67641 .10634 Mdot .7162 .05543 Mdot .75599 .06586 Mdot .79577 .1257 Mdot .83556 .17771 Mdot .87535 .17194 Mdot .91514 .11167 Mdot .95493 .0504 Mdot .99472 .04262 Mdot .5 Mabswid 0 .13734 m .00123 .13727 L .00236 .13726 L .00365 .13732 L .00488 .13744 L .00601 .1376 L .00706 .13781 L .00942 .13844 L .01181 .13932 L .01444 .14055 L .01975 .14387 L .02493 .14813 L .0305 .15371 L .04057 .16606 L .06325 .20007 L .07443 .21709 L .08481 .2313 L .09483 .24251 L .10043 .24741 L .1055 .25086 L .10841 .25239 L .11109 .2535 L .11261 .25399 L .11403 .25437 L .11551 .25468 L .1163 .2548 L .11714 .2549 L .11787 .25497 L .11855 .25501 L .11929 .25503 L .12008 .25503 L .12141 .25496 L .12216 .25489 L .12287 .2548 L .12418 .25458 L .12538 .25431 L .12809 .25347 L .1307 .25237 L .13307 .25112 L .13839 .24747 L .14331 .24308 L .14796 .2381 L .1697 .20604 L .18893 .17096 L .19986 .15117 L .21009 .13429 L .22099 .11921 L .23096 .1089 L .23357 .10683 L .23642 .10489 L Mistroke .23915 .10335 L .24161 .10222 L .24449 .10124 L .24606 .10085 L .24679 .10071 L .24755 .10059 L .24893 .10043 L .2502 .10035 L .25156 .10035 L .25233 .10039 L .25305 .10044 L .25438 .1006 L .25506 .10071 L .2558 .10086 L .25833 .10152 L .25975 .10202 L .26128 .10264 L .26406 .10402 L .26949 .1076 L .27447 .11185 L .28397 .12224 L .29396 .13576 L .31497 .16867 L .33403 .19704 L .34465 .20951 L .34952 .21405 L .35466 .21793 L .35729 .21952 L .36017 .22095 L .36279 .22195 L .36396 .22231 L .36523 .22263 L .366 .22279 L .36671 .22292 L .36738 .22302 L .36812 .22311 L .3694 .2232 L .37075 .22323 L .37152 .22321 L .37222 .22317 L .37299 .2231 L .37381 .223 L .37518 .22276 L .37668 .22241 L .37933 .22155 L .38221 .22028 L .38745 .21709 L .39254 .21292 L .39724 .20819 L .41923 .17695 L Mistroke .44058 .13881 L .46057 .10517 L .47039 .09191 L .47972 .08227 L .48487 .07838 L .48759 .07676 L .49048 .07537 L .49186 .07483 L .49316 .07441 L .49442 .07406 L .49558 .0738 L .49669 .07361 L .49786 .07346 L .49903 .07338 L .50032 .07335 L .50169 .0734 L .50296 .07352 L .50362 .07361 L .50435 .07373 L .50583 .07405 L .50852 .07486 L .51104 .0759 L .51575 .07853 L .52082 .08234 L .53228 .09426 L .54278 .10846 L .56435 .14277 L .57486 .15948 L .58441 .17325 L .5937 .18452 L .59859 .18939 L .60377 .19362 L .60641 .19539 L .60931 .19701 L .61197 .19821 L .61441 .19904 L .61569 .19937 L .61691 .19963 L .61819 .19983 L .61892 .19991 L .61959 .19997 L .62083 .20002 L .62194 .20001 L .62323 .19993 L .62445 .1998 L .62567 .19959 L .62699 .1993 L .62935 .19859 L .63217 .19743 L .63473 .19608 L Mistroke .64031 .19222 L .64558 .18746 L .65489 .17663 L .665 .16196 L .68541 .12682 L .7079 .08844 L .71817 .07412 L .72383 .0676 L .72922 .06246 L .73421 .05869 L .73684 .05711 L .73965 .05574 L .74225 .05478 L .7446 .05415 L .74574 .05393 L .74697 .05376 L .74813 .05366 L .7492 .05362 L .75053 .05364 L .75178 .05372 L .75291 .05386 L .75413 .05407 L .75546 .05437 L .75688 .05477 L .75942 .0557 L .76395 .05804 L .76896 .06156 L .78014 .07272 L .79065 .08663 L .83126 .1511 L .84215 .16525 L .85229 .17517 L .85741 .17875 L .86027 .1803 L .86287 .18142 L .86426 .1819 L .86573 .18232 L .86649 .18251 L .86734 .18268 L .86811 .18281 L .86882 .18291 L .86954 .18299 L .87022 .18304 L .87094 .18308 L .87173 .18309 L .87241 .18308 L .87314 .18304 L .87444 .18292 L .87575 .18273 L .87719 .18243 L Mistroke .87868 .18203 L .88008 .18156 L .88526 .17912 L .88822 .17723 L .89091 .17521 L .89696 .16964 L .90784 .15644 L .91815 .14093 L .95843 .07244 L .96884 .05831 L .9786 .04819 L .9835 .04446 L .98602 .04291 L .98878 .04152 L .99163 .04042 L .99319 .03996 L .99463 .03964 L .996 .03941 L .99723 .03927 L .99791 .03923 L .99866 .0392 L 1 .03921 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{640.688, 396}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol2001lOol00`00Oomoo`0WOol00`00Oomoo`1YOol00`00Oomoo`0LOol000eoo`03 001oogoo0?moobaoo`80019oo`03001oogoo07Uoo`03001oogoo02Qoo`03001oogoo06Qoo`03001o ogoo01aoo`003Goo00<007ooOol0ogoo:Woo0P005Goo0P00N7oo00<007ooOol0:Woo00<007ooOol0 IWoo00<007ooOol07Goo000=Ool00`00Oomoo`3oOolYOol00`00Oomoo`0GOol00`00Oomoo`1dOol0 0`00Oomoo`0/Ool00`00Oomoo`1TOol00`00Oomoo`0NOol000eoo`03001oogoo0?moobQoo`03001o ogoo01Uoo`03001oogoo079oo`03001oogoo02eoo`03001oogoo06Aoo`03001oogoo01ioo`003Goo 00<007ooOol0ogoo9goo00<007ooOol06goo00<007ooOol0LGoo00<007ooOol0;Woo00<007ooOol0 HWoo00<007ooOol07goo000=Ool00`00Oomoo`3oOolVOol00`00Oomoo`0MOol00`00Oomoo`1_Ool0 0`00Oomoo`0`Ool00`00Oomoo`1QOol00`00Oomoo`0OOol000eoo`03001oogoo0?moobEoo`03001o ogoo01ioo`03001oogoo06ioo`03001oogoo039oo`03001oogoo05moo`03001oogoo021oo`003Goo 00<007ooOol0ogoo97oo00<007ooOol087oo00<007ooOol0K7oo00<007ooOol0Woo00<007ooOol0Egoo00<007ooOol097oo000=Ool00`00Oomoo`2AOol2 0003Ool9001oOol00`00Oomoo`0NOol2000;Ool00`00Oomoo`1ROol00`00Oomoo`0lOol00`00Oomo o`1EOol00`00Oomoo`0UOol000eoo`@009=oo`8000Uoo`8007aoo`03001oogoo02eoo`03001oogoo 065oo`03001oogoo03aoo`03001oogoo05Eoo`03001oogoo02Eoo`003Goo00<007ooOol0TWoo0P00 3Goo0P00NWoo00<007ooOol0;Woo00<007ooOol0DGoo0P0037oo00<007ooOol0?Woo00<007ooOol0 Dgoo00<007ooOol09Woo000=Ool00`00Oomoo`2AOol00`00Oomoo`0?Ool2001XOol2000=Ool00`00 Oomoo`0_Ool00`00Oomoo`1@Ool4000;Ool00`00Oomoo`0nOol00`00Oomoo`1BOol00`00Oomoo`0W Ool000eoo`03001oogoo08moo`8001Aoo`03001oogoo06Aoo`@000]oo`03001oogoo035oo`03001o ogoo04moo`@000Yoo`03001oogoo041oo`03001oogoo055oo`03001oogoo02Moo`003Goo00<007oo Ool0SGoo0P005goo00<007ooOol0Hgoo10002goo00<007ooOol0Woo00<007ooOol0E7oo00<007oo Ool0B7oo00<007ooOol0B7oo00<007ooOol0;7oo000=Ool00`00Oomoo`1hOol2000Goo 00<007ooOol0goo00<007ooOol0DGoo00<007ooOol0Cgoo00<007ooOol0?goo00<007ooOol0GWoo00<007oo Ool0Ool00`00Oomoo`1ZOol00`00Oomoo`0lOol0 0`00Oomoo`1@Ool00`00Oomoo`1@Ool00`00Oomoo`0mOol00`00Oomoo`1OOol00`00Oomoo`0aOol0 0`00Oomoo`0gOol000eoo`03001oogoo00ioo`03001oogoo06Uoo`03001oogoo03eoo`03001oogoo 04moo`03001oogoo055oo`03001oogoo03eoo`03001oogoo061oo`03001oogoo02moo`03001oogoo 03Qoo`003Goo00<007ooOol03goo00<007ooOol0Igoo00<007ooOol0?goo00<007ooOol0CGoo00<0 07ooOol0Dgoo00<007ooOol0>goo00<007ooOol0HWoo00<007ooOol0;Woo00<007ooOol0>7oo000= Ool00`00Oomoo`0@Ool00`00Oomoo`1VOol00`00Oomoo`10Ool00`00Oomoo`17oo 00<007ooOol0IGoo00<007ooOol0:Woo00<007ooOol0>Woo000=Ool00`00Oomoo`0COol00`00Oomo o`1QOol00`00Oomoo`14Ool00`00Oomoo`18Ool00`00Oomoo`1GOol00`00Oomoo`0gOol00`00Oomo o`1WOol00`00Oomoo`0XOol00`00Oomoo`0kOol000eoo`03001oogoo01=oo`03001oogoo065oo`03 001oogoo04Aoo`03001oogoo04Qoo`03001oogoo04]oo`8000]oo`03001oogoo03Ioo`03001oogoo 06Qoo`03001oogoo02Ioo`03001oogoo03aoo`003Goo00<007ooOol057oo00<007ooOol0Ggoo00<0 07ooOol0AWoo00<007ooOol0AWoo00<007ooOol0Bgoo10002goo00<007ooOol0=7oo00<007ooOol0 JGoo00<007ooOol09Goo00<007ooOol0?Goo000=Ool00`00Oomoo`0EOol00`00Oomoo`1NOol00`00 Oomoo`17Ool00`00Oomoo`15Ool00`00Oomoo`1;Ool4000;Ool00`00Oomoo`0cOol00`00Oomoo`1[ Ool00`00Oomoo`0SOol00`00Oomoo`0nOol000eoo`03001oogoo01Ioo`03001oogoo05aoo`03001o ogoo04Uoo`03001oogoo04=oo`03001oogoo04eoo`8000eoo`03001oogoo035oo`03001oogoo06eo o`03001oogoo029oo`03001oogoo03ioo`003Goo00<007ooOol05goo00<007ooOol0Fgoo00<007oo Ool0BGoo00<007ooOol0@Woo00<007ooOol0GGoo00<007ooOol0Ool00`00Oomoo`1EOol00`00 Oomoo`0fOol00`00Oomoo`1[Ool00`00Oomoo`0QOol00`00Oomoo`23Ool:001>Ool000eoo`03001o ogoo01ioo`03001oogoo03moo`8000eoo`03001oogoo05Ioo`03001oogoo03Aoo`03001oogoo06eo o`03001oogoo01moo`03001oogoo0=aoo`003Goo00<007ooOol07goo00<007ooOol0?Goo10002goo 00<007ooOol0F7oo00<007ooOol0Aoo`003Goo00<007ooOol08goo 00<007ooOol0AGoo00<007ooOol0GWoo00<007ooOol0;7oo00<007ooOol0NGoo00<007ooOol03goo 00<007ooOol0i7oo000=Ool00`00Oomoo`0TOol00`00Oomoo`13Ool00`00Oomoo`1POol00`00Oomo o`0ZOol00`00Oomoo`1cOol20006Ool2000=Ool2003WOol000eoo`03001oogoo02Aoo`03001oogoo 04=oo`03001oogoo05Ioo`8000Uoo`03001oogoo02Qoo`03001oogoo07=oo`@000Moo`@000Ioo`<0 0>Uoo`003Goo00<007ooOol09Goo00<007ooOol0@Goo00<007ooOol0EWoo10002Goo00<007ooOol0 9Woo00<007ooOol0M7oo10002goo1P00k7oo000=Ool00`00Oomoo`0VOol00`00Oomoo`10Ool00`00 Oomoo`1FOol4000:Ool00`00Oomoo`0UOol00`00Oomoo`1eOol2003nOol00003Ool00000008000Qo o`03001oogoo02Ioo`03001oogoo03moo`03001oogoo05Qoo`8000]oo`03001oogoo02Aoo`03001o ogoo0?moogMoo`0000=oo`00Ool02Woo00<007ooOol09goo00<007ooOol0?Woo00<007ooOol0IWoo 00<007ooOol08Woo00<007ooOol0ogooN7oo0002Ool00`00Oomoo`08Ool5000VOol00`00Oomoo`0l Ool00`00Oomoo`1XOol00`00Oomoo`0QOol00`00Oomoo`3oOomhOol000=oo`03001oogoo00Moo`03 001oogoo02Qoo`03001oogoo03]oo`03001oogoo06Yoo`03001oogoo01moo`03001oogoo0?moogUo o`0000Eoo`00Oomoo`0000Qoo`03001oogoo02Uoo`03001oogoo03Yoo`03001oogoo06]oo`03001o ogoo01eoo`03001oogoo0?moogYoo`000Woo0P002Goo00<007ooOol0:Woo00<007ooOol0>7oo00<0 07ooOol0KGoo00<007ooOol06goo00<007ooOol0ogooNgoo000=Ool00`00Oomoo`0[Ool00`00Oomo o`0fOol00`00Oomoo`1_Ool00`00Oomoo`0IOol00`00Oomoo`3oOomlOol000eoo`03001oogoo02]o o`03001oogoo03Ioo`03001oogoo071oo`03001oogoo01Ioo`800?moogmoo`003Goo00<007ooOol0 ;7oo00<007ooOol0=7oo00<007ooOol0LWoo0P005Goo00<007ooOol0ogooOgoo000=Ool00`00Oomo o`0]Ool00`00Oomoo`0bOol00`00Oomoo`1eOol00`00Oomoo`0@Ool2003oOon2Ool000eoo`03001o ogoo02eoo`03001oogoo039oo`03001oogoo07Ioo`<000eoo`800?moohAoo`003Goo00<007ooOol0 ;Woo00<007ooOol0<7oo00<007ooOol0Mgoo1@002Woo00<007ooOol0ogooQ7oo000=Ool00`00Oomo o`0_Ool00`00Oomoo`0^Ool00`00Oomoo`1hOol400000goo00000008003oOon7Ool000eoo`03001o ogoo02moo`03001oogoo02ioo`03001oogoo07Uoo`800?mooi=oo`003Goo00<007ooOol0<7oo00<0 07ooOol0;7oo00<007ooOol0ogooogoo47oo000=Ool00`00Oomoo`0aOol00`00Oomoo`0[Ool00`00 Oomoo`3oOoooOol@Ool000eoo`03001oogoo035oo`03001oogoo02Yoo`03001oogoo0?mooomooa5o o`003Goo00<007ooOol0"], ImageRangeCache->{{{0, 639.688}, {395, 0}} -> {-0.556265, -1.06257, \ 0.0403297, 0.0233677}}] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Second-order forward difference approximation", "Subsection"], Cell[CellGroupData[{ Cell["Find the update formula", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[k, a, \[Lambda]]; \ order = 2;\), "\[IndentingNewLine]", \(s1forward[q[t], order]\)}], "Input"], Cell[BoxData[ \(\(-\(\(3\ q[t] - 4\ q[k + t] + q[2\ k + t]\)\/\(2\ k\)\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(fmodel[t, q, \[Lambda], a]\)], "Input"], Cell[BoxData[ \(\(-\[Lambda]\)\ q[t] + a\ Sin[t]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(fq\)\(\ \)\(=\)\(\ \)\(First[\ Solve[s1forward[q[t], order]\ \[Equal] \ fmodel[t, q, \[Lambda], a], q[order\ k + t]]\ ]\)\(\ \)\)\)], "Input"], Cell[BoxData[ \({q[2\ k + t] \[Rule] \(-3\)\ q[t] + 2\ k\ \[Lambda]\ q[t] + 4\ q[k + t] - 2\ a\ k\ Sin[t]}\)], "Output"] }, Open ]], Cell["Define the update rule for the sequence of values Qf2", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[Q, n];\)\), "\[IndentingNewLine]", \(UpdateRule\ = \ Simplify[\ \(\(fq\ \ /. \ {t \[Rule] n\ k}\)\ /. \ q[t_] \[Rule] Q[t/k]\)\ /. \ n \[Rule] n - order\ ]\)}], "Input"], Cell[BoxData[ \({Q[n] \[Rule] \((\(-3\) + 2\ k\ \[Lambda])\)\ Q[\(-2\) + n] + 4\ Q[\(-1\) + n] - 2\ a\ k\ Sin[k\ \((\(-2\) + n)\)]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[Qupdate];\)\), "\[IndentingNewLine]", \(Qupdate[n_]\ = \ Q[n]\ /. \ UpdateRule\)}], "Input"], Cell[BoxData[ \(\((\(-3\) + 2\ k\ \[Lambda])\)\ Q[\(-2\) + n] + 4\ Q[\(-1\) + n] - 2\ a\ k\ Sin[k\ \((\(-2\) + n)\)]\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Carry out the computation for specific numerical parameters", \ "Subsubsection"], Cell[BoxData[{ \(\(Clear[Q, k, \[Lambda], a];\)\), "\[IndentingNewLine]", \(k = 0.1; \ \[Lambda] = 0.05; \ a = 1. ;\), "\[IndentingNewLine]", \(Q[n_]\ := \ \(Q[n]\ = \ \ Qupdate[n]\)\)}], "Input"], Cell["Set initial conditions", "Text"], Cell[BoxData[ \(For[n = 0, n < order, \(n++\), Q[n] = qmodel[n\ k, \[Lambda], a]]\)], "Input"], Cell["Carry out computations and compare with exact solution", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Qvalues = Table[{k\ n, Q[n]}, {n, 0, T/k}];\)\), "\[IndentingNewLine]", \(\(Qplot\ = \ ListPlot[Qvalues, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(qplot\ = \ Plot[qmodel[t, \[Lambda], a], {t, 0, T}, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(Show[Qplot, qplot, DisplayFunction \[Rule] $DisplayFunction, PlotRange \[Rule] {{0, T}, {qmin, qmax}}];\)\)}], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -1.14383e-017 0.0397887 0.0686704 0.0686704 [ [.19894 .05617 -3 -9 ] [.19894 .05617 3 0 ] [.39789 .05617 -6 -9 ] [.39789 .05617 6 0 ] [.59683 .05617 -6 -9 ] [.59683 .05617 6 0 ] [.79577 .05617 -6 -9 ] [.79577 .05617 6 0 ] [.99472 .05617 -6 -9 ] [.99472 .05617 6 0 ] [-0.0125 .20601 -6 -4.5 ] [-0.0125 .20601 0 4.5 ] [-0.0125 .34335 -6 -4.5 ] [-0.0125 .34335 0 4.5 ] [-0.0125 .48069 -6 -4.5 ] [-0.0125 .48069 0 4.5 ] [-0.0125 .61803 -6 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .19894 .06867 m .19894 .07492 L s [(5)] .19894 .05617 0 1 Mshowa .39789 .06867 m .39789 .07492 L s [(10)] .39789 .05617 0 1 Mshowa .59683 .06867 m .59683 .07492 L s [(15)] .59683 .05617 0 1 Mshowa .79577 .06867 m .79577 .07492 L s [(20)] .79577 .05617 0 1 Mshowa .99472 .06867 m .99472 .07492 L s [(25)] .99472 .05617 0 1 Mshowa .125 Mabswid .03979 .06867 m .03979 .07242 L s .07958 .06867 m .07958 .07242 L s .11937 .06867 m .11937 .07242 L s .15915 .06867 m .15915 .07242 L s .23873 .06867 m .23873 .07242 L s .27852 .06867 m .27852 .07242 L s .31831 .06867 m .31831 .07242 L s .3581 .06867 m .3581 .07242 L s .43768 .06867 m .43768 .07242 L s .47746 .06867 m .47746 .07242 L s .51725 .06867 m .51725 .07242 L s .55704 .06867 m .55704 .07242 L s .63662 .06867 m .63662 .07242 L s .67641 .06867 m .67641 .07242 L s .7162 .06867 m .7162 .07242 L s .75599 .06867 m .75599 .07242 L s .83556 .06867 m .83556 .07242 L s .87535 .06867 m .87535 .07242 L s .91514 .06867 m .91514 .07242 L s .95493 .06867 m .95493 .07242 L s .25 Mabswid 0 .06867 m 1 .06867 L s 0 .20601 m .00625 .20601 L s [(2)] -0.0125 .20601 1 0 Mshowa 0 .34335 m .00625 .34335 L s [(4)] -0.0125 .34335 1 0 Mshowa 0 .48069 m .00625 .48069 L s [(6)] -0.0125 .48069 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(8)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .10301 m .00375 .10301 L s 0 .13734 m .00375 .13734 L s 0 .17168 m .00375 .17168 L s 0 .24035 m .00375 .24035 L s 0 .27468 m .00375 .27468 L s 0 .30902 m .00375 .30902 L s 0 .37769 m .00375 .37769 L s 0 .41202 m .00375 .41202 L s 0 .44636 m .00375 .44636 L s 0 .51503 m .00375 .51503 L s 0 .54936 m .00375 .54936 L s 0 .5837 m .00375 .5837 L s 0 .03434 m .00375 .03434 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .008 w 0 .13734 Mdot .00398 .13734 Mdot .00796 .13803 Mdot .01194 .1394 Mdot .01592 .1415 Mdot .01989 .14442 Mdot .02387 .14858 Mdot .02785 .15521 Mdot .03183 .16817 Mdot .03581 .19906 Mdot .03979 .28286 Mdot .04377 .52483 Mdot .5 Mabswid 0 .13734 m .00123 .13727 L .00236 .13726 L .00365 .13732 L .00488 .13744 L .00601 .1376 L .00706 .13781 L .00942 .13844 L .01181 .13932 L .01444 .14055 L .01975 .14387 L .02493 .14813 L .0305 .15371 L .04057 .16606 L .06325 .20007 L .07443 .21709 L .08481 .2313 L .09483 .24251 L .10043 .24741 L .1055 .25086 L .10841 .25239 L .11109 .2535 L .11261 .25399 L .11403 .25437 L .11551 .25468 L .1163 .2548 L .11714 .2549 L .11787 .25497 L .11855 .25501 L .11929 .25503 L .12008 .25503 L .12141 .25496 L .12216 .25489 L .12287 .2548 L .12418 .25458 L .12538 .25431 L .12809 .25347 L .1307 .25237 L .13307 .25112 L .13839 .24747 L .14331 .24308 L .14796 .2381 L .1697 .20604 L .18893 .17096 L .19986 .15117 L .21009 .13429 L .22099 .11921 L .23096 .1089 L .23357 .10683 L .23642 .10489 L Mistroke .23915 .10335 L .24161 .10222 L .24449 .10124 L .24606 .10085 L .24679 .10071 L .24755 .10059 L .24893 .10043 L .2502 .10035 L .25156 .10035 L .25233 .10039 L .25305 .10044 L .25438 .1006 L .25506 .10071 L .2558 .10086 L .25833 .10152 L .25975 .10202 L .26128 .10264 L .26406 .10402 L .26949 .1076 L .27447 .11185 L .28397 .12224 L .29396 .13576 L .31497 .16867 L .33403 .19704 L .34465 .20951 L .34952 .21405 L .35466 .21793 L .35729 .21952 L .36017 .22095 L .36279 .22195 L .36396 .22231 L .36523 .22263 L .366 .22279 L .36671 .22292 L .36738 .22302 L .36812 .22311 L .3694 .2232 L .37075 .22323 L .37152 .22321 L .37222 .22317 L .37299 .2231 L .37381 .223 L .37518 .22276 L .37668 .22241 L .37933 .22155 L .38221 .22028 L .38745 .21709 L .39254 .21292 L .39724 .20819 L .41923 .17695 L Mistroke .44058 .13881 L .46057 .10517 L .47039 .09191 L .47972 .08227 L .48487 .07838 L .48759 .07676 L .49048 .07537 L .49186 .07483 L .49316 .07441 L .49442 .07406 L .49558 .0738 L .49669 .07361 L .49786 .07346 L .49903 .07338 L .50032 .07335 L .50169 .0734 L .50296 .07352 L .50362 .07361 L .50435 .07373 L .50583 .07405 L .50852 .07486 L .51104 .0759 L .51575 .07853 L .52082 .08234 L .53228 .09426 L .54278 .10846 L .56435 .14277 L .57486 .15948 L .58441 .17325 L .5937 .18452 L .59859 .18939 L .60377 .19362 L .60641 .19539 L .60931 .19701 L .61197 .19821 L .61441 .19904 L .61569 .19937 L .61691 .19963 L .61819 .19983 L .61892 .19991 L .61959 .19997 L .62083 .20002 L .62194 .20001 L .62323 .19993 L .62445 .1998 L .62567 .19959 L .62699 .1993 L .62935 .19859 L .63217 .19743 L .63473 .19608 L Mistroke .64031 .19222 L .64558 .18746 L .65489 .17663 L .665 .16196 L .68541 .12682 L .7079 .08844 L .71817 .07412 L .72383 .0676 L .72922 .06246 L .73421 .05869 L .73684 .05711 L .73965 .05574 L .74225 .05478 L .7446 .05415 L .74574 .05393 L .74697 .05376 L .74813 .05366 L .7492 .05362 L .75053 .05364 L .75178 .05372 L .75291 .05386 L .75413 .05407 L .75546 .05437 L .75688 .05477 L .75942 .0557 L .76395 .05804 L .76896 .06156 L .78014 .07272 L .79065 .08663 L .83126 .1511 L .84215 .16525 L .85229 .17517 L .85741 .17875 L .86027 .1803 L .86287 .18142 L .86426 .1819 L .86573 .18232 L .86649 .18251 L .86734 .18268 L .86811 .18281 L .86882 .18291 L .86954 .18299 L .87022 .18304 L .87094 .18308 L .87173 .18309 L .87241 .18308 L .87314 .18304 L .87444 .18292 L .87575 .18273 L .87719 .18243 L Mistroke .87868 .18203 L .88008 .18156 L .88526 .17912 L .88822 .17723 L .89091 .17521 L .89696 .16964 L .90784 .15644 L .91815 .14093 L .95843 .07244 L .96884 .05831 L .9786 .04819 L .9835 .04446 L .98602 .04291 L .98878 .04152 L .99163 .04042 L .99319 .03996 L .99463 .03964 L .996 .03941 L .99723 .03927 L .99791 .03923 L .99866 .0392 L 1 .03921 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{640.438, 395.813}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgWoo00<007ooOol0EWoo00<0 07ooOol09Goo000=Ool00`00Oomoo`2FOol80020Ool00`00Oomoo`0[Ool00`00Oomoo`1ROol00`00 Oomoo`0lOol00`00Oomoo`1EOol00`00Oomoo`0UOol000eoo`03001oogoo09Aoo`8000Qoo`<007ao o`03001oogoo02aoo`03001oogoo069oo`03001oogoo03aoo`03001oogoo05Aoo`03001oogoo02Io o`003Goo1000TGoo0P003Goo0P00NWoo00<007ooOol0;Goo00<007ooOol0H7oo00<007ooOol0?Woo 00<007ooOol0DWoo00<007ooOol09goo000=Ool00`00Oomoo`2AOol00`00Oomoo`0?Ool2001gOol0 0`00Oomoo`0_Ool00`00Oomoo`1OOol00`00Oomoo`0nOol00`00Oomoo`1BOol00`00Oomoo`0WOol0 00eoo`03001oogoo08moo`8001Aoo`03001oogoo07=oo`03001oogoo031oo`03001oogoo05ioo`03 001oogoo041oo`03001oogoo051oo`03001oogoo02Qoo`003Goo00<007ooOol0SGoo0P005goo00<0 07ooOol0LWoo00<007ooOol0Ool00`00Oomoo`0YOol000eoo`03001o ogoo08]oo`03001oogoo01Yoo`03001oogoo06moo`03001oogoo03=oo`03001oogoo05Yoo`03001o ogoo04Aoo`03001oogoo04aoo`03001oogoo02Yoo`003Goo00<007ooOol0RWoo00<007ooOol077oo 00<007ooOol0KGoo00<007ooOol0=Goo00<007ooOol0FGoo00<007ooOol0AGoo00<007ooOol0Bgoo 00<007ooOol0:Woo000=Ool00`00Oomoo`2:Ool00`00Oomoo`0MOol00`00Oomoo`1/Ool00`00Oomo o`0fOol00`00Oomoo`1GOol00`00Oomoo`16Ool00`00Oomoo`1:Ool00`00Oomoo`0[Ool000eoo`03 001oogoo08Uoo`03001oogoo01moo`03001oogoo06Yoo`03001oogoo03Moo`03001oogoo05Ioo`03 001oogoo04Qoo`03001oogoo04Uoo`03001oogoo02]oo`003Goo00<007ooOol0R7oo00<007ooOol0 8Goo00<007ooOol0J7oo00<007ooOol0>Goo00<007ooOol0EGoo00<007ooOol0B7oo00<007ooOol0 B7oo00<007ooOol0;7oo000=Ool00`00Oomoo`27Ool00`00Oomoo`0SOol00`00Oomoo`1WOol00`00 Oomoo`0jOol00`00Oomoo`1COol00`00Oomoo`1:Ool00`00Oomoo`17Ool00`00Oomoo`0/Ool000eo o`03001oogoo08Ioo`03001oogoo02Eoo`03001oogoo06Eoo`03001oogoo03]oo`03001oogoo05=o o`03001oogoo04]oo`03001oogoo04Eoo`03001oogoo02eoo`003Goo00<007ooOol0QGoo00<007oo Ool09goo00<007ooOol0I7oo00<007ooOol0?7oo00<007ooOol0DGoo00<007ooOol0C7oo00<007oo Ool0A7oo00<007ooOol0;Woo000=Ool00`00Oomoo`25Ool00`00Oomoo`0WOol00`00Oomoo`1SOol0 0`00Oomoo`0nOol00`00Oomoo`1@Ool00`00Oomoo`1=Ool00`00Oomoo`13Ool00`00Oomoo`0^Ool0 00eoo`03001oogoo08Aoo`03001oogoo02Uoo`03001oogoo069oo`03001oogoo03ioo`03001oogoo 04moo`03001oogoo04moo`03001oogoo045oo`03001oogoo02moo`003Goo00<007ooOol0Pgoo00<0 07ooOol0:goo00<007ooOol0H7oo00<007ooOol0@7oo00<007ooOol0CGoo00<007ooOol0D7oo00<0 07ooOol0@Goo00<007ooOol0;goo000=Ool00`00Oomoo`22Ool00`00Oomoo`0/Ool00`00Oomoo`1O Ool00`00Oomoo`12Ool00`00Oomoo`1Goo00<007ooOol0Woo00<007ooOol0DGoo00<007ooOol0 Cgoo00<007ooOol0?goo00<007ooOol0GWoo00<007ooOol07oo000=Ool00`00Oomoo`0>Ool3001XOol00`00Oomoo`0oOol0 0`00Oomoo`1=Ool00`00Oomoo`1COol00`00Oomoo`0kOol00`00Oomoo`1SOol00`00Oomoo`0/Ool0 0`00Oomoo`0iOol000eoo`03001oogoo00eoo`D006Moo`03001oogoo041oo`03001oogoo04]oo`03 001oogoo05Aoo`03001oogoo03Yoo`03001oogoo06Aoo`03001oogoo02aoo`03001oogoo03Uoo`00 3Goo00<007ooOol03Goo100000=oo`00Ool0I7oo00<007ooOol0@Goo00<007ooOol0Bgoo00<007oo Ool0EGoo00<007ooOol0>Goo00<007ooOol0IGoo00<007ooOol0:Woo00<007ooOol0>Woo000=Ool0 0`00Oomoo`0>Ool20003Ool00`00Oomoo`1ROol00`00Oomoo`12Ool00`00Oomoo`19Ool00`00Oomo o`1FOol00`00Oomoo`0hOol00`00Oomoo`1WOol00`00Oomoo`0XOol00`00Oomoo`0kOol000eoo`03 001oogoo01=oo`03001oogoo065oo`03001oogoo04Aoo`03001oogoo04Moo`03001oogoo05Qoo`03 001oogoo03Moo`03001oogoo06Qoo`03001oogoo02Ioo`03001oogoo03aoo`003Goo00<007ooOol0 57oo00<007ooOol0Ggoo00<007ooOol0AGoo00<007ooOol0Agoo00<007ooOol0FGoo00<007ooOol0 =Goo00<007ooOol0JGoo00<007ooOol09Goo00<007ooOol0?Goo000=Ool00`00Oomoo`0EOol00`00 Oomoo`1NOol00`00Oomoo`16Ool00`00Oomoo`15Ool00`00Oomoo`1JOol00`00Oomoo`0dOol00`00 Oomoo`1[Ool00`00Oomoo`0SOol00`00Oomoo`0nOol000eoo`03001oogoo01Ioo`03001oogoo05ao o`03001oogoo04Qoo`03001oogoo04Aoo`03001oogoo05]oo`03001oogoo039oo`03001oogoo06eo o`03001oogoo029oo`03001oogoo03ioo`003Goo00<007ooOol047oo0P001Goo00<007ooOol0Fgoo 00<007ooOol0B7oo00<007ooOol0@goo00<007ooOol0GGoo00<007ooOol0Ool00`00 Oomoo`0nOol00`00Oomoo`1ROol00`00Oomoo`0[Ool00`00Oomoo`1eOol00`00Oomoo`0GOol20015 Ool000eoo`03001oogoo01Yoo`03001oogoo05Eoo`03001oogoo04moo`03001oogoo03aoo`03001o ogoo06Aoo`03001oogoo02Uoo`03001oogoo07Moo`03001oogoo01Eoo`03001oogoo04Eoo`003Goo 00<007ooOol06goo00<007ooOol0Dgoo00<007ooOol0D7oo00<007ooOol0?7oo00<007ooOol0IGoo 00<007ooOol09goo00<007ooOol0NGoo00<007ooOol04goo00<007ooOol0AWoo000=Ool00`00Oomo o`0LOol00`00Oomoo`1BOol00`00Oomoo`1AOol00`00Oomoo`0jOol00`00Oomoo`1WOol00`00Oomo o`0UOol00`00Oomoo`1kOol00`00Oomoo`0BOol00`00Oomoo`16Ool000eoo`03001oogoo01aoo`03 001oogoo055oo`03001oogoo05=oo`03001oogoo03Qoo`03001oogoo06Qoo`03001oogoo02Aoo`03 001oogoo07eoo`03001oogoo00ioo`<004Uoo`003Goo00<007ooOol07Goo00<007ooOol0D7oo00<0 07ooOol0Dgoo00<007ooOol0>7oo00<007ooOol0JGoo00<007ooOol08Woo00<007ooOol0Ogoo1000 2Woo0P00C7oo000=Ool00`00Oomoo`0NOol00`00Oomoo`1>Ool00`00Oomoo`1EOol00`00Oomoo`0f Ool00`00Oomoo`1[Ool00`00Oomoo`0QOol00`00Oomoo`23Ool:001>Ool000eoo`03001oogoo01io o`03001oogoo04ioo`03001oogoo05Ioo`03001oogoo03Aoo`03001oogoo06eoo`03001oogoo01mo o`03001oogoo0=aoo`003Goo00<007ooOol07goo00<007ooOol0C7oo00<007ooOol0Egoo00<007oo Ool0=7oo00<007ooOol0KWoo00<007ooOol07Goo00<007ooOol0gGoo000=Ool00`00Oomoo`0POol0 0`00Oomoo`1;Ool00`00Oomoo`1HOol00`00Oomoo`0bOol00`00Oomoo`1`Ool00`00Oomoo`0KOol0 0`00Oomoo`3NOol000eoo`03001oogoo021oo`03001oogoo04Yoo`03001oogoo05Yoo`03001oogoo 031oo`03001oogoo079oo`03001oogoo01Uoo`03001oogoo0=moo`003Goo00<007ooOol08Goo00<0 07ooOol0B7oo00<007ooOol0Fgoo00<007ooOol0<7oo00<007ooOol0Lgoo00<007ooOol05goo00<0 07ooOol0h7oo000=Ool00`00Oomoo`0ROol00`00Oomoo`17Ool00`00Oomoo`1LOol00`00Oomoo`0^ Ool00`00Oomoo`1eOol00`00Oomoo`0EOol00`00Oomoo`3QOol000eoo`03001oogoo029oo`03001o ogoo04Ioo`03001oogoo05ioo`03001oogoo02aoo`03001oogoo07Moo`03001oogoo019oo`800>Ao o`003Goo00<007ooOol08goo00<007ooOol0AGoo00<007ooOol0GWoo00<007ooOol0:goo00<007oo Ool0NGoo0P004Goo00<007ooOol0i7oo000=Ool00`00Oomoo`0COol2000?Ool00`00Oomoo`13Ool0 0`00Oomoo`1POol00`00Oomoo`0ZOol00`00Oomoo`1kOol2000=Ool2003WOol000eoo`03001oogoo 019oo`@000ioo`03001oogoo04=oo`03001oogoo065oo`03001oogoo02Qoo`03001oogoo07ioo`<0 00Moo`<00>Uoo`003Goo00<007ooOol04Woo10003goo00<007ooOol0@Goo00<007ooOol0Hgoo00<0 07ooOol09Woo00<007ooOol0PWoo1`00k7oo000=Ool00`00Oomoo`0COol2000AOol00`00Oomoo`10 Ool00`00Oomoo`1TOol00`00Oomoo`0UOol00`00Oomoo`3oOomfOol00003Ool00000008000Qoo`03 001oogoo02Ioo`03001oogoo03moo`03001oogoo06Eoo`03001oogoo02Aoo`03001oogoo0?moogMo o`0000=oo`00Ool02Woo00<007ooOol09goo00<007ooOol0?Woo00<007ooOol0IWoo00<007ooOol0 8Woo00<007ooOol0ogooN7oo0002Ool00`00Oomoo`08Ool5000VOol00`00Oomoo`0lOol00`00Oomo o`1XOol00`00Oomoo`0QOol00`00Oomoo`3oOomhOol000=oo`03001oogoo00Moo`03001oogoo02Qo o`03001oogoo03]oo`03001oogoo06Yoo`03001oogoo01moo`03001oogoo0?moogUoo`0000Eoo`00 Oomoo`0000Qoo`03001oogoo02Uoo`03001oogoo03Yoo`03001oogoo06]oo`03001oogoo01eoo`03 001oogoo0?moogYoo`000Woo0P002Goo00<007ooOol0:Goo00<007ooOol0>Goo00<007ooOol0KGoo 00<007ooOol06goo00<007ooOol0ogooNgoo000=Ool00`00Oomoo`0ZOol00`00Oomoo`0gOol00`00 Oomoo`1_Ool00`00Oomoo`0IOol00`00Oomoo`3oOomlOol000eoo`03001oogoo02]oo`03001oogoo 03Eoo`03001oogoo075oo`03001oogoo01Moo`03001oogoo0?moogeoo`003Goo00<007ooOol0:goo 00<007ooOol0=Goo00<007ooOol0LWoo0P005Woo00<007ooOol0ogooOWoo000=Ool00`00Oomoo`0/ Ool00`00Oomoo`0cOol00`00Oomoo`1eOol00`00Oomoo`0BOol00`00Oomoo`3oOomoOol000eoo`03 001oogoo02eoo`03001oogoo035oo`03001oogoo07Moo`80011oo`800?mooh9oo`003Goo00<007oo Ool0;Goo00<007ooOol0Woo0P00 6Woo00<007ooOol0ogooogoo6Goo000=Ool00`00Oomoo`0lOol00`00Oomoo`0EOol2003oOoooOolL Ool000eoo`03001oogoo03eoo`03001oogoo01=oo`03001oogoo0?mooomooaaoo`003Goo00<007oo Ool0?Woo00<007ooOol047oo0P00ogooogoo7goo000=Ool00`00Oomoo`0oOol00`00Oomoo`0>Ool0 0`00Oomoo`3oOoooOolOOol000eoo`03001oogoo041oo`<000Yoo`<00?mooomoob9oo`003Goo00<0 07ooOol0@goo2P00ogooogoo9Goo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03 001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOomb Ool000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomo o`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo1000ogooogooLGoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`0EOol2003oOoooOomKOol000eoo`03001oogoo01Aoo`@00?mooomooeYoo`00 3Goo00<007ooOol057oo1000ogooogooFWoo000=Ool00`00Oomoo`0EOol2003oOoooOomKOol000eo o`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0 ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo 00<007ooOol0ogooogooLWoo000=Ool4003oOoooOomaOol000eoo`03001oogoo0?mooomoog9oo`00 3Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03 001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOomb Ool000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomo o`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000= Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogoo ogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000=oo`<000Moo`03001oogoo0?mooomoog9oo`00 17oo00<007ooOol01Woo00<007ooOol0ogooogooLWoo00000goo000000030007Ool5003oOoooOom` Ool00005Ool007ooOol00008Ool00`00Oomoo`3oOoooOombOol0009oo`03001oo`0000Qoo`03001o ogoo0?mooomoog9oo`000goo0P0027oo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0 ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo 00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomo og9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001o ogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool4003oOoooOomaOol000eoo`03 001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOomb Ool000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomo o`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000= Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogoo ogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<0 07ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9o o`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo 0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool4003oOoooOomaOol000eoo`03001o ogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol0 00eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3o OoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool0 0`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogoo LWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007oo Ool0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`00 3Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool4003oOoooOomaOol000eoo`03001oogoo 0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eo o`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0 ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo 00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomo og9oo`0000Aoo`00000000Uoo`03001oogoo0?mooomoog9oo`0000Eoo`00Oomoo`0000Qoo`03001o ogoo0?mooomoog9oo`0000Eoo`00Oomoo`0000Qoo`D00?mooomoog1oo`0000Aoo`00000000Uoo`03 001oogoo0?mooomoog9oo`000Woo00<007ooOol027oo00<007ooOol0ogooogooLWoo0002Ool30008 Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogoo ogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<0 07ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9o o`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo 0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eo o`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool4003o OoooOomaOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool0 0`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol067oo0P00 ogooogooF7oo000=Ool00`00Oomoo`0GOol4003oOoooOomGOol000eoo`03001oogoo01Moo`@00?mo oomooeMoo`003Goo00<007ooOol067oo0P00ogooogooF7oo000=Ool00`00Oomoo`3oOoooOombOol0 00eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3o OoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool0 0`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogoo LWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007oo Ool0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`@00?mooomoog5oo`003Goo00<0 07ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9o o`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo 0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eo o`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0 ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`@00?mooomoog5oo`003Goo00<007oo Ool0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`00 3Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03 001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOomb Ool000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomo o`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000= Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`000Woo0P002Goo00<007oo Ool0ogooogooLWoo00001Goo001oogoo000027oo00<007ooOol0ogooogooLWoo00001Goo001oogoo 000027oo1@00ogooogooL7oo0002Ool2003oOoooOomnOol00005Ool007ooOol0003oOoooOommOol0 009oo`800?mooomoogioo`00ogooogooPWoo003oOoooOon2Ool00?mooomooh9oo`00ogooogooPWoo 0000\ \>"], ImageRangeCache->{{{0, 639.438}, {394.813, 0}} -> {-0.556361, -1.06219, \ 0.0403458, 0.023377}}] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Second-order backward difference approximation", "Subsection"], Cell[CellGroupData[{ Cell["Find the update formula", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[k, a, \[Lambda]]; \ order = 2;\), "\[IndentingNewLine]", \(s1backward[q[t], order]\)}], "Input"], Cell[BoxData[ \(\(3\ q[t] + q[\(-2\)\ k + t] - 4\ q[\(-k\) + t]\)\/\(2\ k\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(fmodel[t, q, \[Lambda], a]\)], "Input"], Cell[BoxData[ \(\(-\[Lambda]\)\ q[t] + a\ Sin[t]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(fq\)\(\ \)\(=\)\(\ \)\(First[\ Solve[s1backward[q[t], order]\ \[Equal] \ fmodel[t, q, \[Lambda], a], q[t]]\ ]\)\(\ \)\)\)], "Input"], Cell[BoxData[ \({q[t] \[Rule] \(\(-q[\(-2\)\ k + t]\) + 4\ q[\(-k\) + t] + 2\ a\ k\ \ Sin[t]\)\/\(3 + 2\ k\ \[Lambda]\)}\)], "Output"] }, Open ]], Cell["Define the update rule for the sequence of values Qf2", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[Q, n]\), "\[IndentingNewLine]", \(UpdateRule\ = \ Simplify[\ \(fq\ \ /. \ {t \[Rule] n\ k}\)\ /. \ q[t_] \[Rule] Q[t/k]\ \ ]\)}], "Input"], Cell[BoxData[ \({Q[n] \[Rule] \(\(-Q[\(-2\) + n]\) + 4\ Q[\(-1\) + n] + 2\ a\ k\ Sin[k\ \ n]\)\/\(3 + 2\ k\ \[Lambda]\)}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Clear[Qupdate];\)\), "\[IndentingNewLine]", \(Qupdate[n_]\ = \ Q[n]\ /. \ UpdateRule\)}], "Input"], Cell[BoxData[ \(\(\(-Q[\(-2\) + n]\) + 4\ Q[\(-1\) + n] + 2\ a\ k\ Sin[k\ n]\)\/\(3 + 2\ \ k\ \[Lambda]\)\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Carry out the computation for specific numerical parameters", \ "Subsubsection"], Cell[BoxData[{ \(\(Clear[Q, k, \[Lambda], a];\)\), "\[IndentingNewLine]", \(k = 1. ; \ \[Lambda] = 0.05; \ a = 1. ;\), "\[IndentingNewLine]", \(Q[n_]\ := \ \(Q[n]\ = \ \ Qupdate[n]\)\)}], "Input"], Cell["Set initial conditions", "Text"], Cell[BoxData[ \(For[n = 0, n < order, \(n++\), Q[n] = qmodel[n\ k, \[Lambda], a]]\)], "Input"], Cell["Carry out computations and compare with exact solution", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(Qvalues = Table[{k\ n, Q[n]}, {n, 0, T/k}];\)\), "\[IndentingNewLine]", \(\(Qplot\ = \ ListPlot[Qvalues, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(qplot\ = \ Plot[qmodel[t, \[Lambda], a], {t, 0, T}, DisplayFunction \[Rule] Identity];\)\), "\[IndentingNewLine]", \(\(Show[Qplot, qplot, DisplayFunction \[Rule] $DisplayFunction, PlotRange \[Rule] {{0, T}, {qmin, qmax}}];\)\)}], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -1.14383e-017 0.0397887 0.0686704 0.0686704 [ [.19894 .05617 -3 -9 ] [.19894 .05617 3 0 ] [.39789 .05617 -6 -9 ] [.39789 .05617 6 0 ] [.59683 .05617 -6 -9 ] [.59683 .05617 6 0 ] [.79577 .05617 -6 -9 ] [.79577 .05617 6 0 ] [.99472 .05617 -6 -9 ] [.99472 .05617 6 0 ] [-0.0125 .20601 -6 -4.5 ] [-0.0125 .20601 0 4.5 ] [-0.0125 .34335 -6 -4.5 ] [-0.0125 .34335 0 4.5 ] [-0.0125 .48069 -6 -4.5 ] [-0.0125 .48069 0 4.5 ] [-0.0125 .61803 -6 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .19894 .06867 m .19894 .07492 L s [(5)] .19894 .05617 0 1 Mshowa .39789 .06867 m .39789 .07492 L s [(10)] .39789 .05617 0 1 Mshowa .59683 .06867 m .59683 .07492 L s [(15)] .59683 .05617 0 1 Mshowa .79577 .06867 m .79577 .07492 L s [(20)] .79577 .05617 0 1 Mshowa .99472 .06867 m .99472 .07492 L s [(25)] .99472 .05617 0 1 Mshowa .125 Mabswid .03979 .06867 m .03979 .07242 L s .07958 .06867 m .07958 .07242 L s .11937 .06867 m .11937 .07242 L s .15915 .06867 m .15915 .07242 L s .23873 .06867 m .23873 .07242 L s .27852 .06867 m .27852 .07242 L s .31831 .06867 m .31831 .07242 L s .3581 .06867 m .3581 .07242 L s .43768 .06867 m .43768 .07242 L s .47746 .06867 m .47746 .07242 L s .51725 .06867 m .51725 .07242 L s .55704 .06867 m .55704 .07242 L s .63662 .06867 m .63662 .07242 L s .67641 .06867 m .67641 .07242 L s .7162 .06867 m .7162 .07242 L s .75599 .06867 m .75599 .07242 L s .83556 .06867 m .83556 .07242 L s .87535 .06867 m .87535 .07242 L s .91514 .06867 m .91514 .07242 L s .95493 .06867 m .95493 .07242 L s .25 Mabswid 0 .06867 m 1 .06867 L s 0 .20601 m .00625 .20601 L s [(2)] -0.0125 .20601 1 0 Mshowa 0 .34335 m .00625 .34335 L s [(4)] -0.0125 .34335 1 0 Mshowa 0 .48069 m .00625 .48069 L s [(6)] -0.0125 .48069 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(8)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .10301 m .00375 .10301 L s 0 .13734 m .00375 .13734 L s 0 .17168 m .00375 .17168 L s 0 .24035 m .00375 .24035 L s 0 .27468 m .00375 .27468 L s 0 .30902 m .00375 .30902 L s 0 .37769 m .00375 .37769 L s 0 .41202 m .00375 .41202 L s 0 .44636 m .00375 .44636 L s 0 .51503 m .00375 .51503 L s 0 .54936 m .00375 .54936 L s 0 .5837 m .00375 .5837 L s 0 .03434 m .00375 .03434 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .008 w 0 .13734 Mdot .03979 .16502 Mdot .07958 .21113 Mdot .11937 .22766 Mdot .15915 .19433 Mdot .19894 .13705 Mdot .23873 .10398 Mdot .27852 .12128 Mdot .31831 .169 Mdot .3581 .19941 Mdot .39789 .1809 Mdot .43768 .12701 Mdot .47746 .08397 Mdot .51725 .08821 Mdot .55704 .13283 Mdot .59683 .17397 Mdot .63662 .17109 Mdot .67641 .12426 Mdot .7162 .07409 Mdot .75599 .06437 Mdot .79577 .10182 Mdot .83556 .1499 Mdot .87535 .1624 Mdot .91514 .12591 Mdot .95493 .07218 Mdot .99472 .04887 Mdot .5 Mabswid 0 .13734 m .00123 .13727 L .00236 .13726 L .00365 .13732 L .00488 .13744 L .00601 .1376 L .00706 .13781 L .00942 .13844 L .01181 .13932 L .01444 .14055 L .01975 .14387 L .02493 .14813 L .0305 .15371 L .04057 .16606 L .06325 .20007 L .07443 .21709 L .08481 .2313 L .09483 .24251 L .10043 .24741 L .1055 .25086 L .10841 .25239 L .11109 .2535 L .11261 .25399 L .11403 .25437 L .11551 .25468 L .1163 .2548 L .11714 .2549 L .11787 .25497 L .11855 .25501 L .11929 .25503 L .12008 .25503 L .12141 .25496 L .12216 .25489 L .12287 .2548 L .12418 .25458 L .12538 .25431 L .12809 .25347 L .1307 .25237 L .13307 .25112 L .13839 .24747 L .14331 .24308 L .14796 .2381 L .1697 .20604 L .18893 .17096 L .19986 .15117 L .21009 .13429 L .22099 .11921 L .23096 .1089 L .23357 .10683 L .23642 .10489 L Mistroke .23915 .10335 L .24161 .10222 L .24449 .10124 L .24606 .10085 L .24679 .10071 L .24755 .10059 L .24893 .10043 L .2502 .10035 L .25156 .10035 L .25233 .10039 L .25305 .10044 L .25438 .1006 L .25506 .10071 L .2558 .10086 L .25833 .10152 L .25975 .10202 L .26128 .10264 L .26406 .10402 L .26949 .1076 L .27447 .11185 L .28397 .12224 L .29396 .13576 L .31497 .16867 L .33403 .19704 L .34465 .20951 L .34952 .21405 L .35466 .21793 L .35729 .21952 L .36017 .22095 L .36279 .22195 L .36396 .22231 L .36523 .22263 L .366 .22279 L .36671 .22292 L .36738 .22302 L .36812 .22311 L .3694 .2232 L .37075 .22323 L .37152 .22321 L .37222 .22317 L .37299 .2231 L .37381 .223 L .37518 .22276 L .37668 .22241 L .37933 .22155 L .38221 .22028 L .38745 .21709 L .39254 .21292 L .39724 .20819 L .41923 .17695 L Mistroke .44058 .13881 L .46057 .10517 L .47039 .09191 L .47972 .08227 L .48487 .07838 L .48759 .07676 L .49048 .07537 L .49186 .07483 L .49316 .07441 L .49442 .07406 L .49558 .0738 L .49669 .07361 L .49786 .07346 L .49903 .07338 L .50032 .07335 L .50169 .0734 L .50296 .07352 L .50362 .07361 L .50435 .07373 L .50583 .07405 L .50852 .07486 L .51104 .0759 L .51575 .07853 L .52082 .08234 L .53228 .09426 L .54278 .10846 L .56435 .14277 L .57486 .15948 L .58441 .17325 L .5937 .18452 L .59859 .18939 L .60377 .19362 L .60641 .19539 L .60931 .19701 L .61197 .19821 L .61441 .19904 L .61569 .19937 L .61691 .19963 L .61819 .19983 L .61892 .19991 L .61959 .19997 L .62083 .20002 L .62194 .20001 L .62323 .19993 L .62445 .1998 L .62567 .19959 L .62699 .1993 L .62935 .19859 L .63217 .19743 L .63473 .19608 L Mistroke .64031 .19222 L .64558 .18746 L .65489 .17663 L .665 .16196 L .68541 .12682 L .7079 .08844 L .71817 .07412 L .72383 .0676 L .72922 .06246 L .73421 .05869 L .73684 .05711 L .73965 .05574 L .74225 .05478 L .7446 .05415 L .74574 .05393 L .74697 .05376 L .74813 .05366 L .7492 .05362 L .75053 .05364 L .75178 .05372 L .75291 .05386 L .75413 .05407 L .75546 .05437 L .75688 .05477 L .75942 .0557 L .76395 .05804 L .76896 .06156 L .78014 .07272 L .79065 .08663 L .83126 .1511 L .84215 .16525 L .85229 .17517 L .85741 .17875 L .86027 .1803 L .86287 .18142 L .86426 .1819 L .86573 .18232 L .86649 .18251 L .86734 .18268 L .86811 .18281 L .86882 .18291 L .86954 .18299 L .87022 .18304 L .87094 .18308 L .87173 .18309 L .87241 .18308 L .87314 .18304 L .87444 .18292 L .87575 .18273 L .87719 .18243 L Mistroke .87868 .18203 L .88008 .18156 L .88526 .17912 L .88822 .17723 L .89091 .17521 L .89696 .16964 L .90784 .15644 L .91815 .14093 L .95843 .07244 L .96884 .05831 L .9786 .04819 L .9835 .04446 L .98602 .04291 L .98878 .04152 L .99163 .04042 L .99319 .03996 L .99463 .03964 L .996 .03941 L .99723 .03927 L .99791 .03923 L .99866 .0392 L 1 .03921 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{640.438, 395.813}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg7oo00<007ooOol0F7oo00<007ooOol097oo000=Ool00`00Oomoo`3oOolQOol00`00Oomoo`0WOol0 0`00Oomoo`1UOol00`00Oomoo`0iOol00`00Oomoo`1HOol00`00Oomoo`0TOol000eoo`03001oogoo 0?moob1oo`03001oogoo02Uoo`03001oogoo06Aoo`03001oogoo03Yoo`03001oogoo05Ioo`03001o ogoo02Eoo`003Goo00<007ooOol0UWoo2000P7oo00<007ooOol0:goo00<007ooOol0HWoo00<007oo Ool0>7oo0P000Woo00<007ooOol0EGoo00<007ooOol09Goo000=Ool00`00Oomoo`2AOol200000goo 00000008Ool3001lOol00`00Oomoo`0/Ool00`00Oomoo`1ROol00`00Oomoo`0gOol400000goo001o o`1EOol00`00Oomoo`0VOol000eoo`@008moo`@000eoo`8007Yoo`03001oogoo02eoo`03001oogoo 061oo`03001oogoo03Qoo`@0009oo`03001oogoo059oo`03001oogoo02Moo`003Goo00<007ooOol0 T7oo10003goo0P00Mgoo00<007ooOol0;goo00<007ooOol0Ggoo00<007ooOol0>Goo0P000goo00<0 07ooOol0DWoo00<007ooOol09goo000=Ool00`00Oomoo`2?Ool4000BOol00`00Oomoo`1cOol00`00 Oomoo`0`Ool00`00Oomoo`1NOol00`00Oomoo`10Ool00`00Oomoo`1@Ool00`00Oomoo`0XOol000eo o`03001oogoo08eoo`8001Moo`03001oogoo079oo`03001oogoo035oo`03001oogoo05aoo`03001o ogoo049oo`03001oogoo04moo`03001oogoo02Qoo`003Goo00<007ooOol0S7oo00<007ooOol067oo 00<007ooOol0L7oo00<007ooOol0Woo00<007ooOol0Dgoo00<007ooOol0BWoo00<0 07ooOol0Agoo00<007ooOol0;7oo000=Ool00`00Oomoo`26Ool00`00Oomoo`0POol400000goo001o o`1VOol00`00Oomoo`0kOol00`00Oomoo`1COol00`00Oomoo`1;Ool00`00Oomoo`15Ool00`00Oomo o`0]Ool000eoo`03001oogoo08Eoo`03001oogoo025oo`@0009oo`03001oogoo06Aoo`03001oogoo 03aoo`03001oogoo04Yoo`8000Eoo`03001oogoo04aoo`03001oogoo04Aoo`03001oogoo02ioo`00 3Goo00<007ooOol0QGoo00<007ooOol08Woo0P000goo00<007ooOol0Fgoo0P001Woo00<007ooOol0 ?Woo00<007ooOol0B7oo100017oo00<007ooOol0CGoo00<007ooOol0>Woo0P001goo00<007ooOol0 ;Woo000=Ool00`00Oomoo`24Ool00`00Oomoo`0YOol00`00Oomoo`1IOol40005Ool00`00Oomoo`0n Ool00`00Oomoo`18Ool40003Ool00`00Oomoo`1?Ool00`00Oomoo`0hOol40005Ool00`00Oomoo`0_ Ool000eoo`03001oogoo08=oo`03001oogoo02]oo`03001oogoo05Qoo`@000Aoo`03001oogoo041o o`03001oogoo04Qoo`8000=oo`03001oogoo051oo`03001oogoo03Qoo`@000Eoo`03001oogoo02mo o`003Goo00<007ooOol0PWoo00<007ooOol0;7oo00<007ooOol0FGoo0P0017oo00<007ooOol0@Woo 00<007ooOol0C7oo00<007ooOol0DGoo00<007ooOol0>7oo0P001Goo00<007ooOol0<7oo000=Ool0 0`00Oomoo`22Ool00`00Oomoo`0]Ool00`00Oomoo`1NOol00`00Oomoo`11Ool2001=Ool00`00Oomo o`1BOol00`00Oomoo`0nOol00`00Oomoo`0aOol000eoo`03001oogoo085oo`03001oogoo02moo`03 001oogoo05aoo`03001oogoo045oo`@004aoo`03001oogoo05=oo`03001oogoo03eoo`03001oogoo 035oo`003Goo00<007ooOol0P7oo00<007ooOol0<7oo00<007ooOol0G7oo00<007ooOol0@Goo1@00 BWoo00<007ooOol0EGoo00<007ooOol0>goo00<007ooOol07oo00<007ooOol0goo00<0 07ooOol0D7oo00<007ooOol0D7oo00<007ooOol0?Woo00<007ooOol0Ggoo00<007ooOol00P00;goo 00<007ooOol0=goo000=Ool00`00Oomoo`0>Ool00`00Oomoo`1ZOol00`00Oomoo`0lOol00`00Oomo o`1?Ool00`00Oomoo`1AOol00`00Oomoo`0mOol00`00Oomoo`1POol5000]Ool00`00Oomoo`0hOol0 00eoo`03001oogoo00moo`03001oogoo06Qoo`03001oogoo03eoo`03001oogoo04ioo`03001oogoo 059oo`03001oogoo03aoo`03001oogoo069oo`@002eoo`03001oogoo03Qoo`003Goo00<007ooOol0 47oo00<007ooOol0IWoo00<007ooOol0?goo00<007ooOol0CGoo00<007ooOol0Dgoo00<007ooOol0 >goo00<007ooOol0Hgoo0P00;Goo00<007ooOol0>Goo000=Ool00`00Oomoo`0AOol00`00Oomoo`1U Ool00`00Oomoo`10Ool00`00Oomoo`1;Ool00`00Oomoo`1DOol00`00Oomoo`0jOol00`00Oomoo`1T Ool00`00Oomoo`0/Ool00`00Oomoo`0iOol000eoo`03001oogoo019oo`03001oogoo06=oo`03001o ogoo045oo`03001oogoo04]oo`03001oogoo05Eoo`03001oogoo03Uoo`03001oogoo06Eoo`03001o ogoo02Yoo`03001oogoo03Yoo`003Goo00<007ooOol04goo00<007ooOol0HWoo00<007ooOol0@Woo 00<007ooOol0BGoo00<007ooOol0EWoo00<007ooOol0>7oo00<007ooOol0Igoo00<007ooOol0:7oo 00<007ooOol0>goo000=Ool00`00Oomoo`0COol00`00Oomoo`1QOol00`00Oomoo`14Ool00`00Oomo o`17Ool00`00Oomoo`1HOol00`00Oomoo`0gOol00`00Oomoo`1XOol00`00Oomoo`0COol2000AOol0 0`00Oomoo`0lOol000eoo`03001oogoo01Aoo`03001oogoo05moo`03001oogoo04Eoo`03001oogoo 04Moo`03001oogoo05Uoo`03001oogoo03Eoo`03001oogoo06Uoo`03001oogoo019oo`@000moo`03 001oogoo03eoo`003Goo00<007ooOol05Goo0P00Ggoo00<007ooOol0AWoo00<007ooOol0AGoo00<0 07ooOol0FWoo00<007ooOol0=7oo00<007ooOol0Jgoo00<007ooOol04Goo10003Woo00<007ooOol0 ?Woo000=Ool00`00Oomoo`0DOol4001MOol00`00Oomoo`18Ool00`00Oomoo`14Ool00`00Oomoo`1K Ool00`00Oomoo`0bOol00`00Oomoo`1]Ool00`00Oomoo`0AOol2000?Ool00`00Oomoo`0nOol000eo o`03001oogoo01Aoo`@005eoo`03001oogoo04Qoo`03001oogoo04=oo`03001oogoo05eoo`03001o ogoo035oo`03001oogoo06ioo`03001oogoo021oo`03001oogoo03moo`003Goo00<007ooOol05Goo 0P0000=oo`00Ool0FWoo00<007ooOol0BWoo00@007oo0000041oo`03001oogoo05moo`03001oogoo 02moo`03001oogoo071oo`03001oogoo01ioo`03001oogoo041oo`003Goo00<007ooOol067oo00<0 07ooOol0FGoo00<007ooOol0Bgoo1000?goo00<007ooOol0Ggoo00<007ooOol07goo0P003Goo00<0 07ooOol0LWoo00<007ooOol077oo00<007ooOol0@Goo000=Ool4000HOol00`00Oomoo`1GOol00`00 Oomoo`1Ool3 0019Ool000eoo`03001oogoo01eoo`03001oogoo051oo`03001oogoo05=oo`03001oogoo02Yoo`@0 00Yoo`03001oogoo06Uoo`03001oogoo029oo`03001oogoo07moo`@000Yoo`8004aoo`003Goo00<0 07ooOol07Woo00<007ooOol0CWoo00<007ooOol0EGoo00<007ooOol0:Goo10002Goo00<007ooOol0 Jgoo00<007ooOol08Goo00<007ooOol0Pgoo2P00CWoo000=Ool00`00Oomoo`0NOol00`00Oomoo`1> Ool00`00Oomoo`1FOol00`00Oomoo`0YOol20009Ool00`00Oomoo`1]Ool00`00Oomoo`0OOol00`00 Oomoo`3LOol000eoo`03001oogoo01moo`03001oogoo04aoo`03001oogoo05Moo`03001oogoo03Ao o`03001oogoo06ioo`03001oogoo01eoo`03001oogoo0=eoo`003Goo00<007ooOol087oo00<007oo Ool0Bgoo00<007ooOol0F7oo00<007ooOol01oo`003Goo00<007ooOol08Woo00<007ooOol0>Woo0P002goo00<007ooOol0G7oo00<007ooOol0 ;Woo00<007ooOol0MGoo00<007ooOol05Goo00<007ooOol0hGoo000=Ool00`00Oomoo`0ROol00`00 Oomoo`0iOol40009Ool00`00Oomoo`1NOol00`00Oomoo`0/Ool00`00Oomoo`1gOol00`00Oomoo`0B Ool2003TOol000eoo`03001oogoo02=oo`03001oogoo03Qoo`@000Uoo`03001oogoo05ioo`03001o ogoo02]oo`03001oogoo07Uoo`80015oo`03001oogoo0>Aoo`003Goo00<007ooOol097oo00<007oo Ool0>7oo0P002Goo00<007ooOol0H7oo00<007ooOol0:Woo00<007ooOol0Ngoo0P003Goo0P00igoo 000=Ool00`00Oomoo`0TOol00`00Oomoo`13Ool00`00Oomoo`1QOol00`00Oomoo`0:Ool2000LOol0 0`00Oomoo`1nOol30007Ool3003YOol000eoo`03001oogoo02Eoo`03001oogoo045oo`03001oogoo 06=oo`03001oogoo00Qoo`@001Yoo`03001oogoo089oo`L00>aoo`003Goo00<007ooOol09Woo00<0 07ooOol0@7oo00<007ooOol0I7oo00<007ooOol01goo10006Woo00<007ooOol0ogooMWoo00000goo 000000020008Ool00`00Oomoo`0VOol00`00Oomoo`0oOol00`00Oomoo`1UOol00`00Oomoo`08Ool2 000JOol00`00Oomoo`3oOomgOol00003Ool007oo00Yoo`03001oogoo02Moo`03001oogoo03ioo`03 001oogoo06Ioo`03001oogoo029oo`03001oogoo0?moogQoo`000Woo00<007ooOol027oo1@009Woo 00<007ooOol0?7oo00<007ooOol0J7oo00<007ooOol08Goo00<007ooOol0ogooN7oo0003Ool00`00 Oomoo`07Ool00`00Oomoo`0XOol00`00Oomoo`0kOol00`00Oomoo`1ZOol00`00Oomoo`0OOol00`00 Oomoo`3oOomiOol00005Ool007ooOol00008Ool00`00Oomoo`0YOol00`00Oomoo`02Ool2000fOol0 0`00Oomoo`1[Ool00`00Oomoo`0MOol00`00Oomoo`3oOomjOol0009oo`8000Uoo`03001oogoo02Uo o`04001oogooOol4000dOol00`00Oomoo`1]Ool00`00Oomoo`0KOol00`00Oomoo`3oOomkOol000eo o`03001oogoo02Yoo`03001oogoo00@003=oo`03001oogoo06moo`03001oogoo01Uoo`03001oogoo 0?moogaoo`003Goo00<007ooOol0:goo00<007ooOol00P00Goo00<007ooOol06goo00<007ooOol0ogooogoo67oo000=Ool00`00Oomoo`0jOol2000J Ool00`00Oomoo`3oOoooOolIOol000eoo`03001oogoo03aoo`03001oogoo01Eoo`800?mooomooaao o`003Goo00<007ooOol0?Goo00<007ooOol04goo00<007ooOol0ogooogoo77oo000=Ool00`00Oomo o`0nOol00`00Oomoo`0@Ool2003oOoooOolOOol000eoo`03001oogoo03moo`03001oogoo00ioo`03 001oogoo0?mooomooamoo`003Goo00<007ooOol0@7oo0`002Woo0`00ogooogoo8Woo000=Ool00`00 Oomoo`13Ool:003oOoooOolUOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogoo ogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<0 07ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9o o`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo 0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool4003oOoooOomaOol000eoo`03001o ogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol0 00eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3o OoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool0 0`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogoo LWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007oo Ool0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`00 3Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`@0 0?mooomoog5oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eo o`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0 ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo 00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomo og9oo`000goo0`001goo00<007ooOol0ogooogooLWoo0004Ool00`00Oomoo`06Ool00`00Oomoo`3o OoooOombOol00003Ool0000000<000Moo`D00?mooomoog1oo`0000Eoo`00Oomoo`0000Qoo`03001o ogoo0?mooomoog9oo`000Woo00<007oo000027oo00<007ooOol0ogooogooLWoo0003Ool20008Ool0 0`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogoo LWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007oo Ool0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`00 3Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03 001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOomb Ool000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomo o`3oOoooOombOol000eoo`@00?mooomoog5oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0 ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo 00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomo og9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001o ogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol0 00eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3o OoooOombOol000eoo`@00?mooomoog5oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomo o`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000= Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogoo ogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<0 07ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9o o`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo 0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eo o`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`@00?mooomoog5oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3o OoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool0 0`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogoo LWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007oo Ool0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`00 3Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03 001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000017oo000000002Goo00<007oo Ool0ogooogooLWoo00001Goo001oogoo000027oo00<007ooOol0ogooogooLWoo00001Goo001oogoo 000027oo1@00ogooogooL7oo000017oo000000002Goo00<007ooOol0ogooogooLWoo0002Ool00`00 Oomoo`08Ool00`00Oomoo`3oOoooOombOol0009oo`<000Qoo`03001oogoo0?mooomoog9oo`003Goo 00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomo og9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001o ogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol0 00eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3o OoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool0 0`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogoo LWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`@00?mooomoog5oo`003Goo00<007ooOol0ogoo ogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<0 07ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9o o`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo 0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eo o`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOooo OombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00 Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo1000ogooogoo LGoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007oo Ool0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`00 3Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mo oomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03 001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOomb Ool000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomo o`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000= Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo1000ogooogooLGoo 000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0 ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo 00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001oogoo0?mooomo og9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol000eoo`03001o ogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3oOoooOombOol0 00eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo000=Ool00`00Oomoo`3o OoooOombOol000eoo`03001oogoo0?mooomoog9oo`003Goo00<007ooOol0ogooogooLWoo0002Ool2 0009Ool00`00Oomoo`3oOoooOombOol00005Ool007ooOol00008Ool00`00Oomoo`3oOoooOombOol0 0005Ool007ooOol00008Ool5003oOoooOom`Ool0009oo`800?mooomoogioo`0000Eoo`00Oomoo`00 0?mooomoogeoo`000Woo0P00ogooogooOWoo003oOoooOon2Ool00?mooomooh9oo`00ogooogooPWoo 003oOoooOon2Ool00001\ \>"], ImageRangeCache->{{{0, 639.438}, {394.813, 0}} -> {-0.556361, -1.06219, \ 0.0403458, 0.023377}}] }, Open ]] }, Open ]] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell["2. Stability of the above algorithms", "Section"], Cell[CellGroupData[{ Cell["First-order forward difference approximation", "Subsection"], Cell[CellGroupData[{ Cell["Find the update formula", "Subsubsection"], Cell[CellGroupData[{ Cell[BoxData[{ \(Clear[k, a, \[Lambda], Q]; \ order = 1; \ a = 0;\), "\[IndentingNewLine]", \(fd = s1forward[q[t], order]\ - fmodel[t, q, \[Lambda], a]\)}], "Input"], Cell[BoxData[ \(\[Lambda]\ q[t] + \(\(-q[t]\) + q[k + t]\)\/k\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(CharPoly\ = \ Simplify[\ \(\(k\ fd\ \ /. \ {t \[Rule] n\ k}\)\ /. \ q[t_] \[Rule] Q[t/k]\)\ /. \ \[Lambda] \[Rule] z/k\ ]\)], "Input"], Cell[BoxData[ \(\((\(-1\) + z)\)\ Q[n] + Q[1 + n]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(rhs\ = \ Coefficient[CharPoly, z]\)], "Input"], Cell[BoxData[ \(Q[n]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(lhs\ = \ Simplify[CharPoly\ - \ z*rhs]\)], "Input"], Cell[BoxData[ \(\(-Q[n]\) + Q[1 + n]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Q[n_]\ = \ \[Zeta]^n\)], "Input"], Cell[BoxData[ \(\[Zeta]\^n\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\[Rho][\[Zeta]_] = \ lhs\ /. \ n \[Rule] 0\)], "Input"], Cell[BoxData[ \(\(-1\) + \[Zeta]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\[Sigma][\[Zeta]_] = rhs\ /. \ n \[Rule] 0\)], "Input"], Cell[BoxData[ \(1\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\[Pi]1[\[Zeta]_, z_] = \[Rho][\[Zeta]]\ - \ z\ \[Sigma][\[Zeta]]\)], "Input"], Cell[BoxData[ \(\(-1\) - z + \[Zeta]\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[\[Pi]1[\[ExponentialE]\^\[ImaginaryI]\[Theta], z] \[Equal] 0, z]\)], "Input"], Cell[BoxData[ \({{z \[Rule] \(-1\) + \[ExponentialE]\^\[ImaginaryI]\[Theta]}}\)], \ "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(z[\[Theta]_]\ = \ z\ /. \ First[ Solve[\[Pi]1[\[ExponentialE]\^\(\[ImaginaryI]\ \[Theta]\), z] \[Equal] 0, z]]\)], "Input"], Cell[BoxData[ \(\(-1\) + \[ExponentialE]\^\(\[ImaginaryI]\ \[Theta]\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(ParametricPlot[{Re[z[\[Theta]]], Im[z[\[Theta]]]}, {\[Theta], 0, 2 \[Pi]}, AspectRatio \[Rule] Automatic];\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: 1 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.97619 0.476192 0.500001 0.476192 [ [.02381 .4875 -6 -9 ] [.02381 .4875 6 0 ] [.2619 .4875 -12 -9 ] [.2619 .4875 12 0 ] [.5 .4875 -6 -9 ] [.5 .4875 6 0 ] [.73809 .4875 -12 -9 ] [.73809 .4875 12 0 ] [.96369 .02381 -12 -4.5 ] [.96369 .02381 0 4.5 ] [.96369 .26191 -24 -4.5 ] [.96369 .26191 0 4.5 ] [.96369 .7381 -18 -4.5 ] [.96369 .7381 0 4.5 ] [.96369 .97619 -6 -4.5 ] [.96369 .97619 0 4.5 ] [ 0 0 0 0 ] [ 1 1 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .02381 .5 m .02381 .50625 L s [(-2)] .02381 .4875 0 1 Mshowa .2619 .5 m .2619 .50625 L s [(-1.5)] .2619 .4875 0 1 Mshowa .5 .5 m .5 .50625 L s [(-1)] .5 .4875 0 1 Mshowa .73809 .5 m .73809 .50625 L s [(-0.5)] .73809 .4875 0 1 Mshowa .125 Mabswid .07143 .5 m .07143 .50375 L s .11904 .5 m .11904 .50375 L s .16666 .5 m .16666 .50375 L s .21428 .5 m .21428 .50375 L s .30952 .5 m .30952 .50375 L s .35714 .5 m .35714 .50375 L s .40476 .5 m .40476 .50375 L s .45238 .5 m .45238 .50375 L s .54762 .5 m .54762 .50375 L s .59524 .5 m .59524 .50375 L s .64286 .5 m .64286 .50375 L s .69048 .5 m .69048 .50375 L s .78571 .5 m .78571 .50375 L s .83333 .5 m .83333 .50375 L s .88095 .5 m .88095 .50375 L s .92857 .5 m .92857 .50375 L s .25 Mabswid 0 .5 m 1 .5 L s .97619 .02381 m .98244 .02381 L s [(-1)] .96369 .02381 1 0 Mshowa .97619 .26191 m .98244 .26191 L s [(-0.5)] .96369 .26191 1 0 Mshowa .97619 .7381 m .98244 .7381 L s [(0.5)] .96369 .7381 1 0 Mshowa .97619 .97619 m .98244 .97619 L s [(1)] .96369 .97619 1 0 Mshowa .125 Mabswid .97619 .07143 m .97994 .07143 L s .97619 .11905 m .97994 .11905 L s .97619 .16667 m .97994 .16667 L s .97619 .21429 m .97994 .21429 L s .97619 .30952 m .97994 .30952 L s .97619 .35714 m .97994 .35714 L s .97619 .40476 m .97994 .40476 L s .97619 .45238 m .97994 .45238 L s .97619 .54762 m .97994 .54762 L s .97619 .59524 m .97994 .59524 L s .97619 .64286 m .97994 .64286 L s .97619 .69048 m .97994 .69048 L s .97619 .78572 m .97994 .78572 L s .97619 .83334 m .97994 .83334 L s .97619 .88096 m .97994 .88096 L s .97619 .92857 m .97994 .92857 L s .25 Mabswid .97619 0 m .97619 1 L s 0 0 m 1 0 L 1 1 L 0 1 L closepath clip newpath .5 Mabswid .97619 .5 m .97618 .50369 L .97614 .50705 L .97606 .51093 L .97597 .5146 L .97572 .52112 L .97536 .52818 L .97483 .53596 L .97414 .5442 L .97253 .55894 L .97039 .57409 L .96795 .58822 L .96081 .62007 L .95129 .65197 L .94065 .68052 L .91017 .74191 L .87673 .79126 L .83352 .83988 L .78673 .88019 L .73155 .91611 L .7028 .93085 L .67066 .94456 L .64214 .95449 L .61017 .96327 L .5949 .96664 L .57842 .96969 L .5641 .97186 L .54845 .97372 L .53926 .97457 L .53074 .9752 L .52259 .97566 L .51847 .97584 L .51392 .97599 L .50983 .97609 L .50609 .97615 L .50197 .97619 L .4999 .97619 L .49764 .97619 L .49381 .97615 L .48973 .97608 L .48587 .97598 L .48236 .97587 L .47512 .97554 L .46715 .97506 L .45109 .97367 L .43626 .97191 L .42244 .96983 L .39155 .96368 L .3589 .95481 L .32884 .94437 L Mistroke .29695 .93073 L .24178 .9001 L .19207 .86324 L .14489 .81726 L .10369 .764 L .07215 .70905 L .05955 .681 L .04825 .65059 L .03956 .62147 L .03328 .59452 L .02828 .56509 L .02631 .54879 L .0256 .54128 L .02497 .53331 L .02454 .52649 L .02436 .52287 L .02418 .51897 L .02406 .5154 L .02396 .51213 L .02389 .509 L .02384 .50571 L .02381 .50211 L .02381 .49827 L .02384 .49465 L .02388 .49136 L .02397 .48749 L .0241 .4833 L .02443 .47571 L .02486 .46828 L .02546 .46033 L .02688 .446 L .02913 .42902 L .03171 .41361 L .03905 .38047 L .04818 .3496 L .07174 .29179 L .10012 .24143 L .13486 .19433 L .17963 .14769 L .22702 .10982 L .28309 .07608 L .3408 .05121 L .37094 .04163 L .38817 .03713 L .40459 .03347 L .41978 .03062 L .43635 .02808 L .45149 .02629 L .4655 .02506 L .47334 .02456 L Mistroke .48068 .0242 L .48489 .02405 L .48874 .02394 L .49087 .0239 L .49319 .02386 L .49538 .02383 L .49741 .02382 L .49954 .02381 L .50184 .02381 L .50401 .02383 L .50601 .02385 L .5098 .02391 L .5138 .02401 L .5181 .02415 L .52212 .02432 L .53124 .02484 L .5394 .02544 L .54693 .02613 L .56393 .02812 L .57886 .03038 L .59461 .0333 L .62161 .0396 L .65091 .04836 L .70722 .07126 L .76408 .10374 L .81405 .14205 L .86069 .18909 L .89804 .23861 L .92672 .28866 L .94037 .3188 L .95129 .34804 L .95954 .37519 L .9665 .4044 L .9707 .42788 L .97256 .44131 L .97395 .45391 L .97494 .4655 L .97559 .47611 L .97585 .48197 L .97596 .48507 L .97605 .48836 L .97612 .49153 L .97616 .4944 L .97618 .49729 L .97619 .5 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 288}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg001jOol00`00Oomoo`04Ool20006Ool0081oo`T0 00ioo`P006Yoo`@000Aoo`03001oogoo00Aoo`03001oogoo00Eoo`00Mgoo2@007goo2@00JGoo00<0 07ooOol017oo00<007ooOol01Goo001bOol5000aOol4001TOol20006Ool00`00Oomoo`05Ool006eo o`D003Yoo`@006Qoo`03001oogoo00Eoo`00JGoo1000@goo1000I7oo00<007ooOol01Goo001VOol3 001;Ool4001POol00`00Oomoo`05Ool006=oo`<0059oo`@005aoo`03001oogoo00Eoo`00H7oo0`00 FGoo0P00FWoo00<007ooOol01Goo001NOol2001NOol2001HOol00`00Oomoo`05Ool005aoo`80069o o`8005Ioo`03001oogoo00Eoo`00FWoo0P00IWoo0`00Dgoo00<007ooOol01Goo001GOol3001[Ool2 001AOol00`00Oomoo`05Ool005Eoo`80071oo`8004moo`03001oogoo00Eoo`00Dgoo0P00M7oo0P00 CGoo0P001Woo001AOol2001hOol2001;Ool00`00Oomoo`05Ool004moo`8007aoo`8004Uoo`03001o ogoo00Eoo`00CWoo00<007ooOol0OWoo0P00Agoo00<007ooOol01Goo001Ool2000oOol00`00Oomoo`05Ool004Ioo`03001oogoo091o o`8003eoo`03001oogoo00Eoo`00A7oo0P00UGoo00<007ooOol0>Woo00<007ooOol01Goo0012Ool2 002HOol00`00Oomoo`0iOol00`00Oomoo`05Ool0045oo`03001oogoo09Uoo`8003Uoo`03001oogoo 00Eoo`00@7oo00<007ooOol0W7oo00<007ooOol0=Woo00<007ooOol01Goo000nOol2002POol00`00 Oomoo`0eOol00`00Oomoo`05Ool003eoo`03001oogoo0:5oo`03001oogoo03Aoo`8000Ioo`00?7oo 00<007ooOol0Xgoo0P00=7oo00<007ooOol01Goo000kOol00`00Oomoo`2VOol00`00Oomoo`0aOol0 0`00Oomoo`05Ool003Uoo`800:Yoo`03001oogoo031oo`03001oogoo00Eoo`00>7oo00<007ooOol0 Zgoo0P00<7oo00<007ooOol01Goo000gOol00`00Oomoo`2^Ool00`00Oomoo`0]Ool00`00Oomoo`05 Ool003Eoo`800;9oo`03001oogoo02aoo`03001oogoo00Eoo`00=7oo00<007ooOol0/goo00<007oo Ool0:goo00<007ooOol01Goo000cOol00`00Oomoo`2eOol00`00Oomoo`0ZOol00`00Oomoo`05Ool0 039oo`03001oogoo0;Moo`03001oogoo02Uoo`03001oogoo00Eoo`001oo`03001oogoo01Eoo`03001oogoo00Eoo`007Goo00<0 07ooOol0hGoo00<007ooOol057oo00<007ooOol01Goo000LOol00`00Oomoo`3SOol00`00Oomoo`0C Ool00`00Oomoo`05Ool001]oo`03001oogoo0>Aoo`03001oogoo01=oo`03001oogoo00Eoo`006goo 00<007ooOol0iGoo00<007ooOol04Woo00<007ooOol01Goo000JOol00`00Oomoo`3VOol00`00Oomo o`0BOol00`00Oomoo`05Ool001Yoo`03001oogoo0>Moo`8000Aoo`8000Aoo`8000Ioo`03001oogoo 00Eoo`006Goo00<007ooOol0igoo0P0000=oo`00Ool01goo00@007ooOol000Eoo`03001oogoo00Eo o`006Goo00<007ooOol0igoo00@007oo000000]oo`03001oogoo00=oo`8000Ioo`0067oo00<007oo Ool0hWoo10000Woo00@007ooOol000Uoo`8000Ioo`03001oogoo00Eoo`005goo00<007ooOol0jGoo 00@007ooOol000Uoo`03001oogoo00Eoo`03001oogoo00Eoo`005goo00<007ooOol0jWoo0P0000=o o`00Ool01goo0`001Goo00<007ooOol01Goo000FOol00`00Oomoo`3^Ool00`00Oomoo`0>Ool00`00 Oomoo`05Ool001Ioo`03001oogoo0>moo`03001oogoo00eoo`03001oogoo00Eoo`005Goo00<007oo Ool0l7oo00<007ooOol03Goo00<007ooOol01Goo000EOol00`00Oomoo`3aOol00`00Oomoo`0Ool00`00Oomoo`3oOol00`00 Oomoo`05Ool00`00Oomoo`05Ool000ioo`03001oogoo0?moo`03001oogoo00Eoo`03001oogoo00Eo o`003Goo00<007ooOol0ogoo0Woo00<007ooOol017oo00<007ooOol01Goo000=Ool00`00Oomoo`3o Ool2Ool00`00Oomoo`04Ool00`00Oomoo`05Ool000eoo`03001oogoo0?moo`9oo`03001oogoo00Ao o`03001oogoo00Eoo`003Goo00<007ooOol0ogoo0goo00<007ooOol00goo0P001Woo000Ool20007Ool000Moo`03001oogoo0?moo`io o`8000Moo`001goo00<007ooOol0ogoo3Woo0`001Woo0007Ool00`00Oomoo`3oOol>Ool20007Ool0 00Moo`@003]oo`D0009oo`8000Aoo`8003]oo`D003Yoo`8000Aoo`8000Aoo`8003Yoo`8000Moo`00 1goo00<007ooOol0?Woo00<007ooOol01goo00@007ooOol003aoo`03001oogoo03Uoo`04001oogoo 0008Ool01000Oomoo`00>Goo0P001goo0007Ool2000oOol00`00Oomoo`0:Ool00`00Oomoo`0jOol0 0`00Oomoo`0iOol01000Oomoo`002goo00<007ooOol0=goo0P001goo00000goo000000020002Ool0 0`00Ool0000fOol40004Ool00`00Oomoo`08Ool2000eOol40004Ool00`00Oomoo`0cOol40002Ool0 1000Oomoo`002Goo0P00>Woo0P001goo0007Ool01000Oomoo`00?Goo00<007ooOol027oo00<007oo Ool0?7oo00<007ooOol0>Goo00@007ooOol000Uoo`03001oogoo03Uoo`8000Moo`001goo0`00?Goo 0P002Woo0`00>goo0P00?7oo0P002Woo0`00>Goo0P001goo0006Ool00`00Oomoo`3oOol?Ool20007 Ool000Ioo`03001oogoo0?mooa1oo`03001oogoo00Eoo`001Woo00<007ooOol0ogoo47oo00<007oo Ool01Goo0006Ool00`00Oomoo`3oOol@Ool00`00Oomoo`05Ool000Ioo`03001oogoo0?mooa1oo`03 001oogoo00Eoo`001Woo00<007ooOol0ogoo47oo00<007ooOol01Goo003o000Q000000Ioo`03001o ogoo00]oo`03001oogoo00]oo`03001oogoo00Yoo`03001oogoo00]oo`03001oogoo00]oo`03001o ogoo00Yoo`03001oogoo00]oo`03001oogoo00]oo`03001oogoo00Yoo`03001oogoo00]oo`03001o ogoo00]oo`03001oogoo00Yoo`03001oogoo00]oo`03001oogoo00]oo`03001oogoo00Yoo`03001o ogoo00]oo`03001oogoo00]oo`03001oogoo00Yoo`03001oogoo00]oo`03001oogoo00]oo`03001o ogoo00Eoo`001Woo00<007ooOol0@Woo00<007ooOol0@Goo00<007ooOol0@Goo00<007ooOol0@Woo 00<007ooOol01Goo0006Ool00`00Oomoo`3oOol@Ool00`00Oomoo`05Ool000Ioo`03001oogoo0?mo oa1oo`03001oogoo00Eoo`001Woo00<007ooOol0ogoo47oo00<007ooOol01Goo0006Ool00`00Oomo o`3oOol@Ool00`00Oomoo`05Ool000Ioo`03001oogoo0?moo`moo`8000Moo`001goo00<007ooOol0 ogoo3Woo0P001goo0007Ool00`00Oomoo`3oOol>Ool20007Ool000Moo`03001oogoo0?moo`ioo`80 00Moo`001goo00<007ooOol0ogoo3Woo0P001goo0007Ool00`00Oomoo`3oOol>Ool20007Ool000Mo o`03001oogoo0?moo`ioo`8000Moo`001goo00<007ooOol0ogoo3Woo0`001Woo0007Ool00`00Oomo o`3oOol>Ool20007Ool000Moo`03001oogoo0?moo`ioo`8000Moo`0027oo00<007ooOol0ogoo3Goo 0P001goo0008Ool00`00Oomoo`3oOol=Ool20007Ool000Qoo`03001oogoo0?moo`aoo`03001oo`00 00Moo`0027oo00<007ooOol0ogoo37oo00<007oo00001goo0008Ool00`00Oomoo`3oOolmoo`03001oogoo00eoo`03 001oogoo00Eoo`005Woo00<007ooOol0kWoo00<007ooOol03Woo00<007ooOol01Goo000GOol00`00 Oomoo`3]Ool00`00Oomoo`0>Ool00`00Oomoo`05Ool001Moo`03001oogoo0>Yoo`<000=oo`8000Ao o`8000Ioo`03001oogoo00Eoo`0067oo00<007ooOol0j7oo00@007ooOol000Qoo`04001oogoo0005 Ool00`00Oomoo`05Ool001Qoo`03001oogoo0>Qoo`04001oo`00000;Ool00`00Oomoo`03Ool20006 Ool001Uoo`03001oogoo0>Moo`04001oo`000009Ool20006Ool00`00Oomoo`05Ool001Yoo`03001o ogoo0>Ioo`800003Ool007oo00Qoo`03001oogoo00Eoo`03001oogoo00Eoo`006Woo00<007ooOol0 igoo0P002Woo0`001Goo00<007ooOol01Goo000KOol00`00Oomoo`3UOol00`00Oomoo`0BOol00`00 Oomoo`05Ool001]oo`03001oogoo0>Aoo`03001oogoo01=oo`03001oogoo00Eoo`0077oo00<007oo Ool0hgoo00<007ooOol04goo00<007ooOol01Goo000LOol00`00Oomoo`3ROol00`00Oomoo`0DOol0 0`00Oomoo`05Ool001eoo`03001oogoo0>1oo`03001oogoo01Eoo`03001oogoo00Eoo`007Woo00<0 07ooOol0ggoo00<007ooOol05Goo00<007ooOol01Goo000OOol00`00Oomoo`3MOol00`00Oomoo`0F Ool00`00Oomoo`05Ool001moo`03001oogoo0=aoo`03001oogoo01Moo`03001oogoo00Eoo`0087oo 00<007ooOol0fgoo00<007ooOol05goo00<007ooOol01Goo000QOol00`00Oomoo`3IOol00`00Oomo o`0HOol00`00Oomoo`05Ool0029oo`03001oogoo0=Moo`03001oogoo01Uoo`8000Ioo`008goo00<0 07ooOol0eWoo00<007ooOol06Goo00<007ooOol01Goo000SOol00`00Oomoo`3EOol00`00Oomoo`0J Ool00`00Oomoo`05Ool002Aoo`03001oogoo0==oo`03001oogoo01]oo`03001oogoo00Eoo`009Goo 00<007ooOol0dGoo00<007ooOol077oo00<007ooOol01Goo000VOol00`00Oomoo`3?Ool00`00Oomo o`0MOol00`00Oomoo`05Ool002Moo`03001oogoo0goo00<007ooOol0 YGoo0P00"], ImageRangeCache->{{{0, 287}, {287, 0}} -> {-2.05, -1.05001, 0.00731712, \ 0.00731712}}] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, FrontEndVersion->"5.0 for Microsoft Windows", ScreenRectangle->{{29, 1280}, {0, 937}}, AutoGeneratedPackage->None, WindowSize->{860, 782}, WindowMargins->{{Automatic, 26}, {Automatic, 0}}, CellLabelAutoDelete->True, StyleDefinitions -> "Classroom.nb" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[1754, 51, 53, 0, 41, "Subtitle"], Cell[1810, 53, 202, 4, 48, "Text"], Cell[2015, 59, 86, 2, 50, "Input", InitializationCell->True], Cell[CellGroupData[{ Cell[2126, 65, 66, 1, 62, "Section", InitializationCell->True], Cell[CellGroupData[{ Cell[2217, 70, 77, 1, 52, "Subsection", InitializationCell->True], Cell[2297, 73, 1369, 27, 270, "Input", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell[3703, 105, 55, 1, 52, "Subsection", InitializationCell->True], Cell[CellGroupData[{ Cell[3783, 110, 184, 4, 90, "Input", InitializationCell->True], Cell[3970, 116, 66, 1, 49, "Output"], Cell[4039, 119, 76, 1, 49, "Output"], Cell[4118, 122, 73, 1, 49, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[4240, 129, 100, 1, 52, "Subsection", InitializationCell->True], Cell[4343, 132, 619, 13, 150, "Input", InitializationCell->True] }, Open ]], Cell[CellGroupData[{ Cell[4999, 150, 55, 1, 52, "Subsection", InitializationCell->True], Cell[CellGroupData[{ Cell[5079, 155, 84, 2, 50, "Input", InitializationCell->True], Cell[5166, 159, 147, 2, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5350, 166, 82, 2, 50, "Input", InitializationCell->True], Cell[5435, 170, 708, 21, 62, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6180, 196, 85, 2, 50, "Input", InitializationCell->True], Cell[6268, 200, 147, 2, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6452, 207, 82, 2, 50, "Input", InitializationCell->True], Cell[6537, 211, 708, 21, 62, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7282, 237, 85, 2, 50, "Input", InitializationCell->True], Cell[7370, 241, 144, 2, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7551, 248, 82, 2, 50, "Input", InitializationCell->True], Cell[7636, 252, 712, 21, 62, "Output"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[8409, 280, 110, 3, 62, "Section"], Cell[CellGroupData[{ Cell[8544, 287, 51, 0, 52, "Subsection"], Cell[CellGroupData[{ Cell[8620, 291, 489, 9, 130, "Input"], Cell[9112, 302, 42, 1, 49, "Output"], Cell[9157, 305, 235, 3, 65, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9429, 313, 236, 5, 70, "Input"], Cell[9668, 320, 41388, 1289, 356, 14302, 950, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[51071, 1612, 59, 1, 50, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[51167, 1618, 66, 0, 36, "Subsection"], Cell[CellGroupData[{ Cell[51258, 1622, 48, 0, 45, "Subsubsection"], Cell[CellGroupData[{ Cell[51331, 1626, 127, 2, 70, "Input"], Cell[51461, 1630, 61, 1, 62, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[51559, 1636, 59, 1, 50, "Input"], Cell[51621, 1639, 66, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[51724, 1645, 187, 3, 50, "Input"], Cell[51914, 1650, 93, 1, 49, "Output"] }, Open ]], Cell[52022, 1654, 69, 0, 29, "Text"], Cell[CellGroupData[{ Cell[52116, 1658, 222, 4, 70, "Input"], Cell[52341, 1664, 130, 2, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[52508, 1671, 128, 2, 70, "Input"], Cell[52639, 1675, 111, 2, 49, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[52799, 1683, 86, 1, 45, "Subsubsection"], Cell[52888, 1686, 212, 3, 90, "Input"], Cell[53103, 1691, 38, 0, 29, "Text"], Cell[53144, 1693, 105, 2, 50, "Input"], Cell[53252, 1697, 70, 0, 29, "Text"], Cell[CellGroupData[{ Cell[53347, 1701, 508, 10, 110, "Input"], Cell[53858, 1713, 31506, 766, 405, 6815, 457, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[85425, 2486, 67, 0, 36, "Subsection"], Cell[CellGroupData[{ Cell[85517, 2490, 48, 0, 45, "Subsubsection"], Cell[CellGroupData[{ Cell[85590, 2494, 128, 2, 70, "Input"], Cell[85721, 2498, 61, 1, 62, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[85819, 2504, 59, 1, 50, "Input"], Cell[85881, 2507, 66, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[85984, 2513, 179, 3, 50, "Input"], Cell[86166, 2518, 105, 2, 63, "Output"] }, Open ]], Cell[86286, 2523, 69, 0, 29, "Text"], Cell[CellGroupData[{ Cell[86380, 2527, 191, 4, 70, "Input"], Cell[86574, 2533, 108, 2, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[86719, 2540, 128, 2, 70, "Input"], Cell[86850, 2544, 91, 1, 63, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[86990, 2551, 86, 1, 45, "Subsubsection"], Cell[87079, 2554, 212, 3, 90, "Input"], Cell[87294, 2559, 38, 0, 29, "Text"], Cell[87335, 2561, 105, 2, 50, "Input"], Cell[87443, 2565, 70, 0, 29, "Text"], Cell[CellGroupData[{ Cell[87538, 2569, 508, 10, 110, "Input"], Cell[88049, 2581, 31448, 766, 404, 6813, 457, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[119558, 3354, 67, 0, 36, "Subsection"], Cell[CellGroupData[{ Cell[119650, 3358, 48, 0, 45, "Subsubsection"], Cell[CellGroupData[{ Cell[119723, 3362, 127, 2, 70, "Input"], Cell[119853, 3366, 92, 1, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[119982, 3372, 59, 1, 50, "Input"], Cell[120044, 3375, 66, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[120147, 3381, 187, 3, 50, "Input"], Cell[120337, 3386, 135, 2, 49, "Output"] }, Open ]], Cell[120487, 3391, 69, 0, 29, "Text"], Cell[CellGroupData[{ Cell[120581, 3395, 227, 4, 70, "Input"], Cell[120811, 3401, 160, 2, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[121008, 3408, 128, 2, 70, "Input"], Cell[121139, 3412, 141, 2, 49, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[121329, 3420, 86, 1, 45, "Subsubsection"], Cell[121418, 3423, 212, 3, 90, "Input"], Cell[121633, 3428, 38, 0, 29, "Text"], Cell[121674, 3430, 105, 2, 50, "Input"], Cell[121782, 3434, 70, 0, 29, "Text"], Cell[CellGroupData[{ Cell[121877, 3438, 508, 10, 110, "Input"], Cell[122388, 3450, 30619, 745, 404, 6552, 443, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[153068, 4202, 68, 0, 36, "Subsection"], Cell[CellGroupData[{ Cell[153161, 4206, 48, 0, 45, "Subsubsection"], Cell[CellGroupData[{ Cell[153234, 4210, 128, 2, 70, "Input"], Cell[153365, 4214, 93, 1, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[153495, 4220, 59, 1, 50, "Input"], Cell[153557, 4223, 66, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[153660, 4229, 179, 3, 50, "Input"], Cell[153842, 4234, 138, 2, 63, "Output"] }, Open ]], Cell[153995, 4239, 69, 0, 29, "Text"], Cell[CellGroupData[{ Cell[154089, 4243, 191, 4, 70, "Input"], Cell[154283, 4249, 138, 2, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[154458, 4256, 128, 2, 70, "Input"], Cell[154589, 4260, 123, 2, 63, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[154761, 4268, 86, 1, 45, "Subsubsection"], Cell[154850, 4271, 212, 3, 90, "Input"], Cell[155065, 4276, 38, 0, 29, "Text"], Cell[155106, 4278, 105, 2, 50, "Input"], Cell[155214, 4282, 70, 0, 29, "Text"], Cell[CellGroupData[{ Cell[155309, 4286, 508, 10, 110, "Input"], Cell[155820, 4298, 31139, 762, 404, 6813, 457, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell[187032, 5068, 55, 0, 62, "Section"], Cell[CellGroupData[{ Cell[187112, 5072, 66, 0, 52, "Subsection"], Cell[CellGroupData[{ Cell[187203, 5076, 48, 0, 45, "Subsubsection"], Cell[CellGroupData[{ Cell[187276, 5080, 180, 3, 70, "Input"], Cell[187459, 5085, 79, 1, 62, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[187575, 5091, 178, 3, 50, "Input"], Cell[187756, 5096, 67, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[187860, 5102, 67, 1, 50, "Input"], Cell[187930, 5105, 38, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[188005, 5111, 73, 1, 50, "Input"], Cell[188081, 5114, 54, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[188172, 5120, 54, 1, 50, "Input"], Cell[188229, 5123, 44, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[188310, 5129, 76, 1, 50, "Input"], Cell[188389, 5132, 50, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[188476, 5138, 76, 1, 50, "Input"], Cell[188555, 5141, 35, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[188627, 5147, 107, 2, 50, "Input"], Cell[188737, 5151, 54, 1, 49, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[188828, 5157, 110, 2, 52, "Input"], Cell[188941, 5161, 97, 2, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[189075, 5168, 179, 4, 52, "Input"], Cell[189257, 5174, 87, 1, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[189381, 5180, 155, 2, 50, "Input"], Cell[189539, 5184, 20009, 533, 296, 4882, 342, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)