(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 4.0, MathReader 4.0, or any compatible application. 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[ 44296, 1108]*) (*NotebookOutlinePosition[ 45443, 1147]*) (* CellTagsIndexPosition[ 45346, 1140]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["The package OscillatorPlot - 2", "Title"], Cell["\<\ Bernd Thaller Institute for Mathematics University of Graz\ \>", "Subsubsection"], Cell[CellGroupData[{ Cell["Defining a package context", "Section"], Cell[TextData[{ "As a next step, we put the commands exported by the package in a \"package \ context\". We name this package context \"Graphics`OscillatorPlot`\". Using \ the BeginPackage - EndPackage commands ensures that the package context is \ put on the context search path after loading the package. If you place the \ file OscillatorPlot2.m in the Graphics-subdirectory of AddOns/Applications, \ then it is automatically on ", StyleBox["Mathematica", FontSlant->"Italic"], "'s search path. You can then load it with the command \ <{"CompatibilityType"->0}], "OscillatorFunction[n_,x_] :=\n \ HermiteH[n,x]*Exp[-x^2/2]/Sqrt[2^n*n!]/Pi^(1/4)\n", StyleBox["\n", FontVariations->{"CompatibilityType"->0}], "OscillatorPlot[f_,{x_Symbol,x1_:-Infinity,x2_:Infinity},\n \ n1_:0,n2_:10]:=\n Module[{coeff,coloredLine,lineTable,up,upper,lower,func},\ \n func=Function[{x},f];\n Do[coeff[n]=\n \ NIntegrate[OscillatorFunction[n,x]*func[x],{x,x1,x2}],{n,n1,n2}];\n \ coloredLine[n_]:={Thickness[0.01],Hue[Arg[coeff[n]]/(2*Pi)],\n \ Line[{{n+1/2,Abs[coeff[n]]},{n+1/2,0.}}]};\n \ lineTable=Graphics[Table[coloredLine[n],{n,n1,n2}]];\n \ up=Max[Table[Abs[coeff[n]],{n,n1,n2}]];\n upper=up+up/10; \ lower=-(up/10);\n Show[lineTable, Frame\[Rule]True, \ PlotRange\[Rule]{{n1,n2+1}, {lower,upper}}]\n ]\n \nEnd[]\n\n\ EndPackage[]" }], "Program"], Cell[TextData[{ "\n", StyleBox["Start a new kernel session here!", FontWeight->"Bold"], "\n\nNow the package can be loaded with" }], "Text"], Cell[BoxData[ \(<< Graphics`OscillatorPlot2`\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \($ContextPath\)], "Input"], Cell[BoxData[ \({"Graphics`OscillatorPlot`", "Global`", "System`"}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(?OscillatorPlot\)\)], "Input"], Cell[BoxData[ \("OscillatorPlot[f[x],{x,a,b},n1,n2] plots the energy representation of \ a given function in the basis of harmonic oscillator eigenfunctions. \ Numerical region for the determination of the expansion coefficients is the \ interval (a,b). The graph shows the expansion coefficients c_n for n between \ n1 and n2. Default values are: a=-Infinity, b=Infinity, n1=0, n2=10."\)], \ "Print"] }, Open ]], Cell["\<\ Presently, the symbol OscillatorFunction is in the private context \ of the package. This context is not searched during the main session. \ Therefore this symbol would not conflict with any user-defined \ OscillatorFunction. The symbol OscillatorPlot is in the context \ Graphics`OscillatorPlot` which is on the context path. Hence it was necessary \ to start a new session before loading the package. Otherwise the earlier \ defined symbol OscillatorPlot (in the global context) would have shadowed the \ symbol in the package.\ \>", "Text"], Cell[TextData[{ "Note: You do not need to store the package files in the ", StyleBox["Mathematica", FontSlant->"Italic"], " folder. During development, it is better to place the files somewhere. \ Instead of worrying whether \"somewhere\" is on the search path, we keep the \ file in an open notebook window and load the package just by activating the \ window and evaluating the notebook." }], "Text"] }, Closed]], Cell[CellGroupData[{ Cell["Default values", "Section"], Cell["\<\ The package OscillatorPlot2 already defines defaults in order to \ facilitate the input. Since the defaults are simple numbers, we can put it \ just into the input-pattern\ \>", "Text"], Cell[BoxData[ \(n1_: 0\)], "Input"], Cell[TextData[{ "This pattern matches any expression, later referred to as \"n1\". If no \ input is given for n1, then the value 0 is automatically assumed. ", StyleBox["Mathematica", FontSlant->"Italic"], " gives you several possibilities to define default values for arguments of \ functions. Let us quote here the text from the ", StyleBox["Mathematica", FontSlant->"Italic"], " book:" }], "Text"], Cell[TextData[{ "When the pattern objects ", StyleBox["x", "TI"], StyleBox["_:", "MR"], StyleBox["v", "TI"], " and ", StyleBox["x", "TI"], StyleBox["_.", "MR"], " appear as arguments of functions, they represent arguments which may be \ omitted. When the argument corresponding to ", StyleBox["x", "TI"], StyleBox["_:", "MR"], StyleBox["v", "TI"], " is omitted, ", StyleBox["x", "TI"], " is taken to have value ", StyleBox["v", "TI"], ". When the argument corresponding to ", StyleBox["x", "TI"], StyleBox["_.", "MR"], " is omitted, ", StyleBox["x", "TI"], " is taken to have a ", StyleBox["default value", "TI"], " that is associated with the function in which it appears. You can specify \ this default value by making assignments for ", StyleBox["Default[", "MR"], StyleBox["f", "TI"], StyleBox["]", "MR"], " and so on." }], "Text"], Cell[BoxData[ FormBox[ StyleBox[ FrameBox[ StyleBox[GridBox[{ {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Default[\\\\\\\"\\\", \ \\\"MR\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"f\\\\\\\"\\\", \ \\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"]\\\\\\\"\\\", \\\"MR\\\"]\\) \ \"\>", "\<\"default value for \\!\\(\\*StyleBox[\\\"\\\\\\\"x\\\\\\\"\\\", \\\ \"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"_.\\\\\\\"\\\", \\\"MR\\\"]\\) when \ it appears as any argument of the function \ \\!\\(\\*StyleBox[\\\"\\\\\\\"f\\\\\\\"\\\", \\\"TI\\\"]\\) \"\>"}, {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Default[\\\\\\\"\\\", \ \\\"MR\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"f\\\\\\\"\\\", \ \\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\ \\(\\*StyleBox[\\\"\\\\\\\"n\\\\\\\"\\\", \ \\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"]\\\\\\\"\\\", \\\"MR\\\"]\\) \ \"\>", RowBox[{"\<\"default value for \\!\\(\\*StyleBox[\\\"\\\\\\\ \"x\\\\\\\"\\\", \\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"_.\\\\\\\"\\\", \ \\\"MR\\\"]\\) when it appears as the \\!\\(\\*StyleBox[\\\"\\\\\\\"n\\\\\\\"\ \\\", \\\"TI\\\"]\\)\"\>", StyleBox[\(\[Null]\^"th"\), "InlineFormula"], "\<\" argument (negative \ \\!\\(\\*StyleBox[\\\"\\\\\\\"n\\\\\\\"\\\", \\\"TI\\\"]\\) count from the \ end) \"\>"}]}, {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Default[\\\\\\\"\\\", \ \\\"MR\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"f\\\\\\\"\\\", \ \\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\ \\(\\*StyleBox[\\\"\\\\\\\"n\\\\\\\"\\\", \ \\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\ \\(\\*StyleBox[\\\"\\\\\\\"tot\\\\\\\"\\\", \ \\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"]\\\\\\\"\\\", \\\"MR\\\"]\\) \ \"\>", RowBox[{"\<\"default value for the \\!\\(\\*StyleBox[\\\"\\\ \\\\\"n\\\\\\\"\\\", \\\"TI\\\"]\\)\"\>", StyleBox[\(\[Null]\^"th"\), "InlineFormula"], "\<\" argument when there are a total \ of \\!\\(\\*StyleBox[\\\"\\\\\\\"tot\\\\\\\"\\\", \\\"TI\\\"]\\) arguments \"\ \>"}]} }, ColumnAlignments->{Right, Left}], GridBoxOptions->{RowLines->False}]], "2ColumnBox"], TraditionalForm]], "Text", CellTags->{"SA.5.1", "1.5"}], Cell["For example:", "Text"], Cell[BoxData[ \(\(Default[f]\ = \ 2;\)\)], "Input"], Cell[BoxData[ \(f[x_. ] := 3 x^2 + 1\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(f[]\)], "Input"], Cell[BoxData[ \(13\)], "Output"] }, Open ]], Cell[BoxData[ \(g[x_: 2] := 3\ x^2\ + \ 1\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(g[]\)], "Input"], Cell[BoxData[ \(13\)], "Output"] }, Open ]], Cell["\<\ In the package OscillatorPlot, we use default values for the \ numerical region and for the range of indices. Thus you could simply type \ OscillatorPlot[Exp[-x^2],{x}] instead of \ OscillatorPlot[[Exp[-x^2],{x,-Infinity, Infinity}, 0,10]. If you specify \ other values, you have to type the input using the longer form.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(OscillatorPlot[Exp[\(-x^2\) - 2 I\ x], {x, \(-4\), 4}, 8, 12]\)], "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.6 0.2 0.0515028 4.5485 [ [.2 -0.0125 -3 -9 ] [.2 -0.0125 3 0 ] [.4 -0.0125 -6 -9 ] [.4 -0.0125 6 0 ] [.6 -0.0125 -6 -9 ] [.6 -0.0125 6 0 ] [.8 -0.0125 -6 -9 ] [.8 -0.0125 6 0 ] [1 -0.0125 -6 -9 ] [1 -0.0125 6 0 ] [ 0 0 -0.125 0 ] [-0.0125 .0515 -6 -4.5 ] [-0.0125 .0515 0 4.5 ] [-0.0125 .14247 -24 -4.5 ] [-0.0125 .14247 0 4.5 ] [-0.0125 .23344 -24 -4.5 ] [-0.0125 .23344 0 4.5 ] [-0.0125 .32441 -24 -4.5 ] [-0.0125 .32441 0 4.5 ] [-0.0125 .41538 -24 -4.5 ] [-0.0125 .41538 0 4.5 ] [-0.0125 .50635 -18 -4.5 ] [-0.0125 .50635 0 4.5 ] [-0.0125 .59732 -24 -4.5 ] [-0.0125 .59732 0 4.5 ] [ 0 0 -0.125 0 ] [ 0 .61803 .125 0 ] [ 1 0 .125 0 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .2 0 m .2 .00625 L s [(9)] .2 -0.0125 0 1 Mshowa .4 0 m .4 .00625 L s [(10)] .4 -0.0125 0 1 Mshowa .6 0 m .6 .00625 L s [(11)] .6 -0.0125 0 1 Mshowa .8 0 m .8 .00625 L s [(12)] .8 -0.0125 0 1 Mshowa 1 0 m 1 .00625 L s [(13)] 1 -0.0125 0 1 Mshowa .125 Mabswid .04 0 m .04 .00375 L s .08 0 m .08 .00375 L s .12 0 m .12 .00375 L s .16 0 m .16 .00375 L s .24 0 m .24 .00375 L s .28 0 m .28 .00375 L s .32 0 m .32 .00375 L s .36 0 m .36 .00375 L s .44 0 m .44 .00375 L s .48 0 m .48 .00375 L s .52 0 m .52 .00375 L s .56 0 m .56 .00375 L s .64 0 m .64 .00375 L s .68 0 m .68 .00375 L s .72 0 m .72 .00375 L s .76 0 m .76 .00375 L s .84 0 m .84 .00375 L s .88 0 m .88 .00375 L s .92 0 m .92 .00375 L s .96 0 m .96 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .0515 m .00625 .0515 L s [(0)] -0.0125 .0515 1 0 Mshowa 0 .14247 m .00625 .14247 L s [(0.02)] -0.0125 .14247 1 0 Mshowa 0 .23344 m .00625 .23344 L s [(0.04)] -0.0125 .23344 1 0 Mshowa 0 .32441 m .00625 .32441 L s [(0.06)] -0.0125 .32441 1 0 Mshowa 0 .41538 m .00625 .41538 L s [(0.08)] -0.0125 .41538 1 0 Mshowa 0 .50635 m .00625 .50635 L s [(0.1)] -0.0125 .50635 1 0 Mshowa 0 .59732 m .00625 .59732 L s [(0.12)] -0.0125 .59732 1 0 Mshowa .125 Mabswid 0 .07425 m .00375 .07425 L s 0 .09699 m .00375 .09699 L s 0 .11973 m .00375 .11973 L s 0 .16522 m .00375 .16522 L s 0 .18796 m .00375 .18796 L s 0 .2107 m .00375 .2107 L s 0 .25619 m .00375 .25619 L s 0 .27893 m .00375 .27893 L s 0 .30167 m .00375 .30167 L s 0 .34716 m .00375 .34716 L s 0 .3699 m .00375 .3699 L s 0 .39264 m .00375 .39264 L s 0 .43813 m .00375 .43813 L s 0 .46087 m .00375 .46087 L s 0 .48361 m .00375 .48361 L s 0 .5291 m .00375 .5291 L s 0 .55184 m .00375 .55184 L s 0 .57458 m .00375 .57458 L s 0 .02876 m .00375 .02876 L s 0 .00602 m .00375 .00602 L s .25 Mabswid 0 0 m 0 .61803 L s 0 .61178 m 0 .61803 L s .2 .61178 m .2 .61803 L s .4 .61178 m .4 .61803 L s .6 .61178 m .6 .61803 L s .8 .61178 m .8 .61803 L s .125 Mabswid .04 .61428 m .04 .61803 L s .08 .61428 m .08 .61803 L s .12 .61428 m .12 .61803 L s .16 .61428 m .16 .61803 L s .24 .61428 m .24 .61803 L s .28 .61428 m .28 .61803 L s .32 .61428 m .32 .61803 L s .36 .61428 m .36 .61803 L s .44 .61428 m .44 .61803 L s .48 .61428 m .48 .61803 L s .52 .61428 m .52 .61803 L s .56 .61428 m .56 .61803 L s .64 .61428 m .64 .61803 L s .68 .61428 m .68 .61803 L s .72 .61428 m .72 .61803 L s .76 .61428 m .76 .61803 L s .84 .61428 m .84 .61803 L s .88 .61428 m .88 .61803 L s .92 .61428 m .92 .61803 L s .96 .61428 m .96 .61803 L s .25 Mabswid 0 .61803 m 1 .61803 L s .99375 .0515 m 1 .0515 L s .99375 .14247 m 1 .14247 L s .99375 .23344 m 1 .23344 L s .99375 .32441 m 1 .32441 L s .99375 .41538 m 1 .41538 L s .99375 .50635 m 1 .50635 L s .99375 .59732 m 1 .59732 L s .125 Mabswid .99625 .07425 m 1 .07425 L s .99625 .09699 m 1 .09699 L s .99625 .11973 m 1 .11973 L s .99625 .16522 m 1 .16522 L s .99625 .18796 m 1 .18796 L s .99625 .2107 m 1 .2107 L s .99625 .25619 m 1 .25619 L s .99625 .27893 m 1 .27893 L s .99625 .30167 m 1 .30167 L s .99625 .34716 m 1 .34716 L s .99625 .3699 m 1 .3699 L s .99625 .39264 m 1 .39264 L s .99625 .43813 m 1 .43813 L s .99625 .46087 m 1 .46087 L s .99625 .48361 m 1 .48361 L s .99625 .5291 m 1 .5291 L s .99625 .55184 m 1 .55184 L s .99625 .57458 m 1 .57458 L s .99625 .02876 m 1 .02876 L s .99625 .00602 m 1 .00602 L s .25 Mabswid 1 0 m 1 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 1 0 0 r .01 w .1 .56653 m .1 .0515 L s .5 0 1 r .3 .43047 m .3 .0515 L s 0 1 1 r .5 .32736 m .5 .0515 L s .5 1 0 r .7 .25036 m .7 .0515 L s 1 0 0 r .9 .19366 m .9 .0515 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgk^kS@0oooo10000000100@4100oooo0?ooo`3oool400000003 010@403oool0oooo0300oooo100000000`0@4100oooo0?ooo`05000003@0oooo10000000200@4100 oooo0?ooo`2HV9P0410@0000000b@3oool0 0`000000oooo0?ooo`020?ooo`0507MgM`2HV9P0oooo09RHV01gMgL0=P3oool00`000000oooo0?oo o`050?ooo`030000003oool0oooo03@0oooo00<000000?ooo`3oool00P3oool01@2HV9P08B4Q0>k^ kP3oool0A4A403H0oooo00<000000?ooo`3oool00P3oool01@0@4100^[Zj0?ooo`1gMgL0@3oool00`000000 oooo0?ooo`020?ooo`050000003oool0oooo0?ooo`000000=P3oool00`000000oooo0?ooo`050?oo o`030000003oool0oooo03@0oooo00<000000?ooo`3oool0103oool00`2HV9P08B4Q0>k^kP0f0?oo o`030000003oool0oooo00D0oooo00<0ZJVY010@403oool00P3oool005h0oooo00D08B4Q08N7Q`3o ool0V9RH0000000i0?ooo`030000003oool0oooo0080oooo00D000000?ooo`3oool0oooo0000000f 0?ooo`030000003oool0oooo00D0oooo00<000000?ooo`3oool0=03oool00`000000oooo0?ooo`05 0?ooo`0307MgM`1VIVH0oooo03D0oooo00<000000?ooo`3oool0103oool00`0@410000000;Zj^P03 0?ooo`00GP3oool01@000000oooo0?ooo`3oool0000003T0oooo00<000000?ooo`3oool00P3oool0 1@0@4100oooo0?ooo`3oool0410@03H0oooo00<000000?ooo`3oool01@3oool00`000000oooo0?oo o`0d0?ooo`030000003oool0oooo00H0oooo00<000000?ooo`3oool0=03oool00`000000oooo0?oo o`050?ooo`030;Zj^P0Q8B40oooo0080oooo001N0?ooo`05024Q8@2HV9P0oooo07MgM`1EEED0>@3o ool00`000000oooo0?ooo`020?ooo`0507MgM`2HV9P0oooo09RHV01gMgL0=P3oool00`000000oooo 0?ooo`050?ooo`030000003oool0oooo03@0oooo00<000000?ooo`3oool00P3oool01@000000^[Zj 0?ooo`2HV9P0410@03H0oooo00<000000?ooo`3oool00P3oool01@0@4100V9RH0?ooo`2j^[X00000 00<0oooo001N0?ooo`050k^kP0d0?ooo`80EEEE00<000000?ooo`3oool00`3oool2 05EEE@030000003oool0oooo0380oooo0P1EEED00`000000oooo0?ooo`020?ooo`0509RHV00@4100 0000024Q8@3k^kP0bk^kP0bk^kPD0oooo0P0Q8B440?ooo`050>k^kP0bk^kP06 0?ooo`030000003oool0oooo01`0oooo0`3o000m0?ooo`<0O`3o?03oool3003oocd0oooo0`1oo`0l 0?ooo`<0o`007P3oool00`000000oooo0?ooo`060?ooo`000`3oool00`000000oooo0?ooo`090?oo o`050000003oool0oooo0?ooo`0000001@3oool00`2HV9P08B4Q0>k^kP050?ooo`<00000703oool3 0?l003d0oooo0`1o0?ll0?ooo`<00?oo?@3oool307oo03`0oooo0`3o000L0?ooo`<00000203oool0 00<0oooo00<000000?ooo`3oool02@3oool01@000000oooo0?ooo`3oool0000000H0oooo00<0MgMg 06IVIP3oool0103oool00`000000oooo0?ooo`0L0?ooo`<0o`00?@3oool307l0oc`0oooo0`00oolm 0?ooo`<0Ool0?03oool30?l001h0oooo00<000000?ooo`3oool01P3oool000<0oooo00<0410@0?oo o`3oool02@3oool01@0@4100oooo0?ooo`3oool0410@00L0oooo00<000000?ooo`3oool00`3oool0 0`000000oooo0?ooo`0L0?ooo`<0o`00?@3oool307l0oc`0oooo0`00oolm0?ooo`<0Ool0?03oool3 0?l001h0oooo00<000000?ooo`3oool01P3oool0000409RHV03oool0V9RH07MgM`/0oooo00D0MgMg 09RHV03oool0V9RH07MgM`030?ooo`050000002j^[X0oooo09RHV00@41001@3oool00`000000oooo 0?ooo`0L0?ooo`<0o`00?@3oool307l0oc`0oooo0`00oolm0?ooo`<0Ool0?03oool30?l001h0oooo 00<000000?ooo`3oool01P3oool00004038bk^kP/0oooo00D0k^k^038bk^kP030?ooo`0509RHV00@41000000024Q8@3k^kP0bk^kPD0 oooo0P0Q8B440?ooo`050>k^kP0bk^kP/0oooo00D0k^k^038bk^kP040?ooo`040k^kP/0oooo00D0k^k^038bk^kP030?ooo`05 0:VYZ@0@41000000010@402YZJT01@3oool00`000000oooo0?ooo`0L0?ooo`<0o`00?@3oool307l0 om`0oooo00<000000?ooo`3oool01P3oool00240oooo00<000000?ooo`3oool0703oool30?l003d0 oooo0`1o0?oL0?ooo`030000003oool0oooo00H0oooo000Q0?ooo`<00000703oool30?l00?l0oooo 703oool2000000P0oooo000Q0?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<0 00000?ooo`3oool01P3oool00240oooo00<000000?ooo`3oool0703oool30?l00?l0oooo7@3oool0 0`000000oooo0?ooo`060?ooo`008@3oool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?oo o`030000003oool0oooo00H0oooo000Q0?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0 oooo00<000000?ooo`3oool01P3oool00240oooo00<000000?ooo`3oool0703oool30?l00?l0oooo 7@3oool00`000000oooo0?ooo`060?ooo`008@3oool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3o oolM0?ooo`030000003oool0oooo00H0oooo000Q0?ooo`030000003oool0oooo01`0oooo0`3o003o 0?oooad0oooo00<000000?ooo`3oool01P3oool00240oooo0`00000L0?ooo`<0o`00o`3ooolL0?oo o`800000203oool00240oooo00<000000?ooo`3oool0703oool30?l00?l0oooo7@3oool00`000000 oooo0?ooo`060?ooo`008@3oool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`030000 003oool0oooo00H0oooo000Q0?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<0 00000?ooo`3oool01P3oool00240oooo00<000000?ooo`3oool0703oool30?l00?l0oooo7@3oool0 0`000000oooo0?ooo`060?ooo`008@3oool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?oo o`030000003oool0oooo00H0oooo000Q0?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0 oooo00<000000?ooo`3oool01P3oool00240oooo0`00000L0?ooo`<0o`00o`3ooolL0?ooo`800000 203oool00240oooo00<000000?ooo`3oool0703oool30?l00?l0oooo7@3oool00`000000oooo0?oo o`060?ooo`008@3oool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`030000003oool0 oooo00H0oooo000Q0?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<000000?oo o`3oool01P3oool00240oooo00<000000?ooo`3oool0703oool30?l00?l0oooo7@3oool00`000000 oooo0?ooo`060?ooo`001`3oool01@3^k^h00?ooo`030000003oool0oooo00@0oooo00<000000?ooo`3oool0703oool30?l00?l0 oooo7@3oool00`000000oooo0?ooo`060?ooo`001`3oool01@0@4100oooo0?ooo`3oool0410@00h0 oooo00<000000?ooo`3oool0103oool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`03 0000003oool0oooo00H0oooo00070?ooo`0507MgM`2HV9P0oooo09RHV01gMgL03P3oool00`000000 oooo0?ooo`040?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<000000?ooo`3o ool01P3oool000L0oooo00D0k^k^038bk^kP0<0?ooo`80EEEE00<000000?ooo`3o ool0103oool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`030000003oool0oooo00H0 oooo000Q0?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<000000?ooo`3oool0 1P3oool00240oooo0`00000L0?ooo`<0o`00o`3ooolL0?ooo`800000203oool00240oooo00<00000 0?ooo`3oool0703oool30?l00?l0oooo7@3oool00`000000oooo0?ooo`060?ooo`008@3oool00`00 0000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`030000003oool0oooo00H0oooo000Q0?ooo`03 0000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<000000?ooo`3oool01P3oool00240oooo 00<000000?ooo`3oool0703oool30?l00?l0oooo7@3oool00`000000oooo0?ooo`060?ooo`008@3o ool00`000000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`030000003oool0oooo00H0oooo000Q 0?ooo`030000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<000000?ooo`3oool01P3oool0 0240oooo00<000000?ooo`3oool0703oool30?l00?l0oooo7@3oool00`000000oooo0?ooo`060?oo o`008@3oool3000001`0oooo0`3o003o0?oooa`0oooo0P0000080?ooo`008@3oool00`000000oooo 0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`030000003oool0oooo00H0oooo000Q0?ooo`030000003o ool0oooo01`0oooo0`3o003o0?oooad0oooo00<000000?ooo`3oool01P3oool00240oooo00<00000 0?ooo`3oool0703oool30?l00?l0oooo7@3oool00`000000oooo0?ooo`060?ooo`008@3oool00`00 0000oooo0?ooo`0L0?ooo`<0o`00o`3ooolM0?ooo`030000003oool0oooo00H0oooo000Q0?ooo`03 0000003oool0oooo01`0oooo0`3o003o0?oooad0oooo00<000000?ooo`3oool01P3oool00240oooo 00<000000?ooo`3oool0o`3oooll0?ooo`030000003oool0oooo00H0oooo000Q0?ooo`<00000o`3o oolk0?ooo`800000203oool00240oooo00<000000?ooo`3oool0o`3oooll0?ooo`030000003oool0 oooo00H0oooo000Q0?ooo`030000003oool0oooo0?l0oooo?03oool00`000000oooo0?ooo`060?oo o`008@3oool00`000000oooo0?ooo`3o0?oooc`0oooo00<000000?ooo`3oool01P3oool00004038b k^kPD0oooo0P0Q8B450?ooo`@0000000<0410@0?ooo`3oool01@0000050?ooo`03 0000003oool0oooo0?l0oooo?03oool00`000000oooo0?ooo`060?ooo`0000@0V9RH0?ooo`2HV9P0 MgMg1@3oool2024Q8@L0oooo00<000000?ooo`3oool00P3oool01@2HV9P08B4Q0>k^kP3oool0A4A4 00D0oooo00<000000?ooo`3oool0o`3oooll0?ooo`030000003oool0oooo00H0oooo00030?ooo`03 010@403oool0oooo00`0oooo00<000000?ooo`3oool00`3oool00`2HV9P08B4Q0>k^kP060?ooo`03 0000003oool0oooo0?l0oooo?03oool00`000000oooo0?ooo`060?ooo`000`3oool00`000000oooo 0?ooo`0<0?ooo`030000003oool0oooo00@0oooo00<0V9RH024Q8@3^k^h01@3oool300000?l0oooo >P3oool3000000P0oooo00030?ooo`030000003oool0oooo00`0oooo00<000000?ooo`3oool01@3o ool00`1gMgL0IVIV0?ooo`040?ooo`030000003oool0oooo0?l0oooo?03oool00`000000oooo0?oo o`060?ooo`000`3oool00`0@4100oooo0?ooo`0<0?ooo`030000003oool0oooo00H0oooo00<00000 0?ooo`3oool00`3oool00`000000oooo0?ooo`3o0?oooc`0oooo00<000000?ooo`3oool01P3oool0 000409RHV03oool0V9RH07MgM`h0oooo00<000000?ooo`3oool00P3oool01@000000^[Zj0?ooo`2H V9P0410@00D0oooo00<000000?ooo`3oool0o`3oooll0?ooo`030000003oool0oooo00H0oooo0000 100b"], ImageRangeCache->{{{0, 359}, {221.375, 0}} -> {7.46508, -0.0252928, \ 0.0196833, 0.000865486}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]], Cell["Options", "Section"] }, Open ]] }, FrontEndVersion->"4.0 for Macintosh", ScreenRectangle->{{0, 1280}, {0, 1004}}, WindowSize->{793, 806}, WindowMargins->{{Automatic, 153}, {21, Automatic}}, Magnification->1.25, StyleDefinitions -> "ArticleClassic.nb", MacintoshSystemPageSetup->"\<\ 00<0001804P000000`d26_oQon@3:`8g0dL5N`?P0080001804P000000]P2:001 0000I00000400`<30?l00BL?00400@00000000000000060801T1T00000000000 00000000004000000000000000000000\>" ] (*********************************************************************** 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->{ "SA.5.1"->{ Cell[7450, 215, 2476, 48, 88, "Text", CellTags->{"SA.5.1", "1.5"}]}, "1.5"->{ Cell[7450, 215, 2476, 48, 88, "Text", CellTags->{"SA.5.1", "1.5"}]} } *) (*CellTagsIndex CellTagsIndex->{ {"SA.5.1", 45152, 1130}, {"1.5", 45242, 1133} } *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1739, 51, 47, 0, 84, "Title"], Cell[1789, 53, 91, 4, 71, "Subsubsection"], Cell[CellGroupData[{ Cell[1905, 61, 45, 0, 67, "Section"], Cell[1953, 63, 589, 11, 107, "Text"], Cell[2545, 76, 1424, 25, 556, "Program"], Cell[3972, 103, 150, 5, 88, "Text"], Cell[4125, 110, 61, 1, 34, "Input"], Cell[CellGroupData[{ Cell[4211, 115, 45, 1, 34, "Input"], Cell[4259, 118, 84, 1, 34, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4380, 124, 52, 1, 34, "Input"], Cell[4435, 127, 405, 6, 123, "Print"] }, Open ]], Cell[4855, 136, 553, 9, 107, "Text"], Cell[5411, 147, 412, 8, 69, "Text"] }, Closed]], Cell[CellGroupData[{ Cell[5860, 160, 33, 0, 34, "Section"], Cell[5896, 162, 195, 4, 50, "Text"], Cell[6094, 168, 40, 1, 34, "Input"], Cell[6137, 171, 418, 10, 69, "Text"], Cell[6558, 183, 889, 30, 88, "Text"], Cell[7450, 215, 2476, 48, 88, "Text", CellTags->{"SA.5.1", "1.5"}], Cell[9929, 265, 28, 0, 31, "Text"], Cell[9960, 267, 56, 1, 34, "Input"], Cell[10019, 270, 54, 1, 34, "Input"], Cell[CellGroupData[{ Cell[10098, 275, 36, 1, 34, "Input"], Cell[10137, 278, 36, 1, 34, "Output"] }, Open ]], Cell[10188, 282, 60, 1, 34, "Input"], Cell[CellGroupData[{ Cell[10273, 287, 36, 1, 34, "Input"], Cell[10312, 290, 36, 1, 34, "Output"] }, Open ]], Cell[10363, 294, 344, 6, 69, "Text"], Cell[CellGroupData[{ Cell[10732, 304, 102, 2, 34, "Input"], Cell[10837, 308, 33257, 788, 233, 4760, 432, "GraphicsData", "PostScript", \ "Graphics"], Cell[44097, 1098, 130, 3, 34, "Output"] }, Open ]] }, Closed]], Cell[44254, 1105, 26, 0, 34, "Section"] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)