(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.1' 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[ 80730, 2093]*) (*NotebookOutlinePosition[ 81403, 2116]*) (* CellTagsIndexPosition[ 81359, 2112]*) (*WindowFrame->Normal*) Notebook[{ Cell["Homework 1 \[LongDash] Solution", "Title"], Cell["\<\ 3.016 Mathematical Methods for Materials Scientists and Engineers S. M. Allen September 30, 2005\ \>", "Subsubtitle"], Cell[CellGroupData[{ Cell["Individual Exercise I-1: ", "Subsection"], Cell["\<\ First, recognize that you will need to use ListPlot to make the \ final plot of your data. ListPlot produces a plot from (x,y) pairs that are \ stored in a Table. The program you need to write has to produce output that \ will be written to a Table with the proper format. The function Random[Integer] is used here; it returns either 0 or 1 with \ probability 1/2. Quoting Charles Cantrell: \"All of the problem is contained in one line of \ input. First I create a table of data members that contains the trial number \ and the number of required steps to reach the endpoint. Then I plot the data \ using the ListPlot command. The While loop runs until the position is either \ 100 or -100, and the Table command is run 100 times--giving us 100 data \ points.\"\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(data = Table[{i, location = 0; timestep = 0; While[Abs[location] < 100, If[Random[Integer] \[Equal] 0, location = location + 1, location = location - 1]; timestep = timestep + 1]; timestep}, {i, 100}];\)\), "\n", \(ListPlot[data, AxesLabel \[Rule] {"\", "\"}]\)}], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.00952381 0.0147151 1.12846e-05 [ [.21429 .00222 -6 -9 ] [.21429 .00222 6 0 ] [.40476 .00222 -6 -9 ] [.40476 .00222 6 0 ] [.59524 .00222 -6 -9 ] [.59524 .00222 6 0 ] [.78571 .00222 -6 -9 ] [.78571 .00222 6 0 ] [.97619 .00222 -9 -9 ] [.97619 .00222 9 0 ] [1.025 .01472 0 -6 ] [1.025 .01472 76 6 ] [.01131 .12756 -30 -4.5 ] [.01131 .12756 0 4.5 ] [.01131 .24041 -30 -4.5 ] [.01131 .24041 0 4.5 ] [.01131 .35325 -30 -4.5 ] [.01131 .35325 0 4.5 ] [.01131 .4661 -30 -4.5 ] [.01131 .4661 0 4.5 ] [.01131 .57894 -30 -4.5 ] [.01131 .57894 0 4.5 ] [.02381 .64303 -29 0 ] [.02381 .64303 29 12 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21429 .01472 m .21429 .02097 L s [(20)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(40)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(60)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(80)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(100)] .97619 .00222 0 1 Mshowa .125 Mabswid .07143 .01472 m .07143 .01847 L s .11905 .01472 m .11905 .01847 L s .16667 .01472 m .16667 .01847 L s .2619 .01472 m .2619 .01847 L s .30952 .01472 m .30952 .01847 L s .35714 .01472 m .35714 .01847 L s .45238 .01472 m .45238 .01847 L s .5 .01472 m .5 .01847 L s .54762 .01472 m .54762 .01847 L s .64286 .01472 m .64286 .01847 L s .69048 .01472 m .69048 .01847 L s .7381 .01472 m .7381 .01847 L s .83333 .01472 m .83333 .01847 L s .88095 .01472 m .88095 .01847 L s .92857 .01472 m .92857 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s gsave 1.025 .01472 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 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 13.000 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 13.000 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 (Trial) show 99.000 13.000 moveto (number) show 135.000 13.000 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 .02381 .12756 m .03006 .12756 L s [(10000)] .01131 .12756 1 0 Mshowa .02381 .24041 m .03006 .24041 L s [(20000)] .01131 .24041 1 0 Mshowa .02381 .35325 m .03006 .35325 L s [(30000)] .01131 .35325 1 0 Mshowa .02381 .4661 m .03006 .4661 L s [(40000)] .01131 .4661 1 0 Mshowa .02381 .57894 m .03006 .57894 L s [(50000)] .01131 .57894 1 0 Mshowa .125 Mabswid .02381 .03728 m .02756 .03728 L s .02381 .05985 m .02756 .05985 L s .02381 .08242 m .02756 .08242 L s .02381 .10499 m .02756 .10499 L s .02381 .15013 m .02756 .15013 L s .02381 .1727 m .02756 .1727 L s .02381 .19527 m .02756 .19527 L s .02381 .21784 m .02756 .21784 L s .02381 .26298 m .02756 .26298 L s .02381 .28555 m .02756 .28555 L s .02381 .30811 m .02756 .30811 L s .02381 .33068 m .02756 .33068 L s .02381 .37582 m .02756 .37582 L s .02381 .39839 m .02756 .39839 L s .02381 .42096 m .02756 .42096 L s .02381 .44353 m .02756 .44353 L s .02381 .48867 m .02756 .48867 L s .02381 .51124 m .02756 .51124 L s .02381 .53381 m .02756 .53381 L s .02381 .55638 m .02756 .55638 L s .02381 .60151 m .02756 .60151 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .02381 .64303 -90 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 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 13.000 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 13.000 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 (TimeSteps) show 117.000 13.000 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 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .008 w .03333 .22858 Mdot .04286 .23106 Mdot .05238 .10479 Mdot .0619 .0685 Mdot .07143 .04816 Mdot .08095 .04085 Mdot .09048 .21247 Mdot .1 .08396 Mdot .10952 .07418 Mdot .11905 .07599 Mdot .12857 .1523 Mdot .1381 .03182 Mdot .14762 .0996 Mdot .15714 .0629 Mdot .16667 .08021 Mdot .17619 .10249 Mdot .18571 .0636 Mdot .19524 .25634 Mdot .20476 .2809 Mdot .21429 .13235 Mdot .22381 .18055 Mdot .23333 .60332 Mdot .24286 .1265 Mdot .25238 .10276 Mdot .2619 .31252 Mdot .27143 .27491 Mdot .28095 .05525 Mdot .29048 .36264 Mdot .3 .09829 Mdot .30952 .05466 Mdot .31905 .08098 Mdot .32857 .29484 Mdot .3381 .07048 Mdot .34762 .05942 Mdot .35714 .13821 Mdot .36667 .23109 Mdot .37619 .04437 Mdot .38571 .08992 Mdot .39524 .04886 Mdot .40476 .06936 Mdot .41429 .08475 Mdot .42381 .02726 Mdot .43333 .03979 Mdot .44286 .06617 Mdot .45238 .09567 Mdot .4619 .36025 Mdot .47143 .25745 Mdot .48095 .22644 Mdot .49048 .39453 Mdot .5 .46576 Mdot .50952 .20136 Mdot .51905 .07488 Mdot .52857 .03421 Mdot .5381 .12467 Mdot .54762 .03679 Mdot .55714 .06699 Mdot .56667 .34529 Mdot .57619 .11023 Mdot .58571 .05633 Mdot .59524 .14356 Mdot .60476 .15094 Mdot .61429 .09382 Mdot .62381 .03965 Mdot .63333 .13625 Mdot .64286 .05951 Mdot .65238 .28654 Mdot .6619 .19233 Mdot .67143 .11361 Mdot .68095 .35233 Mdot .69048 .12321 Mdot .7 .03178 Mdot .70952 .13404 Mdot .71905 .19213 Mdot .72857 .0429 Mdot .7381 .06107 Mdot .74762 .12973 Mdot .75714 .11014 Mdot .76667 .21256 Mdot .77619 .0404 Mdot .78571 .16807 Mdot .79524 .13063 Mdot .80476 .11501 Mdot .81429 .12544 Mdot .82381 .20168 Mdot .83333 .04153 Mdot .84286 .03598 Mdot .85238 .07358 Mdot .8619 .04317 Mdot .87143 .15214 Mdot .88095 .068 Mdot .89048 .299 Mdot .9 .16182 Mdot .90952 .30236 Mdot .91905 .54211 Mdot .92857 .0494 Mdot .9381 .14374 Mdot .94762 .03291 Mdot .95714 .10876 Mdot .96667 .04929 Mdot .97619 .14002 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{472.438, 291.938}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOooool5o`0000;ooooo00Co0000oooooooooooo0000>_ooool01Ol0003ooooooooooooooooo 000000;ooooo00Co0000oooooooooooo0000>_ooool01Ol0003ooooooooooooooooo000000;ooooo 00Co0000oooooooooooo0000>?ooool00ol0003oooooooooo`02ooooo`04o`000?ooooooooooo`00 00;ooooo00Co0000oooooooooooo0000Eoooool006Cooooo00Go0000ooooooooooooooooo`000002 ooooo`03o`000?oooooooooo03Oooooo00Co0000oooooooooooo00000oooool01?l0003ooooooooo ool0000jooooo`Co00000oooool01?l0003oooooooooool0000kooooo`?o00000oooool01?l0003o ooooooooool0000hooooo`03o`000?oooooooooo00;ooooo00Co0000oooooooooooo00000_ooool0 1?l0003oooooooooool0001Gooooo`00HOooool01Ol0003ooooooooooooooooo000000;ooooo00Co 0000oooooooooooo0000>_ooool00ol0003oooooo`000003ooooo`04o`000?ooooooooooo`0003[o oooo00?o0000ooooooooool01?ooool01?l0003oooooooooool0000jooooo`05o`000?oooooooooo ooooool000000_ooool01?l0003oooooooooool0000hooooo`03o`000?oooooooooo00;ooooo00Co 0000oooooooooooo00000_ooool01?l0003oooooooooool0001Gooooo`00HOooool01Ol0003ooooo oooooooooooo000000;ooooo00Co0000oooooooooooo0000>oooool2o`0000?ooooo00Co0000oooo oooooooo0000>oooool00ol0003oooooooooo`03ooooo`04o`000?ooooooooooo`0003[ooooo00Go 0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0003Oooooo0_l00004ooooo`04 o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo0000Eoooool006;ooooo0ol00004oooo o`;o0000?Oooool00ol0003oooooooooo`02ooooo`;o0000?Oooool2o`0000Cooooo0_l0000loooo o`?o00001?ooool2o`0003Wooooo00?o0000ooooooooool00oooool2o`0000Cooooo0_l0001Hoooo o`008_ooool00ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo ]?ooool002;ooooo00?o0000ooooooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3o oooookCooooo000Rooooo`03o`000?oooooooooo0?ooooooJ_ooool3o`0000;ooooo0ol00003oooo o`Go00000_ooool4o`000003ooooool0003o000000?o00001oooool3o`000003ooooool0003o0000 00;ooooo0_l000002?oooooo0000o`000?oooooo0000ooooool0003ooooo1?l00003ooooo`?o0000 0_ooool3o`0000Kooooo000Rooooo`03o`000?oooooooooo0?ooooooJoooool00ol0003ooooooooo o`02ooooo`03o`000?oooooooooo00Cooooo00Go0000ooooooooooooooooo`000002ooooo`03o`00 0?oooooooooo00;ooooo00?o0000ooooooooool02?ooool01?l0003oooooooooool00002ooooo`0: o`000?ooooooooooo`000?oooooo0000ooooool0003oooooo`0000;ooooo00Ko0000oooooooooooo 0000ooooool00006ooooo`03o`000?oooooooooo00Gooooo000Ioooooooo0000Jol0000ooooo`03 o`000?oooooooooo00kooooo00?o0000ooooooooool03oooool00ol0003oooooooooo`0>ooooo`03 o`000?oooooooooo00kooooo00?o0000ooooooooool03_ooool00ol0003oooooooooo`0>ooooo`03 o`000?oooooooooo00kooooo00?o0000ooooooooool03oooool00ol0003oooooooooo`0>ooooo`03 o`000?oooooooooo00kooooo00?o0000ooooooooool03_ooool00ol0003oooooooooo`0>ooooo`03 o`000?oooooooooo00kooooo00?o0000ooooooooool03oooool00ol0003oooooooooo`0>ooooo`03 o`000?oooooooooo00kooooo00?o0000ooooooooool03_ooool00ol0003oooooooooo`0>ooooo`03 o`000?oooooooooo00kooooo00?o0000ooooooooool04oooool00ol0003oooooooooo`02ooooo`;o 00001Oooool01Ol0003ooooooooooooooooo000000;ooooo00?o0000ooooooooool00_ooool00ol0 003oooooooooo`08ooooo`04o`000?ooooooooooo`0000;ooooo00[o0000oooooooooooo0000oooo ool0003oooooo`000?oooooo00000_ooool01_l0003oooooooooool0003oooooo`0000?ooooo00?o 0000ooooooooool00_l00006ooooo`008_ooool00ol0003oooooooooo`12ooooo`03o`000?oooooo oooo047ooooo00?o0000ooooooooool0@_ooool00ol0003oooooooooo`12ooooo`03o`000?oooooo oooo047ooooo00?o0000ooooooooool04Oooool01_l0003oooooo`000?oooooo0000ooooo`;o0000 00Coooooo`000?l0003ooooo0ol00004ooooo`;o00001Oooool00ol0003oooooooooo`07ooooo`Co 00000_ooool2o`000004ooooool0003o0000ooooo`;o000000?oooooo`000?ooool00_ooool3o`00 00?ooooo0ol00002ooooo`;o000000?oooooo`000?l000001?ooool002;ooooo00?o0000oooooooo ool0oooooomYooooo`05o`000?oooooo0000ooooool000005Oooool00ol0003oooooooooo`0Joooo o`03o`000?oooooooooo017ooooo000Rooooo`03o`000?oooooooooo08cooooo0_l0003Jooooo`Go 00002Oooool00ol0003oooooooooo`07ooooo`?o00006oooool2o`0001?ooooo000Rooooo`03o`00 0?oooooooooo08cooooo0_l0003oooooobKooooo000Rooooo`03o`000?oooooooooo02Gooooo0_l0 002;ooooo`;o0000??ooool2o`0005Oooooo0_l0001Xooooo`008_ooool00ol0003oooooooooo`0U ooooo`;o0000Roooool2o`0000Gooooo0_l0000eooooo`;o0000ooooo`008_ooool00ol0003oooooooooo`2Toooo o`;o0000ooooool>ooooo`008_ooool2o`000?oooooo]Oooool002;ooooo00?o0000ooooooooool0 oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooo oooo0?oooooo]?ooool002;ooooo00?o0000ooooooooool0oooooondooooo`008_ooool00ol0003o ooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o 0000ooooooooool0oooooondooooo`008_ooool2o`000?oooooo]Oooool002;ooooo00?o0000oooo ooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`00 0?oooooooooo0?oooooo]?ooool002;ooooo00?o0000ooooooooool0oooooondooooo`008_ooool0 0ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;o oooo00?o0000ooooooooool0oooooondooooo`008_ooool2o`000?oooooo]Oooool002;ooooo00?o 0000ooooooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Roooo o`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o0000ooooooooool0oooooondooooo`00 8_ooool00ol0003oooooooooo`3ooooookCooooo0002ooooo`?o00000oooool2o`0000Cooooo0_l0 0004ooooo`;o00001?ooool2o`0000Kooooo00?o0000ooooooooool0oooooondooooo`000oooool0 1Ol0003ooooooooooooooooo000000;ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooo ooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00?o0000ooooooooool00oooool00ol0 003oooooooooo`3ooooookCooooo0005o`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0 003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo0000 1Oooool3o`000:Sooooo0_l0003oooooo`[ooooo00001?l0003oooooooooool00003ooooo`04o`00 0?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002 ooooo`04o`000?ooooooooooo`0000Gooooo00?o0000ooooooooool0Z?ooool2o`000?oooooo2_oo ool00004ooooool0003oooooo`0000?ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooo ooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00001Oooool0 0ol0003oooooooooo`3ooooookCooooo0002ooooo`;o00000oooool01?l0003oooooooooool00002 ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooo ooooool00005ooooo`03o`000?oooooooooo0?oooooo]?ooool000?ooooo00?o0000ooooooooool0 0_ooool2o`0000Cooooo0_l00004ooooo`;o00001?ooool2o`0000Kooooo00?o0000ooooooooool0 oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooo oooo0?oooooo]?ooool002;ooooo00?o0000ooooooooool0oooooondooooo`008_ooool2o`000?oo oooo]Oooool002;ooooo00?o0000ooooooooool0oooooondooooo`008_ooool00ol0003ooooooooo o`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o0000oooo ooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`00 0?oooooooooo0?oooooo]?ooool002;ooooo00?o0000ooooooooool0oooooondooooo`008_ooool0 0ol0003oooooooooo`3ooooookCooooo000Rooooo`;o0000ooooooneooooo`008_ooool00ol0003o ooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o 0000ooooooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Roooo o`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o0000ooooooooool0oooooondooooo`00 8_ooool00ol0003oooooooooo`3ooooookCooooo000Rooooo`;o0000ooooooneooooo`008_ooool0 0ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;o oooo00?o0000ooooooooool0oooooom0ooooo`;o0000L_ooool002;ooooo00?o0000ooooooooool0 oooooom0ooooo`;o0000L_ooool002;ooooo00?o0000ooooooooool0oooooondooooo`008_ooool0 0ol0003oooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;o oooo0_l0003ooooookGooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o 0000ooooooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Roooo o`03o`000?oooooooooo0?oooooo]?ooool00004ooooool0003o0000o`0000Cooooo0_l00004oooo o`;o00001?ooool2o`0000Cooooo0_l00006ooooo`03o`000?oooooooooo0?oooooo]?ooool00005 o`000?ooooooooooooooool000000_ooool01?l0003oooooooooool00002ooooo`04o`000?oooooo ooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00005ooooo`03 o`000?oooooooooo0?oooooo]?ooool000Cooooo00Co0000oooooooooooo00000_ooool01?l0003o ooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_oo ool00ol0003oooooooooo`03ooooo`03o`000?oooooooooo0?oooooo]?ooool000Cooooo00Co0000 oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;o oooo00Co0000oooooooooooo00000_ooool00ol0003oooooooooo`03ooooo`?o0000oooooondoooo o`001?l00003ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool0 1?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000Gooooo00?o0000ooooooooool0 oooooondooooo`0000?o0000ooooooooool01?ooool01?l0003oooooooooool00002ooooo`04o`00 0?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00005 ooooo`03o`000?oooooooooo0?oooooo]?ooool000Go00000oooool2o`0000Cooooo0_l00004oooo o`;o00001?ooool2o`0000Kooooo00?o0000ooooooooool0oooooondooooo`008_ooool00ol0003o ooooooooo`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o 0000ooooooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo000Roooo o`;o0000ooooooneooooo`008_ooool00ol0003oooooooooo`18ooooo`;o0000oooooomZooooo`00 8_ooool00ol0003oooooooooo`18ooooo`;o0000oooooomZooooo`008_ooool00ol0003ooooooooo o`3ooooookCooooo000Rooooo`03o`000?oooooooooo0?oooooo]?ooool002;ooooo00?o0000oooo ooooool0oooooondooooo`008_ooool00ol0003oooooooooo`3ooooookCooooo003oooooomWooooo 003oooooomWooooo003oooooomWooooo003oooooomWooooo003oooooomWooooo003oooooomWooooo 003oooooomWooooo003oooooomWooooo003oooooomWooooo003oooooomWooooo000aooooo`?o0000 oooooonUooooo`00<_ooool00ol0003oooooooooo`3oooooojCooooo0008ooooo`?o00000_ooool5 o`000006ooooool0003oooooo`000?oooooo00000_ooool3o`0000?ooooo0ol00004ooooo`;o0000 0oooool3o`0000?ooooo0ol00002ooooo`Co0000oooooonNooooo`002Oooool00ol0003ooooooooo o`03ooooo`0;o`000?ooooooooooooooool0003oooooo`000?oooooo0000ooooool000001_ooool0 1?l0003oooooooooool00002ooooo`06o`000?ooooooooooo`000?oooooo00001_ooool01?l0003o ooooooooool00005ooooo`03o`000?oooooooooo0?ooooooVoooool000Wooooo00?o0000oooooooo ool00oooool02_l0003ooooooooooooooooo0000ooooool0003oooooo`000?ooool5o`0000Gooooo 00Co0000oooooooooooo00001?ooool5o`0000;ooooo00Co0000oooooooooooo00000_ooool3o`00 0?ooooooW_ooool000Wooooo00?o0000ooooooooool00oooool02ol0003ooooooooooooooooo0000 ooooool0003oooooo`000?oooooo000000?ooooo00Co0000oooooooooooooooo0_l00003ooooo`03 o`000?oooooooooo00;ooooo00Go0000ooooooooooooooooo`000002ooooo`06o`000?oooooooooo o`000?oooooo0000oooooonQooooo`001oooool01_l0003oooooo`000?oooooo0000ooooo`?o0000 0oooool2o`000003ooooool0003ooooo00;ooooo0ol00003ooooo`03o`000?oooooooooo00;ooooo 1?l00003ooooo`?o00000_ooool4o`0000?ooooo1?l0003ooooooigooooo0007ooooo`05o`000?oo oooo0000ooooool000005?ooool01?l0003oooooooooool00002ooooo`03o`000?oooooooooo0?oo oooo/?ooool000Oooooo1Ol00003ooooo`03o`000?oooooooooo00oooooo0ol00002ooooo`03o`00 0?oooooooooo0?oooooo/?ooool00?oooooofOooool00?oooooofOooool00?oooooofOooool00?oo oooofOooool00?oooooofOooool00?oooooofOooool00?oooooofOooool00?oooooofOooool00?oo oooofOooool00?oooooofOooool00?oooooofOooool00?oooooofOooool00?oooooofOooool00?oo oooofOooool00?oooooofOooool00?oooooofOooool00?oooooofOooool00?oooooofOooool00?oo oooofOooool00?oooooofOooool00?oooooofOooool00001\ \>"], ImageRangeCache->{{{0, 471.438}, {290.938, 0}} -> {-10.0509, -8103.98, \ 0.291263, 245.816}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Individual Exercise I-2:", "Subsection"], Cell["\<\ Note: Some features of Omar Fabian's work were incorporated in the \ solution to this exercise.\ \>", "Text"], Cell[CellGroupData[{ Cell["1.", "Subsubsection"], Cell[TextData[StyleBox["Input the LJ potential, differentiate and set \ derivative == 0 to find minima.", "Text"]], "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(LJ = a\/r^12 - b\/r^6\)], "Input"], Cell[BoxData[ \(a\/r\^12 - b\/r\^6\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(MinimaOfLJ\ = \ Solve[D[LJ, r] == \ 0, r]\)], "Input"], Cell[BoxData[ \({{r \[Rule] \(-\(\(2\^\(1/6\)\ a\^\(1/6\)\)\/b\^\(1/6\)\)\)}, {r \ \[Rule] \(2\^\(1/6\)\ a\^\(1/6\)\)\/b\^\(1/6\)}, {r \[Rule] \(-\(\(\((\(-1\))\ \)\^\(1/3\)\ 2\^\(1/6\)\ a\^\(1/6\)\)\/b\^\(1/6\)\)\)}, {r \[Rule] \ \(\((\(-1\))\)\^\(1/3\)\ 2\^\(1/6\)\ a\^\(1/6\)\)\/b\^\(1/6\)}, {r \[Rule] \ \(-\(\(\((\(-1\))\)\^\(2/3\)\ 2\^\(1/6\)\ a\^\(1/6\)\)\/b\^\(1/6\)\)\)}, {r \ \[Rule] \(\((\(-1\))\)\^\(2/3\)\ 2\^\(1/6\)\ a\^\(1/6\)\)\/b\^\(1/6\)}}\)], \ "Output"] }, Open ]], Cell["\<\ Note that there are six roots but only one of them is \"physical\" \ \[LongDash] i.e., real and positive (the second root). Use this to set rmin. \ Note how Replace is used to extract the value of the second root.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Minr = r /. \ MinimaOfLJ[\([2]\)]\)], "Input"], Cell[BoxData[ \(\(2\^\(1/6\)\ a\^\(1/6\)\)\/b\^\(1/6\)\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["2. ", "Subsubsection"], Cell["Replace r with rmin in LJ to get EMin = LJ(rmin).", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(MinE = \ \ LJ /. r -> Minr\)], "Input"], Cell[BoxData[ \(\(-\(b\^2\/\(4\ a\)\)\)\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["3.", "Subsubsection"], Cell["\<\ The Solve function can be used to eliminate specific variables from \ a set of equations. Here we have the equation for rMin in terms of (a,b) and \ the equation for EMin in terms of (a,b). Solve is used to find expressions \ for a and b in terms of rMin and EMin.\ \>", "Text"], Cell[BoxData[ \(Clear[rMin, EMin]\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(NewVars = Solve[{rMin \[Equal] Minr, EMin \[Equal] MinE}, {a, b}]\)], "Input"], Cell[BoxData[ \({{a \[Rule] \(-EMin\)\ rMin\^12, b \[Rule] \(-2\)\ EMin\ rMin\^6}}\)], "Output"] }, Open ]], Cell["\<\ All that remains is to Replace a and b with these new \ variables.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(LJModified = LJ /. NewVars[\([1]\)]\)], "Input"], Cell[BoxData[ \(\(2\ EMin\ rMin\^6\)\/r\^6 - \(EMin\ rMin\^12\)\/r\^12\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["4.", "Subsubsection"], Cell["The force is obtained by differentiation", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(Force\ := \ \(-\[PartialD]\_r\ LJModified\)\), "\[IndentingNewLine]", \(Force\)}], "Input"], Cell[BoxData[ \(\(12\ EMin\ rMin\^6\)\/r\^7 - \(12\ EMin\ rMin\^12\)\/r\^13\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["5.", "Subsubsection"], Cell["\<\ Define new \"normalized\" variables rNorm = r / rMin and LJNorm = \ LJ / EMin and use them to define a normalized LJ function LJNorm(rNorm). \ Recongnizing that I will eventually want to plot these functions, I define \ them with delayed assignment :=\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(LJNorm := \ Simplify[\(-\((LJModified/EMin)\)\) /. rMin \[Rule] r/rNorm, Assumptions \[Rule] {a > 0\ && b > 0}]\), "\[IndentingNewLine]", \(LJNorm\)}], "Input"], Cell[BoxData[ \(1\/rNorm\^12 - 2\/rNorm\^6\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(ForceNorm := Simplify[\(-\((Force\ rMin/EMin)\)\) /. rMin \[Rule] r/rNorm, Assumptions \[Rule] {a > 0\ && b > 0}]\), "\[IndentingNewLine]", \(ForceNorm\)}], "Input"], Cell[BoxData[ \(\(-\(\(12\ \((\(-1\) + rNorm\^6)\)\)\/rNorm\^13\)\)\)], "Output"] }, Open ]], Cell["To demonstrate that LJNorm and ForceNorm are dimensionless:", "Text"], Cell[BoxData[{ \(\(EMinUnits = Mass\ Length^2/Time^2;\)\), "\[IndentingNewLine]", \(\(ForceUnits = Mass\ Length^2/\((Time^2\ Length)\);\)\), "\[IndentingNewLine]", \(\(rUnits\ = Length;\)\), "\[IndentingNewLine]", \(\(rMinUnits\ = \ Length;\)\)}], "Input"], Cell[CellGroupData[{ Cell[BoxData[{ \(\(rBarUnits\ = \ rUnits/rMinUnits;\)\), "\[IndentingNewLine]", \(rBarUnits\ \[Equal] 1\)}], "Input"], Cell[BoxData[ \(True\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ \(\(FBarUnits = ForceUnits\ rMinUnits/EMinUnits;\)\), "\[IndentingNewLine]", \(FBarUnits\ \[Equal] 1\)}], "Input"], Cell[BoxData[ \(True\)], "Output"] }, Open ]], Cell[TextData[{ "So, the variables ", Cell[BoxData[ \(TraditionalForm\`r\&_\)]], " and ", Cell[BoxData[ \(TraditionalForm\`F\&_\)]], " are dimensionless." }], "Text"] }, Open ]], Cell[CellGroupData[{ Cell["6.", "Subsubsection"], Cell["\<\ Now for the plot. Note in the use of the DisplayFunction option to \ control whether or not the graphics are rendered to the screen.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(NewPlot\ = \ Plot[{LJNorm, ForceNorm}, \ {rNorm, \ 0, \ 2}, \ PlotRange\ \[Rule] \ {\(-4\), 4}, \ \ \ PlotStyle\ \[Rule] \ {{Thickness[ .01], \ Hue[0.3]}, {Thickness[ .01], \ Hue[1]}}, AxesLabel -> {"\", "\"}, DisplayFunction \[Rule] Identity]\)], "Input"], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Show[NewPlot, Graphics[Text["\", {1.5\ , \(-2\)}]], Graphics[Text["\", { .6, 2}]], DisplayFunction \[Rule] $DisplayFunction]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.47619 0.309017 0.0772542 [ [.2619 .29652 -9 -9 ] [.2619 .29652 9 0 ] [.5 .29652 -3 -9 ] [.5 .29652 3 0 ] [.7381 .29652 -9 -9 ] [.7381 .29652 9 0 ] [.97619 .29652 -3 -9 ] [.97619 .29652 3 0 ] [1.025 .30902 0 -6 ] [1.025 .30902 40 6 ] [.01131 0 -12 -4.5 ] [.01131 0 0 4.5 ] [.01131 .07725 -12 -4.5 ] [.01131 .07725 0 4.5 ] [.01131 .15451 -12 -4.5 ] [.01131 .15451 0 4.5 ] [.01131 .23176 -12 -4.5 ] [.01131 .23176 0 4.5 ] [.01131 .38627 -6 -4.5 ] [.01131 .38627 0 4.5 ] [.01131 .46353 -6 -4.5 ] [.01131 .46353 0 4.5 ] [.01131 .54078 -6 -4.5 ] [.01131 .54078 0 4.5 ] [.01131 .61803 -6 -4.5 ] [.01131 .61803 0 4.5 ] [.02381 .64303 -59 0 ] [.02381 .64303 59 12 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .2619 .30902 m .2619 .31527 L s [(0.5)] .2619 .29652 0 1 Mshowa .5 .30902 m .5 .31527 L s [(1)] .5 .29652 0 1 Mshowa .7381 .30902 m .7381 .31527 L s [(1.5)] .7381 .29652 0 1 Mshowa .97619 .30902 m .97619 .31527 L s [(2)] .97619 .29652 0 1 Mshowa .125 Mabswid .07143 .30902 m .07143 .31277 L s .11905 .30902 m .11905 .31277 L s .16667 .30902 m .16667 .31277 L s .21429 .30902 m .21429 .31277 L s .30952 .30902 m .30952 .31277 L s .35714 .30902 m .35714 .31277 L s .40476 .30902 m .40476 .31277 L s .45238 .30902 m .45238 .31277 L s .54762 .30902 m .54762 .31277 L s .59524 .30902 m .59524 .31277 L s .64286 .30902 m .64286 .31277 L s .69048 .30902 m .69048 .31277 L s .78571 .30902 m .78571 .31277 L s .83333 .30902 m .83333 .31277 L s .88095 .30902 m .88095 .31277 L s .92857 .30902 m .92857 .31277 L s .25 Mabswid 0 .30902 m 1 .30902 L s gsave 1.025 .30902 -61 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 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 13.000 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 13.000 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 (r) 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 13.000 moveto (\\220) show 75.000 13.000 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 (rMin) show 99.000 13.000 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 .02381 0 m .03006 0 L s [(-4)] .01131 0 1 0 Mshowa .02381 .07725 m .03006 .07725 L s [(-3)] .01131 .07725 1 0 Mshowa .02381 .15451 m .03006 .15451 L s [(-2)] .01131 .15451 1 0 Mshowa .02381 .23176 m .03006 .23176 L s [(-1)] .01131 .23176 1 0 Mshowa .02381 .38627 m .03006 .38627 L s [(1)] .01131 .38627 1 0 Mshowa .02381 .46353 m .03006 .46353 L s [(2)] .01131 .46353 1 0 Mshowa .02381 .54078 m .03006 .54078 L s [(3)] .01131 .54078 1 0 Mshowa .02381 .61803 m .03006 .61803 L s [(4)] .01131 .61803 1 0 Mshowa .125 Mabswid .02381 .01545 m .02756 .01545 L s .02381 .0309 m .02756 .0309 L s .02381 .04635 m .02756 .04635 L s .02381 .0618 m .02756 .0618 L s .02381 .09271 m .02756 .09271 L s .02381 .10816 m .02756 .10816 L s .02381 .12361 m .02756 .12361 L s .02381 .13906 m .02756 .13906 L s .02381 .16996 m .02756 .16996 L s .02381 .18541 m .02756 .18541 L s .02381 .20086 m .02756 .20086 L s .02381 .21631 m .02756 .21631 L s .02381 .24721 m .02756 .24721 L s .02381 .26266 m .02756 .26266 L s .02381 .27812 m .02756 .27812 L s .02381 .29357 m .02756 .29357 L s .02381 .32447 m .02756 .32447 L s .02381 .33992 m .02756 .33992 L s .02381 .35537 m .02756 .35537 L s .02381 .37082 m .02756 .37082 L s .02381 .40172 m .02756 .40172 L s .02381 .41717 m .02756 .41717 L s .02381 .43262 m .02756 .43262 L s .02381 .44807 m .02756 .44807 L s .02381 .47898 m .02756 .47898 L s .02381 .49443 m .02756 .49443 L s .02381 .50988 m .02756 .50988 L s .02381 .52533 m .02756 .52533 L s .02381 .55623 m .02756 .55623 L s .02381 .57168 m .02756 .57168 L s .02381 .58713 m .02756 .58713 L s .02381 .60258 m .02756 .60258 L s .25 Mabswid .02381 0 m .02381 .61803 L s gsave .02381 .64303 -120 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 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 13.000 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 13.000 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 (Normalized) show 129.000 13.000 moveto (LJ) show 147.000 13.000 moveto (and) show 171.000 13.000 moveto (F) show 177.000 13.000 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 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .2 1 0 r .01 w .41563 .61803 m .42234 .51359 L .42754 .45307 L .43247 .40632 L .44174 .34078 L .44692 .31408 L .45172 .29421 L .457 .27692 L .46262 .26273 L .46797 .25254 L .4729 .24548 L .47831 .23981 L .48103 .23766 L .48402 .23575 L .4866 .23446 L .48905 .23349 L .49041 .23307 L .49188 .23268 L .4931 .23241 L .49445 .23217 L .49571 .232 L .4969 .23189 L .49816 .23181 L .49886 .23178 L .49951 .23177 L .5006 .23177 L .50179 .2318 L .50303 .23187 L .50419 .23197 L .50546 .2321 L .50681 .23228 L .50923 .23268 L .51471 .23391 L .52464 .23703 L .54322 .24451 L .5811 .26058 L .62142 .27453 L .66024 .28428 L .6814 .28834 L .7015 .29154 L .73973 .29622 L .75944 .29807 L .78042 .29971 L .82204 .30221 L .86215 .30391 L .90074 .30511 L .94178 .30604 L .97619 .30662 L s 1 0 0 r .48244 .61803 m .48336 .5919 L .48844 .48394 L .49397 .38958 L .49915 .31911 L .50402 .26602 L .50892 .22323 L .51341 .19184 L .5184 .16435 L .52379 .1419 L .52867 .12685 L .53402 .11514 L .53671 .11082 L .5396 .10722 L .54232 .10468 L .54359 .10376 L .54479 .10303 L .546 .10244 L .5471 .102 L .54832 .10163 L .54963 .10137 L .55089 .10124 L .55161 .10121 L .55227 .10122 L .55348 .10131 L .55477 .10151 L .55596 .10178 L .55725 .10216 L .55959 .10307 L .56397 .10546 L .56924 .10932 L .5748 .11428 L .58476 .12472 L .62277 .1699 L .66324 .21134 L .70218 .24048 L .72345 .25253 L .74358 .2619 L .78195 .27546 L .80176 .28074 L .82278 .28536 L .86057 .2917 L .87994 .2942 L .90081 .29645 L .94199 .29984 L .97619 .30189 L s 0 g gsave .7381 .15451 -78 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 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 13.000 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 13.000 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 (Force) show 93.000 13.000 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 gsave .30952 .46353 -90 -10 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 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 13.000 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 13.000 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 (Potential) show 117.000 13.000 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 % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{513.812, 317.5}, ImageMargins->{{0, 0}, {0, 6}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgooooo`00<_ooool00ol0003oooooooooo`3o ooooolgooooo000Yooooo`Go000000?oooooo`000?l000000ol0003oooooolkooooo000_ooooo`04 o`000?ooooooooooo`0000Sooooo0ol0003oooooolCooooo000`ooooo`03o`000?oooooo000000So oooo00?o0000ooooooooool0ooooooo4ooooo`00oooool00ol0003oooooooooo`3ooooo olCooooo000kooooo`;o0000ooooooo5ooooo`00>oooool00ol0003oooooooooo`3oooooolCooooo 000kooooo`03o`000?oooooooooo0?ooooooa?ooool003_ooooo00?o0000ooooooooool0ooooooo4 ooooo`00>oooool00ol0003oooooooooo`3oooooolCooooo000kooooo`03o`000?oooooooooo0?oo ooooa?ooool003_ooooo0_l0003oooooolGooooo000kooooo`03o`000?oooooooooo0?ooooooa?oo ool003_ooooo00?o0000ooooooooool0ooooooo4ooooo`00>oooool00ol0003oooooooooo`3ooooo olCooooo000kooooo`03o`000?oooooooooo0?ooooooa?ooool003_ooooo00?o0000ooooooooool0 ooooooo4ooooo`00>oooool00ol0003oooooooooo`3oooooolCooooo000kooooo`;o0000ooooooo5 ooooo`00>oooool00ol0003oooooooooo`3oooooolCooooo000kooooo`03o`000?oooooooooo0?oo ooooa?ooool003_ooooo00?o0000ooooooooool0ooooooo4ooooo`00>oooool00ol0003ooooooooo o`3oooooolCooooo000kooooo`03o`000?oooooooooo0?ooooooa?ooool003_ooooo0_l0003ooooo olGooooo000kooooo`03o`000?oooooooooo0?ooooooa?ooool003_ooooo00?o0000ooooooooool0 ooooooo4ooooo`00>oooool00ol0003oooooooooo`3oooooolCooooo000`ooooo`?o00002?ooool0 0ol0003oooooooooo`3oooooolCooooo000_ooooo`05o`000?ooooooooooooooool000001oooool0 0ol0003oooooooooo`3oooooolCooooo000Yooooo`Go00001Oooool00ol0003oooooooooo`05oooo o`03o`000?oooooooooo0?ooooooa?ooool0037ooooo0_l00008ooooo`?o0000ooooooo4ooooo`00 oooool00ol0003oooooooooo`3oooooolCooooo000kooooo`;o0000ooooooo5ooooo`00 >oooool00ol0003oooooooooo`3oooooolCooooo000kooooo`03o`000?oooooooooo0=Cooooo1Ooo 003Zooooo`00>oooool00ol0003oooooooooo`3Aooooo`coo`00i_ooool003_ooooo00?o0000oooo ooooool0dOooool=ool00>Gooooo000kooooo`03o`000?oooooooooo0oooool2o`000oooool00ol0003oooooooooo`3>ooooo`Ooo`001Ooo ool7ool00>;ooooo000kooooo`03o`000?oooooooooo0oooool00ol0003o ooooooooo`3oooool2o`000oooool00ol0003oooooooooo`3:ooooo`Ko o`004_ooool6ool00=_ooooo000kooooo`03o`000?oooooooooo0<[ooooo1Ooo000Cooooo`Ooo`00 f_ooool003_ooooo00?o0000ooooooooool0b_ooool5ool001Cooooo1ooo003Iooooo`00>oooool0 0ol0003oooooooooo`39ooooo`Goo`005_ooool7ool00=Sooooo000kooooo`;o0000b_ooool5ool0 01Oooooo1_oo003Hooooo`00>oooool00ol0003oooooooooo`39ooooo`Goo`006?ooool6ool00=Oo oooo000kooooo`03o`000?oooooooooo0oooool2o`000oooool00ol0003oooooooooo`36ooooo`Goo`009oooool7ool00<[ooooo000koooo o`03o`000?oooooooooo0oooool00ol0003oooooooooo`36ooooo`Coo`00 :oooool7ool00oooool0 0ol0003oooooooooo`35ooooo`Goo`00;Oooool7ool00oooool00ol0003oooooooooo`35ooooo`Coo`00oooool2o`000oooool00ol0003oooooooooo`34ooooo`Coo`00=oooool7ool00;gooooo000kooooo`03o`000?oo oooooooo0oooool00ol0003oooooooooo`34ooooo`Coo`00>_ooool8ool0 0;Wooooo000kooooo`;o0000a?ooool5ool003_ooooo2Ooo002gooooo`00>oooool00ol0003ooooo ooooo`33ooooo`Goo`00?Oooool8ool00;Kooooo000kooooo`03o`000?oooooooooo0ooo ool00ol0003oooooooooo`2dooooo`KoOooool:ool0 0:;ooooo000kooooo`03o`000?oooooooooo0;;ooooo1_lco`0:ooooo`Coo`003Oooool>oc?o03So oooo2_oo002Pooooo`00>oooool00ol0003oooooooooo`2aooooo`Kooc?o03Gooooo3?oo002Iooooo`00>oooool2o`000;7ooooo1Olco`0oc?o03Cooooo3Ooo002Fooooo`00>oooool00ol0003oooooooooo`2_ooooo`Ko oooool0 0ol0003oooooooooo`2^ooooo`Gooooool003_ooooo00?o0000 ooooooooool0Eoooool01?l0003oooooooooool0000;ooooo`03o`000?oooooooooo04Cooooo1Olc o`0>ooooo`Goo`00oooool00ol0003oooooooooo`1G ooooo`04o`000?ooooooooooo`0000_ooooo00?o0000ooooooooool0@oooool6oc?o00kooooo1Ooo 000doooooa?oOoooolFoc?o00[ooooo00?o 0000ooooooooool01oooool4o`00017ooooo6?oo000_ooooo`05o`000?ooooooooooooooool00000 >oooool003_ooooo00?o0000ooooooooool0Eoooool01?l0003oooooooooool00007ooooo`03o`00 0?oooooooooo04Oooooo1Olco`0>ooooo`03o`000?oo003oo`0000;oo`00?OoooolJoc?o0003oooo ool0003o000000Wooooo00?o0000ooooooooool05OoooolMool002Oooooo00Go0000oooooooooooo ooooo`00000kooooo`00>oooool00ol0003oooooooooo`1Hooooo`;o00002?ooool5o`0004Cooooo 1Olco`0@ooooo`Coo`00@OoooolMoc?o00Gooooo1Ol0000Hooooob7oo`007oooool3o`0003cooooo 000kooooo`;o0000[?ooool5oc?o013ooooo1?oo0015ooooob;ooooool00ol0003oooooooooo`2Zooooo`Gooooool003_ooooo00?o0000ooooooooool0Z_ooool5oc?o017ooooo1?oo001Koooood;o oooool00ol0003oooooooooo`0@ ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool04Oooool00ol0003oooooooooo`0@ ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool04Oooool00ol0003oooooooooo`0@ ooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool03?ooool5oc?o0003o`000?oooooo oooo00kooooo1Ooo000Aooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool04?ooool0 0ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool04?ooool0 0ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo:_lco`0Eooooo`;o00001_ooool0 1?l0003oooooooooooooool2o`0000?ooooo00Go0000ooooooooooooooooo`000003ooooo`03o`00 0?oooooooooo00;ooooo00Co0000oooooooooooo00001?ooool003_ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000ooooooooool0 4?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo00cooooo1Olco`000ol0003ooooo ooooo`0>ooooo`Coo`004_ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo 00?o0000ooooooooool04Oooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo 00?o0000ooooooooool04Oooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo013ooooo 00?o0000ooooooooool04Oooool00ol0003oooooooooo`0Cooooo`;o000000?oooooo`000?l00000 1?ooool00ol0003oooooooooo`02o`00000:ooooool0003o0000ooooool0003oooooo`000?oooooo 0000ooooo`?o00000oooool4o`0000Gooooo000kooooo`03o`000?oooooooooo05oooooo00?o0000 ooooooooool0Aoooool4oc?o01;ooooo1?oo001Qooooo`03o`000?oooooooooo05oooooo00?o0000 ooooooooool07Oooool00ol0003oooooooooo`05ooooo`;o000000?oooooo`000?l000004?ooool0 03_ooooo00?o0000ooooooooool0ZOooool4oc?o01;ooooo1?oo003[ooooo`04o`000?oooooooooo ooooo`;o00000_ooool00ol0003oooooooooo`0:ooooo`00>oooool00ol0003oooooooooo`2Xoooo o`Gooooool0 0ol0003oooooooooo`2Xooooo`Cooooool00ol0003oooooooooo`2Wooooo`Gooooool00ol0003o ooooooooo`2Vooooo`Gooooool2o`000:Oooooo1?lco`0Dooooo`Coo`00ooooool2 ooooo`00>oooool00ol0003oooooooooo`2Vooooo`Cooooool00ol0003ooooooooo o`2Uooooo`Gooooool2o`000:Kooooo1?lco`0Eooooo`Coo`00ooooool2ooooo`00 >oooool00ol0003oooooooooo`2Uooooo`Cooooool00ol0003oooooooooo`2Sooooo`Gooooool00ol0003o ooooooooo`2Sooooo`Cooooool00ol0003oooooooooo`2Sooooo`Cooooool00ol0003ooooooooo o`2Rooooo`Gooooool2o`000:?ooooo1?lco`0Fooooo`Coo`00ooooool4ooooo`00 >oooool00ol0003oooooooooo`2Rooooo`Cooooool00ol0003oooooooooo`2Qoooo o`Gooooool2o`000:;ooooo1?lco`0Gooooo`Coo`00ooooool4ooooo`00>oooool0 0ol0003oooooooooo`2Qooooo`Cooooool00ol0003oooooooooo`2Pooooo`Cooooool00ol0003ooooo ooooo`2Pooooo`Cooooool2o`000:7ooooo1?lco`0Gooooo`Coo`00ooooool5oooo o`00>oooool00ol0003oooooooooo`2Pooooo`Cooooool00ol0003oooooooooo`2O ooooo`Goooo ool00ol0003oooooooooo`2Oooooo`Cooooool00ol0003oooooooooo`2Oooooo`Co oooool00ol0 003oooooooooo`2Nooooo`Gooooool00ol0003oooooooooo`2Nooooo`Cooooool00ol0003ooooo ooooo`2Mooooo`Gooooool2o`0009kooooo1?lco`0Hooooo`Coo`00ooooool7oooo o`00>oooool00ol0003oooooooooo`2Mooooo`Cooooool00ol0003oooooooooo`2M ooooo`Coooo ool00ol0003oooooooooo`2Mooooo`Cooooool00ol0003oooooooooo`2Mooooo`Co oooool00ol0 003oooooooooo`2Looooo`Co"], ImageRangeCache->{{{0, 512.812}, {316.5, 0}} -> {-0.29978, -4.67687, \ 0.00508086, 0.0313181}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "7.", StyleBox["\n", "Subsection"] }], "Subsubsection"], Cell["\<\ I will use the same approach as above, using as my constants the \ equilibrium position, xMin, and minimum potential energy, \ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(PotEnergy\ = \ A\ + B\ x\ + \ k\ x^2/2\)], "Input"], Cell[BoxData[ \(A + B\ x - \(36\ EMin\ x\^2\)\/rMin\^2\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(SpringForceEq = \ \((D[PotEnergy, x] /. x -> \ xMin\ )\) == \ 0\)\(\[IndentingNewLine]\) \)\)], "Input"], Cell[BoxData[ \(B - \(72\ EMin\ xMin\)\/rMin\^2 \[Equal] 0\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(MinEnergyEq = \ \ \((\ PotEnergy /. x \[Rule] xMin)\) == \ \(-PotMin\)\)], "Input"], Cell[BoxData[ \(A + B\ xMin - \(36\ EMin\ xMin\^2\)\/rMin\^2 \[Equal] \(-PotMin\)\)], \ "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(SolutionPot = \ Solve[{SpringForceEq, MinEnergyEq}, {A, B}]\)\(\[IndentingNewLine]\) \)\)], "Input"], Cell[BoxData[ \({{A \[Rule] \(-PotMin\) - \(36\ EMin\ xMin\^2\)\/rMin\^2, B \[Rule] \(72\ EMin\ xMin\)\/rMin\^2}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(ModifiedPotEnergy\ = \ PotEnergy\ /. \ SolutionPot[\([1]\)]\)\)\)], "Input"], Cell[BoxData[ \(\(-PotMin\) - \(36\ EMin\ x\^2\)\/rMin\^2 + \(72\ EMin\ x\ xMin\)\/rMin\ \^2 - \(36\ EMin\ xMin\^2\)\/rMin\^2\)], "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["8.", "Subsubsection"], Cell["\<\ Use function Series to get the Taylor series representation, \ include Normal to eliminate the O[r - rMin]^3 term.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\[IndentingNewLine]\)\(HarmonicApprox = Normal[Series[LJModified, {r, rMin, 2}]]\)\)\)], "Input"], Cell[BoxData[ \(EMin - \(36\ EMin\ \((r - rMin)\)\^2\)\/rMin\^2\)], "Output"] }, Open ]], Cell["\<\ Note that this expression is a parabolic fit to the bottom of the \ potential well for the LJ function.\ \>", "Text"] }, Open ]], Cell[CellGroupData[{ Cell["9.", "Subsubsection"], Cell["\<\ The coefficient of the second-order term of HarmonicApprox is equal \ to k/2, where k is the spring constant. We can use the relation between \ vibrational frequency and spring constant to get the desired \ expression.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(k\ = \ 2\ SeriesCoefficient[Series[LJModified, {r, rMin, 2}], 2]\)], "Input"], Cell[BoxData[ \(\(-\(\(72\ EMin\)\/rMin\^2\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(VibFrequency\ = \ \(1\/\(2 \[Pi]\)\) \@\(k\/Mass\)\)], "Input"], Cell[BoxData[ \(\(3\ \@2\ \@\(-\(EMin\/\(Mass\ rMin\^2\)\)\)\)\/\[Pi]\)], "Output"] }, Open ]] }, Open ]] }, Open ]] }, FrontEndVersion->"5.1 for Macintosh", ScreenRectangle->{{0, 1148}, {0, 746}}, WindowSize->{694, 599}, WindowMargins->{{4, Automatic}, {Automatic, 4}}, StyleDefinitions -> "3016_Carter.nb" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[1754, 51, 48, 0, 64, "Title"], Cell[1805, 53, 127, 4, 110, "Subsubtitle"], Cell[CellGroupData[{ Cell[1957, 61, 47, 0, 65, "Subsection"], Cell[2007, 63, 791, 13, 250, "Text"], Cell[CellGroupData[{ Cell[2823, 80, 402, 8, 150, "Input"], Cell[3228, 90, 30559, 707, 320, 7666, 420, "GraphicsData", "PostScript", \ "Graphics"], Cell[33790, 799, 130, 3, 65, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[33969, 808, 46, 0, 65, "Subsection"], Cell[34018, 810, 119, 3, 66, "Text"], Cell[CellGroupData[{ Cell[34162, 817, 27, 0, 43, "Subsubsection"], Cell[34192, 819, 123, 1, 46, "Text"], Cell[CellGroupData[{ Cell[34340, 824, 54, 1, 74, "Input"], Cell[34397, 827, 52, 1, 90, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34486, 833, 75, 1, 50, "Input"], Cell[34564, 836, 477, 7, 202, "Output"] }, Open ]], Cell[35056, 846, 237, 4, 86, "Text"], Cell[CellGroupData[{ Cell[35318, 854, 66, 1, 50, "Input"], Cell[35387, 857, 72, 1, 92, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[35508, 864, 28, 0, 43, "Subsubsection"], Cell[35539, 866, 65, 0, 46, "Text"], Cell[CellGroupData[{ Cell[35629, 870, 59, 1, 50, "Input"], Cell[35691, 873, 57, 1, 92, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[35797, 880, 27, 0, 43, "Subsubsection"], Cell[35827, 882, 290, 5, 106, "Text"], Cell[36120, 889, 50, 1, 50, "Input"], Cell[CellGroupData[{ Cell[36195, 894, 105, 2, 50, "Input"], Cell[36303, 898, 109, 2, 67, "Output"] }, Open ]], Cell[36427, 903, 90, 3, 46, "Text"], Cell[CellGroupData[{ Cell[36542, 910, 68, 1, 50, "Input"], Cell[36613, 913, 88, 1, 92, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[36750, 920, 27, 0, 43, "Subsubsection"], Cell[36780, 922, 56, 0, 46, "Text"], Cell[CellGroupData[{ Cell[36861, 926, 118, 2, 70, "Input"], Cell[36982, 930, 93, 1, 92, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[37124, 937, 27, 0, 43, "Subsubsection"], Cell[37154, 939, 275, 5, 86, "Text"], Cell[CellGroupData[{ Cell[37454, 948, 203, 4, 90, "Input"], Cell[37660, 954, 60, 1, 89, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[37757, 960, 208, 4, 110, "Input"], Cell[37968, 966, 85, 1, 92, "Output"] }, Open ]], Cell[38068, 970, 75, 0, 46, "Text"], Cell[38146, 972, 286, 5, 110, "Input"], Cell[CellGroupData[{ Cell[38457, 981, 128, 2, 70, "Input"], Cell[38588, 985, 38, 1, 65, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38663, 991, 148, 3, 70, "Input"], Cell[38814, 996, 38, 1, 65, "Output"] }, Open ]], Cell[38867, 1000, 188, 8, 48, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[39092, 1013, 27, 0, 43, "Subsubsection"], Cell[39122, 1015, 156, 3, 66, "Text"], Cell[CellGroupData[{ Cell[39303, 1022, 372, 7, 150, "Input"], Cell[39678, 1031, 130, 3, 65, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[39845, 1039, 191, 3, 90, "Input"], Cell[40039, 1044, 37603, 915, 352, 10939, 581, "GraphicsData", "PostScript", \ "Graphics"], Cell[77645, 1961, 130, 3, 65, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[77824, 1970, 75, 3, 64, "Subsubsection"], Cell[77902, 1975, 149, 3, 66, "Text"], Cell[CellGroupData[{ Cell[78076, 1982, 74, 1, 50, "Input"], Cell[78153, 1985, 72, 1, 94, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[78262, 1991, 141, 3, 70, "Input"], Cell[78406, 1996, 76, 1, 92, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[78519, 2002, 114, 2, 50, "Input"], Cell[78636, 2006, 110, 3, 94, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[78783, 2014, 135, 3, 70, "Input"], Cell[78921, 2019, 140, 2, 94, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[79098, 2026, 132, 2, 70, "Input"], Cell[79233, 2030, 143, 2, 94, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[79425, 2038, 27, 0, 43, "Subsubsection"], Cell[79455, 2040, 138, 3, 66, "Text"], Cell[CellGroupData[{ Cell[79618, 2047, 128, 2, 70, "Input"], Cell[79749, 2051, 81, 1, 94, "Output"] }, Open ]], Cell[79845, 2055, 127, 3, 66, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[80009, 2063, 27, 0, 43, "Subsubsection"], Cell[80039, 2065, 242, 5, 86, "Text"], Cell[CellGroupData[{ Cell[80306, 2074, 105, 2, 70, "Input"], Cell[80414, 2078, 64, 1, 92, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[80515, 2084, 85, 1, 83, "Input"], Cell[80603, 2087, 87, 1, 128, "Output"] }, Open ]] }, Open ]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)