(************** 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[ 53540, 1691]*) (*NotebookOutlinePosition[ 54271, 1716]*) (* CellTagsIndexPosition[ 54227, 1712]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{\(Remove["\"];\), "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"programMain", ":=", RowBox[{"(", "\[IndentingNewLine]", RowBox[{\(A = 10^\(-17\)\), ";", \(Ar = 10^\(-13\)\), ";", \(Ah = 10^\(-15\)\), ";", "\[IndentingNewLine]", \(N1 = 10^15\), ";", \(Nh = 3*10^14\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", \(solveDiffeq[n10_, nh0_, nv0_, nc0_, R_, tfinal_] := Module[{t}, sol = \[IndentingNewLine]NDSolve[{\(n1'\)[t] \[Equal] A*\((N1 - n1[t])\)*nc[t], \[IndentingNewLine]\(nh'\)[ t] \[Equal] \(-Ar\)*nh[t]*nc[t] + Ah*nv[t]*\((Nh - nh[t])\), \ \[IndentingNewLine]\(nv'\)[ t] \[Equal] R - Ah*nv[t]*\((Nh - nh[t])\), \ \(nc'\)[ t] \[Equal] \(nh'\)[t] + \(nv'\)[t] - \(n1'\)[t], n1[0] \[Equal] n10, nc[0] \[Equal] nc0, \ nv[0] \[Equal] nv0, nh[0] \[Equal] nh0}, {n1, nh, nv, nc}, {t, 0, tfinal}, MaxSteps \[Rule] 50000]]\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", \(initValues[b_, d_] := Module[{}, n10 = Last[n1[d] /. b]; nh0 = Last[nh[d] /. b]; nv0 = Last[nv[d] /. b]; nc0 = Last[nc[d] /. b]]\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", StyleBox[\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(---\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\), FontColor->RGBColor[0, 1, 0]], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", \( (*\ irradiation\ \ *) \), "\[IndentingNewLine]", \(R = 10^14\), ";", \(tfinal9 = irrTime\), ";", "\[IndentingNewLine]", \(solveDiffeq[n10, nh0, nv0, nc0, R, tfinal9]\), ";", "\[IndentingNewLine]", \(sol9 = sol\), ";", "\[IndentingNewLine]", "\[IndentingNewLine]", \( (*\ Relaxation\ stage\ *) \), "\[IndentingNewLine]", \(initValues[sol9, tfinal9]\), ";", "\[IndentingNewLine]", \(R = 0\), ";", \(tfinal10 = 60\), ";", "\[IndentingNewLine]", \(solveDiffeq[n10, nh0, nv0, nc0, R, tfinal10]\), ";", "\[IndentingNewLine]", \(sol10 = sol\), ";"}], "\[IndentingNewLine]", ")"}]}], "\[IndentingNewLine]", RowBox[{"(*", StyleBox[\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(-\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(---\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\), FontColor->RGBColor[0, 1, 0]], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", \(tstart = 0; tend = 20; tstep = 1;\), "\[IndentingNewLine]", \(n1List1 = {}; n1List2 = {};\), "\[IndentingNewLine]", \(For[tloop = tstart, tloop \[LessEqual] tend, tloop += tstep, \[IndentingNewLine]n10 = 0; nh0 = 0; nv0 = 0; nc0 = 0; \[IndentingNewLine]irrTime = tloop; \[IndentingNewLine]programMain; \[IndentingNewLine]initValues[ sol9, tfinal9]; (*find\ n \((t)\) at\ the\ end\ of\ irradiation\ stage*) ; \ \[IndentingNewLine]AppendTo[ n1List1, {tloop, n10}]; \[IndentingNewLine]initValues[sol10, tfinal10]; (*find\ n \((t + T)\)\ at\ the\ end\ of\ relaxation\ stage*) ; \ \[IndentingNewLine]AppendTo[ n1List2, {tloop, n10}];\[IndentingNewLine]]\), "\[IndentingNewLine]", \(<< Graphics`MultipleListPlot`\), "\[IndentingNewLine]", \(gr1 = ListPlot[n1List1, PlotRange \[Rule] All, PlotJoined \[Rule] True, PlotLabel \[Rule] "\", ImageSize \[Rule] 72\ 3];\), "\[IndentingNewLine]", \(\(gr2 = ListPlot[n1List2, PlotRange \[Rule] All, PlotJoined \[Rule] True, PlotLabel \[Rule] "\", ImageSize \[Rule] 72\ 3]\)\(;\)\(Show[{gr1, gr2}, PlotLabel \[Rule] "\"]\)\(;\)\(\ \[IndentingNewLine]\)\)}], "Input"], 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.0238095 0.047619 0.0147151 1.15897e-014 [ [.2619 .00222 -3 -9 ] [.2619 .00222 3 0 ] [.5 .00222 -6 -9 ] [.5 .00222 6 0 ] [.7381 .00222 -6 -9 ] [.7381 .00222 6 0 ] [.97619 .00222 -6 -9 ] [.97619 .00222 6 0 ] [.01131 .13061 -33 -6.25 ] [.01131 .13061 0 6.25 ] [.01131 .24651 -33 -6.25 ] [.01131 .24651 0 6.25 ] [.01131 .36241 -33 -6.25 ] [.01131 .36241 0 6.25 ] [.01131 .4783 -33 -6.25 ] [.01131 .4783 0 6.25 ] [.01131 .5942 -33 -6.25 ] [.01131 .5942 0 6.25 ] [.5 .63053 -32 0 ] [.5 .63053 32 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 .2619 .01472 m .2619 .02097 L s [(5)] .2619 .00222 0 1 Mshowa .5 .01472 m .5 .02097 L s [(10)] .5 .00222 0 1 Mshowa .7381 .01472 m .7381 .02097 L s [(15)] .7381 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(20)] .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 .21429 .01472 m .21429 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .40476 .01472 m .40476 .01847 L s .45238 .01472 m .45238 .01847 L s .54762 .01472 m .54762 .01847 L s .59524 .01472 m .59524 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .78571 .01472 m .78571 .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 .13061 m .03006 .13061 L s gsave .01131 .13061 -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 (13) show grestore .02381 .24651 m .03006 .24651 L s gsave .01131 .24651 -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 (13) show grestore .02381 .36241 m .03006 .36241 L s gsave .01131 .36241 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (3) 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 (13) show grestore .02381 .4783 m .03006 .4783 L s gsave .01131 .4783 -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 (13) show grestore .02381 .5942 m .03006 .5942 L s gsave .01131 .5942 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (5) 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 (13) show grestore .125 Mabswid .02381 .03789 m .02756 .03789 L s .02381 .06107 m .02756 .06107 L s .02381 .08425 m .02756 .08425 L s .02381 .10743 m .02756 .10743 L s .02381 .15379 m .02756 .15379 L s .02381 .17697 m .02756 .17697 L s .02381 .20015 m .02756 .20015 L s .02381 .22333 m .02756 .22333 L s .02381 .26969 m .02756 .26969 L s .02381 .29287 m .02756 .29287 L s .02381 .31605 m .02756 .31605 L s .02381 .33923 m .02756 .33923 L s .02381 .38558 m .02756 .38558 L s .02381 .40876 m .02756 .40876 L s .02381 .43194 m .02756 .43194 L s .02381 .45512 m .02756 .45512 L s .02381 .50148 m .02756 .50148 L s .02381 .52466 m .02756 .52466 L s .02381 .54784 m .02756 .54784 L s .02381 .57102 m .02756 .57102 L s .02381 .61738 m .02756 .61738 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .5 .63053 -93 -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 (n) 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 69.000 12.813 moveto (H) show 75.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 81.000 12.813 moveto (L) show 93.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 111.000 12.813 moveto (t) show 123.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 .07143 .02016 L .11905 .03429 L .16667 .05467 L .21429 .07955 L .2619 .10765 L .30952 .13801 L .35714 .16993 L .40476 .20288 L .45238 .23647 L .5 .27044 L .54762 .30455 L .59524 .33868 L .64286 .37269 L .69048 .40652 L .7381 .44011 L .78571 .47341 L .83333 .5064 L .88095 .53905 L .92857 .57136 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]HGAYHf4PAg9QL6QYHg7oo0P00MGoo000VOol00`00Oomoo`0jOol2001cOol0 02Ioo`03001oogoo03aoo`80075oo`009Woo00<007ooOol0?Woo0P00Kgoo000VOol00`00Oomoo`10 Ool00`00Oomoo`1/Ool002Ioo`03001oogoo045oo`03001oogoo06]oo`009Woo00<007ooOol0@Woo 0P00Jgoo000VOol00`00Oomoo`14Ool00`00Oomoo`1XOol002Ioo`03001oogoo04Eoo`03001oogoo 06Moo`000goo10000Woo00D007ooOomoo`0000=oo`D0009oo`8000aoo`03001oogoo04Ioo`8006Mo o`000goo00<007ooOol017oo00<007oo00001Woo00D007ooOomoo`00009oo`03001oogoo00Uoo`03 001oogoo04Qoo`03001oogoo06Aoo`0017oo00<007ooOol017oo00<007ooOol01Goo00D007ooOomo o`00009oo`03001oogoo00Uoo`03001oogoo04Uoo`8006Aoo`001Goo00<007ooOol00Woo00<007oo 00001Woo00D007ooOomoo`00009oo`03001oogoo00Uoo`8004aoo`03001oogoo065oo`000goo00@0 07ooOol0009oo`05001oogooOol00005Ool01@00Oomoogoo00000Woo00<007ooOol02Goo00<007oo Ool0C7oo0P00HGoo0004Ool2000Ool002Ioo`03001oogoo065oo`03 001oogoo04]oo`0017oo0P000goo00D007ooOomoo`0000=oo`D0009oo`8000aoo`03001oogoo069o o`8004]oo`000goo00@007ooOol000=oo`03001oo`0000Ioo`05001oogooOol00002Ool00`00Oomo o`09Ool00`00Oomoo`1TOol00`00Oomoo`18Ool000Eoo`03001oogoo00=oo`03001oogoo00Eoo`05 001oogooOol00002Ool00`00Oomoo`09Ool00`00Oomoo`1UOol20018Ool000Ioo`07001oogooOol0 07oo00001Woo00D007ooOomoo`00009oo`03001oogoo00Uoo`03001oogoo06Moo`03001oogoo04Eo o`000goo00@007ooOol0009oo`05001oogooOol00005Ool01@00Oomoogoo00000Woo00<007ooOol0 2Goo0P00JGoo0P00AGoo0004Ool2000Woo000VOol00`00Oomoo`1cOol2000j Ool002Ioo`03001oogoo07Eoo`03001oogoo03Moo`009Woo00<007ooOol0MWoo0P00=goo000VOol0 0`00Oomoo`1hOol00`00Oomoo`0dOol002Ioo`03001oogoo07Uoo`8003Aoo`009Woo00<007ooOol0 Ngoo00<007ooOol0Ool00`00Oomoo`0NOol002Ioo`03001oogoo08moo`8001io o`009Woo00<007ooOol0TGoo00<007ooOol06goo000VOol00`00Oomoo`2BOol2000KOol002Ioo`03 001oogoo09Aoo`03001oogoo01Qoo`009Woo00<007ooOol0UGoo0P0067oo000VOol00`00Oomoo`2G Ool00`00Oomoo`0EOol002Ioo`03001oogoo09Qoo`8001Eoo`009Woo00<007ooOol0VWoo0P004goo 000VOol00`00Oomoo`2LOol00`00Oomoo`0@Ool000Aoo`8000=oo`05001oogooOol00003Ool50002 Ool2000"], ImageRangeCache->{{{0, 215}, {132.438, 0}} -> {-4.53082, -5.47902*^12, \ 0.119005, 4.88959*^11}}], 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.047619 0.0147151 1.03328e-014 [ [.2619 .00222 -3 -9 ] [.2619 .00222 3 0 ] [.5 .00222 -6 -9 ] [.5 .00222 6 0 ] [.7381 .00222 -6 -9 ] [.7381 .00222 6 0 ] [.97619 .00222 -6 -9 ] [.97619 .00222 6 0 ] [.01131 .11804 -33 -6.25 ] [.01131 .11804 0 6.25 ] [.01131 .22137 -33 -6.25 ] [.01131 .22137 0 6.25 ] [.01131 .3247 -33 -6.25 ] [.01131 .3247 0 6.25 ] [.01131 .42803 -33 -6.25 ] [.01131 .42803 0 6.25 ] [.01131 .53135 -33 -6.25 ] [.01131 .53135 0 6.25 ] [.5 .63053 -38 0 ] [.5 .63053 38 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 .2619 .01472 m .2619 .02097 L s [(5)] .2619 .00222 0 1 Mshowa .5 .01472 m .5 .02097 L s [(10)] .5 .00222 0 1 Mshowa .7381 .01472 m .7381 .02097 L s [(15)] .7381 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(20)] .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 .21429 .01472 m .21429 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .40476 .01472 m .40476 .01847 L s .45238 .01472 m .45238 .01847 L s .54762 .01472 m .54762 .01847 L s .59524 .01472 m .59524 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .78571 .01472 m .78571 .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 .11804 m .03006 .11804 L s gsave .01131 .11804 -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 (13) show grestore .02381 .22137 m .03006 .22137 L s gsave .01131 .22137 -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 (13) show grestore .02381 .3247 m .03006 .3247 L s gsave .01131 .3247 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (3) 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 (13) show grestore .02381 .42803 m .03006 .42803 L s gsave .01131 .42803 -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 (13) show grestore .02381 .53135 m .03006 .53135 L s gsave .01131 .53135 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (5) 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 (13) show grestore .125 Mabswid .02381 .03538 m .02756 .03538 L s .02381 .05605 m .02756 .05605 L s .02381 .07671 m .02756 .07671 L s .02381 .09738 m .02756 .09738 L s .02381 .13871 m .02756 .13871 L s .02381 .15937 m .02756 .15937 L s .02381 .18004 m .02756 .18004 L s .02381 .2007 m .02756 .2007 L s .02381 .24204 m .02756 .24204 L s .02381 .2627 m .02756 .2627 L s .02381 .28337 m .02756 .28337 L s .02381 .30403 m .02756 .30403 L s .02381 .34536 m .02756 .34536 L s .02381 .36603 m .02756 .36603 L s .02381 .38669 m .02756 .38669 L s .02381 .40736 m .02756 .40736 L s .02381 .44869 m .02756 .44869 L s .02381 .46936 m .02756 .46936 L s .02381 .49002 m .02756 .49002 L s .02381 .51069 m .02756 .51069 L s .02381 .55202 m .02756 .55202 L s .02381 .57268 m .02756 .57268 L s .02381 .59335 m .02756 .59335 L s .02381 .61401 m .02756 .61401 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .5 .63053 -99 -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 (n) 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 69.000 12.813 moveto (H) show 75.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 81.000 12.813 moveto %%IncludeResource: font Mathematica1Mono %%IncludeFont: Mathematica1Mono /Mathematica1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (+) show 87.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 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 (t) show 135.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 .07143 .04986 L .11905 .08109 L .16667 .11207 L .21429 .14285 L .2619 .17343 L .30952 .20381 L .35714 .23398 L .40476 .26393 L .45238 .29365 L .5 .32312 L .54762 .35234 L .59524 .38131 L .64286 .41001 L .69048 .43845 L .7381 .46661 L .78571 .4945 L .83333 .52212 L .88095 .54946 L .92857 .57653 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]HGAYHf4PAg9QL6QYHgWoo0P00MGoo000NOol00`00Oomoo`02Ool01000Oomoo`00?7oo0P00Lgoo000NOol0 0`00Oomoo`03Ool00`00Ool0000nOol2001aOol001eoo`8000=oo`<00003Ool007oo03moo`8006mo o`009Woo00<007ooOol0@7oo00<007ooOol0K7oo000VOol00`00Oomoo`11Ool2001/Ool002Ioo`03 001oogoo04=oo`03001oogoo06Uoo`009Woo00<007ooOol0A7oo0P00JGoo000VOol00`00Oomoo`16 Ool2001WOol002Ioo`03001oogoo04Qoo`03001oogoo06Aoo`009Woo00<007ooOol0BGoo0P00I7oo 000VOol00`00Oomoo`1;Ool00`00Oomoo`1QOol000Aoo`8000=oo`05001oogooOol00003Ool50002 Ool2000Woo000NOol00`00Oomoo`02Ool0 1000Oomoo`00Mgoo0P00>7oo000NOol00`00Oomoo`03Ool00`00Ool0001iOol2000fOol001eoo`80 00=oo`<00003Ool007oo07Yoo`8003Aoo`009Woo00<007ooOol0Ngoo00<007ooOol0Ool002Ioo`03001oogoo0:5oo`8000aoo`00 9Woo00<007ooOol0Xgoo0P002Woo000VOol00`00Oomoo`2UOol00`00Oomoo`07Ool002Ioo`03001o ogoo0:moo`009Woo00<007ooOol0[goo000VOol00`00Oomoo`2_Ool00=Qoo`00f7oo003HOol00=Qo o`00H7oo00<007ooOol04Woo00<007ooOol0H7oo001OOol00`00Oomoo`09Ool00`00Oomoo`08Ool0 0`00Oomoo`1OOol005Ioo`80009oo`80009oo`03001oogoo00=oo`8000Eoo`04001oogooOol30005 Ool00`00Oomoo`08Ool20003Ool3000:Ool20012Ool005Moo`04001oogoo0003Ool00`00Oomoo`02 Ool01000Oomoo`000Woo1@000Woo00<007ooOol017oo00<007ooOol027oo0P001Woo00<007ooOol0 1Woo00@007ooOol0045oo`00Egoo00@007ooOol000=oo`03001oogoo009oo`03001oogoo00Eoo`03 001oogoo009oo`03001oogoo00Aoo`03001oogoo00Moo`04001oogoo0002Ool30009Ool00`00Oomo o`12Ool005Ioo`@000Aoo`04001oogooOol50004Ool00`00Oomoo`02Ool00`00Oomoo`04Ool00`00 Oomoo`06Ool20002Ool200000goo000000020007Ool50011Ool005moo`05001oogooOol0000:Ool0 1@00Ool007oo00000goo00<007ooOol06Woo00<007ooOol0@Woo001POol00`00Oomoo`0;Ool50002 Ool00`00Oomoo`1POol00=Qoo`00f7oo0000\ \>"], ImageRangeCache->{{{0, 215}, {132.438, 0}} -> {-4.53082, -6.1455*^12, \ 0.119005, 5.48437*^11}}], 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.047619 0.0147151 1.03328e-014 [ [.2619 .00222 -3 -9 ] [.2619 .00222 3 0 ] [.5 .00222 -6 -9 ] [.5 .00222 6 0 ] [.7381 .00222 -6 -9 ] [.7381 .00222 6 0 ] [.97619 .00222 -6 -9 ] [.97619 .00222 6 0 ] [.01131 .11804 -33 -6.25 ] [.01131 .11804 0 6.25 ] [.01131 .22137 -33 -6.25 ] [.01131 .22137 0 6.25 ] [.01131 .3247 -33 -6.25 ] [.01131 .3247 0 6.25 ] [.01131 .42803 -33 -6.25 ] [.01131 .42803 0 6.25 ] [.01131 .53135 -33 -6.25 ] [.01131 .53135 0 6.25 ] [.5 .63053 -74 0 ] [.5 .63053 74 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 .2619 .01472 m .2619 .02097 L s [(5)] .2619 .00222 0 1 Mshowa .5 .01472 m .5 .02097 L s [(10)] .5 .00222 0 1 Mshowa .7381 .01472 m .7381 .02097 L s [(15)] .7381 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(20)] .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 .21429 .01472 m .21429 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .40476 .01472 m .40476 .01847 L s .45238 .01472 m .45238 .01847 L s .54762 .01472 m .54762 .01847 L s .59524 .01472 m .59524 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .78571 .01472 m .78571 .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 .11804 m .03006 .11804 L s gsave .01131 .11804 -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 (13) show grestore .02381 .22137 m .03006 .22137 L s gsave .01131 .22137 -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 (13) show grestore .02381 .3247 m .03006 .3247 L s gsave .01131 .3247 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (3) 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 (13) show grestore .02381 .42803 m .03006 .42803 L s gsave .01131 .42803 -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 (13) show grestore .02381 .53135 m .03006 .53135 L s gsave .01131 .53135 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (5) 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 (13) show grestore .125 Mabswid .02381 .03538 m .02756 .03538 L s .02381 .05605 m .02756 .05605 L s .02381 .07671 m .02756 .07671 L s .02381 .09738 m .02756 .09738 L s .02381 .13871 m .02756 .13871 L s .02381 .15937 m .02756 .15937 L s .02381 .18004 m .02756 .18004 L s .02381 .2007 m .02756 .2007 L s .02381 .24204 m .02756 .24204 L s .02381 .2627 m .02756 .2627 L s .02381 .28337 m .02756 .28337 L s .02381 .30403 m .02756 .30403 L s .02381 .34536 m .02756 .34536 L s .02381 .36603 m .02756 .36603 L s .02381 .38669 m .02756 .38669 L s .02381 .40736 m .02756 .40736 L s .02381 .44869 m .02756 .44869 L s .02381 .46936 m .02756 .46936 L s .02381 .49002 m .02756 .49002 L s .02381 .51069 m .02756 .51069 L s .02381 .55202 m .02756 .55202 L s .02381 .57268 m .02756 .57268 L s .02381 .59335 m .02756 .59335 L s .02381 .61401 m .02756 .61401 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .5 .63053 -135 -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 (n) 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 69.000 12.813 moveto (H) show 75.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 81.000 12.813 moveto %%IncludeResource: font Mathematica1Mono %%IncludeFont: Mathematica1Mono /Mathematica1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (+) show 87.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 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 (and) show 129.000 12.813 moveto (n) 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 159.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 177.000 12.813 moveto (time) show 207.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 .07143 .01957 L .11905 .03217 L .16667 .05034 L .21429 .07252 L .2619 .09757 L .30952 .12464 L .35714 .15309 L .40476 .18247 L .45238 .21242 L .5 .2427 L .54762 .27312 L .59524 .30354 L .64286 .33387 L .69048 .36403 L .7381 .39397 L .78571 .42366 L .83333 .45307 L .88095 .48219 L .92857 .51099 L .97619 .53948 L s .02381 .01472 m .07143 .04986 L .11905 .08109 L .16667 .11207 L .21429 .14285 L .2619 .17343 L .30952 .20381 L .35714 .23398 L .40476 .26393 L .45238 .29365 L .5 .32312 L .54762 .35234 L .59524 .38131 L .64286 .41001 L .69048 .43845 L .7381 .46661 L .78571 .4945 L .83333 .52212 L .88095 .54946 L .92857 .57653 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]HGAYHf4PAg9QL6QYHgOol3002BOol002Ioo`03001oogoo 00Yoo`03001oogoo011oo`80091oo`009Woo00<007ooOol02goo0P004Woo0P00SWoo000VOol00`00 Oomoo`0=Ool00`00Oomoo`0AOol2002Woo0P005goo00<007ooOol0Fgoo 000NOol00`00Oomoo`02Ool01000Oomoo`00?7oo0P005Woo00<007ooOol0FWoo000NOol00`00Oomo o`03Ool00`00Ool0000nOol2000EOol2001JOol001eoo`8000=oo`<00003Ool007oo03moo`8001Eo o`03001oogoo05Moo`009Woo00<007ooOol0@7oo00<007ooOol04goo00<007ooOol0EWoo000VOol0 0`00Oomoo`11Ool2000DOol2001FOol002Ioo`03001oogoo04=oo`03001oogoo01=oo`03001oogoo 05=oo`009Woo00<007ooOol0A7oo0P0057oo0P00Dgoo000VOol00`00Oomoo`16Ool2000DOol00`00 Oomoo`1@Ool002Ioo`03001oogoo04Qoo`03001oogoo019oo`80051oo`009Woo00<007ooOol0BGoo 0P0057oo0P00CWoo000VOol00`00Oomoo`1;Ool00`00Oomoo`0COol00`00Oomoo`1;Ool000Aoo`80 00=oo`05001oogooOol00003Ool50002Ool20007oo000VOol00`00Oomoo`1QOol00`00Oomoo`0COol2000fOol002Io o`03001oogoo069oo`8001Eoo`8003Aoo`009Woo00<007ooOol0I7oo0P005Goo00<007ooOol0Ool002Ioo`03001oogoo 0:5oo`8000aoo`009Woo00<007ooOol0Xgoo0P002Woo000VOol00`00Oomoo`2UOol00`00Oomoo`07 Ool002Ioo`03001oogoo0:moo`009Woo00<007ooOol0[goo000VOol00`00Oomoo`2_Ool00=Qoo`00 f7oo003HOol00=Qoo`00?7oo00<007ooOol04Woo00<007ooOol0:Woo00<007ooOol01Woo00<007oo Ool0CWoo000kOol00`00Oomoo`09Ool00`00Oomoo`08Ool00`00Oomoo`0XOol00`00Oomoo`08Ool0 0`00Oomoo`1=Ool0039oo`80009oo`80009oo`03001oogoo00=oo`8000Eoo`04001oogooOol30005 Ool00`00Oomoo`07Ool70002Ool20002Ool40006Ool20002Ool20002Ool00`00Oomoo`03Ool20005 Ool00`00Oomoo`08Ool20003Ool3000:Ool20004Ool300001goo000007oo001oo`00009oo`<001eo o`00"], ImageRangeCache->{{{0, 215}, {132.438, 0}} -> {-4.53082, -6.1455*^12, \ 0.119005, 5.48437*^11}}] }, Open ]] }, Open ]] }, FrontEndVersion->"5.0 for Microsoft Windows", ScreenRectangle->{{0, 800}, {0, 527}}, WindowSize->{792, 500}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, CellLabelAutoDelete->True ] (******************************************************************* 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, 4414, 76, 950, "Input"], Cell[CellGroupData[{ Cell[6215, 133, 14464, 477, 142, 6740, 378, "GraphicsData", "PostScript", \ "Graphics"], Cell[20682, 612, 14984, 500, 142, 7237, 401, "GraphicsData", "PostScript", \ "Graphics"], Cell[35669, 1114, 17843, 573, 142, 8395, 453, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)