(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' 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[ 56566, 1750]*) (*NotebookOutlinePosition[ 57210, 1772]*) (* CellTagsIndexPosition[ 57166, 1768]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["The Nusselt Problem", "Subsubtitle", FontSize->18, FontVariations->{"CompatibilityType"->0}], StyleBox["\n", FontWeight->"Plain", FontVariations->{"CompatibilityType"->0}], StyleBox["Author's Data", FontSize->14, FontWeight->"Bold"], StyleBox[": ", FontSize->14], StyleBox["Housam BINOUS\nDepartment of Chemical Engineering\nNational \ Institute of Applied Sciences and Technology\nTunis, TUNISIA\nEmail: \ binoushousam@yahoo.com ", FontSize->14, FontWeight->"Plain"] }], "Title", TextAlignment->Center, Background->RGBColor[0.605478, 0.996109, 0.605478]], Cell[TextData[{ StyleBox["Acknowledgement :", FontWeight->"Bold"], "\n\nGoverning equations and analytical solution of the nusselt problem can \ be found in: Rice and Do, Applied Mathematics and Modeling for Chemical \ Engineers, Wiley, 1995.\n" }], "Subsubtitle", Background->RGBColor[0.773449, 0.996109, 0.996109]], Cell[BoxData[ \(Off[General::"\"]\)], "Input"], Cell[CellGroupData[{ Cell[TextData[StyleBox["Solution using NDSolve", FontSize->18]], "Subsubtitle", TextAlignment->Left, Background->RGBColor[1, 1, 0.682353]], Cell[CellGroupData[{ Cell[BoxData[ \(sol = NDSolve[{D[\[Theta][\[Tau], \[Xi]], {\[Tau], 1}] \[Equal] \ \(1\/\[Xi]\) D[\[Xi]\ D[\[Theta][\[Tau], \[Xi]], {\[Xi], 1}], {\[Xi], 1}], \[IndentingNewLine]\[Theta][0, \[Xi]] \[Equal] Exp[\(-1000\)\ \((1 - \[Xi])\)], \ \[Theta][\[Tau], 1] \[Equal] \ 1, \((D[\[Theta][\[Tau], \[Xi]], {\[Xi], 1}] /. \[Xi] \[Rule] 10^\(-8\))\) \[Equal] \ 0}, \[IndentingNewLine]{\[Theta][\[Tau], \[Xi]]}, {\[Tau], 0, 10}, {\[Xi], 10^\(-8\), 1}]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{\(\[Theta][\[Tau], \[Xi]]\), "\[Rule]", RowBox[{ TagBox[\(InterpolatingFunction[{{0.`, 10.`}, {1.`*^-8, 1.`}}, "<>"]\), False, Editable->False], "[", \(\[Tau], \[Xi]\), "]"}]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[BoxData[ \(\(\[Tau] = 1;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(plt1 = Plot[First[\[Theta][\[Tau], \[Xi]] /. sol], {\[Xi], 10^\(-8\), 1}, PlotRange \[Rule] \ All, PlotStyle \[Rule] \ RGBColor[1, 0, 0], DisplayFunction \[Rule] \ Identity]\)], "Input"], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[BoxData[ \(\(\[Tau] = 0.1;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(plt2 = Plot[First[\[Theta][\[Tau], \[Xi]] /. sol], {\[Xi], 10^\(-8\), 1}, PlotRange \[Rule] \ All, PlotStyle \[Rule] \ RGBColor[0, 0, 1], DisplayFunction \[Rule] \ Identity]\)], "Input"], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[BoxData[ \(\(\[Tau] = 0.01;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(plt3 = Plot[First[\[Theta][\[Tau], \[Xi]] /. sol], {\[Xi], 10^\(-8\), 1}, PlotRange \[Rule] \ All, PlotStyle \[Rule] \ RGBColor[0, 1, 0], DisplayFunction \[Rule] \ Identity]\)], "Input"], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Show[plt1, plt2, plt3, DisplayFunction \[Rule] \ $DisplayFunction]\)], "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.952381 0.0147205 0.588598 [ [.21429 .00222 -9 -9 ] [.21429 .00222 9 0 ] [.40476 .00222 -9 -9 ] [.40476 .00222 9 0 ] [.59524 .00222 -9 -9 ] [.59524 .00222 9 0 ] [.78571 .00222 -9 -9 ] [.78571 .00222 9 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [.01131 .13244 -18 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .25016 -18 -4.5 ] [.01131 .25016 0 4.5 ] [.01131 .36788 -18 -4.5 ] [.01131 .36788 0 4.5 ] [.01131 .4856 -18 -4.5 ] [.01131 .4856 0 4.5 ] [.01131 .60332 -6 -4.5 ] [.01131 .60332 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 .21429 .01472 m .21429 .02097 L s [(0.2)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(0.4)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(0.6)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(0.8)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(1)] .97619 .00222 0 1 Mshowa .125 Mabswid .07143 .01472 m .07143 .01847 L s .11905 .01472 m .11905 .01847 L s .16667 .01472 m .16667 .01847 L s .2619 .01472 m .2619 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .45238 .01472 m .45238 .01847 L s .5 .01472 m .5 .01847 L s .54762 .01472 m .54762 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .7381 .01472 m .7381 .01847 L s .83333 .01472 m .83333 .01847 L s .88095 .01472 m .88095 .01847 L s .92857 .01472 m .92857 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s .02381 .13244 m .03006 .13244 L s [(0.2)] .01131 .13244 1 0 Mshowa .02381 .25016 m .03006 .25016 L s [(0.4)] .01131 .25016 1 0 Mshowa .02381 .36788 m .03006 .36788 L s [(0.6)] .01131 .36788 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(0.8)] .01131 .4856 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(1)] .01131 .60332 1 0 Mshowa .125 Mabswid .02381 .04415 m .02756 .04415 L s .02381 .07358 m .02756 .07358 L s .02381 .10301 m .02756 .10301 L s .02381 .16187 m .02756 .16187 L s .02381 .1913 m .02756 .1913 L s .02381 .22073 m .02756 .22073 L s .02381 .27959 m .02756 .27959 L s .02381 .30902 m .02756 .30902 L s .02381 .33845 m .02756 .33845 L s .02381 .39731 m .02756 .39731 L s .02381 .42674 m .02756 .42674 L s .02381 .45617 m .02756 .45617 L s .02381 .51503 m .02756 .51503 L s .02381 .54446 m .02756 .54446 L s .02381 .57389 m .02756 .57389 L s .25 Mabswid .02381 0 m .02381 .61803 L s 1 0 0 r .5 Mabswid .02381 .60044 m .02499 .60044 L .02605 .60044 L .02729 .60044 L .02846 .60044 L .03053 .60044 L .03279 .60045 L .03527 .60045 L .0379 .60045 L .04262 .60045 L .04749 .60045 L .05205 .60045 L .06244 .60045 L .07305 .60046 L .08274 .60046 L .10458 .60047 L .12357 .60049 L .14429 .60051 L .18493 .60056 L .22406 .60063 L .26565 .60071 L .30571 .6008 L .34426 .6009 L .38527 .60101 L .42475 .60114 L .46273 .60126 L .50315 .6014 L .54206 .60155 L .58342 .60171 L .62326 .60187 L .66159 .60203 L .70238 .6022 L .74164 .60236 L .77939 .60252 L .8196 .60269 L .85828 .60285 L .89942 .60302 L .93905 .60318 L .97619 .60332 L s 0 0 1 r .02381 .10394 m .02499 .10394 L .02605 .10394 L .02729 .10395 L .02846 .10395 L .03053 .10397 L .03279 .10399 L .03527 .10402 L .0379 .10406 L .04262 .10415 L .04749 .10428 L .05205 .10442 L .06244 .10484 L .07305 .10541 L .08274 .10604 L .10458 .10787 L .12357 .10994 L .14429 .11269 L .18493 .11961 L .22406 .1282 L .26565 .13946 L .30571 .15239 L .34426 .16676 L .38527 .18416 L .42475 .20295 L .46273 .22287 L .50315 .24601 L .54206 .27009 L .58342 .2975 L .62326 .3255 L .66159 .35375 L .70238 .38499 L .74164 .41597 L .77939 .44633 L .8196 .479 L .85828 .51046 L .89942 .54359 L .93905 .57486 L .97619 .60332 L s 0 1 0 r .02381 .01472 m .06244 .01472 L .08255 .01472 L .08806 .01472 L .09396 .01472 L .09948 .01472 L .10458 .01472 L .10969 .01472 L .11238 .01472 L .11523 .01472 L .11642 .01472 L .1177 .01472 L .11881 .01472 L .12002 .01472 L .12076 .01472 L .12144 .01472 L .12209 .01472 L .12278 .01472 L .1241 .01472 L .1253 .01472 L .12648 .01472 L .12773 .01472 L .12879 .01472 L .12996 .01472 L .13257 .01472 L .13501 .01472 L .14415 .01472 L .14916 .01472 L .15053 .01472 L .15199 .01472 L .15323 .01472 L .15394 .01472 L .1546 .01472 L .15581 .01472 L .15709 .01472 L .15818 .01472 L .15937 .01472 L .16065 .01472 L .16203 .01472 L .16333 .01472 L .1645 .01472 L .16689 .01472 L .16908 .01472 L .1715 .01472 L .17405 .01472 L .17834 .01472 L .18306 .01472 L .18827 .01472 L .19322 .01472 L .20251 .01472 L Mistroke .2125 .01472 L .2234 .01472 L .2344 .01472 L .23703 .01472 L .23835 .01472 L .2398 .01472 L .24219 .01472 L .24476 .01472 L .24957 .01472 L .25405 .01472 L .25931 .01472 L .2641 .01472 L .27375 .01472 L .28423 .01472 L .30328 .01472 L .30442 .01472 L .30564 .01472 L .30679 .01472 L .30784 .01472 L .30902 .01472 L .31013 .01472 L .31133 .01472 L .31263 .01472 L .31537 .01472 L .31789 .01472 L .3227 .01472 L .32807 .01472 L .33375 .01472 L .34411 .01472 L .35315 .01472 L .35813 .01472 L .36058 .01472 L .36166 .01472 L .36285 .01472 L .36414 .01472 L .36485 .01472 L .36551 .01472 L .36841 .01472 L .37349 .01472 L .38343 .01472 L .39325 .01472 L .40359 .01472 L .41242 .01473 L .41743 .01473 L .42209 .01473 L .42713 .01473 L .42978 .01474 L .43123 .01474 L .43259 .01474 L .43505 .01474 L Mistroke .43776 .01475 L .44252 .01475 L .45141 .01477 L .4611 .0148 L .47149 .01484 L .48136 .01488 L .49097 .01494 L .49989 .01501 L .51047 .01512 L .52043 .01526 L .53904 .01562 L .54882 .01589 L .55914 .01624 L .57761 .0171 L .5871 .01769 L .59727 .01847 L .6156 .02034 L .62579 .0217 L .63686 .0235 L .65699 .02772 L .67677 .03355 L .69779 .04206 L .71727 .05253 L .73556 .06517 L .75576 .08293 L .77427 .10323 L .81391 .16157 L .85204 .2389 L .89262 .34383 L .93168 .46205 L .9732 .59407 L .97619 .60332 L Mfstroke 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg 003o00D0oooo00<000000?ooo`3oool02P3oool00`000000oooo0?ooo`0:0?ooo`030000003oool0 oooo00X0oooo00<000000?ooo`3oool02P3oool00`000000oooo0?ooo`0:0?ooo`030000003oool0 oooo00X0oooo00<000000?ooo`3oool02P3oool00`000000oooo0?ooo`050?ooo`005@3oool00`00 0000oooo0?ooo`0a0?ooo`030000003oool0oooo0340oooo00<000000?ooo`3oool0<03oool00`00 0000oooo0?ooo`050?ooo`L00?l09@3oool00`000000oooo0?ooo`0a0?ooo`030000003oool0oooo 00D0oooo000E0?ooo`030000003oool0oooo0:L0oooo1000o`1M0?ooo`005@3oool00`000000oooo 0?ooo`2[0?ooo`<00?l0FP3oool001D0oooo00<000000?ooo`3oool0[P3oool3003o05L0oooo000E 0?ooo`030000003oool0oooo0;40oooo0P00o`1E0?ooo`005@3oool00`000000oooo0?ooo`2c0?oo o`800?l0D`3oool001D0oooo0P00002f0?ooo`800?l0D@3oool001D0oooo00<000000?ooo`3oool0 ]`3oool2003o04l0oooo000E0?ooo`030000003oool0oooo0;T0oooo0P00o`1=0?ooo`005@3oool0 0`000000oooo0?ooo`2k0?ooo`03003o003oool0oooo04X0oooo000E0?ooo`030000003oool0oooo 0;`0oooo0P00o`1:0?ooo`005@3oool00`000000oooo0?ooo`2n0?ooo`03003o003oool0oooo04L0 oooo000E0?ooo`030000003oool0oooo0;l0oooo00<00?l00?ooo`3oool0AP3oool001D0oooo00<0 00000?ooo`3oool0`03oool00`00o`00oooo0?ooo`150?ooo`005@3oool200000<80oooo00<00?l0 0?ooo`3oool0A03oool001D0oooo00<000000?ooo`3oool0`P3oool00`00o`00oooo0?ooo`130?oo o`005@3oool00`000000oooo0?ooo`330?ooo`03003o003oool0oooo0480oooo000E0?ooo`030000 003oool0oooo0<@0oooo00<00?l00?ooo`3oool0@@3oool001D0oooo00<000000?ooo`3oool0a@3o ool00`00o`00oooo0?ooo`100?ooo`005@3oool00`000000oooo0?ooo`360?ooo`03003o003oool0 oooo03l0oooo000E0?ooo`030000003oool0oooo0`3oool001D0oooo00<000000?ooo`3oool05P3oool50000ojl0oooo00<00?l0 0?ooo`3oool0>`3oool001D0oooo00<000000?ooo`3oool06`3oool90000ojL0oooo00<00?l00?oo o`3oool0>P3oool001D0oooo00<000000?ooo`3oool0903oool70000oj40oooo00<00?l00?ooo`3o ool0>@3oool001D0oooo00<000000?ooo`3oool0:`3oool30000oih0oooo00<00?l00?ooo`3oool0 >@3oool00080oooo0P0000040?ooo`8000000`3oool4000000@0oooo00<000000?ooo`3oool0;P3o ool40000oi/0oooo00<00?l00?ooo`3oool0>03oool000050?ooo`000000oooo0?ooo`000000203o ool00`000000oooo0?ooo`050?ooo`030000003oool0oooo0380oooo0`000?nI0?ooo`03003o003o ool0oooo03L0oooo00001@3oool000000?ooo`3oool0000000T0oooo00<000000?ooo`3oool0103o ool3000003D0oooo10000?nE0?ooo`03003o003oool0oooo03L0oooo00001@3oool000000?ooo`3o ool0000000X0oooo00<000000?ooo`3oool00`3oool00`000000oooo0?ooo`0i0?ooo`@0003oTP3o ool00`00o`00oooo0?ooo`0f0?ooo`0000D0oooo0000003oool0oooo000000080?ooo`040000003o ool0oooo000000@0oooo00<000000?ooo`3oool0?@3oool40000ohl0oooo00<00?l00?ooo`3oool0 =@3oool00080oooo0P00000:0?ooo`8000001@3oool00`000000oooo0?ooo`110?ooo`@0003oS03o ool00`00o`00oooo0?ooo`0d0?ooo`005@3oool00`000000oooo0?ooo`150?ooo`@0003oR03oool0 0`00o`00oooo0?ooo`0d0?ooo`005@3oool00`000000oooo0?ooo`190?ooo`<0003oQP3oool00`00 o`00oooo0?ooo`0c0?ooo`005@3oool00`000000oooo0?ooo`1<0?ooo`80003oQ@3oool00`00o`00 oooo0?ooo`0b0?ooo`005@3oool00`000000oooo0?ooo`1>0?ooo`<0003oPP3oool00`00o`00oooo 0?ooo`0b0?ooo`005@3oool200000580oooo0P000?n10?ooo`03003o003oool0oooo0340oooo000E 0?ooo`030000003oool0oooo05<0oooo0`000?mn0?ooo`03003o003oool0oooo0340oooo000E0?oo o`030000003oool0oooo05H0oooo0P000?mm0?ooo`03003o003oool0oooo0300oooo000E0?ooo`03 0000003oool0oooo05P0oooo0P000?mk0?ooo`03003o003oool0oooo0300oooo000E0?ooo`030000 003oool0oooo05X0oooo0P000?mj0?ooo`03003o003oool0oooo02l0oooo000E0?ooo`030000003o ool0oooo05`0oooo0P000?mh0?ooo`03003o003oool0oooo02l0oooo000E0?ooo`030000003oool0 oooo05h0oooo0`000?mf0?ooo`03003o003oool0oooo02h0oooo000E0?ooo`030000003oool0oooo 0640oooo0P000?md0?ooo`03003o003oool0oooo02h0oooo000E0?ooo`800000I03oool20000og<0 oooo00<00?l00?ooo`3oool0;@3oool001D0oooo00<000000?ooo`3oool0I@3oool20000og40oooo 00<00?l00?ooo`3oool0;@3oool001D0oooo00<000000?ooo`3oool0I`3oool20000og00oooo00<0 0?l00?ooo`3oool0;03oool001D0oooo00<000000?ooo`3oool0J@3oool20000ofh0oooo00<00?l0 0?ooo`3oool0;03oool001D0oooo00<000000?ooo`3oool0J`3oool20000ofd0oooo00<00?l00?oo o`3oool0:`3oool001D0oooo00<000000?ooo`3oool0K@3oool20000of/0oooo00<00?l00?ooo`3o ool0:`3oool001D0oooo00<000000?ooo`3oool0K`3oool00`000?l0oooo0?ooo`1Y0?ooo`03003o 003oool0oooo02X0oooo000E0?ooo`030000003oool0oooo0700oooo0P000?mY0?ooo`03003o003o ool0oooo02X0oooo000E0?ooo`800000L`3oool20000ofP0oooo00<00?l00?ooo`3oool0:@3oool0 01D0oooo00<000000?ooo`3oool0M03oool00`000?l0oooo0?ooo`1U0?ooo`03003o003oool0oooo 02T0oooo000E0?ooo`030000003oool0oooo07D0oooo0P000?mV0?ooo`03003o003oool0oooo02P0 oooo000E0?ooo`030000003oool0oooo07L0oooo0P000?mT0?ooo`03003o003oool0oooo02P0oooo 000E0?ooo`030000003oool0oooo07T0oooo0P000?mS0?ooo`03003o003oool0oooo02L0oooo000E 0?ooo`030000003oool0oooo07/0oooo0P000?mQ0?ooo`03003o003oool0oooo02L0oooo00020?oo o`800000103oool2000000D0oooo0`0000030?ooo`030000003oool0oooo07d0oooo0P000?mO0?oo o`03003o003oool0oooo02L0oooo00001@3oool000000?ooo`3oool0000000/0oooo00<000000?oo o`3oool00P3oool00`000000oooo0?ooo`1o0?ooo`030000o`3oool0oooo05d0oooo00<00?l00?oo o`3oool09P3oool000050?ooo`000000oooo0?ooo`000000203oool5000000<0oooo0`0000200?oo o`80003oG@3oool00`00o`00oooo0?ooo`0V0?ooo`0000D0oooo0000003oool0oooo000000080?oo o`040000003oool0oooo000000@0oooo00<000000?ooo`3oool0PP3oool20000oe`0oooo00<00?l0 0?ooo`3oool09@3oool000050?ooo`000000oooo0?ooo`0000002@3oool00`000000oooo00000004 0?ooo`030000003oool0oooo08@0oooo0P000?mJ0?ooo`03003o003oool0oooo02D0oooo00020?oo o`8000002`3oool2000000@0oooo00<000000?ooo`3oool0QP3oool20000oeP0oooo00<00?l00?oo o`3oool09@3oool001D0oooo00<000000?ooo`3oool0R03oool20000oeL0oooo00<00?l00?ooo`3o ool0903oool001D0oooo00<000000?ooo`3oool0RP3oool00`000?l0oooo0?ooo`1D0?ooo`03003o 003oool0oooo02@0oooo000E0?ooo`030000003oool0oooo08/0oooo0P000?mE0?ooo`03003o003o ool0oooo02<0oooo000E0?ooo`030000003oool0oooo08d0oooo00<0003o0?ooo`3oool0DP3oool0 0`00o`00oooo0?ooo`0S0?ooo`005@3oool2000008l0oooo00<0003o0?ooo`3oool0D@3oool00`00 o`00oooo0?ooo`0S0?ooo`005@3oool00`000000oooo0?ooo`2?0?ooo`80003oDP3oool00`00o`00 oooo0?ooo`0R0?ooo`005@3oool00`000000oooo0?ooo`2A0?ooo`030000o`3oool0oooo04l0oooo 00<00?l00?ooo`3oool08P3oool001D0oooo00<000000?ooo`3oool0TP3oool00`000?l0oooo0?oo o`1>0?ooo`03003o003oool0oooo0280oooo000E0?ooo`030000003oool0oooo09<0oooo0P000?m? 0?ooo`03003o003oool0oooo0240oooo000E0?ooo`030000003oool0oooo09D0oooo00<0003o0?oo o`3oool0C03oool00`00o`00oooo0?ooo`0Q0?ooo`005@3oool00`000000oooo0?ooo`2F0?ooo`80 003oC@3oool00`00o`00oooo0?ooo`0P0?ooo`005@3oool00`000000oooo0?ooo`2H0?ooo`030000 o`3oool0oooo04X0oooo00<00?l00?ooo`3oool0803oool001D0oooo0P00002J0?ooo`80003oBP3o ool00`00o`00oooo0?ooo`0P0?ooo`005@3oool00`000000oooo0?ooo`2K0?ooo`80003oB@3oool0 0`00o`00oooo0?ooo`0O0?ooo`005@3oool00`000000oooo0?ooo`2M0?ooo`030000o`3oool0oooo 04H0oooo00<00?l00?ooo`3oool07`3oool001D0oooo00<000000?ooo`3oool0WP3oool20000odL0 oooo00<00?l00?ooo`3oool07P3oool001D0oooo00<000000?ooo`3oool0X03oool00`000?l0oooo 0?ooo`140?ooo`03003o003oool0oooo01h0oooo000E0?ooo`030000003oool0oooo0:40oooo00<0 003o0?ooo`3oool0@`3oool00`00o`00oooo0?ooo`0N0?ooo`005@3oool00`000000oooo0?ooo`2R 0?ooo`80003oA03oool00`00o`00oooo0?ooo`0M0?ooo`005@3oool00`000000oooo0?ooo`2T0?oo o`030000o`3oool0oooo0440oooo00<00?l00?ooo`3oool07@3oool001D0oooo0P00002V0?ooo`03 0000o`3oool0oooo0440oooo00<00?l00?ooo`3oool0703oool001D0oooo00<000000?ooo`3oool0 YP3oool00`000?l0oooo0?ooo`100?ooo`03003o003oool0oooo01`0oooo000E0?ooo`030000003o ool0oooo0:L0oooo0P000?m00?ooo`03003o003oool0oooo01`0oooo000E0?ooo`030000003oool0 oooo0:T0oooo00<0003o0?ooo`3oool0?P3oool00`00o`00oooo0?ooo`0K0?ooo`005@3oool00`00 0000oooo0?ooo`2Z0?ooo`030000o`3oool0oooo03d0oooo00<00?l00?ooo`3oool06`3oool001D0 oooo00<000000?ooo`3oool0Z`3oool00`000?l0oooo0?ooo`0l0?ooo`03003o003oool0oooo01/0 oooo00020?ooo`800000103oool2000000<0oooo0`0000050?ooo`030000003oool0oooo0:`0oooo 0P000?lm0?ooo`03003o003oool0oooo01X0oooo00001@3oool000000?ooo`3oool0000000P0oooo 00@000000?ooo`3oool00000103oool00`000000oooo0?ooo`2^0?ooo`030000o`3oool0oooo03X0 oooo00<00?l00?ooo`3oool06P3oool000050?ooo`000000oooo0?ooo`000000203oool010000000 oooo0?ooo`0000040?ooo`<00000[`3oool00`000?l0oooo0?ooo`0i0?ooo`03003o003oool0oooo 01X0oooo00001@3oool000000?ooo`3oool0000000P0oooo0`0000050?ooo`030000003oool0oooo 0;00oooo00<0003o0?ooo`3oool0>@3oool00`00o`00oooo0?ooo`0I0?ooo`0000D0oooo0000003o ool0oooo000000090?ooo`030000003oool0oooo00@0oooo00<000000?ooo`3oool0/@3oool00`00 0?l0oooo0?ooo`0h0?ooo`03003o003oool0oooo01T0oooo00020?ooo`8000002P3oool3000000@0 oooo00<000000?ooo`3oool0/P3oool20000ocP0oooo00<00?l00?ooo`3oool06@3oool001D0oooo 00<000000?ooo`3oool0]03oool00`000?l0oooo0?ooo`0f0?ooo`03003o003oool0oooo01P0oooo 000E0?ooo`030000003oool0oooo0;D0oooo00<0003o0?ooo`3oool0=@3oool00`00o`00oooo0?oo o`0H0?ooo`005@3oool00`000000oooo0?ooo`2f0?ooo`80003o=@3oool00`00o`00oooo0?ooo`0H 0?ooo`005@3oool00`000000oooo0?ooo`2h0?ooo`030000o`3oool0oooo03<0oooo00<00?l00?oo o`3oool05`3oool001D0oooo0P00002j0?ooo`030000o`3oool0oooo0380oooo00<00?l00?ooo`3o ool05`3oool001D0oooo00<000000?ooo`3oool0^P3oool20000oc80oooo00<00?l00?ooo`3oool0 5`3oool001D0oooo00<000000?ooo`3oool0_03oool00`000?l0oooo0?ooo`0`0?ooo`03003o003o ool0oooo01H0oooo000E0?ooo`030000003oool0oooo0;d0oooo00<0003o0?ooo`3oool0;`3oool0 0`00o`00oooo0?ooo`0F0?ooo`005@3oool00`000000oooo0?ooo`2n0?ooo`80003o;`3oool00`00 o`00oooo0?ooo`0F0?ooo`005@3oool00`000000oooo0?ooo`300?ooo`030000o`3oool0oooo02d0 oooo00<00?l00?ooo`3oool05@3oool001D0oooo00<000000?ooo`3oool0`@3oool00`000?l0oooo 0?ooo`0/0?ooo`03003o003oool0oooo01D0oooo000E0?ooo`030000003oool0oooo0<80oooo0P00 0?l/0?ooo`03003o003oool0oooo01D0oooo000E0?ooo`800000a@3oool00`000?l0oooo0?ooo`0Z 0?ooo`03003o003oool0oooo01@0oooo000E0?ooo`030000003oool0oooo00?ooo`005@3oool00`000000oooo0?ooo`3M0?ooo`80003o603o ool00`00o`00oooo0?ooo`0>0?ooo`005@3oool00`000000oooo0?ooo`3O0?ooo`030000o`3oool0 oooo01H0oooo00<00?l00?ooo`3oool03@3oool001D0oooo00<000000?ooo`3oool0h03oool00`00 0?l0oooo0?ooo`0E0?ooo`03003o003oool0oooo00d0oooo000E0?ooo`800000hP3oool00`000?l0 oooo0?ooo`0D0?ooo`03003o003oool0oooo00d0oooo000E0?ooo`030000003oool0oooo0>80oooo 0P000?lE0?ooo`03003o003oool0oooo00`0oooo000E0?ooo`030000003oool0oooo0>@0oooo00<0 003o0?ooo`3oool04P3oool00`00o`00oooo0?ooo`0<0?ooo`005@3oool00`000000oooo0?ooo`3U 0?ooo`030000o`3oool0oooo0140oooo00<00?l00?ooo`3oool0303oool001D0oooo00<000000?oo o`3oool0iP3oool00`000?l0oooo0?ooo`0A0?ooo`03003o003oool0oooo00/0oooo000E0?ooo`03 0000003oool0oooo0>L0oooo00<0003o0?ooo`3oool0403oool00`00o`00oooo0?ooo`0;0?ooo`00 5@3oool00`000000oooo0?ooo`3X0?ooo`80003o403oool00`00o`00oooo0?ooo`0;0?ooo`005@3o ool00`000000oooo0?ooo`3Z0?ooo`030000o`3oool0oooo00h0oooo00<00?l00?ooo`3oool02P3o ool001D0oooo0P00003/0?ooo`030000o`3oool0oooo00d0oooo00<00?l00?ooo`3oool02P3oool0 01D0oooo00<000000?ooo`3oool0k03oool20000o`d0oooo00<00?l00?ooo`3oool02P3oool001D0 oooo00<000000?ooo`3oool0kP3oool00`000?l0oooo0?ooo`0;0?ooo`03003o003oool0oooo00T0 oooo000E0?ooo`030000003oool0oooo0>l0oooo00<0003o0?ooo`3oool02P3oool00`00o`00oooo 0?ooo`090?ooo`005@3oool00`000000oooo0?ooo`3`0?ooo`80003o2P3oool00`00o`00oooo0?oo o`090?ooo`005@3oool00`000000oooo0?ooo`3b0?ooo`030000o`3oool0oooo00L0oooo00<00?l0 0?ooo`3oool02@3oool001D0oooo00<000000?ooo`3oool0l`3oool00`000?l0oooo0?ooo`070?oo o`03003o003oool0oooo00P0oooo000E0?ooo`030000003oool0oooo0?@0oooo0P000?l70?ooo`03 003o003oool0oooo00P0oooo000E0?ooo`800000m`3oool00`000?l0oooo0?ooo`040?ooo`03003o 003oool0oooo00P0oooo000E0?ooo`030000003oool0oooo0?L0oooo00<0003o0?ooo`3oool0103o ool00`00o`00oooo0?ooo`070?ooo`005@3oool00`000000oooo0?ooo`3h0?ooo`80003o103oool0 0`00o`00oooo0?ooo`070?ooo`005@3oool00`000000oooo0?ooo`3j0?ooo`050000o`3oool0oooo 0?ooo`00o`002@3oool001D0oooo00<000000?ooo`3oool0n`3oool01@000?l0oooo0?ooo`3oool0 0?l000P0oooo000E0?ooo`030000003oool0oooo0?`0oooo00@0003o0?ooo`3oool00?l0203oool0 00d0oooo1@0000030?ooo`030000003oool0oooo0?d0oooo0P000?l00`00o`00oooo0?ooo`060?oo o`003`3oool00`000000oooo0?ooo`030?ooohP0o`00NP3oool00`000?l00?l00?ooo`060?ooo`00 3`3oool00`000000oooo0?ooo`030?ooo`<00000Q@3ooomk0?l00003003o003oool0oooo00D0oooo 000?0?ooo`030000003oool0oooo00<0oooo00<000000?ooo`3oool0o`3oool90?ooo`003`3oool0 0`000000oooo0?ooo`030?ooo`030000003oool0oooo0?l0oooo2@3oool000h0oooo0P0000050?oo o`030000003oool0oooo0?l0oooo2@3oool001D0oooo00<000000?ooo`3oool0o`3oool90?ooo`00 \ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.0829852, -0.07748, \ 0.00386179, 0.00624857}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[TextData[StyleBox["Analytical Solution gives same result as the one \ obtained using NDSolve", FontSize->18]], "Subsubtitle", TextAlignment->Left, Background->RGBColor[1, 1, 0.682353]], Cell[BoxData[ \(\(n = 1;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(While[ n < 10, {\[Lambda]\_n = \(FindRoot[ BesselJ[0, x] \[Equal] \ 0, {x, n\ Pi}]\)[\([1, 2]\)], Print[\[Lambda]\_n], \(n++\)}]\)], "Input"], Cell[BoxData[ \(2.4048255576957724`\)], "Print"], Cell[BoxData[ \(5.5200781102863115`\)], "Print"], Cell[BoxData[ \(8.65372791291101`\)], "Print"], Cell[BoxData[ \(11.79153443901428`\)], "Print"], Cell[BoxData[ \(14.930917708487783`\)], "Print"], Cell[BoxData[ \(18.07106396791092`\)], "Print"], Cell[BoxData[ \(21.21163662987926`\)], "Print"], Cell[BoxData[ \(24.352471530749302`\)], "Print"], Cell[BoxData[ \(27.493479132040253`\)], "Print"] }, Open ]], Cell[BoxData[ \(\[Theta]anal[\[Tau]_, \[Xi]_] := 1 - 2\ \(\[Sum]\+\(n = 1\)\%9\( BesselJ[ 0, \[Lambda]\_n\ \[Xi]]\/\(\(\[Lambda]\_n\) BesselJ[1, \[Lambda]\_n]\)\) Exp[\(-\[Lambda]\_n^2\)\ \[Tau]]\)\)], "Input"], Cell[BoxData[ \(\(\[Tau] = 1;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(plt4 = Plot[First[\[Theta]anal[\[Tau], \[Xi]] /. sol], {\[Xi], 10^\(-8\), 1}, PlotRange \[Rule] \ All, PlotStyle \[Rule] \ RGBColor[1, 0, 0], DisplayFunction \[Rule] \ Identity]\)], "Input"], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[BoxData[ \(\(\[Tau] = 0.1;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(plt5 = Plot[First[\[Theta]anal[\[Tau], \[Xi]] /. sol], {\[Xi], 10^\(-8\), 1}, PlotRange \[Rule] \ All, PlotStyle \[Rule] \ RGBColor[0, 0, 1], DisplayFunction \[Rule] \ Identity]\)], "Input"], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[BoxData[ \(\(\[Tau] = 0.01;\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(plt6 = Plot[First[\[Theta]anal[\[Tau], \[Xi]] /. sol], {\[Xi], 10^\(-8\), 1}, PlotRange \[Rule] \ All, PlotStyle \[Rule] \ RGBColor[0, 1, 0], DisplayFunction \[Rule] \ Identity]\)], "Input"], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Show[plt1, plt2, plt3, plt4, plt5, plt6, DisplayFunction \[Rule] \ $DisplayFunction]\)], "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.952381 0.0147349 0.588584 [ [.21429 .00223 -9 -9 ] [.21429 .00223 9 0 ] [.40476 .00223 -9 -9 ] [.40476 .00223 9 0 ] [.59524 .00223 -9 -9 ] [.59524 .00223 9 0 ] [.78571 .00223 -9 -9 ] [.78571 .00223 9 0 ] [.97619 .00223 -3 -9 ] [.97619 .00223 3 0 ] [.01131 .13245 -18 -4.5 ] [.01131 .13245 0 4.5 ] [.01131 .25017 -18 -4.5 ] [.01131 .25017 0 4.5 ] [.01131 .36789 -18 -4.5 ] [.01131 .36789 0 4.5 ] [.01131 .4856 -18 -4.5 ] [.01131 .4856 0 4.5 ] [.01131 .60332 -6 -4.5 ] [.01131 .60332 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 .21429 .01473 m .21429 .02098 L s [(0.2)] .21429 .00223 0 1 Mshowa .40476 .01473 m .40476 .02098 L s [(0.4)] .40476 .00223 0 1 Mshowa .59524 .01473 m .59524 .02098 L s [(0.6)] .59524 .00223 0 1 Mshowa .78571 .01473 m .78571 .02098 L s [(0.8)] .78571 .00223 0 1 Mshowa .97619 .01473 m .97619 .02098 L s [(1)] .97619 .00223 0 1 Mshowa .125 Mabswid .07143 .01473 m .07143 .01848 L s .11905 .01473 m .11905 .01848 L s .16667 .01473 m .16667 .01848 L s .2619 .01473 m .2619 .01848 L s .30952 .01473 m .30952 .01848 L s .35714 .01473 m .35714 .01848 L s .45238 .01473 m .45238 .01848 L s .5 .01473 m .5 .01848 L s .54762 .01473 m .54762 .01848 L s .64286 .01473 m .64286 .01848 L s .69048 .01473 m .69048 .01848 L s .7381 .01473 m .7381 .01848 L s .83333 .01473 m .83333 .01848 L s .88095 .01473 m .88095 .01848 L s .92857 .01473 m .92857 .01848 L s .25 Mabswid 0 .01473 m 1 .01473 L s .02381 .13245 m .03006 .13245 L s [(0.2)] .01131 .13245 1 0 Mshowa .02381 .25017 m .03006 .25017 L s [(0.4)] .01131 .25017 1 0 Mshowa .02381 .36789 m .03006 .36789 L s [(0.6)] .01131 .36789 1 0 Mshowa .02381 .4856 m .03006 .4856 L s [(0.8)] .01131 .4856 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(1)] .01131 .60332 1 0 Mshowa .125 Mabswid .02381 .04416 m .02756 .04416 L s .02381 .07359 m .02756 .07359 L s .02381 .10302 m .02756 .10302 L s .02381 .16188 m .02756 .16188 L s .02381 .19131 m .02756 .19131 L s .02381 .22074 m .02756 .22074 L s .02381 .2796 m .02756 .2796 L s .02381 .30903 m .02756 .30903 L s .02381 .33846 m .02756 .33846 L s .02381 .39731 m .02756 .39731 L s .02381 .42674 m .02756 .42674 L s .02381 .45617 m .02756 .45617 L s .02381 .51503 m .02756 .51503 L s .02381 .54446 m .02756 .54446 L s .02381 .57389 m .02756 .57389 L s .25 Mabswid .02381 0 m .02381 .61803 L s 1 0 0 r .5 Mabswid .02381 .60044 m .02499 .60044 L .02605 .60044 L .02729 .60044 L .02846 .60044 L .03053 .60045 L .03279 .60045 L .03527 .60045 L .0379 .60045 L .04262 .60045 L .04749 .60045 L .05205 .60045 L .06244 .60045 L .07305 .60046 L .08274 .60046 L .10458 .60047 L .12357 .60049 L .14429 .60051 L .18493 .60056 L .22406 .60063 L .26565 .60071 L .30571 .6008 L .34426 .6009 L .38527 .60101 L .42475 .60114 L .46273 .60126 L .50315 .6014 L .54206 .60155 L .58342 .60171 L .62326 .60187 L .66159 .60203 L .70238 .6022 L .74164 .60236 L .77939 .60252 L .8196 .60269 L .85828 .60285 L .89942 .60302 L .93905 .60318 L .97619 .60332 L s 0 0 1 r .02381 .10395 m .02499 .10395 L .02605 .10395 L .02729 .10396 L .02846 .10396 L .03053 .10398 L .03279 .104 L .03527 .10403 L .0379 .10407 L .04262 .10417 L .04749 .10429 L .05205 .10444 L .06244 .10486 L .07305 .10542 L .08274 .10605 L .10458 .10789 L .12357 .10995 L .14429 .1127 L .18493 .11962 L .22406 .12821 L .26565 .13947 L .30571 .1524 L .34426 .16677 L .38527 .18417 L .42475 .20296 L .46273 .22288 L .50315 .24601 L .54206 .2701 L .58342 .29751 L .62326 .32551 L .66159 .35376 L .70238 .385 L .74164 .41598 L .77939 .44634 L .8196 .479 L .85828 .51046 L .89942 .54359 L .93905 .57487 L .97619 .60332 L s 0 1 0 r .02381 .01474 m .06244 .01474 L .08255 .01473 L .08806 .01473 L .09396 .01473 L .09948 .01473 L .10458 .01473 L .10969 .01473 L .11238 .01473 L .11523 .01473 L .11642 .01473 L .1177 .01473 L .11881 .01473 L .12002 .01473 L .12076 .01473 L .12144 .01473 L .12209 .01473 L .12278 .01473 L .1241 .01473 L .1253 .01473 L .12648 .01473 L .12773 .01473 L .12879 .01473 L .12996 .01473 L .13257 .01473 L .13501 .01473 L .14415 .01473 L .14916 .01473 L .15053 .01473 L .15199 .01473 L .15323 .01473 L .15394 .01473 L .1546 .01473 L .15581 .01473 L .15709 .01473 L .15818 .01473 L .15937 .01473 L .16065 .01473 L .16203 .01473 L .16333 .01473 L .1645 .01473 L .16689 .01473 L .16908 .01473 L .1715 .01473 L .17405 .01473 L .17834 .01473 L .18306 .01473 L .18827 .01473 L .19322 .01473 L .20251 .01473 L Mistroke .2125 .01473 L .2234 .01473 L .2344 .01473 L .23703 .01473 L .23835 .01473 L .2398 .01473 L .24219 .01473 L .24476 .01473 L .24957 .01473 L .25405 .01473 L .25931 .01473 L .2641 .01473 L .27375 .01473 L .28423 .01473 L .30328 .01473 L .30442 .01473 L .30564 .01473 L .30679 .01473 L .30784 .01473 L .30902 .01473 L .31013 .01473 L .31133 .01473 L .31263 .01473 L .31537 .01473 L .31789 .01473 L .3227 .01473 L .32807 .01473 L .33375 .01473 L .34411 .01473 L .35315 .01473 L .35813 .01473 L .36058 .01473 L .36166 .01473 L .36285 .01473 L .36414 .01473 L .36485 .01473 L .36551 .01473 L .36841 .01473 L .37349 .01473 L .38343 .01473 L .39325 .01473 L .40359 .01474 L .41242 .01474 L .41743 .01474 L .42209 .01474 L .42713 .01475 L .42978 .01475 L .43123 .01475 L .43259 .01475 L .43505 .01476 L Mistroke .43776 .01476 L .44252 .01477 L .45141 .01479 L .4611 .01481 L .47149 .01485 L .48136 .01489 L .49097 .01495 L .49989 .01502 L .51047 .01513 L .52043 .01527 L .53904 .01564 L .54882 .0159 L .55914 .01625 L .57761 .01711 L .5871 .0177 L .59727 .01848 L .6156 .02036 L .62579 .02172 L .63686 .02351 L .65699 .02774 L .67677 .03356 L .69779 .04207 L .71727 .05255 L .73556 .06518 L .75576 .08294 L .77427 .10324 L .81391 .16158 L .85204 .23891 L .89262 .34384 L .93168 .46206 L .9732 .59407 L .97619 .60332 L Mfstroke 1 0 0 r .02381 .60042 m .02499 .60042 L .02605 .60042 L .02729 .60042 L .02846 .60042 L .03053 .60042 L .03279 .60042 L .03527 .60042 L .0379 .60042 L .04262 .60042 L .04749 .60042 L .05205 .60042 L .06244 .60042 L .07305 .60043 L .08274 .60043 L .10458 .60045 L .12357 .60046 L .14429 .60048 L .18493 .60053 L .22406 .6006 L .26565 .60068 L .30571 .60077 L .34426 .60087 L .38527 .60099 L .42475 .60111 L .46273 .60124 L .50315 .60139 L .54206 .60153 L .58342 .60169 L .62326 .60186 L .66159 .60201 L .70238 .60219 L .74164 .60235 L .77939 .60251 L .8196 .60268 L .85828 .60285 L .89942 .60302 L .93905 .60317 L .97619 .60332 L s 0 0 1 r .02381 .10399 m .02499 .10399 L .02605 .10399 L .02729 .104 L .02846 .104 L .03053 .10402 L .03279 .10404 L .03527 .10407 L .0379 .10411 L .04262 .1042 L .04749 .10433 L .05205 .10447 L .06244 .10488 L .07305 .10544 L .08274 .10607 L .10458 .1079 L .12357 .10996 L .14429 .11271 L .18493 .11963 L .22406 .12823 L .26565 .13948 L .30571 .1524 L .34426 .16678 L .38527 .18417 L .42475 .20295 L .46273 .22287 L .50315 .24601 L .54206 .27009 L .58342 .29749 L .62326 .3255 L .66159 .35374 L .70238 .38498 L .74164 .41596 L .77939 .44632 L .8196 .47899 L .85828 .51045 L .89942 .54359 L .93905 .57486 L .97619 .60332 L s 0 1 0 r .02381 .01472 m .02499 .01472 L .02605 .01472 L .02729 .01472 L .02846 .01472 L .03053 .01472 L .03279 .01472 L .03527 .01472 L .0379 .01472 L .04262 .01472 L .05205 .01472 L .06244 .01472 L .08274 .01473 L .09394 .01473 L .10458 .01474 L .11464 .01474 L .11995 .01474 L .12562 .01474 L .12831 .01474 L .13085 .01474 L .13333 .01474 L .1356 .01474 L .13777 .01474 L .13887 .01474 L .14005 .01474 L .14124 .01474 L .14235 .01474 L .14367 .01474 L .14487 .01474 L .14605 .01474 L .14715 .01474 L .14837 .01474 L .14967 .01474 L .15234 .01474 L .15478 .01474 L .15925 .01474 L .16395 .01474 L .17387 .01474 L .18471 .01474 L .20596 .01473 L .21601 .01473 L .22527 .01473 L .22995 .01473 L .23234 .01473 L .23491 .01473 L .23625 .01473 L .23748 .01473 L .23873 .01473 L .23989 .01473 L .24123 .01473 L Mistroke .24194 .01473 L .24269 .01473 L .24341 .01473 L .24409 .01473 L .24536 .01473 L .24646 .01473 L .24767 .01473 L .24893 .01473 L .25011 .01473 L .25222 .01473 L .25447 .01473 L .25963 .01473 L .26432 .01473 L .27426 .01473 L .28344 .01473 L .30422 .01474 L .31399 .01474 L .32451 .01474 L .34353 .01474 L .3616 .01474 L .37115 .01475 L .37649 .01475 L .38134 .01475 L .3867 .01475 L .39175 .01475 L .39639 .01475 L .40142 .01475 L .40676 .01475 L .41241 .01476 L .42253 .01477 L .43199 .01478 L .44214 .0148 L .45173 .01482 L .46047 .01484 L .47037 .01488 L .47953 .01493 L .48959 .01499 L .50028 .01508 L .50964 .01518 L .51973 .01531 L .538 .01566 L .54756 .0159 L .55802 .01624 L .56906 .01669 L .57911 .01721 L .5987 .01859 L .60946 .01961 L .61964 .02079 L .63963 .02388 L .65807 .02789 L Mistroke .67681 .03343 L .69744 .04175 L .71687 .05222 L .73775 .06699 L .7589 .08638 L .77899 .10957 L .81721 .16839 L .85787 .25389 L .89702 .35728 L .93862 .48361 L .97619 .60332 L Mfstroke 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg 003o00D0oooo00<000000?ooo`3oool02P3oool00`000000oooo0?ooo`0:0?ooo`030000003oool0 oooo00X0oooo00<000000?ooo`3oool02P3oool00`000000oooo0?ooo`0:0?ooo`030000003oool0 oooo00X0oooo00<000000?ooo`3oool02P3oool00`000000oooo0?ooo`050?ooo`005@3oool00`00 0000oooo0?ooo`0a0?ooo`030000003oool0oooo0340oooo00<000000?ooo`3oool0<03oool00`00 0000oooo0?ooo`030?ooo`T00?l09@3oool00`000000oooo0?ooo`0a0?ooo`030000003oool0oooo 00D0oooo000E0?ooo`030000003oool0oooo0:H0oooo1@00o`1M0?ooo`005@3oool00`000000oooo 0?ooo`2X0?ooo`H00?l0FP3oool001D0oooo00<000000?ooo`3oool0[03oool5003o05L0oooo000E 0?ooo`030000003oool0oooo0;40oooo0P00o`1E0?ooo`005@3oool00`000000oooo0?ooo`2c0?oo o`800?l0D`3oool001D0oooo0P00002f0?ooo`800?l0D@3oool001D0oooo00<000000?ooo`3oool0 ]`3oool2003o04l0oooo000E0?ooo`030000003oool0oooo0;T0oooo0P00o`1=0?ooo`005@3oool0 0`000000oooo0?ooo`2j0?ooo`800?l0C03oool001D0oooo00<000000?ooo`3oool0_03oool2003o 04X0oooo000E0?ooo`030000003oool0oooo0;d0oooo0P00o`190?ooo`005@3oool00`000000oooo 0?ooo`2o0?ooo`03003o003oool0oooo04H0oooo000E0?ooo`030000003oool0oooo0<00oooo00<0 0?l00?ooo`3oool0A@3oool001D0oooo0P0000320?ooo`03003o003oool0oooo04@0oooo000E0?oo o`030000003oool0oooo0<80oooo0P00o`140?ooo`005@3oool00`000000oooo0?ooo`330?ooo`80 0?l0@`3oool001D0oooo00<000000?ooo`3oool0a03oool2003o0480oooo000E0?ooo`030000003o ool0oooo00?ooo`<0003oPP3oool00`00o`00oooo 0?ooo`0b0?ooo`005@3oool200000580oooo0P000?n10?ooo`03003o003oool0oooo0340oooo000E 0?ooo`030000003oool0oooo05<0oooo0`000?mn0?ooo`03003o003oool0oooo0340oooo000E0?oo o`030000003oool0oooo05H0oooo0P000?mm0?ooo`03003o003oool0oooo0300oooo000E0?ooo`03 0000003oool0oooo05P0oooo0P000?mk0?ooo`03003o003oool0oooo0300oooo000E0?ooo`030000 003oool0oooo05X0oooo0P000?mj0?ooo`03003o003oool0oooo02l0oooo000E0?ooo`030000003o ool0oooo05`0oooo0P000?mh0?ooo`03003o003oool0oooo02l0oooo000E0?ooo`030000003oool0 oooo05h0oooo0`000?mf0?ooo`03003o003oool0oooo02h0oooo000E0?ooo`030000003oool0oooo 0640oooo0P000?md0?ooo`03003o003oool0oooo02h0oooo000E0?ooo`800000I03oool20000og<0 oooo00<00?l00?ooo`3oool0;@3oool001D0oooo00<000000?ooo`3oool0I@3oool20000og40oooo 00<00?l00?ooo`3oool0;@3oool001D0oooo00<000000?ooo`3oool0I`3oool20000og00oooo00<0 0?l00?ooo`3oool0;03oool001D0oooo00<000000?ooo`3oool0J@3oool20000ofh0oooo00<00?l0 0?ooo`3oool0;03oool001D0oooo00<000000?ooo`3oool0J`3oool20000ofd0oooo00<00?l00?oo o`3oool0:`3oool001D0oooo00<000000?ooo`3oool0K@3oool20000of/0oooo00<00?l00?ooo`3o ool0:`3oool001D0oooo00<000000?ooo`3oool0K`3oool00`000?l0oooo0?ooo`1Y0?ooo`03003o 003oool0oooo02X0oooo000E0?ooo`030000003oool0oooo0700oooo0P000?mY0?ooo`03003o003o ool0oooo02X0oooo000E0?ooo`800000L`3oool20000ofP0oooo00<00?l00?ooo`3oool0:@3oool0 01D0oooo00<000000?ooo`3oool0M03oool00`000?l0oooo0?ooo`1U0?ooo`03003o003oool0oooo 02T0oooo000E0?ooo`030000003oool0oooo07D0oooo0P000?mV0?ooo`03003o003oool0oooo02P0 oooo000E0?ooo`030000003oool0oooo07L0oooo0P000?mT0?ooo`03003o003oool0oooo02P0oooo 000E0?ooo`030000003oool0oooo07T0oooo0P000?mS0?ooo`03003o003oool0oooo02L0oooo000E 0?ooo`030000003oool0oooo07/0oooo0P000?mQ0?ooo`03003o003oool0oooo02L0oooo00020?oo o`800000103oool2000000D0oooo0`0000030?ooo`030000003oool0oooo07d0oooo0P000?mO0?oo o`800?l0:03oool000050?ooo`000000oooo0?ooo`0000002`3oool00`000000oooo0?ooo`020?oo o`030000003oool0oooo07l0oooo00<0003o0?ooo`3oool0G@3oool00`00o`00oooo0?ooo`0V0?oo o`0000D0oooo0000003oool0oooo000000080?ooo`D000000`3oool300000800oooo0P000?mM0?oo o`800?l09`3oool000050?ooo`000000oooo0?ooo`000000203oool010000000oooo0?ooo`000004 0?ooo`030000003oool0oooo0880oooo0P000?mL0?ooo`03003o003oool0oooo02D0oooo00001@3o ool000000?ooo`3oool0000000T0oooo00<000000?ooo`000000103oool00`000000oooo0?ooo`24 0?ooo`80003oFP3oool00`00o`00oooo0?ooo`0U0?ooo`000P3oool2000000/0oooo0P0000040?oo o`030000003oool0oooo08H0oooo0P000?mH0?ooo`800?l09P3oool001D0oooo00<000000?ooo`3o ool0R03oool20000oeL0oooo00<00?l00?ooo`3oool0903oool001D0oooo00<000000?ooo`3oool0 RP3oool00`000?l0oooo0?ooo`1D0?ooo`03003o003oool0oooo02@0oooo000E0?ooo`030000003o ool0oooo08/0oooo0P000?mE0?ooo`03003o003oool0oooo02<0oooo000E0?ooo`030000003oool0 oooo08d0oooo00<0003o0?ooo`3oool0DP3oool00`00o`00oooo0?ooo`0S0?ooo`005@3oool20000 08l0oooo00<0003o0?ooo`3oool0D@3oool00`00o`00oooo0?ooo`0S0?ooo`005@3oool00`000000 oooo0?ooo`2?0?ooo`80003oDP3oool00`00o`00oooo0?ooo`0R0?ooo`005@3oool00`000000oooo 0?ooo`2A0?ooo`030000o`3oool0oooo04l0oooo00<00?l00?ooo`3oool08P3oool001D0oooo00<0 00000?ooo`3oool0TP3oool00`000?l0oooo0?ooo`1>0?ooo`800?l08`3oool001D0oooo00<00000 0?ooo`3oool0T`3oool20000odl0oooo00<00?l00?ooo`3oool08@3oool001D0oooo00<000000?oo o`3oool0U@3oool00`000?l0oooo0?ooo`1<0?ooo`03003o003oool0oooo0240oooo000E0?ooo`03 0000003oool0oooo09H0oooo0P000?m=0?ooo`03003o003oool0oooo0200oooo000E0?ooo`030000 003oool0oooo09P0oooo00<0003o0?ooo`3oool0BP3oool00`00o`00oooo0?ooo`0P0?ooo`005@3o ool2000009X0oooo0P000?m:0?ooo`03003o003oool0oooo0200oooo000E0?ooo`030000003oool0 oooo09/0oooo0P000?m90?ooo`03003o003oool0oooo01l0oooo000E0?ooo`030000003oool0oooo 09d0oooo00<0003o0?ooo`3oool0AP3oool00`00o`00oooo0?ooo`0O0?ooo`005@3oool00`000000 oooo0?ooo`2N0?ooo`80003oAP3oool2003o0200oooo000E0?ooo`030000003oool0oooo0:00oooo 00<0003o0?ooo`3oool0A03oool00`00o`00oooo0?ooo`0N0?ooo`005@3oool00`000000oooo0?oo o`2Q0?ooo`030000o`3oool0oooo04<0oooo00<00?l00?ooo`3oool07P3oool001D0oooo00<00000 0?ooo`3oool0XP3oool20000od<0oooo0P00o`0O0?ooo`005@3oool00`000000oooo0?ooo`2T0?oo o`030000o`3oool0oooo0440oooo00<00?l00?ooo`3oool07@3oool001D0oooo0P00002V0?ooo`03 0000o`3oool0oooo0400oooo0P00o`0N0?ooo`005@3oool00`000000oooo0?ooo`2V0?ooo`030000 o`3oool0oooo0400oooo00<00?l00?ooo`3oool0703oool001D0oooo00<000000?ooo`3oool0Y`3o ool20000od00oooo00<00?l00?ooo`3oool0703oool001D0oooo00<000000?ooo`3oool0Z@3oool0 0`000?l0oooo0?ooo`0m0?ooo`800?l07@3oool001D0oooo00<000000?ooo`3oool0ZP3oool00`00 0?l0oooo0?ooo`0m0?ooo`03003o003oool0oooo01/0oooo000E0?ooo`030000003oool0oooo0:/0 oooo00<0003o0?ooo`3oool0?03oool00`00o`00oooo0?ooo`0K0?ooo`000P3oool2000000@0oooo 0P0000030?ooo`<000001@3oool00`000000oooo0?ooo`2/0?ooo`80003o?03oool2003o01`0oooo 00001@3oool000000?ooo`3oool0000000P0oooo00@000000?ooo`3oool00000103oool00`000000 oooo0?ooo`2^0?ooo`030000o`3oool0oooo03X0oooo00<00?l00?ooo`3oool06P3oool000050?oo o`000000oooo0?ooo`000000203oool010000000oooo0?ooo`0000040?ooo`<00000[`3oool00`00 0?l0oooo0?ooo`0i0?ooo`03003o003oool0oooo01X0oooo00001@3oool000000?ooo`3oool00000 00P0oooo0`0000050?ooo`030000003oool0oooo0;00oooo00<0003o0?ooo`3oool0>03oool2003o 01/0oooo00001@3oool000000?ooo`3oool0000000T0oooo00<000000?ooo`3oool0103oool00`00 0000oooo0?ooo`2a0?ooo`030000o`3oool0oooo03P0oooo00<00?l00?ooo`3oool06@3oool00080 oooo0P00000:0?ooo`<00000103oool00`000000oooo0?ooo`2b0?ooo`80003o>03oool00`00o`00 oooo0?ooo`0I0?ooo`005@3oool00`000000oooo0?ooo`2d0?ooo`030000o`3oool0oooo03D0oooo 0P00o`0J0?ooo`005@3oool00`000000oooo0?ooo`2e0?ooo`030000o`3oool0oooo03D0oooo00<0 0?l00?ooo`3oool0603oool001D0oooo00<000000?ooo`3oool0]P3oool20000ocD0oooo00<00?l0 0?ooo`3oool0603oool001D0oooo00<000000?ooo`3oool0^03oool00`000?l0oooo0?ooo`0c0?oo o`03003o003oool0oooo01L0oooo000E0?ooo`800000^P3oool00`000?l0oooo0?ooo`0b0?ooo`03 003o003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo0;X0oooo0P000?lb0?ooo`03003o 003oool0oooo01L0oooo000E0?ooo`030000003oool0oooo0;`0oooo00<0003o0?ooo`3oool0<03o ool00`00o`00oooo0?ooo`0F0?ooo`005@3oool00`000000oooo0?ooo`2m0?ooo`030000o`3oool0 oooo02l0oooo00<00?l00?ooo`3oool05P3oool001D0oooo00<000000?ooo`3oool0_P3oool20000 obl0oooo00<00?l00?ooo`3oool05P3oool001D0oooo00<000000?ooo`3oool0`03oool00`000?l0 oooo0?ooo`0]0?ooo`03003o003oool0oooo01D0oooo000E0?ooo`030000003oool0oooo0<40oooo 00<0003o0?ooo`3oool0;03oool00`00o`00oooo0?ooo`0E0?ooo`005@3oool00`000000oooo0?oo o`320?ooo`80003o;03oool00`00o`00oooo0?ooo`0E0?ooo`005@3oool20000000oooo00<0003o0?ooo`3oool05@3oool00`00o`00oooo0?ooo`0=0?ooo`005@3oool200000>80 oooo00<0003o0?ooo`3oool0503oool2003o00h0oooo000E0?ooo`030000003oool0oooo0>80oooo 0P000?lE0?ooo`03003o003oool0oooo00`0oooo000E0?ooo`030000003oool0oooo0>@0oooo00<0 003o0?ooo`3oool04P3oool00`00o`00oooo0?ooo`0<0?ooo`005@3oool00`000000oooo0?ooo`3U 0?ooo`030000o`3oool0oooo0140oooo0P00o`0=0?ooo`005@3oool00`000000oooo0?ooo`3V0?oo o`030000o`3oool0oooo0140oooo00<00?l00?ooo`3oool02`3oool001D0oooo00<000000?ooo`3o ool0i`3oool00`000?l0oooo0?ooo`0@0?ooo`03003o003oool0oooo00/0oooo000E0?ooo`030000 003oool0oooo0>P0oooo0P000?l@0?ooo`800?l0303oool001D0oooo00<000000?ooo`3oool0jP3o ool00`000?l0oooo0?ooo`0>0?ooo`03003o003oool0oooo00X0oooo000E0?ooo`800000k03oool0 0`000?l0oooo0?ooo`0=0?ooo`03003o003oool0oooo00X0oooo000E0?ooo`030000003oool0oooo 0>`0oooo0P000?l=0?ooo`03003o003oool0oooo00X0oooo000E0?ooo`030000003oool0oooo0>h0 oooo00<0003o0?ooo`3oool02`3oool00`00o`00oooo0?ooo`090?ooo`005@3oool00`000000oooo 0?ooo`3_0?ooo`030000o`3oool0oooo00X0oooo00<00?l00?ooo`3oool02@3oool001D0oooo00<0 00000?ooo`3oool0l03oool20000o`X0oooo00<00?l00?ooo`3oool02@3oool001D0oooo00<00000 0?ooo`3oool0lP3oool00`000?l0oooo0?ooo`070?ooo`800?l02P3oool001D0oooo00<000000?oo o`3oool0l`3oool00`000?l0oooo0?ooo`070?ooo`03003o003oool0oooo00P0oooo000E0?ooo`03 0000003oool0oooo0?@0oooo0P000?l70?ooo`03003o003oool0oooo00P0oooo000E0?ooo`800000 m`3oool00`000?l0oooo0?ooo`040?ooo`800?l02@3oool001D0oooo00<000000?ooo`3oool0m`3o ool00`000?l0oooo0?ooo`040?ooo`03003o003oool0oooo00L0oooo000E0?ooo`030000003oool0 oooo0?P0oooo0P000?l40?ooo`03003o003oool0oooo00L0oooo000E0?ooo`030000003oool0oooo 0?X0oooo00D0003o0?ooo`3oool0oooo003o00090?ooo`005@3oool00`000000oooo0?ooo`3k0?oo o`050000o`3oool0oooo0?ooo`00o`00203oool001D0oooo00<000000?ooo`3oool0o03oool01000 0?l0oooo0?ooo`00o`080?ooo`003@3oool5000000<0oooo00<000000?ooo`3oool0o@3oool20000 o`03003o003oool0oooo00H0oooo000?0?ooo`030000003oool0oooo00<0ooooR03o001j0?ooo`03 0000o`00o`00oooo00H0oooo000?0?ooo`030000003oool0oooo00<0oooo0`0000250?ooog/0o`00 00<00?l00?ooo`3oool01@3oool000l0oooo00<000000?ooo`3oool00`3oool00`000000oooo0?oo o`3o0?ooo`T0oooo000?0?ooo`030000003oool0oooo00<0oooo00<000000?ooo`3oool0o`3oool9 0?ooo`003P3oool2000000D0oooo00<000000?ooo`3oool0o`3oool90?ooo`005@3oool00`000000 oooo0?ooo`3o0?ooo`T0oooo0000\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.0829754, -0.0774884, \ 0.00386173, 0.00624861}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]] }, Open ]] }, FrontEndVersion->"5.2 for Microsoft Windows", ScreenRectangle->{{0, 800}, {0, 517}}, WindowSize->{495, 359}, WindowMargins->{{17, Automatic}, {Automatic, 39}} ] (******************************************************************* 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[CellGroupData[{ Cell[1776, 53, 631, 19, 229, "Title"], Cell[2410, 74, 326, 7, 167, "Subsubtitle"], Cell[2739, 83, 59, 1, 30, "Input"], Cell[CellGroupData[{ Cell[2823, 88, 144, 3, 49, "Subsubtitle"], Cell[CellGroupData[{ Cell[2992, 95, 589, 10, 170, "Input"], Cell[3584, 107, 358, 10, 48, "Output"] }, Open ]], Cell[3957, 120, 48, 1, 30, "Input"], Cell[CellGroupData[{ Cell[4030, 125, 231, 4, 70, "Input"], Cell[4264, 131, 130, 3, 29, "Output"] }, Open ]], Cell[4409, 137, 50, 1, 30, "Input"], Cell[CellGroupData[{ Cell[4484, 142, 231, 4, 70, "Input"], Cell[4718, 148, 130, 3, 29, "Output"] }, Open ]], Cell[4863, 154, 51, 1, 30, "Input"], Cell[CellGroupData[{ Cell[4939, 159, 231, 4, 70, "Input"], Cell[5173, 165, 130, 3, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5340, 173, 106, 2, 50, "Input"], Cell[5449, 177, 22700, 626, 186, 6254, 418, "GraphicsData", "PostScript", \ "Graphics"], Cell[28152, 805, 130, 3, 29, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[28331, 814, 195, 4, 74, "Subsubtitle"], Cell[28529, 820, 43, 1, 30, "Input"], Cell[CellGroupData[{ Cell[28597, 825, 190, 4, 70, "Input"], Cell[28790, 831, 52, 1, 25, "Print"], Cell[28845, 834, 52, 1, 25, "Print"], Cell[28900, 837, 50, 1, 25, "Print"], Cell[28953, 840, 51, 1, 25, "Print"], Cell[29007, 843, 52, 1, 25, "Print"], Cell[29062, 846, 51, 1, 25, "Print"], Cell[29116, 849, 51, 1, 25, "Print"], Cell[29170, 852, 52, 1, 25, "Print"], Cell[29225, 855, 52, 1, 25, "Print"] }, Open ]], Cell[29292, 859, 277, 5, 51, "Input"], Cell[29572, 866, 48, 1, 30, "Input"], Cell[CellGroupData[{ Cell[29645, 871, 235, 4, 70, "Input"], Cell[29883, 877, 130, 3, 29, "Output"] }, Open ]], Cell[30028, 883, 50, 1, 30, "Input"], Cell[CellGroupData[{ Cell[30103, 888, 235, 4, 70, "Input"], Cell[30341, 894, 130, 3, 29, "Output"] }, Open ]], Cell[30486, 900, 51, 1, 30, "Input"], Cell[CellGroupData[{ Cell[30562, 905, 235, 4, 70, "Input"], Cell[30800, 911, 130, 3, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30967, 919, 124, 2, 50, "Input"], Cell[31094, 923, 25299, 817, 186, 9309, 615, "GraphicsData", "PostScript", \ "Graphics"], Cell[56396, 1742, 130, 3, 29, "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)