(************** 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[ 45830, 1418]*) (*NotebookOutlinePosition[ 46472, 1440]*) (* CellTagsIndexPosition[ 46428, 1436]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{\(Remove["\"];\), "\[IndentingNewLine]", RowBox[{ RowBox[{"programMain", ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{\(k1 = 8.617*10^\(-5\)\), ";", "\[IndentingNewLine]", \(A1 = 10^\(-22\)\), ";", \(A2 = 30*10^\(-22\)\), ";", \(Ak = 10^\(-22\)\), ";", "\[IndentingNewLine]", \(N1 = 10^23\), ";", \(N2 = 10^23\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", \(solveDiffeq[n10_, n20_, nc0_, R_, tfinal_] := \[IndentingNewLine]Module[{t}, sol = \[IndentingNewLine]NDSolve[{\(n1'\)[t] \[Equal] A1*\((N1 - n1[t])\)*nc[t], \(n2'\)[t] \[Equal] A2*\((N2 - n2[t])\)* nc[t], \[IndentingNewLine]\ \(nc'\)[t] \[Equal] R - \(n1'\)[t] - \(n2'\)[t] - Ak*nc[t]*\((n1[t] + n2[t] + nc[t])\), \ n1[0] \[Equal] n10, n2[0] \[Equal] n20, nc[0] \[Equal] nc0}, {n1, n2, nc}, {t, 0, tfinal}, MaxSteps \[Rule] 50000]]\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", \(initValues[b_, d_] := Module[{}, n10 = Last[n1[d] /. b]; n20 = Last[n2[d] /. b]; nc0 = Last[nc[d] /. b]]\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", StyleBox[\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(---\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\), FontColor->RGBColor[0, 1, 0]], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", \( (*irradiate\ for\ irrTime*) \), "\[IndentingNewLine]", \(R = 10^21\), ";", \(tfinal7 = irrTime\), " ", ";", "\[IndentingNewLine]", \(solveDiffeq[n10, n20, nc0, R, tfinal7]\), ";", "\[IndentingNewLine]", \(sol7 = sol\), ";", "\[IndentingNewLine]", \(graphAllEleven[sol7, tfinal7]\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", \( (*\ Relaxation\ stage\ *) \), "\[IndentingNewLine]", \(initValues[sol7, tfinal7]\), ";", "\[IndentingNewLine]", \(R = 0\), ";", \(tfinal8 = 60\), ";", "\[IndentingNewLine]", \(solveDiffeq[n10, n20, nc0, R, tfinal8]\), ";", "\[IndentingNewLine]", \(sol8 = sol\), ";", "\[IndentingNewLine]", \(graphAllEleven[sol8, tfinal8]\), ";"}], "\[IndentingNewLine]", ")"}]}], "\[IndentingNewLine]", RowBox[{"(*", StyleBox[\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(-\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(---\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\), FontColor->RGBColor[0, 1, 0]], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", \(tstart = 1; tend = 800; tstep = 30;\), "\[IndentingNewLine]", RowBox[{\(n1List1 = {}; n2List1 = {}; ncList1 = {}; logn1List1 = {};\), "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{\(For[tloop = tstart, tloop \[LessEqual] tend, tloop += tstep, \[IndentingNewLine]irrTime = tloop; n10 = 0; n20 = 0; nc0 = 0; \[IndentingNewLine]programMain; \ \[IndentingNewLine]initValues[sol8, tfinal8]; (*find\ n1, n2, nc\ at\ end\ of\ relaxation\ section*) ; \ \[IndentingNewLine]AppendTo[n1List1, {irrTime, n10}]; AppendTo[n2List1, {irrTime, n20}]; \[IndentingNewLine]AppendTo[ ncList1, {irrTime, nc0}]; AppendTo[ logn1List1, {Log[irrTime], Log[n10]}];\[IndentingNewLine]\[IndentingNewLine]]\), "\[IndentingNewLine]"}], "\[IndentingNewLine]", \(ListPlot[n1List1, PlotRange \[Rule] All, PlotJoined \[Rule] True, PlotLabel \[Rule] "\", ImageSize \[Rule] 72\ 3];\), "\[IndentingNewLine]", \(ListPlot[ n2List1, PlotRange \[Rule] All, PlotJoined \[Rule] True, PlotLabel \[Rule] "\", ImageSize \[Rule] 72\ 3];\), "\[IndentingNewLine]", \(ListPlot[ logn1List1, PlotRange \[Rule] All, PlotJoined \[Rule] True, PlotLabel \[Rule] "\", ImageSize \[Rule] 72\ 3];\), "\n", "Null"}], "Input"], Cell[BoxData[ RowBox[{\(Remove::"rmnsm"\), \(\(:\)\(\ \)\), "\<\"There are no symbols \ matching \\\"\\!\\(\\\"Global`*\\\"\\)\\\". \\!\\(\\*ButtonBox[\\\"More\ \[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"Remove::rmnsm\\\"]\\)\"\>"}]], "Message"], Cell[CellGroupData[{ Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 %%ImageSize: 216 133.495 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.0225885 0.001221 0.0145179 6.08696e-024 [ [.26679 .00202 -9 -9 ] [.26679 .00202 9 0 ] [.51099 .00202 -9 -9 ] [.51099 .00202 9 0 ] [.75519 .00202 -9 -9 ] [.75519 .00202 9 0 ] [.99939 .00202 -9 -9 ] [.99939 .00202 9 0 ] [.01009 .13626 -33 -6.25 ] [.01009 .13626 0 6.25 ] [.01009 .258 -33 -6.25 ] [.01009 .258 0 6.25 ] [.01009 .37974 -33 -6.25 ] [.01009 .37974 0 6.25 ] [.01009 .50147 -33 -6.25 ] [.01009 .50147 0 6.25 ] [.5 .63053 -23 0 ] [.5 .63053 23 12.5625 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .26679 .01452 m .26679 .02077 L s [(200)] .26679 .00202 0 1 Mshowa .51099 .01452 m .51099 .02077 L s [(400)] .51099 .00202 0 1 Mshowa .75519 .01452 m .75519 .02077 L s [(600)] .75519 .00202 0 1 Mshowa .99939 .01452 m .99939 .02077 L s [(800)] .99939 .00202 0 1 Mshowa .125 Mabswid .08364 .01452 m .08364 .01827 L s .14469 .01452 m .14469 .01827 L s .20574 .01452 m .20574 .01827 L s .32784 .01452 m .32784 .01827 L s .38889 .01452 m .38889 .01827 L s .44994 .01452 m .44994 .01827 L s .57204 .01452 m .57204 .01827 L s .63309 .01452 m .63309 .01827 L s .69414 .01452 m .69414 .01827 L s .81624 .01452 m .81624 .01827 L s .87729 .01452 m .87729 .01827 L s .93834 .01452 m .93834 .01827 L s .25 Mabswid 0 .01452 m 1 .01452 L s .02259 .13626 m .02884 .13626 L s gsave .01009 .13626 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (2) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .02259 .258 m .02884 .258 L s gsave .01009 .258 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (4) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .02259 .37974 m .02884 .37974 L s gsave .01009 .37974 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (6) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .02259 .50147 m .02884 .50147 L s gsave .01009 .50147 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (8) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .125 Mabswid .02259 .04495 m .02634 .04495 L s .02259 .07539 m .02634 .07539 L s .02259 .10582 m .02634 .10582 L s .02259 .16669 m .02634 .16669 L s .02259 .19713 m .02634 .19713 L s .02259 .22756 m .02634 .22756 L s .02259 .28843 m .02634 .28843 L s .02259 .31887 m .02634 .31887 L s .02259 .3493 m .02634 .3493 L s .02259 .41017 m .02634 .41017 L s .02259 .44061 m .02634 .44061 L s .02259 .47104 m .02634 .47104 L s .02259 .53191 m .02634 .53191 L s .02259 .56234 m .02634 .56234 L s .02259 .59278 m .02634 .59278 L s .25 Mabswid .02259 0 m .02259 .61803 L s gsave .5 .63053 -84 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20.5625 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (n1) show 81.000 12.813 moveto (vs) show 99.000 12.813 moveto (t) show 105.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .02381 .01472 m .06044 .02161 L .09707 .03152 L .1337 .04858 L .17033 .08932 L .20696 .15937 L .24359 .22391 L .28022 .27951 L .31685 .32738 L .35348 .36857 L .39011 .40404 L .42674 .43456 L .46337 .46083 L .5 .48344 L .53663 .50291 L .57326 .51966 L .60989 .53408 L .64652 .5465 L .68315 .55718 L .71978 .56638 L .75641 .57429 L .79304 .5811 L .82967 .58697 L .8663 .59202 L .90293 .59636 L .93956 .6001 L .97619 .60332 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{216, 133.438}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol002=oo`03001oogoo025oo`03001oogoo08ioo`008goo00<007ooOol0 8Woo00<007ooOol0SGoo000SOol00`00Oomoo`0ROol00`00Oomoo`2=Ool002=oo`03001oogoo02=o o`03001oogoo08aoo`008goo00<007ooOol097oo00<007ooOol0Rgoo000SOol00`00Oomoo`0TOol0 0`00Oomoo`2;Ool000=oo`<00003Ool007oo009oo`04001oogooOol50002Ool2000;Ool00`00Oomo o`0UOol00`00Oomoo`2:Ool000Aoo`07001oogooOol007oo00001Woo00D007ooOomoo`00009oo`03 001oogoo00Qoo`03001oogoo02Ioo`03001oogoo08Uoo`0000=oo`0000000`000goo00<007ooOol0 1Goo00D007ooOomoo`00009oo`03001oogoo00Qoo`03001oogoo02Moo`03001oogoo08Qoo`0000Eo o`00Oomoo`0000=oo`03001oo`0000Ioo`05001oogooOol00002Ool00`00Oomoo`08Ool2000XOol0 0`00Oomoo`28Ool0009oo`03001oo`00009oo`05001oogooOol00005Ool01@00Oomoogoo00000Woo 00<007ooOol027oo00<007ooOol0:7oo00<007ooOol0Qgoo0003Ool2000;Ool20004Ool20003Ool3 0002Ool4000[Ool00`00Oomoo`26Ool001]oo`8000=oo`800003Ool007oo02Yoo`03001oogoo08Io o`007Goo00<007ooOol00Woo0P00;7oo00<007ooOol0QGoo000MOol00`00Oomoo`02Ool2000]Ool0 0`00Oomoo`24Ool001]oo`<0009oo`@002ioo`03001oogoo08=oo`008goo00<007ooOol0;7oo00<0 07ooOol0Pgoo000SOol00`00Oomoo`0]Ool00`00Oomoo`22Ool002=oo`03001oogoo02ioo`03001o ogoo085oo`008goo00<007ooOol0;goo00<007ooOol0P7oo000SOol00`00Oomoo`0_Ool00`00Oomo o`20Ool002=oo`03001oogoo031oo`03001oogoo07moo`008goo00<007ooOol07oo00<0 07ooOol0Mgoo0002Ool00`00Oomoo`02Ool01@00Oomoogoo00001Goo00D007ooOomoo`00009oo`03 001oogoo00Qoo`8003Yoo`03001oogoo07Ioo`000Woo0`002goo0P0017oo0P002goo00<007ooOol0 >Woo00<007ooOol0MGoo000KOol30002Ool4000mOol00`00Oomoo`1dOol001]oo`8000=oo`800003 Ool007oo03eoo`03001oogoo07=oo`007Goo00<007ooOol00Woo0P00?goo00<007ooOol0LWoo000M Ool00`00Oomoo`02Ool20010Ool00`00Oomoo`1aOol001]oo`<0009oo`@0045oo`03001oogoo071o o`008goo00<007ooOol0@7oo0P00L7oo000SOol00`00Oomoo`12Ool00`00Oomoo`1]Ool002=oo`03 001oogoo04=oo`03001oogoo06aoo`008goo00<007ooOol0A7oo0P00K7oo000SOol00`00Oomoo`16 Ool00`00Oomoo`1YOol002=oo`03001oogoo04Moo`03001oogoo06Qoo`008goo00<007ooOol0B7oo 0P00J7oo000SOol00`00Oomoo`1:Ool00`00Oomoo`1UOol002=oo`03001oogoo04]oo`8006Eoo`00 8goo00<007ooOol0CGoo00<007ooOol0HWoo000SOol00`00Oomoo`1>Ool2001ROol0009oo`8000=o o`05001oogooOol00003Ool50002Ool2000;Ool00`00Oomoo`1@Ool2001POol00005Ool007ooOol0 0003Ool00`00Ool00006Ool01@00Oomoogoo00000Woo00<007ooOol027oo00<007ooOol0DWoo0P00 GWoo00001Goo001oogoo000017oo00<007ooOol01Goo00D007ooOomoo`00009oo`03001oogoo00Qo o`03001oogoo05Aoo`8005aoo`000Woo0P0017oo00<007oo00001Woo00D007ooOomoo`00009oo`03 001oogoo00Qoo`8005Moo`8005Yoo`0000Eoo`00Oomoo`00009oo`05001oogooOol00005Ool01@00 Oomoogoo00000Woo00<007ooOol027oo00<007ooOol0F7oo0P00F7oo0002Ool2000Ool0069oo`04001oogoo0004 Ool00`00Oomoo`07Ool01000Oomoo`000Woo0`002Goo00<007ooOol0Cgoo001QOol40005Ool00`00 Oomoo`06Ool20002Ool200000goo000000020007Ool5001>Ool006Yoo`03001oogoo01Uoo`03001o ogoo04moo`00JGoo0P00KGoo003HOol00=Qoo`00f7oo0000\ \>"], ImageRangeCache->{{{0, 215}, {132.438, 0}} -> {-164.796, -1.09631*^22, \ 4.6835, 9.39476*^20}}], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 %%ImageSize: 216 133.495 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.0225885 0.001221 0.0089651 5.94354e-024 [ [.26679 -0.00353 -9 -9 ] [.26679 -0.00353 9 0 ] [.51099 -0.00353 -9 -9 ] [.51099 -0.00353 9 0 ] [.75519 -0.00353 -9 -9 ] [.75519 -0.00353 9 0 ] [.99939 -0.00353 -9 -9 ] [.99939 -0.00353 9 0 ] [.01009 .12784 -33 -6.25 ] [.01009 .12784 0 6.25 ] [.01009 .24671 -33 -6.25 ] [.01009 .24671 0 6.25 ] [.01009 .36558 -33 -6.25 ] [.01009 .36558 0 6.25 ] [.01009 .48445 -33 -6.25 ] [.01009 .48445 0 6.25 ] [.01009 .60332 -33 -6.25 ] [.01009 .60332 0 6.25 ] [.5 .63053 -23 0 ] [.5 .63053 23 12.5625 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .26679 .00897 m .26679 .01522 L s [(200)] .26679 -0.00353 0 1 Mshowa .51099 .00897 m .51099 .01522 L s [(400)] .51099 -0.00353 0 1 Mshowa .75519 .00897 m .75519 .01522 L s [(600)] .75519 -0.00353 0 1 Mshowa .99939 .00897 m .99939 .01522 L s [(800)] .99939 -0.00353 0 1 Mshowa .125 Mabswid .08364 .00897 m .08364 .01272 L s .14469 .00897 m .14469 .01272 L s .20574 .00897 m .20574 .01272 L s .32784 .00897 m .32784 .01272 L s .38889 .00897 m .38889 .01272 L s .44994 .00897 m .44994 .01272 L s .57204 .00897 m .57204 .01272 L s .63309 .00897 m .63309 .01272 L s .69414 .00897 m .69414 .01272 L s .81624 .00897 m .81624 .01272 L s .87729 .00897 m .87729 .01272 L s .93834 .00897 m .93834 .01272 L s .25 Mabswid 0 .00897 m 1 .00897 L s .02259 .12784 m .02884 .12784 L s gsave .01009 .12784 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (2) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .02259 .24671 m .02884 .24671 L s gsave .01009 .24671 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (4) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .02259 .36558 m .02884 .36558 L s gsave .01009 .36558 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (6) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .02259 .48445 m .02884 .48445 L s gsave .01009 .48445 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (8) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (22) show grestore .02259 .60332 m .02884 .60332 L s gsave .01009 .60332 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (1) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (23) show grestore .125 Mabswid .02259 .03868 m .02634 .03868 L s .02259 .0684 m .02634 .0684 L s .02259 .09812 m .02634 .09812 L s .02259 .15755 m .02634 .15755 L s .02259 .18727 m .02634 .18727 L s .02259 .21699 m .02634 .21699 L s .02259 .27642 m .02634 .27642 L s .02259 .30614 m .02634 .30614 L s .02259 .33586 m .02634 .33586 L s .02259 .3953 m .02634 .3953 L s .02259 .42501 m .02634 .42501 L s .02259 .45473 m .02634 .45473 L s .02259 .51417 m .02634 .51417 L s .02259 .54388 m .02634 .54388 L s .02259 .5736 m .02634 .5736 L s .25 Mabswid .02259 0 m .02259 .61803 L s gsave .5 .63053 -84 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20.5625 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (n2) show 81.000 12.813 moveto (vs) show 99.000 12.813 moveto (t) show 105.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .02381 .01472 m .06044 .18516 L .09707 .34919 L .1337 .49767 L .17033 .59168 L .20696 .60315 L .24359 .60332 L .28022 .60332 L .31685 .60332 L .35348 .60332 L .39011 .60332 L .42674 .60332 L .46337 .60332 L .5 .60332 L .53663 .60332 L .57326 .60332 L .60989 .60332 L .64652 .60332 L .68315 .60332 L .71978 .60332 L .75641 .60332 L .79304 .60332 L .82967 .60332 L .8663 .60332 L .90293 .60332 L .93956 .60332 L .97619 .60332 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{216, 133.438}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol00`00Oomoo`2SOol002=oo`03001oogoo00aoo`03 001oogoo0:=oo`008goo00<007ooOol03Goo00<007ooOol0XWoo000SOol00`00Oomoo`0=Ool00`00 Oomoo`2ROol002=oo`03001oogoo00eoo`03001oogoo0:9oo`008goo00<007ooOol03Goo00<007oo Ool0XWoo000SOol00`00Oomoo`0=Ool00`00Oomoo`2ROol002=oo`03001oogoo00ioo`03001oogoo 0:5oo`008goo00<007ooOol03Woo00<007ooOol0XGoo000SOol00`00Oomoo`0>Ool00`00Oomoo`2Q Ool002=oo`03001oogoo00ioo`03001oogoo0:5oo`000Woo0P000goo00D007ooOomoo`0000=oo`D0 009oo`8000]oo`03001oogoo00moo`03001oogoo0:1oo`0000Eoo`00Oomoo`0000=oo`03001oo`00 00Ioo`05001oogooOol00002Ool00`00Oomoo`08Ool00`00Oomoo`0?Ool00`00Oomoo`2POol00005 Ool007ooOol00004Ool00`00Oomoo`05Ool01@00Oomoogoo00000Woo00<007ooOol027oo00<007oo Ool03goo00<007ooOol0X7oo0002Ool20004Ool00`00Ool00006Ool01@00Oomoogoo00000Woo00<0 07ooOol027oo0P0047oo00<007ooOol0X7oo00001Goo001oogoo00000Woo00D007ooOomoo`0000Eo o`05001oogooOol00002Ool00`00Oomoo`08Ool00`00Oomoo`0@Ool00`00Oomoo`2OOol0009oo`80 00aoo`8000Aoo`8000=oo`<0009oo`@0019oo`03001oogoo09moo`006goo0P000goo0P0000=oo`00 Ool04Goo00<007ooOol0Wgoo000MOol00`00Oomoo`02Ool2000BOol00`00Oomoo`2OOol001eoo`03 001oogoo009oo`8001=oo`03001oogoo09ioo`006goo0`000Woo10004goo00<007ooOol0WWoo000S Ool00`00Oomoo`0BOol00`00Oomoo`2MOol002=oo`03001oogoo019oo`03001oogoo09eoo`008goo 00<007ooOol04goo00<007ooOol0W7oo000SOol00`00Oomoo`0COol00`00Oomoo`2LOol002=oo`03 001oogoo01Aoo`03001oogoo09]oo`008goo00<007ooOol057oo00<007ooOol0Vgoo000SOol00`00 Oomoo`0DOol00`00Oomoo`2KOol002=oo`03001oogoo01Eoo`03001oogoo09Yoo`008goo00<007oo Ool05Goo00<007ooOol0VWoo000SOol00`00Oomoo`0FOol00`00Oomoo`2IOol002=oo`03001oogoo 01Ioo`03001oogoo09Uoo`0000=oo`0000000`0000=oo`00Ool00Woo00@007ooOomoo`D0009oo`80 00]oo`03001oogoo01Moo`03001oogoo09Qoo`000goo00<007ooOol00Woo00<007oo00001Woo00D0 07ooOomoo`00009oo`03001oogoo00Qoo`03001oogoo01Moo`8009Uoo`000goo00<007ooOol00goo 00<007ooOol01Goo00D007ooOomoo`00009oo`03001oogoo00Qoo`03001oogoo01Uoo`<009Ioo`00 0goo00<007ooOol00Woo00<007oo00001Woo00D007ooOomoo`00009oo`03001oogoo00Qoo`8001eo ohP000ioo`000goo00D007ooOomoo`0000=oo`03001oogoo00=oo`05001oogooOol00002Ool00`00 Oomoo`08Ool00`00Oomoo`2bOol0009oo`8000aoo`8000Aoo`8000=oo`<0009oo`@00;Aoo`006goo 0P001Goo0P00]7oo000MOol01@00Oomoogoo0000]Woo000MOol00`00Oomoo`02Ool00`00Oomoo`2c Ool001]oo`<0009oo`<00;Eoo`00f7oo003HOol00=Qoo`00HGoo0P000Woo0P0000=oo`0000000P00 2Goo0P000goo0`002Woo0P00Cgoo001ROol01000Oomoo`000Woo00<007ooOol02Woo0P001Woo00<0 07ooOol01Woo00@007ooOol004ioo`00HWoo00@007ooOol000=oo`03001oogoo00Qoo`04001oogoo 0002Ool30009Ool00`00Oomoo`1?Ool0065oo`@000Eoo`03001oogoo00Ioo`80009oo`800003Ool0 0000008000Moo`D004ioo`00J7oo00@007ooOol001Yoo`03001oogoo04moo`00JGoo0P00KGoo003H Ool00=Qoo`00\ \>"], ImageRangeCache->{{{0, 215}, {132.438, 0}} -> {-164.815, -1.07624*^22, \ 4.68366, 9.62179*^20}}], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 %%ImageSize: 216 133.495 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.142988 -3.29 0.0735613 [ [.1668 .08132 -3 -9 ] [.1668 .08132 3 0 ] [.30979 .08132 -3 -9 ] [.30979 .08132 3 0 ] [.45277 .08132 -3 -9 ] [.45277 .08132 3 0 ] [.59576 .08132 -3 -9 ] [.59576 .08132 3 0 ] [.73875 .08132 -3 -9 ] [.73875 .08132 3 0 ] [.88174 .08132 -3 -9 ] [.88174 .08132 3 0 ] [.01131 .24094 -12 -4.5 ] [.01131 .24094 0 4.5 ] [.01131 .38806 -12 -4.5 ] [.01131 .38806 0 4.5 ] [.01131 .53518 -12 -4.5 ] [.01131 .53518 0 4.5 ] [.5 .63053 -47 0 ] [.5 .63053 47 12.5625 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .1668 .09382 m .1668 .10007 L s [(1)] .1668 .08132 0 1 Mshowa .30979 .09382 m .30979 .10007 L s [(2)] .30979 .08132 0 1 Mshowa .45277 .09382 m .45277 .10007 L s [(3)] .45277 .08132 0 1 Mshowa .59576 .09382 m .59576 .10007 L s [(4)] .59576 .08132 0 1 Mshowa .73875 .09382 m .73875 .10007 L s [(5)] .73875 .08132 0 1 Mshowa .88174 .09382 m .88174 .10007 L s [(6)] .88174 .08132 0 1 Mshowa .125 Mabswid .05241 .09382 m .05241 .09757 L s .081 .09382 m .081 .09757 L s .1096 .09382 m .1096 .09757 L s .1382 .09382 m .1382 .09757 L s .19539 .09382 m .19539 .09757 L s .22399 .09382 m .22399 .09757 L s .25259 .09382 m .25259 .09757 L s .28119 .09382 m .28119 .09757 L s .33838 .09382 m .33838 .09757 L s .36698 .09382 m .36698 .09757 L s .39558 .09382 m .39558 .09757 L s .42418 .09382 m .42418 .09757 L s .48137 .09382 m .48137 .09757 L s .50997 .09382 m .50997 .09757 L s .53857 .09382 m .53857 .09757 L s .56716 .09382 m .56716 .09757 L s .62436 .09382 m .62436 .09757 L s .65296 .09382 m .65296 .09757 L s .68155 .09382 m .68155 .09757 L s .71015 .09382 m .71015 .09757 L s .76735 .09382 m .76735 .09757 L s .79594 .09382 m .79594 .09757 L s .82454 .09382 m .82454 .09757 L s .85314 .09382 m .85314 .09757 L s .91033 .09382 m .91033 .09757 L s .93893 .09382 m .93893 .09757 L s .96753 .09382 m .96753 .09757 L s .99613 .09382 m .99613 .09757 L s .25 Mabswid 0 .09382 m 1 .09382 L s .02381 .24094 m .03006 .24094 L s [(48)] .01131 .24094 1 0 Mshowa .02381 .38806 m .03006 .38806 L s [(50)] .01131 .38806 1 0 Mshowa .02381 .53518 m .03006 .53518 L s [(52)] .01131 .53518 1 0 Mshowa .125 Mabswid .02381 .1306 m .02756 .1306 L s .02381 .16738 m .02756 .16738 L s .02381 .20416 m .02756 .20416 L s .02381 .27772 m .02756 .27772 L s .02381 .3145 m .02756 .3145 L s .02381 .35128 m .02756 .35128 L s .02381 .42484 m .02756 .42484 L s .02381 .46162 m .02756 .46162 L s .02381 .4984 m .02756 .4984 L s .02381 .05704 m .02756 .05704 L s .02381 .02026 m .02756 .02026 L s .02381 .57197 m .02756 .57197 L s .02381 .60875 m .02756 .60875 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .5 .63053 -108 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20.5625 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (ln) show %%IncludeResource: font Mathematica2Mono %%IncludeFont: Mathematica2Mono /Mathematica2Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 75.000 12.813 moveto (H) show 81.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (n1) show %%IncludeResource: font Mathematica2Mono %%IncludeFont: Mathematica2Mono /Mathematica2Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 93.000 12.813 moveto (L) show 105.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (vs) show 123.000 12.813 moveto (ln) show %%IncludeResource: font Mathematica2Mono %%IncludeFont: Mathematica2Mono /Mathematica2Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 135.000 12.813 moveto (H) show 141.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (t) show %%IncludeResource: font Mathematica2Mono %%IncludeFont: Mathematica2Mono /Mathematica2Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 147.000 12.813 moveto (L) show 153.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .02381 .01472 m .51483 .27825 L .61161 .34254 L .66881 .39367 L .70955 .45154 L .74122 .50016 L .76713 .52726 L .78906 .54459 L .80807 .5568 L .82484 .5659 L .83986 .57293 L .85344 .57847 L .86585 .58294 L .87726 .58657 L .88783 .58956 L .89767 .59205 L .90688 .59412 L .91553 .59585 L .92369 .59732 L .93141 .59855 L .93873 .5996 L .9457 .60049 L .95234 .60125 L .95868 .60189 L .96476 .60244 L .97059 .60292 L .97619 .60332 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{216, 133.438}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHggoo00<007ooOol0PWoo000EOol00`00Oomoo`0lOol20022Ool0 01Eoo`03001oogoo03ioo`80081oo`005Goo0P00@Goo0P00OWoo000EOol00`00Oomoo`12Ool2001l Ool001Eoo`03001oogoo04Aoo`8007Yoo`005Goo00<007ooOol0AWoo0P00N7oo000EOol00`00Oomo o`18Ool2001fOol000Yoo`<0009oo`8000Aoo`03001oogoo04Yoo`8007Aoo`002goo00@007ooOol0 009oo`05001oogooOol0001>Ool00`00Oomoo`1aOol000Qoo`D00005Ool007ooOol00003Ool3001= Ool2001aOol000Qoo`04001oogoo0003Ool20004Ool00`00Oomoo`1?Ool2001_Ool000Uoo`03001o o`00009oo`04001oogoo0003Ool00`00Oomoo`1AOol2001]Ool000Yoo`8000=oo`8000Aoo`03001o ogoo05=oo`8006]oo`005Goo00<007ooOol0EGoo0P00JGoo000EOol00`00Oomoo`1GOol2001WOol0 01Eoo`03001oogoo05Uoo`8006Eoo`005Goo0P00G7oo00<007ooOol0HWoo000EOol00`00Oomoo`1L Ool2001ROol001Eoo`03001oogoo05ioo`03001oogoo05moo`005Goo00<007ooOol0Ggoo0P00Ggoo 000EOol00`00Oomoo`1QOol00`00Oomoo`1LOol001Eoo`03001oogoo069oo`8005aoo`005Goo00<0 07ooOol0I7oo00<007ooOol0FGoo000EOol2001VOol2001IOol001Eoo`03001oogoo06Moo`03001o ogoo05Ioo`005Goo00<007ooOol0J7oo0P00EWoo000EOol00`00Oomoo`1ZOol00`00Oomoo`1COol0 01Eoo`03001oogoo06]oo`8005=oo`005Goo00<007ooOol0KGoo00<007ooOol0D7oo000EOol00`00 Oomoo`1^Ool00`00Oomoo`1?Ool001Eoo`80071oo`03001oogoo04ioo`005Goo00<007ooOol0L7oo 00<007ooOol0CGoo000EOol00`00Oomoo`1aOol00`00Oomoo`1goo000EOol0 0`00Oomoo`23Ool00`00Oomoo`0jOol001Eoo`03001oogoo08=oo`03001oogoo03Yoo`005Goo00<0 07ooOol0Q7oo00<007ooOol0>Goo000EOol00`00Oomoo`25Ool00`00Oomoo`0hOol001Eoo`03001o ogoo08Eoo`03001oogoo03Qoo`005Goo0P00Qgoo00<007ooOol0=goo000EOol00`00Oomoo`27Ool0 0`00Oomoo`0fOol001Eoo`03001oogoo08Qoo`03001oogoo03Eoo`005Goo00<007ooOol0RGoo00<0 07ooOol0=7oo000EOol00`00Oomoo`2:Ool00`00Oomoo`0cOol000Uoo`8000=oo`@000=oo`03001o ogoo08]oo`03001oogoo039oo`0027oo00@007ooOol0009oo`03001oogoo00Aoo`03001oogoo08ao o`03001oogoo035oo`002goo00D007ooOomoo`0000Eoo`<008eoo`03001oogoo031oo`002Goo0P00 1Goo00<007ooOol00Woo00<007ooOol0SWoo00<007ooOol0;goo0009Ool00`00Oomoo`02Ool01000 Oomoo`000goo00<007ooOol0Sgoo00<007ooOol0;Woo0009Ool30003Ool20004Ool00`00Oomoo`2@ Ool2000^Ool001Eoo`03001oogoo099oo`03001oogoo02]oo`005Goo00<007ooOol0Tgoo0P00:goo 000EOol00`00Oomoo`2EOol3000XOol001Eoo`8009Uoo`8002Ioo`005Goo00<007ooOol0VWoo0P00 97oo000EOol00`00Oomoo`2LOol3000QOol001Eoo`03001oogoo09moo`@001eoo`005Goo00<007oo Ool0Xgoo1@0067oo000EOol00`00Oomoo`2XOol8000@Ool001Eoo`03001oogoo0;1oo`<000eoo`00 5Goo0P00`Goo000EOol00`00Oomoo`30Ool001Eoo`03001oogoo0<1oo`00f7oo003HOol00=Qoo`00 f7oo001COol00`00Oomoo`0"], ImageRangeCache->{{{0, 215}, {132.438, 0}} -> {-0.805373, 44.7246, \ 0.0367859, 0.0715043}}] }, Open ]] }, Open ]] }, FrontEndVersion->"5.0 for Microsoft Windows", ScreenRectangle->{{0, 800}, {0, 527}}, WindowSize->{792, 500}, WindowMargins->{{0, Automatic}, {Automatic, 0}} ] (******************************************************************* 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, 4471, 74, 990, "Input"], Cell[6250, 129, 299, 4, 22, "Message"], Cell[CellGroupData[{ Cell[6574, 137, 12597, 402, 142, 5406, 310, "GraphicsData", "PostScript", \ "Graphics"], Cell[19174, 541, 13527, 427, 142, 5804, 328, "GraphicsData", "PostScript", \ "Graphics"], Cell[32704, 970, 13098, 444, 142, 6335, 357, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)