(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 4.0, MathReader 4.0, or any compatible application. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 233950, 7303]*) (*NotebookOutlinePosition[ 234913, 7336]*) (* CellTagsIndexPosition[ 234869, 7332]*) (*WindowFrame->Normal*) Notebook[{ Cell["\<\ Approximate angular scattering functions for He-He collisions\ \>", "Subtitle", Evaluatable->False, CellHorizontalScrolling->False, TextAlignment->Center], Cell["\<\ We wish to generate a set of approximate angular distributions for use in \ Monte Carlo calculations. For that purpose we have used the modified the \ screened Coulomb (or Rutherford) formula as described below. This notebook \ is based on our previous work on electrons in N2 and Ar-Ar elastic collisions \ and we will not repeat the derivations given in the Ar-Ar notebook.\ \>", "Text", Evaluatable->False], Cell[CellGroupData[{ Cell["Setup notebook environment", "Text"], Cell["ClearAll[\"Global`*\"]; ", "Input"], Cell["Remove[\"Global`*\"]; ", "Input"], Cell["startclock = SessionTime[];", "Input"], Cell["Off[General::spell]", "Input"], Cell["Off[General::spell1]", "Input"], Cell["Off[NumberForm::sigz]", "Input"], Cell["<< Graphics`Graphics`", "Input"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(now\ = \ StringForm["\<``/``/`` ``:``:``\>", \(Date[]\)[\([2]\)], \ \(Date[]\)[\([3]\)], \(Date[]\)[\([1]\)], \(Date[]\)[\([4]\)], \ \(Date[]\)[\([5]\)], \(Date[]\)[\([6]\)]]\)], "Input"], Cell[BoxData[ InterpretationBox["\<\"\\!\\(12\\)/\\!\\(9\\)/\\!\\(2000\\) \ \\!\\(21\\):\\!\\(20\\):\\!\\(24\\)\"\>", StringForm[ "``/``/`` ``:``:``", 12, 9, 2000, 21, 20, 24], Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["Differential scattering cross section", "Subsection", Evaluatable->False], Cell["Data from Newman, Smith, and Stebbings (1985).", "Subsubsection", Evaluatable->False], Cell["\<\ From their tabulation of the He-He differential cross sections we have\ \>", "SmallText", Evaluatable->False], Cell["!! c:\\mathbook\\helium\\hehecollisions\\newman85.dat", "Input", PageWidth->Infinity, Evaluatable->False], Cell["\<\ Note that entries of 1E-30 are used instead of a blank entry or zero so that \ our plotting program ProPlot will work for log plots.\ \>", "SmallText", Evaluatable->False], Cell["\<\ Read in the data from this file, transform to relative angle, and prepare \ plots\ \>", "Text"], Cell["\<\ inputfile := \ \"c:\\\\mathbook\\\\helium\\\\hehecollisions\\\\newman85.dat\"; stream = OpenRead[inputfile]; Find[stream,\"%\"]; Find[stream,\"%\"]; data1 = ReadList[stream, Number, RecordLists ->True]; Close[stream]; angle1 = Part[Transpose[data1],1]; relangle1 = angle1*2; q1diff500 = Part[Transpose[data1],2]*10^16; q1diff1500 = Part[Transpose[data1],3]*10^16; q1diff5000 = Part[Transpose[data1],4]*10^16; i1expcm500 = \ Table[Part[q1diff500,j]/4/Cos[Part[angle1*Pi/180,j]],{j,1,Length[angle1]}]; i1expcm1500 = \ Table[Part[q1diff1500,j]/4/Cos[Part[angle1*Pi/180,j]],{j,1,Length[angle1]}]; i1expcm5000 = \ Table[Part[q1diff5000,j]/4/Cos[Part[angle1*Pi/180,j]],{j,1,Length[angle1]}]; data1qdiff500 = Transpose[{relangle1,i1expcm500}]; data1qdiff1500 = Transpose[{relangle1,i1expcm1500}]; data1qdiff5000 = Transpose[{relangle1,i1expcm5000}]; data1q500plot = LogLogListPlot[data1qdiff500, PlotStyle ->{PointSize[0.025],Hue[0.55]}, PlotRange -> {{0.1,20},{1,10000}}, (*PlotRange -> {0.1,1000},*) DisplayFunction ->Identity]; data1q1500plot = LogLogListPlot[data1qdiff1500, PlotStyle ->{PointSize[0.025],Hue[1]}, PlotRange -> {{0.1,20},{1,10000}}, DisplayFunction -> Identity]; data1q5000plot = LogLogListPlot[data1qdiff5000, PlotStyle ->{PointSize[0.025],Hue[0.4]}, PlotRange -> {{0.1,20},{1,10000}}, DisplayFunction -> Identity];\ \>", "Input", PageWidth->Infinity, FontFamily->"Arial", FontSize->11, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["\<\ Show[data1q500plot,data1q1500plot,data1q5000plot, DisplayFunction -> $DisplayFunction, GridLines -> Automatic, AxesLabel -> {\"angle (deg)\",\"I(angle) (10^-20 m^2/str)\"}];\ \>", "Input", PageWidth->Infinity, Evaluatable->False, FontFamily->"Times New Roman", FontSize->12, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["!! c:\\mathbook\\helium\\hehecollisions\\nitz87a.dat", "Input", PageWidth->Infinity, Evaluatable->False], Cell["\<\ inputfile := \ \"c:\\\\mathbook\\\\helium\\\\hehecollisions\\\\nitz87a.dat\"; stream = OpenRead[inputfile]; Find[stream,\"%\"]; Find[stream,\"%\"]; data2 = ReadList[stream, Number, RecordLists ->True]; Close[stream]; angle2 = Part[Transpose[data2],1]; relangle2 = angle2*2; q2diff500 = Part[Transpose[data2],2]; q2diff1500 = Part[Transpose[data2],3]; q2diff5000 = Part[Transpose[data2],4]; i2expcm500 = \ Table[Part[q2diff500,j]/4/Cos[Part[angle2*Pi/180,j]],{j,1,Length[angle2]}]; i2expcm1500 = \ Table[Part[q2diff1500,j]/4/Cos[Part[angle2*Pi/180,j]],{j,1,Length[angle2]}]; i2expcm5000 = \ Table[Part[q2diff5000,j]/4/Cos[Part[angle2*Pi/180,j]],{j,1,Length[angle2]}]; data2qdiff500 = Transpose[{relangle2,i2expcm500}]; data2qdiff1500 = Transpose[{relangle2,i2expcm1500}]; data2qdiff5000 = Transpose[{relangle2,i2expcm5000}]; data2q500plot = LogLogListPlot[data2qdiff500, PlotStyle ->{PointSize[0.025],Hue[0.55]}, PlotRange -> {{0.1,10},{1,100000}}, DisplayFunction ->Identity]; data2q1500plot = LogLogListPlot[data2qdiff1500, PlotStyle ->{PointSize[0.025],Hue[1]}, PlotRange -> {{0.1,10},{1,100000}}, DisplayFunction -> Identity]; data2q5000plot = LogLogListPlot[data2qdiff5000, PlotStyle ->{PointSize[0.025],Hue[0.4]}, PlotRange -> {{0.1,10},{1,100000}}, DisplayFunction -> Identity];\ \>", "Input", PageWidth->Infinity, FontFamily->"Arial", FontSize->11, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["\<\ Show[data2q500plot,data2q1500plot,data2q5000plot,data1q500plot,data1q1500plot,\ data1q5000plot, DisplayFunction -> $DisplayFunction, GridLines -> Automatic, AxesLabel -> {\"rel. angle (deg)\",\"I(angle) (10^-20 \ m^2/str)\"}];\ \>", "Input", PageWidth->Infinity, Evaluatable->False, FontFamily->"Times New Roman", FontSize->12, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]] }, Open ]], Cell["Analysis of 500 eV data", "Subsection"], Cell["We will compare with screened-Coulomb scattering", "Subsubsection"], Cell["\<\ In center-of-mass the screened-Coulomb differential scattering formula is\ \>", "SmallText"], Cell[BoxData[ \(\(\(\ \)\(icoulombcm\ := \ a \((1 - \ Cos[thetacm]\ + \ 2*beta)\)^\(-2\)\)\)\)], "Input"], Cell["\<\ where thetacm is the scattering angle in center-of-mass. The \ total,diffusion, and viscosity cross sections are\ \>", "SmallText"], Cell[CellGroupData[{ Cell[BoxData[ \(qtcoulomb\ = \ 2*Pi*Integrate[icoulombcm*Sin[thetacm], {thetacm, 0, Pi}]\ // \ Simplify\)], "Input"], Cell[BoxData[ \(\(a\ \[Pi]\)\/\(beta + beta\^2\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qdcoulomb\ = \ 2*Pi*Integrate[ icoulombcm*\((1 - Cos[thetacm])\)*Sin[thetacm], {thetacm, 0, Pi}]\ // \ Simplify\)], "Input"], Cell[BoxData[ \(\(-\(\(2\ a\ \[Pi]\ \((1 + \((1 + beta)\)\ Log[\(-2\)\ beta] - \((1 + beta)\)\ Log[\(-2\)\ \((1 + beta)\)])\)\)\/\(1 + beta\)\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qvcoulomb\ = \ 2*Pi*Integrate[ icoulombcm*\((Sin[thetacm])\)^3, {thetacm, 0, Pi}]\ // \ Simplify\)], "Input"], Cell[BoxData[ \(\(-4\)\ a\ \[Pi]\ \((2 + \((1 + 2\ beta)\)\ Log[\(-2\)\ beta] - \((1 + 2\ beta)\)\ Log[\(-2\)\ \((1 + beta)\)])\)\)], "Output"] }, Open ]], Cell["\<\ We need the differential scattering expression transformed to the laboratory \ frame. From McDaniel (1964) we need the solution to\ \>", "SmallText"], Cell[CellGroupData[{ Cell[BoxData[ \(coulombrule\ = \ Solve[{\((icoulombcm\ *\ Sin[thetacm]\ /. \ thetacm\ -> \ 2*theta)\)* dthetacm\ == \ icoulomblab*Sin[theta]* dtheta, \ \n\t\t\ \ \ \ \ \ \ \ \ \ dthetacm\ == \ 2*dtheta}, icoulomblab]\)], "Input"], Cell[BoxData[ \({{icoulomblab \[Rule] \(4\ a\ Cos[theta]\)\/\((\(-1\) - 2\ beta + Cos[2\ \ theta])\)\^2}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(icoulomblab\ = \ Part[icoulomblab\ /. \ coulombrule, 1]\)], "Input"], Cell[BoxData[ \(\(4\ a\ Cos[theta]\)\/\((\(-1\) - 2\ beta + Cos[2\ theta])\)\^2\)], \ "Output"] }, Open ]], Cell["\<\ From Zhkarofsky et al. Eq. (7-11) when we replace the relative energy by the \ laboratory energy m*v^2/2 = mr*v^2 = en \ \>", "SmallText"], Cell[BoxData[ \(a\ := \ \(\(10^20\)\(*\)\(\((z^2*0.529*^-10*13.6/ en)\)^2\)\(\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \)\( (*\ 10^\(-20\) m^2/str*) \)\)\)], "Input"], Cell["\<\ We seem to find the best fit with z = 2 and beta = 0.004, so that for en = \ 500 eV we have\ \>", "SmallText"], Cell[BoxData[ \(\(\(\t\)\(coulomb500plot\ = \ LogLogPlot[\((icoulombcm\ //. {z\ -> 1.8, en\ -> 500, beta\ -> \ 0.002, thetacm\ -> \ relangle* Pi/180})\), \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {relangle, 0.1, 20}, \ PlotRange\ -> \ {{0.1, 20}, {1, 100000}}, \ DisplayFunction\ -> \ Identity];\)\)\)], "Input"], Cell["The Rice University experimental data for 500 eV", "Subsubsection"], Cell[BoxData[ \(\(list500\ = \ Union[Drop[Drop[data1qdiff500, {6}], {10}], data2qdiff500];\)\)], "Input"], Cell[BoxData[ \(\(dataqdiff500\ = \ DeleteCases[ Table[\ If[Part[list500, j, 2]\ >= \ 10^\(-8\), Part[list500, j]], {j, 1, Length[list500]}], Null];\)\)], "Input"], Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \ \)\(dataq500plot\ = \ LogLogListPlot[ dataqdiff500, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.1, 20}, {1, 100000}}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity];\)\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[dataq500plot, coulomb500plot, (*\(AxesLabel -> {"\", "\"}\)\(,\)*) \n\t\tDisplayFunction -> $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.434588 0.434588 0 0.123607 [ [.13082 -0.0125 -11 -12.875 ] [.13082 -0.0125 11 0 ] [.30376 -0.0125 -11 -12.875 ] [.30376 -0.0125 11 0 ] [.43459 -0.0125 -5 -12.875 ] [.43459 -0.0125 5 0 ] [.56541 -0.0125 -5 -12.875 ] [.56541 -0.0125 5 0 ] [.73835 -0.0125 -5 -12.875 ] [.73835 -0.0125 5 0 ] [.86918 -0.0125 -8 -12.875 ] [.86918 -0.0125 8 0 ] [1 -0.0125 -8 -12.875 ] [1 -0.0125 8 0 ] [-0.0125 .12361 -16 -6.4375 ] [-0.0125 .12361 0 6.4375 ] [-0.0125 .24721 -22 -6.4375 ] [-0.0125 .24721 0 6.4375 ] [-0.0125 .37082 -28 -6.4375 ] [-0.0125 .37082 0 6.4375 ] [-0.0125 .49443 -34 -6.4375 ] [-0.0125 .49443 0 6.4375 ] [-0.0125 .61803 -46 -6.4375 ] [-0.0125 .61803 0 6.4375 ] [ -0.0005 -0.0005 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .13082 0 m .13082 .00625 L s gsave .13082 -0.0125 -72 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (0.2) show 81.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 .30376 0 m .30376 .00625 L s gsave .30376 -0.0125 -72 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (0.5) show 81.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 .43459 0 m .43459 .00625 L s gsave .43459 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (1) show 69.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 .56541 0 m .56541 .00625 L s gsave .56541 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (2) show 69.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 .73835 0 m .73835 .00625 L s gsave .73835 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (5) show 69.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 .86918 0 m .86918 .00625 L s gsave .86918 -0.0125 -69 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore 1 0 m 1 .00625 L s gsave 1 -0.0125 -69 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (20) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore .001 w .20735 0 m .20735 .00375 L s .26165 0 m .26165 .00375 L s .33818 0 m .33818 .00375 L s .36727 0 m .36727 .00375 L s .39247 0 m .39247 .00375 L s .4147 0 m .4147 .00375 L s .64194 0 m .64194 .00375 L s .69624 0 m .69624 .00375 L s .77276 0 m .77276 .00375 L s .80186 0 m .80186 .00375 L s .82706 0 m .82706 .00375 L s .84929 0 m .84929 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12361 m .00625 .12361 L s gsave -0.0125 .12361 -77 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 .24721 m .00625 .24721 L s gsave -0.0125 .24721 -83 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (100) show 81.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 .37082 m .00625 .37082 L s gsave -0.0125 .37082 -89 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (1000) show 87.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 .49443 m .00625 .49443 L s gsave -0.0125 .49443 -95 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10000) 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 0 .61803 m .00625 .61803 L s gsave -0.0125 .61803 -107 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (100000.) show 105.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 .001 w 0 .03721 m .00375 .03721 L s 0 .05898 m .00375 .05898 L s 0 .07442 m .00375 .07442 L s 0 .0864 m .00375 .0864 L s 0 .09618 m .00375 .09618 L s 0 .10446 m .00375 .10446 L s 0 .11163 m .00375 .11163 L s 0 .11795 m .00375 .11795 L s 0 .16082 m .00375 .16082 L s 0 .18258 m .00375 .18258 L s 0 .19803 m .00375 .19803 L s 0 .21 m .00375 .21 L s 0 .21979 m .00375 .21979 L s 0 .22807 m .00375 .22807 L s 0 .23523 m .00375 .23523 L s 0 .24156 m .00375 .24156 L s 0 .28442 m .00375 .28442 L s 0 .30619 m .00375 .30619 L s 0 .32163 m .00375 .32163 L s 0 .33361 m .00375 .33361 L s 0 .3434 m .00375 .3434 L s 0 .35167 m .00375 .35167 L s 0 .35884 m .00375 .35884 L s 0 .36516 m .00375 .36516 L s 0 .40803 m .00375 .40803 L s 0 .4298 m .00375 .4298 L s 0 .44524 m .00375 .44524 L s 0 .45722 m .00375 .45722 L s 0 .46701 m .00375 .46701 L s 0 .47528 m .00375 .47528 L s 0 .48245 m .00375 .48245 L s 0 .48877 m .00375 .48877 L s 0 .53164 m .00375 .53164 L s 0 .5534 m .00375 .5534 L s 0 .56885 m .00375 .56885 L s 0 .58082 m .00375 .58082 L s 0 .59061 m .00375 .59061 L s 0 .59889 m .00375 .59889 L s 0 .60606 m .00375 .60606 L s 0 .61238 m .00375 .61238 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .10015 .60787 Mdot .11509 .59366 Mdot .12893 .57635 Mdot .14359 .56051 Mdot .15555 .54029 Mdot .1668 .52881 Mdot .17888 .51465 Mdot .18034 .53438 Mdot .18886 .50651 Mdot .19833 .50488 Mdot .20861 .49831 Mdot .21716 .50134 Mdot .22534 .4947 Mdot .23428 .49549 Mdot .24176 .50075 Mdot .24896 .49692 Mdot .25197 .49348 Mdot .25687 .49389 Mdot .26353 .48507 Mdot .27353 .47952 Mdot .28638 .46855 Mdot .29606 .46678 Mdot .2984 .45984 Mdot .30898 .45303 Mdot .31968 .44591 Mdot .32981 .44354 Mdot .33178 .44354 Mdot .3388 .43921 Mdot .34798 .43691 Mdot .35674 .42751 Mdot .3618 .42798 Mdot .36401 .42263 Mdot .37259 .41357 Mdot .38029 .4102 Mdot .38769 .40975 Mdot .39435 .39993 Mdot .40123 .39672 Mdot .40787 .39762 Mdot .41386 .38993 Mdot .4147 .39169 Mdot .42008 .38593 Mdot .43459 .37993 Mdot .45258 .37122 Mdot .469 .35766 Mdot .48411 .3479 Mdot .49809 .34158 Mdot .51362 .33468 Mdot .52564 .31959 Mdot .53024 .31888 Mdot .53695 .31331 Mdot .54761 .30796 Mdot .56915 .29511 Mdot .59982 .27696 Mdot .6262 .25932 Mdot .65055 .24476 Mdot .67103 .23126 Mdot .69049 .218 Mdot .70723 .20868 Mdot .72343 .19874 Mdot .7376 .18856 Mdot .75147 .18036 Mdot .76374 .17243 Mdot .77588 .16496 Mdot .78671 .15836 Mdot .79749 .15187 Mdot .80717 .14421 Mdot .81638 .13883 Mdot .82564 .13124 Mdot .83401 .12533 Mdot .84246 .11899 Mdot .85013 .11521 Mdot .8579 .10597 Mdot .86498 .09857 Mdot .87217 .09269 Mdot 0 g .5 Mabswid 0 .26361 m .04146 .26359 L .0726 .26357 L .10314 .26353 L .12807 .2635 L .16561 .26342 L .18391 .26337 L .19936 .26332 L .22818 .2632 L .25555 .26304 L .28043 .26286 L .3011 .26267 L .3373 .26221 L .36526 .26171 L .3931 .26105 L .41622 .26034 L .43701 .25954 L .45691 .2586 L .48816 .25668 L .51963 .25403 L .54423 .25133 L .58253 .24565 L .61673 .23865 L .64453 .23126 L .67004 .22292 L .70879 .20696 L .74282 .18938 L .77064 .17243 L .7963 .15477 L .81811 .13833 L .83695 .12315 L .85513 .10771 L .87111 .09355 L .8853 .08056 L .89932 .06738 L .9119 .05529 L .92441 .04304 L .93573 .03178 L .94601 .02143 L .95636 .01089 L .96582 .00117 L s .96582 .00117 m .96695 0 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgoooo`00oooo1@0000;oool5 0;?oDoooo`00oooo2P2coiGoool0 033oool00`000?ooooooo`1=oooo202coiSoool0033oool20000COooo`P0/onIoooo000`oooo00<0 003oooooool0B_ooo`T0/onJoooo000`oooo00<0003oooooool0B?ooo`X0/onKoooo000`oooo00<0 003oooooool0A_ooo`X0/onMoooo000`oooo0P0004Coool<0;?oW_ooo`00oooo`T0/onYoooo000`oooo0P0003_oool90;?oZ_ooo`000;?o[oooo`004Oooo`050000oooooooo oooo000000;oool010000?ooooooo`0000;oool010000?ooooooo`0000;oool010000?ooooooo`00 00;oool00`000?ooooooo`03oooo0`0002cooolA0;?o/?ooo`004Oooo`050000oooooooooooo0000 00;oool010000?ooooooo`0000;oool010000?ooooooo`0000;oool010000?ooooooo`0000;oool0 0`000?ooooooo`03oooo00<0003oooooool0:_oooa80/onaoooo000Aoooo00D0003oooooooooool0 00000_ooo`040000oooooooo00000_ooo`040000oooooooo00000_ooo`040000oooooooo00000_oo o`030000oooooooo00?oool00`000?ooooooo`0Woooo502cok;oool0017oool01@000?oooooooooo o`000002oooo00@0003oooooool00002oooo00@0003oooooool00002oooo00@0003oooooool00002 oooo00<0003oooooool00oooo`030000oooooooo02GooolE0;?o/oooo`004?ooo`800004oooo0P00 00Coool200001?ooo`800004oooo0P0000Koool00`000?ooooooo`0Toooo3@2co`?oool30;?o]_oo o`00"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-1.48874, -0.773956, \ 0.00999921, 0.0351561}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \ \)\(ListPlot[ dataqdiff500, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.3, 0.6}, {0, 10000}} (*\(,\)\(\n\)\(\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \)\(DisplayFunction\ -> Identity\)*) ];\)\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -1 3.33333 -1.6263e-018 6.18034e-005 [ [.16667 -0.0125 -12 -9 ] [.16667 -0.0125 12 0 ] [.33333 -0.0125 -9 -9 ] [.33333 -0.0125 9 0 ] [.5 -0.0125 -12 -9 ] [.5 -0.0125 12 0 ] [.66667 -0.0125 -9 -9 ] [.66667 -0.0125 9 0 ] [.83333 -0.0125 -12 -9 ] [.83333 -0.0125 12 0 ] [-0.0125 .12361 -24 -4.5 ] [-0.0125 .12361 0 4.5 ] [-0.0125 .24721 -24 -4.5 ] [-0.0125 .24721 0 4.5 ] [-0.0125 .37082 -24 -4.5 ] [-0.0125 .37082 0 4.5 ] [-0.0125 .49443 -24 -4.5 ] [-0.0125 .49443 0 4.5 ] [-0.0125 .61803 -30 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .16667 0 m .16667 .00625 L s [(0.35)] .16667 -0.0125 0 1 Mshowa .33333 0 m .33333 .00625 L s [(0.4)] .33333 -0.0125 0 1 Mshowa .5 0 m .5 .00625 L s [(0.45)] .5 -0.0125 0 1 Mshowa .66667 0 m .66667 .00625 L s [(0.5)] .66667 -0.0125 0 1 Mshowa .83333 0 m .83333 .00625 L s [(0.55)] .83333 -0.0125 0 1 Mshowa .125 Mabswid .03333 0 m .03333 .00375 L s .06667 0 m .06667 .00375 L s .1 0 m .1 .00375 L s .13333 0 m .13333 .00375 L s .2 0 m .2 .00375 L s .23333 0 m .23333 .00375 L s .26667 0 m .26667 .00375 L s .3 0 m .3 .00375 L s .36667 0 m .36667 .00375 L s .4 0 m .4 .00375 L s .43333 0 m .43333 .00375 L s .46667 0 m .46667 .00375 L s .53333 0 m .53333 .00375 L s .56667 0 m .56667 .00375 L s .6 0 m .6 .00375 L s .63333 0 m .63333 .00375 L s .7 0 m .7 .00375 L s .73333 0 m .73333 .00375 L s .76667 0 m .76667 .00375 L s .8 0 m .8 .00375 L s .86667 0 m .86667 .00375 L s .9 0 m .9 .00375 L s .93333 0 m .93333 .00375 L s .96667 0 m .96667 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12361 m .00625 .12361 L s [(2000)] -0.0125 .12361 1 0 Mshowa 0 .24721 m .00625 .24721 L s [(4000)] -0.0125 .24721 1 0 Mshowa 0 .37082 m .00625 .37082 L s [(6000)] -0.0125 .37082 1 0 Mshowa 0 .49443 m .00625 .49443 L s [(8000)] -0.0125 .49443 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(10000)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .0309 m .00375 .0309 L s 0 .0618 m .00375 .0618 L s 0 .09271 m .00375 .09271 L s 0 .15451 m .00375 .15451 L s 0 .18541 m .00375 .18541 L s 0 .21631 m .00375 .21631 L s 0 .27812 m .00375 .27812 L s 0 .30902 m .00375 .30902 L s 0 .33992 m .00375 .33992 L s 0 .40172 m .00375 .40172 L s 0 .43262 m .00375 .43262 L s 0 .46353 m .00375 .46353 L s 0 .52533 m .00375 .52533 L s 0 .55623 m .00375 .55623 L s 0 .58713 m .00375 .58713 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .26667 .60722 Mdot .3 .61186 Mdot .34667 .51915 Mdot .42 .46816 Mdot .52 .38164 Mdot .6 .36928 Mdot .62 .32447 Mdot .71333 .28584 Mdot .81333 .25031 Mdot .91333 .23949 Mdot .93333 .23949 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0.260791, -884.13, \ 0.00118192, 63.7461}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qdiff500int\ = \ Interpolation[dataqdiff500]\)], "Input"], Cell[BoxData[ TagBox[\(InterpolatingFunction[{{0.036`, 10.16`}}, "<>"]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(norm\ = \ 2*Pi*NIntegrate[Sin[angle*3.1416/180], \ {angle, 0, 180}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(12.566370614189621`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qt500\ = \ 2*Pi*NIntegrate[ qdiff500int[relangle]*Sin[relangle*3.1416/180], \ {relangle, 0.036, 10.16}]*3.1416/180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.7082968749999999`\)."\)], "Message"], Cell[BoxData[ \(13.806061997804553`\)], "Output"] }, Open ]], Cell["\<\ Presumably the following is closer to the true total cross section.\ \>", "SmallText"], Cell[CellGroupData[{ Cell[BoxData[ \(qt500\ = \ 2*Pi*NIntegrate[ qdiff500int[relangle]*Sin[relangle*3.1416/180], \ {relangle, 0.036, 10.16}]*3.1416/180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.7082968749999999`\)."\)], "Message"], Cell[BoxData[ \(13.806061997804553`\)], "Output"] }, Open ]], Cell["\<\ This value is about 15% larger than that given by Newman et al. (1985), but \ is only 7% larger than the experimental value of Savola (1973) as obtained \ using the detector angle of 0.056 deg. This difference may be the result of \ the improved differential scattering data available at low angles from the \ experiments of Nitz as used in our calculation but not in Newman's \ calculation.\ \>", "SmallText"], Cell[CellGroupData[{ Cell[BoxData[ \(qm500\ = \ 2*Pi*NIntegrate[ qdiff500int[relangle]*\((1 - Cos[relangle*3.1416/180])\)* Sin[relangle*3.1416/ 180], \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {relangle, 0.036, 10.16}]*3.1416/180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.8664843749999999`\)."\)], "Message"], Cell[BoxData[ \(0.011621917815295224`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qv500\ = \ 2*Pi*NIntegrate[ qdiff500int[ relangle]*\((Sin[ relangle*3.1416/ 180])\)^3, \ \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \t{relangle, 0.036, 10.16}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.8664843749999999`\)."\)], "Message"], Cell[BoxData[ \(0.023160590632078017`\)], "Output"] }, Open ]], Cell["\<\ These qm and qv values are very small and suggest we need to look very \ carefully at the large angle scattering, e.g., how does experiment compare \ with screened Coulomb scattering.\ \>", "SmallText"], Cell[BoxData[ \(\(angleqdiff500\ = \ Part[Transpose[dataqdiff500], 1];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(83\)], "Output"] }, Open ]], Cell[BoxData[ \(\(listqdiff500\ = \ \ Part[Transpose[dataqdiff500], 2];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(83\)], "Output"] }, Open ]], Cell[BoxData[ \(\(integrandqv500\ = \ Table[\((Part[listqdiff500, j]*\((Sin[ Part[angleqdiff500, j]* Pi/180])\)^3)\), \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {j, 1, Length[angleqdiff500]}];\)\)], "Input"], Cell[BoxData[ \(\(tableintegrandqv500\ = \ Transpose[{angleqdiff500, integrandqv500}];\)\)], "Input"], Cell[BoxData[ \(\(\ \ \ \ \ \ \ \ \ \ \ \(integrand500plot\ = \ LogLogListPlot[tableintegrandqv500, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (* log\ log\ plot*) \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.1, 20}, {1*^-5, 0.1}}, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity]; \)\)\)], "Input", Evaluatable->False], Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \)\(integrand500plot\ = \ ListPlot[ tableintegrandqv500, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*linear\ plot*) \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0, 90}, {0, 0.04}}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity];\)\(\n\)\)\)], "Input"], Cell["\<\ The integrand for the viscosity integral and screened Coulomb scattering \ is\ \>", "SmallText"], Cell[BoxData[ \(\(\tcoulombint500plot\ = \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (* log\ log\ plot*) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \((LogLogPlot[ \((\((icoulombcm*\((Sin[thetacm])\)^3)\)\ //. {z\ -> 1.8, en\ -> 500, beta\ -> \ 0.002, thetacm\ -> \ relangle*Pi/180})\), \n \t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {relangle, 0.1, 20}, \ \ AxesLabel\ -> \ {"\", \ "\"}, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \t\ \t\tPlotRange\ -> \ {{0.1, 10}, {1*^-5, 0.1}}, \ DisplayFunction\ -> \ Identity]; \)\)\)], "Input", Evaluatable->False], Cell["\<\ coulombint500plot = (*linear plot*) Plot[icoulombcm*Sin[thetacm]^3 //. {z -> 1.8, en ->500, beta -> 0.002, thetacm -> relangle*Pi/180}, \ {relangle, 0, 90}, AxesLabel -> {\"rel. angle (deg)\", \"i\"}, PlotRange -> {{0, 90}, {0, 0.04}}, DisplayFunction -> Identity]; \ \>", "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[integrand500plot, coulombint500plot, \ DisplayFunction\ -> \ $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 -1.9082e-017 0.0111111 9.1073e-018 15.4508 [ [.22222 -0.0125 -6 -9 ] [.22222 -0.0125 6 0 ] [.44444 -0.0125 -6 -9 ] [.44444 -0.0125 6 0 ] [.66667 -0.0125 -6 -9 ] [.66667 -0.0125 6 0 ] [.88889 -0.0125 -6 -9 ] [.88889 -0.0125 6 0 ] [-0.0125 .15451 -24 -4.5 ] [-0.0125 .15451 0 4.5 ] [-0.0125 .30902 -24 -4.5 ] [-0.0125 .30902 0 4.5 ] [-0.0125 .46353 -24 -4.5 ] [-0.0125 .46353 0 4.5 ] [-0.0125 .61803 -24 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .22222 0 m .22222 .00625 L s [(20)] .22222 -0.0125 0 1 Mshowa .44444 0 m .44444 .00625 L s [(40)] .44444 -0.0125 0 1 Mshowa .66667 0 m .66667 .00625 L s [(60)] .66667 -0.0125 0 1 Mshowa .88889 0 m .88889 .00625 L s [(80)] .88889 -0.0125 0 1 Mshowa .125 Mabswid .05556 0 m .05556 .00375 L s .11111 0 m .11111 .00375 L s .16667 0 m .16667 .00375 L s .27778 0 m .27778 .00375 L s .33333 0 m .33333 .00375 L s .38889 0 m .38889 .00375 L s .5 0 m .5 .00375 L s .55556 0 m .55556 .00375 L s .61111 0 m .61111 .00375 L s .72222 0 m .72222 .00375 L s .77778 0 m .77778 .00375 L s .83333 0 m .83333 .00375 L s .94444 0 m .94444 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .15451 m .00625 .15451 L s [(0.01)] -0.0125 .15451 1 0 Mshowa 0 .30902 m .00625 .30902 L s [(0.02)] -0.0125 .30902 1 0 Mshowa 0 .46353 m .00625 .46353 L s [(0.03)] -0.0125 .46353 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(0.04)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .0309 m .00375 .0309 L s 0 .0618 m .00375 .0618 L s 0 .09271 m .00375 .09271 L s 0 .12361 m .00375 .12361 L s 0 .18541 m .00375 .18541 L s 0 .21631 m .00375 .21631 L s 0 .24721 m .00375 .24721 L s 0 .27812 m .00375 .27812 L s 0 .33992 m .00375 .33992 L s 0 .37082 m .00375 .37082 L s 0 .40172 m .00375 .40172 L s 0 .43262 m .00375 .43262 L s 0 .49443 m .00375 .49443 L s 0 .52533 m .00375 .52533 L s 0 .55623 m .00375 .55623 L s 0 .58713 m .00375 .58713 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .0004 .00176 Mdot .00058 .00462 Mdot .00073 .00827 Mdot .00089 .01272 Mdot .00107 .01773 Mdot .00122 .02334 Mdot .00138 .02788 Mdot .00156 .03144 Mdot .00171 .03165 Mdot .00189 .0334 Mdot .00204 .03249 Mdot .0022 .02933 Mdot .00238 .02757 Mdot .00253 .02288 Mdot .00269 .02209 Mdot .00287 .02056 Mdot .00289 .03039 Mdot .00302 .0207 Mdot .00318 .02335 Mdot .00336 .02432 Mdot .00351 .02948 Mdot .00367 .02967 Mdot .00384 .03471 Mdot .004 .04312 Mdot .00416 .04501 Mdot .00422 .04429 Mdot .00433 .04824 Mdot .00449 .0455 Mdot .00473 .0481 Mdot .00507 .0481 Mdot .00533 .05428 Mdot .0054 .0495 Mdot .00571 .05159 Mdot .00604 .05356 Mdot .00638 .0602 Mdot .00644 .0621 Mdot .00669 .06407 Mdot .00702 .07102 Mdot .00736 .06851 Mdot .00756 .0749 Mdot .00764 .07022 Mdot .008 .06799 Mdot .00833 .07216 Mdot .00867 .08049 Mdot .00898 .07453 Mdot .00931 .07831 Mdot .00964 .0885 Mdot .00996 .08434 Mdot .01 .08832 Mdot .01029 .08642 Mdot .01111 .09733 Mdot .01222 .11014 Mdot .01333 .11106 Mdot .01444 .11774 Mdot .01556 .13071 Mdot .01689 .14709 Mdot .018 .13442 Mdot .01844 .14274 Mdot .01911 .14311 Mdot .02022 .15346 Mdot .02267 .17008 Mdot .02667 .19746 Mdot .03067 .21613 Mdot .03489 .2426 Mdot .03889 .26114 Mdot .04311 .27782 Mdot .04711 .30462 Mdot .05133 .32727 Mdot .05533 .33891 Mdot .05956 .36249 Mdot .06356 .37984 Mdot .06778 .40057 Mdot .07178 .42041 Mdot .076 .44189 Mdot .08 .44653 Mdot .084 .46723 Mdot .08822 .46948 Mdot .09222 .47999 Mdot .09644 .48733 Mdot .10044 .51254 Mdot .10467 .4877 Mdot .10867 .47504 Mdot .11289 .47679 Mdot 0 g .5 Mabswid 0 0 m .04057 .23851 L .05077 .32948 L .06168 .40256 L .06795 .43234 L .07365 .45244 L .0764 .46001 L .07939 .4668 L .08219 .47195 L .08481 .47579 L .08721 .47854 L .0898 .48075 L .09127 .48169 L .09267 .48238 L .09393 .48284 L .0953 .48318 L .09605 .48329 L .09676 .48336 L .09809 .48338 L .09883 .48332 L .09952 .48324 L .10027 .48311 L .10106 .48293 L .10238 .48253 L .10361 .48206 L .10641 .48066 L .11145 .47709 L .11604 .47287 L .12636 .46095 L .1682 .39931 L .20844 .34229 L .25126 .29183 L .29249 .25245 L .33213 .22135 L .37435 .19384 L .41497 .17159 L .45816 .15146 L .49977 .13478 L .53978 .12076 L .58237 .10762 L .62337 .09641 L .66694 .0858 L .70892 .07666 L .74931 .06873 L .79228 .0611 L .83365 .05445 L .87343 .04864 L .91579 .04301 L .95655 .03808 L .99989 .03332 L Mistroke 1 .03331 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgoooo00<0003oooooool0gOooo`006oooo`030000oooooooo 00koool00`2co`00002co`040;?o3?ooo`030000oooooooo0=koool001_oool00`000?ooooooo`0> oooo0P2co`0300000;?o0;?o0080/ol3oool00005oooo0000oooooooo000000Soool010000?ooooooo`0000;oool010000?oo ooooo`0000Coool00`000?ooooooo`0Aoooo00<00000/ol0/ol00P2co`Soool00`000?ooooooo`3Q oooo00001Oooo`000?ooooooo`000008oooo00@0003oooooool00004oooo00<0003oooooool00ooo o`<0000@oooo0P2co`0300000;?o0;?o0080/ol6oooo00<0003oooooool0h_ooo`0000Goool0003o ooooool000002?ooo`040000oooooooo00001Oooo`030000oooooooo00;oool00`000?ooooooo`0@ oooo0P2co`0300000;?o0;?o0080/ol00oooo`2co`2co`020;?o00<0003oooooool0hoooo`0000Go ool0003oooooool000002?ooo`040000oooooooo00000_ooo`040000oooooooo00001?ooo`030000 oooooooo013oool20;?o00<00000/ol0/ol01P2co`0300000;?o0;?o0>Coool000;oool200002_oo o`800004oooo0P0000Goool00`000?ooooooo`0@oooo0`2co`0300000;?o0;?o00@0/ol01000002c o`2co`2conCoool001_oool00`000?ooooooo`0Aoooo0`2co`0300000;?o0;?o0080/ol00`00002c o`2co`020;?oi?ooo`006oooo`030000oooooooo01;oool30;?o100000D0/ooToooo000Koooo00<0 003oooooool04_ooo`/0/ooUoooo000Koooo00<0003oooooool04oooo`P0/ooWoooo000Koooo0P00 01Coool80;?oioooo`006oooo`030000oooooooo01Coool60;?oj?ooo`006oooo`030000oooooooo 01Goool40;?ojOooo`006oooo`030000oooooooo01Coool60;?oj?ooo`006oooo`030000oooooooo 01Coool60;?oj?ooo`006oooo`030000oooooooo01Coool60;?oj?ooo`006oooo`030000oooooooo 01Coool60;?oj?ooo`006oooo`030000oooooooo01Goool40;?ojOooo`006oooo`80003ooooo1?oo o`006oooo`030000oooooooo0?ooool3oooo000Koooo00<0003oooooool0ooooo`?oool001_oool0 0`000?ooooooo`3ooooo0oooo`006oooo`030000oooooooo0?ooool3oooo000Koooo00<0003ooooo ool0ooooo`?oool001_oool00`000?ooooooo`3ooooo0oooo`006oooo`030000oooooooo0?ooool3 oooo000Koooo0P000?ooool4oooo000Koooo00<0003oooooool0ooooo`?oool001_oool00`000?oo ooooo`3ooooo0oooo`006oooo`030000oooooooo0?ooool3oooo000Koooo00<0003oooooool0oooo o`?oool001_oool00`000?ooooooo`3ooooo0oooo`006oooo`030000oooooooo0?ooool3oooo000K oooo00<0003oooooool0ooooo`?oool001_oool20000ooooo`Coool001_oool00`000?ooooooo`3o oooo0oooo`006oooo`030000oooooooo0?ooool3oooo000Koooo00<0003oooooool0ooooo`?oool0 01_oool00`000?ooooooo`3ooooo0oooo`006oooo`030000oooooooo0?ooool3oooo0002oooo0P00 00Coool200001?ooo`800005oooo0`0000?oool00`000?ooooooo`3ooooo0oooo`0000Goool0003o ooooool000002?ooo`040000oooooooo00001Oooo`030000oooooooo00;oool00`000?ooooooo`3o oooo0oooo`0000Goool0003oooooool000002?ooo`040000oooooooo00000_ooo`D00003oooo0`00 0?ooool3oooo00001Oooo`000?ooooooo`000008oooo00@0003oooooool00002oooo00@0003ooooo ool0003ooooo2_ooo`0000Goool0003oooooool000002?ooo`040000oooooooo00000oooo`030000 oooo00000?ooool:oooo0002oooo0P0000[oool200001Oooo`80003ooooo2_ooo`00ooooob7oool0 0001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-9.55474, -0.00305647, \ 0.347217, 0.000249694}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqt500\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 1.8, en\ -> 500, beta\ -> \ 0.002, thetacm\ -> \ relangle* Pi/180})\)\ \ *\ \ \n\t\t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ Sin[ relangle*3.1416/180], {relangle, 10.16, 180}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(0.6870487069851887`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ delqm500 = 2*Pi*NIntegrate[(icoulombcm //. {z -> 1.8, en -> 500, beta -> \ 0.002, thetacm -> relangle*Pi/180})* (1 - Cos[relangle*Pi/180])* Sin[relangle*Pi/180], {relangle, \ 10.16, 180}]* Pi/180\ \>", "Input"], Cell[BoxData[ \(0.06038622346977897`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqv500\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 1.8, en\ -> 500, beta\ -> \ 0.002, thetacm\ -> \ relangle* Pi/180})\)*\n\t\t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \((Sin[relangle*3.1416/180])\)^3, \ {relangle, 10.16, 180}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.0941696557793748`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(coulombqv500\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 1.8, en\ -> 500, beta\ -> \ 0.002, thetacm -> \ relangle* Pi/180})\)*\n\t\t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \((Sin[ relangle*3.1416/180])\)^3, \ {relangle, 0, 180}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(0.11584281580244375`\)], "Output"] }, Open ]], Cell["\<\ Something may be wrong with these three integrals - they seem too large.\ \>", "Subsubtitle"], Cell["Analysis of 1500 eV data", "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ \(list1500\ = \ Union[Drop[Drop[data1qdiff1500, {6}], {10}], data2qdiff1500]\)], "Input"], Cell[BoxData[ \({{0.036`, 745000.0367642778`}, {0.052`, 690000.0710428806`}, {0.066`, 512500.0850054065`}, {0.08`, 430000.1047884136`}, {0.096`, 297500.10439851263`}, {0.11`, 183000.0843143265`}, {0.124`, 125750.07362347704`}, {0.14`, 81750.0610110994`}, {0.154`, 53250.04808682784`}, {0.17`, 37750.041541286795`}, {0.184`, 31500.040607982617`}, {0.198`, 25500.038065877827`}, {0.214`, 22775.039714702318`}, {0.228`, 22950.045427396588`}, {0.242`, 20700.046160088485`}, {0.258`, 18000.045622342703`}, {0.26`, 16900.043500983942`}, {0.272`, 16900.047609133846`}, {0.286`, 13475.041968815824`}, {0.302`, 12400.043062731713`}, {0.316`, 11375.043250556`}, {0.33`, 10475.04343583356`}, {0.346`, 9525.043419376216`}, {0.36`, 8775.043303067387`}, {0.374`, 7475.039812649196`}, {0.38`, 8075.044399334685`}, {0.39`, 6950.0402513881845`}, {0.404`, 6975.043348486577`}, {0.412`, 6400.041365806245`}, {0.426`, 2.500017275334127`*^-10}, {0.456`, 2.5000197941705397`*^-10}, {0.48`, 4700.041233315392`}, {0.486`, 2.50002248436104`*^-10}, {0.504`, 4025.038930968351`}, {0.514`, 2.500025149805785`*^-10}, {0.544`, 2.50002817128352`*^-10}, {0.574`, 2.500031364118212`*^-10}, {0.58`, 3100.0397088276204`}, {0.596`, 2925.0395628922474`}, {0.602`, 2.5000344987002783`*^-10}, {0.632`, 2.500038022828243`*^-10}, {0.662`, 2.5000417183165424`*^-10}, {0.68`, 2175.038295540816`}, {0.688`, 2132.538435911745`}, {0.72`, 2.500049348833761`*^-10}, {0.75`, 2.5000535469865643`*^-10}, {0.78`, 1697.5393253069099`}, {0.808`, 2.5000621493334707`*^-10}, {0.838`, 2.500066850163158`*^-10}, {0.868`, 1252.5359329038924`}, {0.896`, 2.5000764243651197`*^-10}, {0.9`, 1200.0370119677887`}, {0.926`, 2.500081627885553`*^-10}, {0.962`, 1047.5369131984269`}, {1.`, 957.5364600543431`}, {1.1`, 755.0347867410987`}, {1.2`, 637.53495644617`}, {1.3`, 517.5333031333787`}, {1.4`, 467.534892287762`}, {1.52`, 380.033432385083`}, {1.62`, 365.03647739440595`}, {1.66`, 335.0351531062393`}, {1.72`, 310.0349240101815`}, {1.82`, 292.53689597944685`}, {1.92`, 245.2844292074305`}, {2.04`, 224.7856190339409`}, {2.4`, 148.7826304895223`}, {2.76`, 99.52886762002373`}, {3.14`, 63.773940880825094`}, {3.5`, 40.268781728792725`}, {3.88`, 26.26505445413261`}, {4.24`, 17.812191700152447`}, {4.62`, 12.660288059282525`}, {4.98`, 9.058552910060195`}, {5.36`, 6.757390846510315`}, {5.72`, 5.406734442782891`}, {6.1`, 4.281064201191851`}, {6.46`, 3.5306087276813973`}, {6.84`, 3.0303969391296968`}, {7.2`, 2.4874083321491103`}, {7.56`, 2.1371492796571525`}, {7.94`, 1.7466912941976835`}, {8.3`, 1.546555052111411`}, {8.68`, 1.3990116021205934`}, {9.04`, 1.1360331931850098`}, {9.42`, 2.5084709290516193`*^-15}, {9.78`, 2.509132768849075`*^-15}, {10.16`, 2.5098586376587673`*^-15}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(dataqdiff1500\ = \ DeleteCases[ Table[\ If[Part[list1500, j, 2]\ >= \ 10^\(-8\), Part[list1500, j]], {j, 1, Length[list1500]}], Null]\)], "Input"], Cell[BoxData[ \({{0.036`, 745000.0367642778`}, {0.052`, 690000.0710428806`}, {0.066`, 512500.0850054065`}, {0.08`, 430000.1047884136`}, {0.096`, 297500.10439851263`}, {0.11`, 183000.0843143265`}, {0.124`, 125750.07362347704`}, {0.14`, 81750.0610110994`}, {0.154`, 53250.04808682784`}, {0.17`, 37750.041541286795`}, {0.184`, 31500.040607982617`}, {0.198`, 25500.038065877827`}, {0.214`, 22775.039714702318`}, {0.228`, 22950.045427396588`}, {0.242`, 20700.046160088485`}, {0.258`, 18000.045622342703`}, {0.26`, 16900.043500983942`}, {0.272`, 16900.047609133846`}, {0.286`, 13475.041968815824`}, {0.302`, 12400.043062731713`}, {0.316`, 11375.043250556`}, {0.33`, 10475.04343583356`}, {0.346`, 9525.043419376216`}, {0.36`, 8775.043303067387`}, {0.374`, 7475.039812649196`}, {0.38`, 8075.044399334685`}, {0.39`, 6950.0402513881845`}, {0.404`, 6975.043348486577`}, {0.412`, 6400.041365806245`}, {0.48`, 4700.041233315392`}, {0.504`, 4025.038930968351`}, {0.58`, 3100.0397088276204`}, {0.596`, 2925.0395628922474`}, {0.68`, 2175.038295540816`}, {0.688`, 2132.538435911745`}, {0.78`, 1697.5393253069099`}, {0.868`, 1252.5359329038924`}, {0.9`, 1200.0370119677887`}, {0.962`, 1047.5369131984269`}, {1.`, 957.5364600543431`}, {1.1`, 755.0347867410987`}, {1.2`, 637.53495644617`}, {1.3`, 517.5333031333787`}, {1.4`, 467.534892287762`}, {1.52`, 380.033432385083`}, {1.62`, 365.03647739440595`}, {1.66`, 335.0351531062393`}, {1.72`, 310.0349240101815`}, {1.82`, 292.53689597944685`}, {1.92`, 245.2844292074305`}, {2.04`, 224.7856190339409`}, {2.4`, 148.7826304895223`}, {2.76`, 99.52886762002373`}, {3.14`, 63.773940880825094`}, {3.5`, 40.268781728792725`}, {3.88`, 26.26505445413261`}, {4.24`, 17.812191700152447`}, {4.62`, 12.660288059282525`}, {4.98`, 9.058552910060195`}, {5.36`, 6.757390846510315`}, {5.72`, 5.406734442782891`}, {6.1`, 4.281064201191851`}, {6.46`, 3.5306087276813973`}, {6.84`, 3.0303969391296968`}, {7.2`, 2.4874083321491103`}, {7.56`, 2.1371492796571525`}, {7.94`, 1.7466912941976835`}, {8.3`, 1.546555052111411`}, {8.68`, 1.3990116021205934`}, {9.04`, 1.1360331931850098`}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(70\)], "Output"] }, Open ]], Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \ \)\(dataq1500plot\ = \ LogLogListPlot[ dataqdiff1500, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.1, 20}, {1, 100000}}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity];\)\)\)], "Input"], Cell[BoxData[ \(\(\(\t\)\(coulomb1500plot\ = \ LogLogPlot[\((icoulombcm\ //. {z\ -> 1.8, en\ -> \ 1500, beta\ -> \ 0.001, thetacm\ -> \ relangle*Pi/180})\), {relangle, 0.1, 20}, \ \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ AxesLabel\ -> \ {"\", \ \ "\"}, PlotRange\ -> \ {{0.1, 10}, {1*^-5, 0.1}}, \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> \ Identity];\)\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[dataq1500plot, coulomb1500plot, DisplayFunction\ -> \ $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.434588 0.434588 0 0.123607 [ [.13082 -0.0125 -11 -12.875 ] [.13082 -0.0125 11 0 ] [.30376 -0.0125 -11 -12.875 ] [.30376 -0.0125 11 0 ] [.43459 -0.0125 -5 -12.875 ] [.43459 -0.0125 5 0 ] [.56541 -0.0125 -5 -12.875 ] [.56541 -0.0125 5 0 ] [.73835 -0.0125 -5 -12.875 ] [.73835 -0.0125 5 0 ] [.86918 -0.0125 -8 -12.875 ] [.86918 -0.0125 8 0 ] [1 -0.0125 -8 -12.875 ] [1 -0.0125 8 0 ] [-0.0125 .12361 -16 -6.4375 ] [-0.0125 .12361 0 6.4375 ] [-0.0125 .24721 -22 -6.4375 ] [-0.0125 .24721 0 6.4375 ] [-0.0125 .37082 -28 -6.4375 ] [-0.0125 .37082 0 6.4375 ] [-0.0125 .49443 -34 -6.4375 ] [-0.0125 .49443 0 6.4375 ] [-0.0125 .61803 -46 -6.4375 ] [-0.0125 .61803 0 6.4375 ] [ -0.0005 -0.0005 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .13082 0 m .13082 .00625 L s gsave .13082 -0.0125 -72 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (0.2) show 81.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 .30376 0 m .30376 .00625 L s gsave .30376 -0.0125 -72 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (0.5) show 81.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 .43459 0 m .43459 .00625 L s gsave .43459 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (1) show 69.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 .56541 0 m .56541 .00625 L s gsave .56541 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (2) show 69.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 .73835 0 m .73835 .00625 L s gsave .73835 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (5) show 69.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 .86918 0 m .86918 .00625 L s gsave .86918 -0.0125 -69 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore 1 0 m 1 .00625 L s gsave 1 -0.0125 -69 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (20) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore .001 w .20735 0 m .20735 .00375 L s .26165 0 m .26165 .00375 L s .33818 0 m .33818 .00375 L s .36727 0 m .36727 .00375 L s .39247 0 m .39247 .00375 L s .4147 0 m .4147 .00375 L s .64194 0 m .64194 .00375 L s .69624 0 m .69624 .00375 L s .77276 0 m .77276 .00375 L s .80186 0 m .80186 .00375 L s .82706 0 m .82706 .00375 L s .84929 0 m .84929 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12361 m .00625 .12361 L s gsave -0.0125 .12361 -77 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 .24721 m .00625 .24721 L s gsave -0.0125 .24721 -83 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (100) show 81.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 .37082 m .00625 .37082 L s gsave -0.0125 .37082 -89 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (1000) show 87.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 .49443 m .00625 .49443 L s gsave -0.0125 .49443 -95 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10000) 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 0 .61803 m .00625 .61803 L s gsave -0.0125 .61803 -107 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (100000.) show 105.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 .001 w 0 .03721 m .00375 .03721 L s 0 .05898 m .00375 .05898 L s 0 .07442 m .00375 .07442 L s 0 .0864 m .00375 .0864 L s 0 .09618 m .00375 .09618 L s 0 .10446 m .00375 .10446 L s 0 .11163 m .00375 .11163 L s 0 .11795 m .00375 .11795 L s 0 .16082 m .00375 .16082 L s 0 .18258 m .00375 .18258 L s 0 .19803 m .00375 .19803 L s 0 .21 m .00375 .21 L s 0 .21979 m .00375 .21979 L s 0 .22807 m .00375 .22807 L s 0 .23523 m .00375 .23523 L s 0 .24156 m .00375 .24156 L s 0 .28442 m .00375 .28442 L s 0 .30619 m .00375 .30619 L s 0 .32163 m .00375 .32163 L s 0 .33361 m .00375 .33361 L s 0 .3434 m .00375 .3434 L s 0 .35167 m .00375 .35167 L s 0 .35884 m .00375 .35884 L s 0 .36516 m .00375 .36516 L s 0 .40803 m .00375 .40803 L s 0 .4298 m .00375 .4298 L s 0 .44524 m .00375 .44524 L s 0 .45722 m .00375 .45722 L s 0 .46701 m .00375 .46701 L s 0 .47528 m .00375 .47528 L s 0 .48245 m .00375 .48245 L s 0 .48877 m .00375 .48877 L s 0 .53164 m .00375 .53164 L s 0 .5534 m .00375 .5534 L s 0 .56885 m .00375 .56885 L s 0 .58082 m .00375 .58082 L s 0 .59061 m .00375 .59061 L s 0 .59889 m .00375 .59889 L s 0 .60606 m .00375 .60606 L s 0 .61238 m .00375 .61238 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .06351 .60722 Mdot .08149 .58421 Mdot .10015 .56574 Mdot .11509 .55602 Mdot .12893 .54468 Mdot .14359 .53861 Mdot .15555 .53902 Mdot .1668 .53348 Mdot .17888 .52598 Mdot .18034 .5226 Mdot .18886 .5226 Mdot .19833 .51044 Mdot .20861 .50597 Mdot .21716 .50134 Mdot .22534 .49692 Mdot .23428 .49182 Mdot .24176 .48741 Mdot .24896 .4788 Mdot .25197 .48295 Mdot .25687 .4749 Mdot .26353 .47509 Mdot .26723 .47047 Mdot .29606 .4539 Mdot .30527 .44557 Mdot .33178 .43156 Mdot .33691 .42844 Mdot .3618 .41253 Mdot .36401 .41147 Mdot .38769 .39923 Mdot .40787 .38291 Mdot .4147 .38061 Mdot .42728 .37331 Mdot .43459 .36849 Mdot .45258 .35574 Mdot .469 .34666 Mdot .48411 .33546 Mdot .49809 .33001 Mdot .51362 .31888 Mdot .52564 .31672 Mdot .53024 .31212 Mdot .53695 .30796 Mdot .54761 .30484 Mdot .55771 .29538 Mdot .56915 .29069 Mdot .59982 .26854 Mdot .6262 .24696 Mdot .65055 .22307 Mdot .67103 .19839 Mdot .69049 .17544 Mdot .70723 .1546 Mdot .72343 .13627 Mdot .7376 .1183 Mdot .75147 .10257 Mdot .76374 .0906 Mdot .77588 .07806 Mdot .78671 .06772 Mdot .79749 .05952 Mdot .80717 .04892 Mdot .81638 .04077 Mdot .82564 .02994 Mdot .83401 .02341 Mdot .84246 .01802 Mdot .85013 .00685 Mdot 0 g .5 Mabswid 0 .22004 m .04146 .22 L .0726 .21995 L .10314 .21988 L .12807 .21981 L .16561 .21965 L .18391 .21955 L .19936 .21945 L .22818 .21921 L .25555 .2189 L .28043 .21854 L .3011 .21815 L .3373 .21725 L .36526 .21627 L .3931 .21498 L .41622 .2136 L .45691 .21024 L .48816 .20659 L .51738 .20207 L .54423 .19674 L .58639 .18551 L .61899 .17395 L .66898 .15051 L .71071 .12525 L .74371 .10168 L .77078 .08017 L .79589 .05867 L .81725 .03934 L .8376 .02016 L .85531 .00293 L s .85531 .00293 m .85826 0 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg_ooo`00oooo000`oooo00<0003oooooool0Roooo`80000=oooo1@2codkoool0 033oool00`000?ooooooo`27oooo100000ooool50;?oC_ooo`00 oooo1P2coiWoool0033oool00`000?ooooooo`1>oooo1P2coiWoool0033oool00`000?ooooooo`18 oooo1@2co`;oool40;?oV_ooo`00"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-1.48874, -0.773956, \ 0.00999921, 0.0351561}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \ \)\(ListPlot[ dataqdiff1500, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.3, 0.6}, {0, 10000}} (*\(,\)\(\n\)\(\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \)\(DisplayFunction\ -> Identity\)*) ];\)\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -1 3.33333 -1.6263e-018 6.18034e-005 [ [.16667 -0.0125 -12 -9 ] [.16667 -0.0125 12 0 ] [.33333 -0.0125 -9 -9 ] [.33333 -0.0125 9 0 ] [.5 -0.0125 -12 -9 ] [.5 -0.0125 12 0 ] [.66667 -0.0125 -9 -9 ] [.66667 -0.0125 9 0 ] [.83333 -0.0125 -12 -9 ] [.83333 -0.0125 12 0 ] [-0.0125 .12361 -24 -4.5 ] [-0.0125 .12361 0 4.5 ] [-0.0125 .24721 -24 -4.5 ] [-0.0125 .24721 0 4.5 ] [-0.0125 .37082 -24 -4.5 ] [-0.0125 .37082 0 4.5 ] [-0.0125 .49443 -24 -4.5 ] [-0.0125 .49443 0 4.5 ] [-0.0125 .61803 -30 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .16667 0 m .16667 .00625 L s [(0.35)] .16667 -0.0125 0 1 Mshowa .33333 0 m .33333 .00625 L s [(0.4)] .33333 -0.0125 0 1 Mshowa .5 0 m .5 .00625 L s [(0.45)] .5 -0.0125 0 1 Mshowa .66667 0 m .66667 .00625 L s [(0.5)] .66667 -0.0125 0 1 Mshowa .83333 0 m .83333 .00625 L s [(0.55)] .83333 -0.0125 0 1 Mshowa .125 Mabswid .03333 0 m .03333 .00375 L s .06667 0 m .06667 .00375 L s .1 0 m .1 .00375 L s .13333 0 m .13333 .00375 L s .2 0 m .2 .00375 L s .23333 0 m .23333 .00375 L s .26667 0 m .26667 .00375 L s .3 0 m .3 .00375 L s .36667 0 m .36667 .00375 L s .4 0 m .4 .00375 L s .43333 0 m .43333 .00375 L s .46667 0 m .46667 .00375 L s .53333 0 m .53333 .00375 L s .56667 0 m .56667 .00375 L s .6 0 m .6 .00375 L s .63333 0 m .63333 .00375 L s .7 0 m .7 .00375 L s .73333 0 m .73333 .00375 L s .76667 0 m .76667 .00375 L s .8 0 m .8 .00375 L s .86667 0 m .86667 .00375 L s .9 0 m .9 .00375 L s .93333 0 m .93333 .00375 L s .96667 0 m .96667 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12361 m .00625 .12361 L s [(2000)] -0.0125 .12361 1 0 Mshowa 0 .24721 m .00625 .24721 L s [(4000)] -0.0125 .24721 1 0 Mshowa 0 .37082 m .00625 .37082 L s [(6000)] -0.0125 .37082 1 0 Mshowa 0 .49443 m .00625 .49443 L s [(8000)] -0.0125 .49443 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(10000)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .0309 m .00375 .0309 L s 0 .0618 m .00375 .0618 L s 0 .09271 m .00375 .09271 L s 0 .15451 m .00375 .15451 L s 0 .18541 m .00375 .18541 L s 0 .21631 m .00375 .21631 L s 0 .27812 m .00375 .27812 L s 0 .30902 m .00375 .30902 L s 0 .33992 m .00375 .33992 L s 0 .40172 m .00375 .40172 L s 0 .43262 m .00375 .43262 L s 0 .46353 m .00375 .46353 L s 0 .52533 m .00375 .52533 L s 0 .55623 m .00375 .55623 L s 0 .58713 m .00375 .58713 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .15333 .58868 Mdot .2 .54233 Mdot .24667 .46198 Mdot .26667 .49907 Mdot .3 .42954 Mdot .34667 .43108 Mdot .37333 .39554 Mdot .6 .29048 Mdot .68 .24876 Mdot .93333 .19159 Mdot .98667 .18078 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgWoool40;?o2?ooo`H0/ol2oooo000Qoooo00<0003ooooo ool0ioooo`H0/ol7oooo1P2co`;oool0027oool00`000?ooooooo`3Woooo1P2co`Soool40;?o0ooo o`008Oooo`030000oooooooo0>Ooool60;?o3oooo`008Oooo`030000oooooooo0>Ooool60;?o3ooo o`008Oooo`030000oooooooo0>Soool40;?o4?ooo`008Oooo`030000oooooooo0?coool0027oool0 0`000?ooooooo`3loooo000Qoooo0P000?goool0027oool00`000?ooooooo`3loooo000Qoooo00<0 003oooooool0o?ooo`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo000Q oooo00<0003oooooool0o?ooo`002Oooo`<00002oooo0P0000Coool200001?ooo`800005oooo00<0 003oooooool0Z?ooo`@0/om@oooo000:oooo00@0003oooooool00002oooo00@0003oooooool00002 oooo00@0003oooooool00002oooo00<0003oooooool00_ooo`030000oooooooo0:Ooool60;?oCooo o`001oooo`D000001Oooo`000?ooooooo`000002oooo00@0003oooooool00002oooo00@0003ooooo ool00004oooo0P000:Soool60;?oCoooo`001oooo`040000oooooooo00000_ooo`040000oooooooo 00000_ooo`040000oooooooo00000_ooo`040000oooooooo00001?ooo`030000oooooooo0:Ooool6 0;?oCoooo`002?ooo`030000oooo000000;oool010000?ooooooo`0000;oool010000?ooooooo`00 00;oool010000?ooooooo`0000Coool00`000?ooooooo`2Woooo1P2codooool000Woool200000ooo o`800004oooo0P0000Coool200001Oooo`030000oooooooo0:Soool40;?oD?ooo`008Oooo`030000 oooooooo0?coool0027oool00`000?ooooooo`3loooo000Qoooo00<0003oooooool0o?ooo`008Ooo o`030000oooooooo0?coool0027oool20000oOooo`008Oooo`030000oooooooo09?oool40;?oIOoo o`008Oooo`030000oooooooo09;oool60;?oI?ooo`008Oooo`030000oooooooo09;oool60;?oI?oo o`008Oooo`030000oooooooo09;oool60;?oI?ooo`008Oooo`030000oooooooo09;oool60;?oI?oo o`008Oooo`030000oooooooo09?oool40;?oIOooo`008Oooo`030000oooooooo0?coool0027oool2 0000oOooo`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo000Qoooo00<0 003oooooool0o?ooo`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo000Q oooo00<0003oooooool0o?ooo`008Oooo`030000oooooooo0?coool0027oool20000oOooo`008Ooo o`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo000Qoooo00<0003oooooool0o?oo o`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo0007oooo0`0000Coool2 00001?ooo`800004oooo0P0000Goool00`000?ooooooo`3loooo0007oooo00@0003oooooool00002 oooo00@0003oooooool00002oooo00@0003oooooool00002oooo00@0003oooooool00004oooo00<0 003oooooool0o?ooo`001oooo`040000oooooooo00000_ooo`040000oooooooo00000_ooo`040000 oooooooo00000_ooo`040000oooooooo00001?ooo`80003moooo0007oooo0`0000?oool010000?oo ooooo`0000;oool010000?ooooooo`0000;oool010000?ooooooo`0000Coool00`000?ooooooo`3l oooo0008oooo00<0003oooooool00_ooo`040000oooooooo00000_ooo`040000oooooooo00000_oo o`040000oooooooo00001?ooo`030000oooooooo0?coool000Soool300000oooo`800004oooo0P00 00Coool200001Oooo`030000oooooooo05[oool40;?oW_ooo`008Oooo`030000oooooooo05Woool6 0;?oWOooo`008Oooo`030000oooooooo05Woool60;?oWOooo`008Oooo`030000oooooooo05Woool6 0;?oWOooo`008Oooo`80001Joooo1P2coigoool0027oool00`000?ooooooo`1Joooo102coikoool0 027oool00`000?ooooooo`3loooo000Qoooo00<0003oooooool0o?ooo`008Oooo`030000oooooooo 0?coool0027oool00`000?ooooooo`17oooo102co`Soool40;?oYOooo`008Oooo`030000oooooooo 04Koool60;?o1_ooo`H0/onToooo000Qoooo00<0003oooooool0A_ooo`H0/ol6oooo1P2cojCoool0 027oool20000Aoooo`H0/ol6oooo1P2cojCoool0027oool00`000?ooooooo`16oooo1P2co`Koool6 0;?oY?ooo`008Oooo`030000oooooooo04Ooool40;?o2?ooo`@0/onUoooo000Qoooo00<0003ooooo ool0o?ooo`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`0joooo102cokkoool0 027oool00`000?ooooooo`0ioooo1P2cokgoool0027oool00`000?ooooooo`0ioooo1P2cokgoool0 027oool20000>_ooo`H0/onmoooo000Qoooo00<0003oooooool0>Oooo`H0/onmoooo000Qoooo00<0 003oooooool0>_ooo`@0/onnoooo000Qoooo00<0003oooooool0o?ooo`008Oooo`030000oooooooo 0?coool0027oool00`000?ooooooo`3loooo0008oooo0P0000Coool200001?ooo`800004oooo0P00 00Goool00`000?ooooooo`3loooo0007oooo00@0003oooooool00002oooo00@0003oooooool00002 oooo00@0003oooooool00002oooo00@0003oooooool00004oooo00<0003oooooool0?oooo`@0/oni oooo0007oooo00@0003oooooool00002oooo00@0003oooooool00002oooo00@0003oooooool00002 oooo00@0003oooooool00004oooo0P0003ooool60;?o^?ooo`002?ooo`800003oooo00@0003ooooo ool00002oooo00@0003oooooool00002oooo00@0003oooooool00004oooo00<0003oooooool0?_oo o`H0/onhoooo0007oooo00@0003oooooool00002oooo00@0003oooooool00002oooo00@0003ooooo ool00002oooo00@0003oooooool00004oooo00<0003oooooool0?_ooo`H0/onhoooo0008oooo0P00 00Coool200001?ooo`800004oooo0P0000Goool00`000?ooooooo`0noooo1P2cokSoool0027oool0 0`000?ooooooo`0ooooo102cokWoool0027oool00`000?ooooooo`3loooo000Qoooo00<0003ooooo ool0o?ooo`008Oooo`030000oooooooo0?coool0027oool20000oOooo`008Oooo`030000oooooooo 0?coool0027oool00`000?ooooooo`0^oooo102col[oool0027oool00`000?ooooooo`0]oooo1P2c olWoool0027oool00`000?ooooooo`0]oooo1P2colWoool0027oool00`000?ooooooo`0]oooo1P2c olWoool0027oool00`000?ooooooo`0]oooo1P2colWoool0027oool00`000?ooooooo`0^oooo102c ol[oool0027oool20000oOooo`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3l oooo000Qoooo00<0003oooooool0o?ooo`008Oooo`030000oooooooo0?coool0027oool00`000?oo ooooo`0Roooo102comKoool0027oool00`000?ooooooo`0Qoooo1P2comGoool0027oool200008_oo o`H0/ooEoooo000Qoooo00<0003oooooool08Oooo`H0/ooEoooo000Qoooo00<0003oooooool08Ooo o`H0/ooEoooo000Qoooo00<0003oooooool08_ooo`@0/ooFoooo000Qoooo00<0003oooooool0o?oo o`008Oooo`030000oooooooo0?coool00003oooo0000000000<00002oooo0P0000Coool200001?oo o`800004oooo0P0000Goool00`000?ooooooo`3loooo0003oooo00D0003oooooooooool000000_oo o`040000oooooooo00000_ooo`040000oooooooo00000_ooo`040000oooooooo00000_ooo`030000 oooooooo00;oool00`000?ooooooo`3loooo0003oooo00D0003oooooooooool000000_ooo`040000 oooooooo00000_ooo`040000oooooooo00000_ooo`040000oooooooo00000_ooo`030000oooooooo 00;oool20000oOooo`000oooo`050000oooooooooooo000000;oool010000?ooooooo`0000;oool0 10000?ooooooo`0000;oool010000?ooooooo`0000;oool00`000?ooooooo`3ooooo0_ooo`000ooo o`050000oooooooooooo000000;oool010000?ooooooo`0000;oool010000?ooooooo`0000;oool0 10000?ooooooo`0000;oool00`000?ooooooo`3ooooo0_ooo`000_ooo`800004oooo0P0000Coool2 00001?ooo`800004oooo0P000?ooool5oooo003ooooo8Oooo`00ooooob7oool00?oooolQoooo0000 \ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0.260791, -884.13, \ 0.00118192, 63.7461}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qdiff1500int\ = \ Interpolation[dataqdiff1500]\)], "Input"], Cell[BoxData[ TagBox[\(InterpolatingFunction[{{0.036`, 9.04`}}, "<>"]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(norm\ = \ 2*Pi*NIntegrate[Sin[relangle*3.1416/180], \ {relangle, 0, 180}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(12.566370614189621`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qt1500\ = \ 2*Pi*NIntegrate[ qdiff1500int[relangle]*Sin[relangle*3.1416/180], \ {relangle, 0.036, 9.04}]*3.1416/180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.2881485954920154`\)."\)], "Message"], Cell[BoxData[ \(11.853159672012122`\)], "Output"] }, Open ]], Cell["\<\ This value is about 15% larger than that given by Newman et al. (1985), but \ is only 7% larger than the experimental value of Savola (1973) as obtained \ using the detector angle of 0.056 deg. This difference may be the result of \ the improved differential scattering data available at low angles from the \ experiments of Nitz as used in our calculation but not in Newman's \ calculation.\ \>", "SmallText"], Cell[CellGroupData[{ Cell[BoxData[ \(qm1500\ = \ 2*Pi*NIntegrate[ qdiff1500int[relangle]*\((1 - Cos[relangle*3.1416/180])\)* Sin[relangle*3.1416/ 180], \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {relangle, 0.036, 9.04}]*3.1416/180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(2.3925156249999993`\)."\)], "Message"], Cell[BoxData[ \(0.0031648735707903827`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qv1500\ = \ 2*Pi*NIntegrate[ qdiff1500int[ relangle]*\((Sin[ relangle*3.1416/ 180])\)^3, \ \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \t{relangle, 0.036, 9.04}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(2.3925156249999993`\)."\)], "Message"], Cell[BoxData[ \(0.00631847953384979`\)], "Output"] }, Open ]], Cell[BoxData[ \(\(\t\(coulombint1500plot\ = \ LogLogPlot[ \((\((icoulombcm*\((Sin[thetacm])\)^3)\)\ //. {z\ -> 1.8, en\ -> \ 1500, beta\ -> \ 0.001, thetacm\ -> \ relangle*Pi/180})\), {relangle, 0.1, 20}, \ \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ AxesLabel\ -> \ {"\", \ "\"}, PlotRange\ -> \ {{0.1, 10}, {1*^-5, 0.01}}, \n \t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> \ Identity]; \)\)\)], "Input", Evaluatable->False], Cell[BoxData[ \(\(\(\t\)\(coulombint1500plot\ = \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*linear\ plot*) \n\t Plot[\((\((icoulombcm*\((Sin[thetacm])\)^3)\)\ //. {z\ -> 1.8, en\ -> \ 1500, beta\ -> \ 0.001, thetacm\ -> \ relangle*Pi/180})\), {relangle, 0, 90}, \ \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ AxesLabel\ -> \ {"\", \ \ "\"}, PlotRange\ -> \ {{0, 90}, {0, 0.01}}, \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> \ Identity];\)\)\)], "Input"], Cell[BoxData[ \(\(angleqdiff1500\ = \ Part[Transpose[dataqdiff1500], 1];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(70\)], "Output"] }, Open ]], Cell[BoxData[ \(\(listqdiff1500\ = \ \ Part[Transpose[dataqdiff1500], 2];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(70\)], "Output"] }, Open ]], Cell[BoxData[ \(\(integrandqv1500\ = \ Table[\((Part[listqdiff1500, j]*\((Sin[ Part[angleqdiff1500, j]* Pi/180])\)^3)\), \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {j, 1, Length[angleqdiff1500]}];\)\)], "Input"], Cell[BoxData[ \(\(tableintegrandqv1500\ = \ Transpose[{angleqdiff1500, integrandqv1500}];\)\)], "Input"], Cell[BoxData[ \(\(\ \ \ \ \ \ \ \ \ \ \ \(integrand1500plot\ = \ LogLogListPlot[tableintegrandqv1500, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.1, 20}, {1*^-5, 0.01}}, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity]; \)\n\)\)], "Input", Evaluatable->False], Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \)\(integrand1500plot\ = \ ListPlot[ tableintegrandqv1500, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*linear\ plot*) \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0, 90}, {0, 0.01}}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity];\)\(\n\)\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[integrand1500plot, coulombint1500plot, \ DisplayFunction\ -> \ $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 -1.9082e-017 0.0111111 9.1073e-018 61.8034 [ [.22222 -0.0125 -6 -9 ] [.22222 -0.0125 6 0 ] [.44444 -0.0125 -6 -9 ] [.44444 -0.0125 6 0 ] [.66667 -0.0125 -6 -9 ] [.66667 -0.0125 6 0 ] [.88889 -0.0125 -6 -9 ] [.88889 -0.0125 6 0 ] [-0.0125 .12361 -30 -4.5 ] [-0.0125 .12361 0 4.5 ] [-0.0125 .24721 -30 -4.5 ] [-0.0125 .24721 0 4.5 ] [-0.0125 .37082 -30 -4.5 ] [-0.0125 .37082 0 4.5 ] [-0.0125 .49443 -30 -4.5 ] [-0.0125 .49443 0 4.5 ] [-0.0125 .61803 -24 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .22222 0 m .22222 .00625 L s [(20)] .22222 -0.0125 0 1 Mshowa .44444 0 m .44444 .00625 L s [(40)] .44444 -0.0125 0 1 Mshowa .66667 0 m .66667 .00625 L s [(60)] .66667 -0.0125 0 1 Mshowa .88889 0 m .88889 .00625 L s [(80)] .88889 -0.0125 0 1 Mshowa .125 Mabswid .05556 0 m .05556 .00375 L s .11111 0 m .11111 .00375 L s .16667 0 m .16667 .00375 L s .27778 0 m .27778 .00375 L s .33333 0 m .33333 .00375 L s .38889 0 m .38889 .00375 L s .5 0 m .5 .00375 L s .55556 0 m .55556 .00375 L s .61111 0 m .61111 .00375 L s .72222 0 m .72222 .00375 L s .77778 0 m .77778 .00375 L s .83333 0 m .83333 .00375 L s .94444 0 m .94444 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12361 m .00625 .12361 L s [(0.002)] -0.0125 .12361 1 0 Mshowa 0 .24721 m .00625 .24721 L s [(0.004)] -0.0125 .24721 1 0 Mshowa 0 .37082 m .00625 .37082 L s [(0.006)] -0.0125 .37082 1 0 Mshowa 0 .49443 m .00625 .49443 L s [(0.008)] -0.0125 .49443 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(0.01)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .0309 m .00375 .0309 L s 0 .0618 m .00375 .0618 L s 0 .09271 m .00375 .09271 L s 0 .15451 m .00375 .15451 L s 0 .18541 m .00375 .18541 L s 0 .21631 m .00375 .21631 L s 0 .27812 m .00375 .27812 L s 0 .30902 m .00375 .30902 L s 0 .33992 m .00375 .33992 L s 0 .40172 m .00375 .40172 L s 0 .43262 m .00375 .43262 L s 0 .46353 m .00375 .46353 L s 0 .52533 m .00375 .52533 L s 0 .55623 m .00375 .55623 L s 0 .58713 m .00375 .58713 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .0004 .01142 Mdot .00058 .03188 Mdot .00073 .04841 Mdot .00089 .07234 Mdot .00107 .08649 Mdot .00122 .08003 Mdot .00138 .07878 Mdot .00156 .07371 Mdot .00171 .0639 Mdot .00189 .06094 Mdot .00204 .06448 Mdot .0022 .06504 Mdot .00238 .07334 Mdot .00253 .08938 Mdot .00269 .0964 Mdot .00287 .10157 Mdot .00289 .0976 Mdot .00302 .11175 Mdot .00318 .10358 Mdot .00336 .11222 Mdot .00351 .11794 Mdot .00367 .12369 Mdot .00384 .12964 Mdot .004 .13452 Mdot .00416 .12849 Mdot .00422 .14559 Mdot .00433 .13546 Mdot .00449 .15112 Mdot .00458 .14706 Mdot .00533 .17079 Mdot .0056 .16931 Mdot .00644 .19873 Mdot .00662 .20347 Mdot .00756 .2247 Mdot .00764 .22818 Mdot .00867 .26467 Mdot .00964 .26912 Mdot .01 .28742 Mdot .01069 .30639 Mdot .01111 .31458 Mdot .01222 .33015 Mdot .01333 .36191 Mdot .01444 .37351 Mdot .01556 .42142 Mdot .01689 .43837 Mdot .018 .50974 Mdot .01844 .50336 Mdot .01911 .51814 Mdot .02022 .57919 Mdot .02133 .57013 Mdot .03889 .56625 Mdot .04311 .50295 Mdot .04711 .44491 Mdot .05133 .40889 Mdot .05533 .36623 Mdot .05956 .34042 Mdot .06356 .33083 Mdot .06778 .31749 Mdot .07178 .31076 Mdot .076 .31639 Mdot .08 .30266 Mdot .084 .30079 Mdot .08822 .28454 Mdot .09222 .28753 Mdot .09644 .29719 Mdot .10044 .27235 Mdot 0 g .5 Mabswid 0 0 m .00123 3e-005 L .00236 .00019 L .00365 .00069 L .00488 .00163 L .00601 .00301 L .00706 .00479 L .00821 .00737 L .00942 .01088 L .01196 .02089 L .01434 .03356 L .01975 .07234 L .03042 .16483 L .03568 .20584 L .04057 .23735 L .04589 .26385 L .0487 .2747 L .0517 .28408 L .05445 .29089 L .05704 .29589 L .05936 .29933 L .06188 .30206 L .06331 .3032 L .06465 .30401 L .0654 .30436 L .06612 .30464 L .06694 .30488 L .06769 .30504 L .06843 .30514 L .06921 .30518 L .06989 .30517 L .07063 .3051 L .07142 .30498 L .07227 .30479 L .07308 .30456 L .07382 .3043 L .07638 .30309 L .07917 .30131 L .08481 .29652 L .1054 .27194 L .1274 .24318 L .17097 .19467 L .21295 .16014 L .25334 .13524 L .29631 .11488 L .33768 .09943 L .37746 .08732 L .41982 .07661 L .46058 .06792 L .50392 .06001 L Mistroke .54567 .05343 L .58583 .04787 L .62856 .04263 L .6697 .03815 L .71342 .0339 L .75554 .03022 L .79608 .02703 L .83918 .02395 L .8807 .02127 L .92063 .01893 L .96313 .01666 L 1 .01487 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg7oool0023oool60;?o1oooo`030000oooooooo00Soool30;?o00<00000/ol0/ol0h_ooo`008?oo o`H0/ol7oooo00<0003oooooool01Oooo`D0/ol01000002co`2co`2con;oool0027oool50;?o1ooo o`030000oooooooo00Coool60;?o00@00000/ol0/ol0/ooRoooo000Qoooo1P2co`Ooool00`000?oo ooooo`03oooo1@2co`0400000;?o0;?o0;?ohoooo`008Oooo`H0/ol7oooo00<0003oooooool00ooo o`@0/ol01000002co`2co`2conCoool0027oool60;?o1oooo`030000oooooooo00?oool30;?o00<0 0000/ol0/ol00`2con?oool0027oool60;?o2?ooo`040000oooooooooooo102co`0300000;?o0;?o 00<0/ooSoooo000Qoooo00<00000/ol0/ol00P2co`Woool2000000?oool0/ol0/ol00P2co`030000 0;?o0;?o00@0/ooSoooo000Qoooo1P2co`Woool200000P2co`8000070;?ohoooo`008Oooo`H0/ol9 oooo00@0/ol00000000000080;?oi?ooo`008Oooo`H0/ol9oooo202conSoool0027oool60;?o2?oo o`P0/ooYoooo000Qoooo1P2co`Soool80;?ojOooo`008Oooo`H0/ol8oooo202conWoool0027oool6 0;?o1oooo`P0/ooZoooo000Qoooo1P2co`Ooool60;?ok?ooo`008Oooo`H0/ol6oooo1`2concoool0 027oool60;?o1_ooo`L0/oo/oooo000Qoooo00<00000/ol0/ol00P2co`Ooool60;?okOooo`008Ooo o`030000oooooooo00Woool60;?okOooo`008Oooo`030000oooooooo00[oool40;?ok_ooo`008Ooo o`0300000;?o0;?o0080/oojoooo000Qoooo1P2co`Koool40;?okoooo`008Oooo`H0/ol5oooo1P2c onkoool000;oool200001?ooo`800004oooo0P0000Coool200000oooo`<00005oooo1P2co`Goool6 0;?ok_ooo`0000Goool0003oooooool000002?ooo`040000oooooooo00000_ooo`040000oooooooo 00000_ooo`040000oooooooo00001?ooo`L0/ol4oooo1P2conkoool00005oooo0000oooooooo0000 00Soool010000?ooooooo`0000;oool010000?ooooooo`0000;oool010000?ooooooo`0000Coool0 0`00002co`2co`040;?o1?ooo`H0/oo^oooo00001Oooo`000?ooooooo`000008oooo00@0003ooooo ool00002oooo00@0003oooooool00002oooo0`0000Goool00`00002co`2co`040;?o1Oooo`@0/oo_ oooo00001Oooo`000?ooooooo`000008oooo00@0003oooooool00002oooo00@0003oooooool00003 oooo00<0003oooooool01?ooo`0300000;?o0;?o00@0/oohoooo0002oooo0P0000[oool200001?oo o`800004oooo0`0000Coool00`000?ooo`2co`030;?onOooo`008Oooo`030000oooooooo0?coool0 027oool00`000?ooooooo`3loooo000Qoooo00<0003oooooool0o?ooo`008Oooo`800009oooo102c oo3oool0027oool00`000?ooooooo`07oooo1P2conooool0027oool00`000?ooooooo`07oooo1P2c onooool0027oool00`000?ooo`2co`030;?o1?ooo`H0/oo_oooo000Qoooo00<00000/ol0/ol0102c o`?oool60;?okoooo`008Oooo`0300000;?o0;?o00@0/ol4oooo102coo3oool0027oool00`00002c o`2co`040;?on?ooo`008Oooo`0300000;?o0;?o00@0/oohoooo000Qoooo00<00000/ol0/ol0102c ooSoool0027oool00`00002co`2co`040;?o0oooo`@0/ooaoooo000Qoooo00<00000/ol0/ol0102c o`;oool60;?ol?ooo`008Oooo`0300000;?o0;?o00@0/ol2oooo1P2coo3oool0027oool00`000?oo o`2co`030;?o0oooo`H0/oo`oooo000Qoooo00<0003oooooool01_ooo`H0/oo`oooo000Qoooo00<0 003oooooool01oooo`@0/ooaoooo000Qoooo00<0003oooooool0o?ooo`008Oooo`80003moooo000Q oooo00<0003oooooool0o?ooo`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3l oooo000Qoooo00<0003oooooool0o?ooo`008Oooo`030000oooooooo0?coool000;oool200001?oo o`800004oooo0P0000Coool200001?ooo`800005oooo00<0003oooooool0o?ooo`0000Goool0003o ooooool000002?ooo`040000oooooooo00000_ooo`040000oooooooo00000_ooo`040000oooooooo 00001?ooo`030000oooooooo0?coool00005oooo0000oooooooo000000Soool010000?ooooooo`00 00;oool010000?ooooooo`0000;oool010000?ooooooo`0000Coool2000000?oool0/ol0/ol00P2c o`;oool40;?ol_ooo`0000Goool0003oooooool000002?ooo`040000oooooooo00000_ooo`040000 oooooooo00000oooo`800005oooo00<0003oool0/ol02`2coo7oool00005oooo0000oooooooo0000 00Soool010000?ooooooo`0000;oool010000?ooooooo`0000;oool010000?ooooooo`0000Coool0 0`000?ooo`2co`0;0;?olOooo`000_ooo`80000:oooo0P0000Coool200001?ooo`800005oooo00<0 003oool0/ol02`2coo7oool0027oool00`000?ooo`2co`0;0;?olOooo`008Oooo`030000oooo0;?o 00D0/ol00oooo`2co`2co`020;?ol_ooo`008Oooo`030000oooo0;?o00D0/oogoooo000Qoooo00<0 003oool0/ol01@2cooOoool0027oool200001P2cooOoool0027oool00`000?ooooooo`040;?on?oo o`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo000Qoooo00<0003ooooo ool0o?ooo`008Oooo`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo000Qoooo00<0 003oooooool0o?ooo`008Oooo`800006oooo102coo?oool0027oool00`000?ooooooo`0:0;?ol_oo o`008Oooo`030000oooo0;?o00X0/ooboooo000Qoooo00<0003oool0/ol02P2coo;oool0027oool0 0`000?ooo`2co`0:0;?ol_ooo`008Oooo`030000oooo0;?o00T0/oocoooo000Qoooo00<0003oool0 /ol01@2cooOoool0027oool200001P2cooOoool0027oool00`000?ooooooo`040;?on?ooo`008Ooo o`030000oooooooo0?coool0027oool00`000?ooooooo`3loooo000Qoooo00<0003oooooool0o?oo o`008Oooo`030000oooooooo0?coool000Soool200001?ooo`800004oooo0P0000?oool500000ooo o`030000oooooooo0?coool000Ooool010000?ooooooo`0000Soool010000?ooooooo`0000Coool0 0`000?ooooooo`03oooo00<0003oooooool0o?ooo`001oooo`040000oooooooo00002?ooo`040000 oooooooo00001?ooo`030000oooooooo00?oool20000oOooo`001oooo`040000oooooooo00002?oo o`040000oooooooo00001?ooo`030000oooooooo0?ooool3oooo0007oooo00@0003oooooool00008 oooo00@0003oooooool00004oooo00<0003oooooool0ooooo`?oool000Soool200002_ooo`800004 oooo0P000?ooool5oooo003ooooo8Oooo`00ooooob7oool00?oooolQoooo0000\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-11.7627, -0.000884129, \ 0.354576, 6.37461*^-05}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqt1500\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 1.8, en\ -> \ 1500, beta\ -> \ 0.001, thetacm\ -> \ relangle* Pi/180})\)\ \ *\ \ \n\t\t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ Sin[relangle*3.1416/180], {relangle, 9.04, 180}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.1044576022898847`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqm1500\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 1.8, en\ -> \ 1500, beta\ -> \ 0.001, thetacm\ -> \ relangle*Pi/180})\)\ \ *\ \ \((1 - Cos[relangle*3.1416/180])\)* Sin[relangle*3.1416/ 180], \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {relangle, 9.04, 180}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.007276385240959193`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqv1500\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 1.8, en\ -> \ 1500, beta\ -> \ 0.001, thetacm\ -> \ relangle*Pi/180})\)*\((Sin[ relangle*3.1416/180])\)^2* Sin[relangle*3.1416/ 180], \ \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \t{relangle, 9.04, 180}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.01156648473936736`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqv1500\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 1.8, en\ -> \ 1500, beta\ -> \ 0.001, thetacm\ -> \ relangle*Pi/180})\)*\((Sin[ relangle*3.1416/180])\)^2* Sin[relangle*3.1416/ 180], \ \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \t{relangle, 0, 180}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.014938332543900024`\)], "Output"] }, Open ]], Cell["Analysis of 5000 eV data", "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ \(list5000\ = \ Union[Drop[Drop[data1qdiff5000, {6}], {10}], data2qdiff5000]\)], "Input"], Cell[BoxData[ \({{0.036`, 1.0475000516920552`*^6}, {0.052`, 670000.0689836667`}, {0.066`, 330000.0547351886`}, {0.08`, 213500.0520286658`}, {0.096`, 103000.03614469513`}, {0.11`, 87500.04031422714`}, {0.124`, 69000.04039777268`}, {0.14`, 61250.045711679995`}, {0.154`, 47500.04289435347`}, {0.17`, 40000.044017257525`}, {0.184`, 33750.04350855281`}, {0.198`, 28500.04254421639`}, {0.214`, 23700.0413277034`}, {0.228`, 20350.040280937716`}, {0.242`, 20075.044766366005`}, {0.258`, 16550.041947209542`}, {0.26`, 14700.037838133961`}, {0.272`, 14850.0418340614`}, {0.286`, 13350.041579494713`}, {0.302`, 12600.043757291902`}, {0.316`, 11875.04515167934`}, {0.33`, 11025.045716473986`}, {0.346`, 8850.040342412547`}, {0.36`, 7975.039355209391`}, {0.374`, 7650.040744717905`}, {0.38`, 7275.04000063899`}, {0.39`, 6575.0380795506935`}, {0.404`, 6800.04226089014`}, {0.412`, 2.500016158518065`*^-10}, {0.426`, 5450.037660228396`}, {0.456`, 2.5000197941705397`*^-10}, {0.48`, 3925.034434204875`}, {0.486`, 4075.036649508495`}, {0.504`, 2.5000241807256835`*^-10}, {0.514`, 2.500025149805785`*^-10}, {0.544`, 2900.032678688883`}, {0.574`, 2.500031364118212`*^-10}, {0.58`, 2325.0297816207153`}, {0.596`, 2.500033814437818`*^-10}, {0.602`, 2187.5301863627433`}, {0.632`, 2.500038022828243`*^-10}, {0.662`, 1520.0253647364577`}, {0.68`, 1362.5239897353388`}, {0.688`, 2.500045059685516`*^-10}, {0.72`, 1410.0278327422413`}, {0.75`, 2.5000535469865643`*^-10}, {0.78`, 2.500057916505022`*^-10}, {0.808`, 2.5000621493334707`*^-10}, {0.838`, 2.500066850163158`*^-10}, {0.868`, 2.5000717223630585`*^-10}, {0.896`, 2.5000764243651197`*^-10}, {0.9`, 472.51457346231683`}, {0.926`, 2.500081627885553`*^-10}, {0.962`, 2.500088098325601`*^-10}, {1.`, 320.01218508343584`}, {1.1`, 207.25954907561942`}, {1.2`, 150.5082524629782`}, {1.3`, 106.75686977678878`}, {1.4`, 88.50660527800414`}, {1.52`, 66.00580667740915`}, {1.62`, 2.500249845167164`*^-15}, {1.66`, 44.50466959172432`}, {1.72`, 2.500281645243399`*^-15}, {1.82`, 2.500315350251683`*^-15}, {1.92`, 2.5003509603203922`*^-15}, {2.04`, 23.953795665685803`}, {2.4`, 13.778021747853241`}, {2.76`, 9.052625647851405`}, {3.14`, 6.152309590856069`}, {3.5`, 4.527111486280424`}, {3.88`, 3.4019499102495567`}, {4.24`, 2.5267294406115126`}, {4.62`, 1.9716021720779902`}, {4.98`, 1.4638821691671862`}, {5.36`, 1.258876887331366`}, {5.72`, 0.976215941058022`}, {6.1`, 0.7936241823262086`}, {6.46`, 0.7061217455362796`}, {6.84`, 0.5184232780164026`}, {7.2`, 0.4634144425454033`}, {7.56`, 0.3783230260588864`}, {7.94`, 2.506013334573434`*^-15}, {8.3`, 2.5065722076359985`*^-15}, {8.68`, 2.5071892511121745`*^-15}, {9.04`, 2.5077995434547677`*^-15}, {9.42`, 2.5084709290516193`*^-15}, {9.78`, 2.509132768849075`*^-15}, {10.16`, 2.5098586376587673`*^-15}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(dataqdiff5000\ = \ DeleteCases[ Table[\ If[Part[list5000, j, 2]\ >= \ 10^\(-8\), Part[list5000, j]], {j, 1, Length[list5000]}], Null]\)], "Input"], Cell[BoxData[ \({{0.036`, 1.0475000516920552`*^6}, {0.052`, 670000.0689836667`}, {0.066`, 330000.0547351886`}, {0.08`, 213500.0520286658`}, {0.096`, 103000.03614469513`}, {0.11`, 87500.04031422714`}, {0.124`, 69000.04039777268`}, {0.14`, 61250.045711679995`}, {0.154`, 47500.04289435347`}, {0.17`, 40000.044017257525`}, {0.184`, 33750.04350855281`}, {0.198`, 28500.04254421639`}, {0.214`, 23700.0413277034`}, {0.228`, 20350.040280937716`}, {0.242`, 20075.044766366005`}, {0.258`, 16550.041947209542`}, {0.26`, 14700.037838133961`}, {0.272`, 14850.0418340614`}, {0.286`, 13350.041579494713`}, {0.302`, 12600.043757291902`}, {0.316`, 11875.04515167934`}, {0.33`, 11025.045716473986`}, {0.346`, 8850.040342412547`}, {0.36`, 7975.039355209391`}, {0.374`, 7650.040744717905`}, {0.38`, 7275.04000063899`}, {0.39`, 6575.0380795506935`}, {0.404`, 6800.04226089014`}, {0.426`, 5450.037660228396`}, {0.48`, 3925.034434204875`}, {0.486`, 4075.036649508495`}, {0.544`, 2900.032678688883`}, {0.58`, 2325.0297816207153`}, {0.602`, 2187.5301863627433`}, {0.662`, 1520.0253647364577`}, {0.68`, 1362.5239897353388`}, {0.72`, 1410.0278327422413`}, {0.9`, 472.51457346231683`}, {1.`, 320.01218508343584`}, {1.1`, 207.25954907561942`}, {1.2`, 150.5082524629782`}, {1.3`, 106.75686977678878`}, {1.4`, 88.50660527800414`}, {1.52`, 66.00580667740915`}, {1.66`, 44.50466959172432`}, {2.04`, 23.953795665685803`}, {2.4`, 13.778021747853241`}, {2.76`, 9.052625647851405`}, {3.14`, 6.152309590856069`}, {3.5`, 4.527111486280424`}, {3.88`, 3.4019499102495567`}, {4.24`, 2.5267294406115126`}, {4.62`, 1.9716021720779902`}, {4.98`, 1.4638821691671862`}, {5.36`, 1.258876887331366`}, {5.72`, 0.976215941058022`}, {6.1`, 0.7936241823262086`}, {6.46`, 0.7061217455362796`}, {6.84`, 0.5184232780164026`}, {7.2`, 0.4634144425454033`}, {7.56`, 0.3783230260588864`}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(61\)], "Output"] }, Open ]], Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \ \)\(dataq5000plot\ = \ LogLogListPlot[ dataqdiff5000, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.1, 20}, {1, 100000}}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity];\)\)\)], "Input"], Cell[BoxData[ \(\(\(\t\)\(coulomb5000plot\ = \ LogLogPlot[\((icoulombcm\ //. {z\ -> 2, en\ -> \ 5000, beta\ -> \ 0.0004, thetacm\ -> \ relangle*Pi/180})\), {relangle, 0.1, 20}, \ \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ AxesLabel\ -> \ {"\", \ \ "\"}, PlotRange\ -> \ {{0.1, 10}, {1*^-5, 0.1}}, \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> \ Identity];\)\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[dataq5000plot, coulomb5000plot, DisplayFunction\ -> \ $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.434588 0.434588 0 0.123607 [ [.13082 -0.0125 -11 -12.875 ] [.13082 -0.0125 11 0 ] [.30376 -0.0125 -11 -12.875 ] [.30376 -0.0125 11 0 ] [.43459 -0.0125 -5 -12.875 ] [.43459 -0.0125 5 0 ] [.56541 -0.0125 -5 -12.875 ] [.56541 -0.0125 5 0 ] [.73835 -0.0125 -5 -12.875 ] [.73835 -0.0125 5 0 ] [.86918 -0.0125 -8 -12.875 ] [.86918 -0.0125 8 0 ] [1 -0.0125 -8 -12.875 ] [1 -0.0125 8 0 ] [-0.0125 .12361 -16 -6.4375 ] [-0.0125 .12361 0 6.4375 ] [-0.0125 .24721 -22 -6.4375 ] [-0.0125 .24721 0 6.4375 ] [-0.0125 .37082 -28 -6.4375 ] [-0.0125 .37082 0 6.4375 ] [-0.0125 .49443 -34 -6.4375 ] [-0.0125 .49443 0 6.4375 ] [-0.0125 .61803 -46 -6.4375 ] [-0.0125 .61803 0 6.4375 ] [ -0.0005 -0.0005 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .13082 0 m .13082 .00625 L s gsave .13082 -0.0125 -72 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (0.2) show 81.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 .30376 0 m .30376 .00625 L s gsave .30376 -0.0125 -72 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (0.5) show 81.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 .43459 0 m .43459 .00625 L s gsave .43459 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (1) show 69.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 .56541 0 m .56541 .00625 L s gsave .56541 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (2) show 69.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 .73835 0 m .73835 .00625 L s gsave .73835 -0.0125 -66 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (5) show 69.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 .86918 0 m .86918 .00625 L s gsave .86918 -0.0125 -69 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore 1 0 m 1 .00625 L s gsave 1 -0.0125 -69 -16.875 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (20) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore .001 w .20735 0 m .20735 .00375 L s .26165 0 m .26165 .00375 L s .33818 0 m .33818 .00375 L s .36727 0 m .36727 .00375 L s .39247 0 m .39247 .00375 L s .4147 0 m .4147 .00375 L s .64194 0 m .64194 .00375 L s .69624 0 m .69624 .00375 L s .77276 0 m .77276 .00375 L s .80186 0 m .80186 .00375 L s .82706 0 m .82706 .00375 L s .84929 0 m .84929 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12361 m .00625 .12361 L s gsave -0.0125 .12361 -77 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10) 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 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 .24721 m .00625 .24721 L s gsave -0.0125 .24721 -83 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (100) show 81.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 .37082 m .00625 .37082 L s gsave -0.0125 .37082 -89 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (1000) show 87.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 .49443 m .00625 .49443 L s gsave -0.0125 .49443 -95 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (10000) 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 0 .61803 m .00625 .61803 L s gsave -0.0125 .61803 -107 -10.4375 Mabsadd m 1 1 Mabs scale currentpoint translate /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding ISOLatin1Encoding def currentdict end newfontname exch definefont pop } def 0 20.875 translate 1 -1 scale 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 (100000.) show 105.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 .001 w 0 .03721 m .00375 .03721 L s 0 .05898 m .00375 .05898 L s 0 .07442 m .00375 .07442 L s 0 .0864 m .00375 .0864 L s 0 .09618 m .00375 .09618 L s 0 .10446 m .00375 .10446 L s 0 .11163 m .00375 .11163 L s 0 .11795 m .00375 .11795 L s 0 .16082 m .00375 .16082 L s 0 .18258 m .00375 .18258 L s 0 .19803 m .00375 .19803 L s 0 .21 m .00375 .21 L s 0 .21979 m .00375 .21979 L s 0 .22807 m .00375 .22807 L s 0 .23523 m .00375 .23523 L s 0 .24156 m .00375 .24156 L s 0 .28442 m .00375 .28442 L s 0 .30619 m .00375 .30619 L s 0 .32163 m .00375 .32163 L s 0 .33361 m .00375 .33361 L s 0 .3434 m .00375 .3434 L s 0 .35167 m .00375 .35167 L s 0 .35884 m .00375 .35884 L s 0 .36516 m .00375 .36516 L s 0 .40803 m .00375 .40803 L s 0 .4298 m .00375 .4298 L s 0 .44524 m .00375 .44524 L s 0 .45722 m .00375 .45722 L s 0 .46701 m .00375 .46701 L s 0 .47528 m .00375 .47528 L s 0 .48245 m .00375 .48245 L s 0 .48877 m .00375 .48877 L s 0 .53164 m .00375 .53164 L s 0 .5534 m .00375 .5534 L s 0 .56885 m .00375 .56885 L s 0 .58082 m .00375 .58082 L s 0 .59061 m .00375 .59061 L s 0 .59889 m .00375 .59889 L s 0 .60606 m .00375 .60606 L s 0 .61238 m .00375 .61238 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .01799 .61087 Mdot .0406 .59811 Mdot .06351 .59172 Mdot .08149 .57807 Mdot .10015 .56885 Mdot .11509 .55973 Mdot .12893 .55065 Mdot .14359 .54075 Mdot .15555 .53257 Mdot .1668 .53184 Mdot .17888 .52147 Mdot .18034 .51511 Mdot .18886 .51565 Mdot .19833 .50994 Mdot .20861 .50683 Mdot .21716 .50365 Mdot .22534 .49967 Mdot .23428 .48787 Mdot .24176 .48228 Mdot .24896 .48005 Mdot .25197 .47735 Mdot .25687 .47192 Mdot .26353 .47372 Mdot .27353 .46184 Mdot .29606 .44422 Mdot .2984 .44624 Mdot .31968 .42798 Mdot .33178 .41611 Mdot .3388 .41284 Mdot .35674 .3933 Mdot .3618 .38743 Mdot .37259 .38927 Mdot .4147 .33058 Mdot .43459 .30966 Mdot .45258 .28634 Mdot .469 .26916 Mdot .48411 .25072 Mdot .49809 .24066 Mdot .51362 .22491 Mdot .53024 .20375 Mdot .56915 .1705 Mdot .59982 .14081 Mdot .6262 .11826 Mdot .65055 .09753 Mdot .67103 .08106 Mdot .69049 .06573 Mdot .70723 .04976 Mdot .72343 .03644 Mdot .7376 .02046 Mdot .75147 .01236 Mdot 0 g .5 Mabswid 0 .21166 m .04146 .21154 L .0726 .21142 L .10314 .21125 L .12807 .21107 L .16561 .21069 L .18391 .21044 L .19936 .21018 L .22818 .20959 L .25555 .20884 L .28043 .20794 L .3011 .207 L .3373 .20481 L .36526 .20248 L .3931 .19944 L .41622 .19623 L .45691 .18866 L .48816 .18077 L .54423 .16084 L .58481 .14098 L .61532 .12281 L .66791 .08513 L .70759 .05196 L .74227 .02025 L s .74227 .02025 m .76318 0 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgoooo000`oooo00<0003oooooool0V_ooo`80/ol0 0`00002co`2co`1>oooo000`oooo00<0003oooooool0V_ooo`030;?o00000;?o0080/om>oooo000` oooo0P0009Ooool30;?o0P0000<0/om?oooo000`oooo00<0003oooooool0UOooo`<0/ol00`00002c ooooo`1Boooo000`oooo0P0009Koool20;?o00<00000/ol0/ol0Doooo`00goool0 033oool20000k_ooo`00goool001Ooool01@000?ooooooooooo`000002oooo00@0003oooooool00002 oooo00@0003oooooool00002oooo00<0003oooooool00oooo`030000oooooooo0>goool001Ooool0 1@000?ooooooooooo`000002oooo00@0003oooooool00002oooo00@0003oooooool00002oooo00<0 003oooooool00oooo`030000oooooooo04ooool30;?oVoooo`005_ooo`800004oooo0P0000Coool2 00001?ooo`800006oooo00<0003oooooool0C_ooo`L0/onHoooo000`oooo00<0003oooooool0COoo o`T0/onGoooo000`oooo00<0003oooooool0COooo`T0/onGoooo000`oooo00<0003oooooool0COoo o`T0/onGoooo000`oooo00<0003oooooool0C_ooo`<0/ol01?ooo`2co`2co`2coiSoool0033oool2 0000Boooo`<0/onPoooo000`oooo00<0003oooooool0B?ooo`H0/onOoooo000`oooo00<0003ooooo ool0Aoooo`L0/onOoooo000`oooo00<0003oooooool0A_ooo`P0/onOoooo000`oooo00<0003ooooo ool0AOooo`P0/onPoooo000`oooo0P0004Koool60;?oX_ooo`00_ooo`D0/on^oooo000`oooo0P0003Soool8 0;?o[_ooo`00"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-1.48874, -0.773956, \ 0.00999921, 0.0351561}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \ \)\(ListPlot[ dataqdiff5000, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.3, 0.6}, {0, 10000}} (*\(,\)\(\n\)\(\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \)\(DisplayFunction\ -> Identity\)*) ];\)\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations -1 3.33333 -1.6263e-018 6.18034e-005 [ [.16667 -0.0125 -12 -9 ] [.16667 -0.0125 12 0 ] [.33333 -0.0125 -9 -9 ] [.33333 -0.0125 9 0 ] [.5 -0.0125 -12 -9 ] [.5 -0.0125 12 0 ] [.66667 -0.0125 -9 -9 ] [.66667 -0.0125 9 0 ] [.83333 -0.0125 -12 -9 ] [.83333 -0.0125 12 0 ] [-0.0125 .12361 -24 -4.5 ] [-0.0125 .12361 0 4.5 ] [-0.0125 .24721 -24 -4.5 ] [-0.0125 .24721 0 4.5 ] [-0.0125 .37082 -24 -4.5 ] [-0.0125 .37082 0 4.5 ] [-0.0125 .49443 -24 -4.5 ] [-0.0125 .49443 0 4.5 ] [-0.0125 .61803 -30 -4.5 ] [-0.0125 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .16667 0 m .16667 .00625 L s [(0.35)] .16667 -0.0125 0 1 Mshowa .33333 0 m .33333 .00625 L s [(0.4)] .33333 -0.0125 0 1 Mshowa .5 0 m .5 .00625 L s [(0.45)] .5 -0.0125 0 1 Mshowa .66667 0 m .66667 .00625 L s [(0.5)] .66667 -0.0125 0 1 Mshowa .83333 0 m .83333 .00625 L s [(0.55)] .83333 -0.0125 0 1 Mshowa .125 Mabswid .03333 0 m .03333 .00375 L s .06667 0 m .06667 .00375 L s .1 0 m .1 .00375 L s .13333 0 m .13333 .00375 L s .2 0 m .2 .00375 L s .23333 0 m .23333 .00375 L s .26667 0 m .26667 .00375 L s .3 0 m .3 .00375 L s .36667 0 m .36667 .00375 L s .4 0 m .4 .00375 L s .43333 0 m .43333 .00375 L s .46667 0 m .46667 .00375 L s .53333 0 m .53333 .00375 L s .56667 0 m .56667 .00375 L s .6 0 m .6 .00375 L s .63333 0 m .63333 .00375 L s .7 0 m .7 .00375 L s .73333 0 m .73333 .00375 L s .76667 0 m .76667 .00375 L s .8 0 m .8 .00375 L s .86667 0 m .86667 .00375 L s .9 0 m .9 .00375 L s .93333 0 m .93333 .00375 L s .96667 0 m .96667 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12361 m .00625 .12361 L s [(2000)] -0.0125 .12361 1 0 Mshowa 0 .24721 m .00625 .24721 L s [(4000)] -0.0125 .24721 1 0 Mshowa 0 .37082 m .00625 .37082 L s [(6000)] -0.0125 .37082 1 0 Mshowa 0 .49443 m .00625 .49443 L s [(8000)] -0.0125 .49443 1 0 Mshowa 0 .61803 m .00625 .61803 L s [(10000)] -0.0125 .61803 1 0 Mshowa .125 Mabswid 0 .0309 m .00375 .0309 L s 0 .0618 m .00375 .0618 L s 0 .09271 m .00375 .09271 L s 0 .15451 m .00375 .15451 L s 0 .18541 m .00375 .18541 L s 0 .21631 m .00375 .21631 L s 0 .27812 m .00375 .27812 L s 0 .30902 m .00375 .30902 L s 0 .33992 m .00375 .33992 L s 0 .40172 m .00375 .40172 L s 0 .43262 m .00375 .43262 L s 0 .46353 m .00375 .46353 L s 0 .52533 m .00375 .52533 L s 0 .55623 m .00375 .55623 L s 0 .58713 m .00375 .58713 L s .25 Mabswid 0 0 m 0 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .15333 .54696 Mdot .2 .49288 Mdot .24667 .4728 Mdot .26667 .44962 Mdot .3 .40636 Mdot .34667 .42027 Mdot .42 .33683 Mdot .6 .24258 Mdot .62 .25185 Mdot .81333 .17923 Mdot .93333 .14369 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0.260791, -884.13, \ 0.00118192, 63.7461}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qdiff5000int\ = \ Interpolation[dataqdiff5000]\)], "Input"], Cell[BoxData[ TagBox[\(InterpolatingFunction[{{0.036`, 7.56`}}, "<>"]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(norm\ = \ 2*Pi*NIntegrate[Sin[relangle*3.1416/180], \ {relangle, 0, 90}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(6.283208386607923`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qt5000\ = \ 2*Pi*NIntegrate[ qdiff5000int[relangle]*Sin[relangle*3.1416/180], \ {relangle, 0.036, 7.56}]*3.1416/180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.5356406249999999`\)."\)], "Message"], Cell[BoxData[ \(7.727604727922762`\)], "Output"] }, Open ]], Cell["\<\ This value is about 15% larger than that given by Newman et al. (1985), but \ is only 7% larger than the experimental value of Savola (1973) as obtained \ using the detector angle of 0.056 deg. This difference may be the result of \ the improved differential scattering data available at low angles from the \ experiments of Nitz as used in our calculation but not in Newman's \ calculation.\ \>", "SmallText"], Cell[CellGroupData[{ Cell[BoxData[ \(qm5000\ = \ 2*Pi*NIntegrate[ qdiff5000int[relangle]*\((1 - Cos[relangle*3.1416/180])\)* Sin[relangle*3.1416/ 180], \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {relangle, 0.036, 7.56}]*3.1416/180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.8883281249999999`\)."\)], "Message"], Cell[BoxData[ \(0.0004927448116516897`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qv5000\ = \ 2*Pi*NIntegrate[ qdiff5000int[ relangle]*\((Sin[ relangle*3.1416/ 180])\)^3, \ \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \t{relangle, 0.036, 7.56}]*3.1416/ 180\)], "Input"], Cell[BoxData[ \(NIntegrate::"ncvb" \(\(:\)\(\ \)\) "NIntegrate failed to converge to prescribed accuracy after \!\(7\) \ recursive bisections in \!\(relangle\) near \!\(relangle\) = \ \!\(0.8883281249999999`\)."\)], "Message"], Cell[BoxData[ \(0.000984317960030413`\)], "Output"] }, Open ]], Cell["\<\ These qm and qv values are very small and suggests we need to look very \ carefully at the large angle scattering, e.g., how does experiment compare \ with screened Coulomb scattering.\ \>", "SmallText"], Cell[BoxData[ \(\(\t\(coulomb5000plot\ = \ LogLogPlot[ \((\((icoulombcm*\((Sin[thetacm])\)^3)\)\ //. {z\ -> 2, en\ -> \ 5000, beta\ -> \ 0.0004, thetacm\ -> \ relangle*Pi/180})\), {relangle, 0.1, 20}, \ \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ AxesLabel\ -> \ {"\", \ "\"}, PlotRange\ -> \ {{0.1, 10}, {1*^-5, 0.001}}, \n \t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> \ Identity]; \)\)\)], "Input", Evaluatable->False], Cell[BoxData[ \(\(\(\t\)\(coulomb5000plot\ = \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*linear\ plot*) \ \ \ Plot[\((\((icoulombcm*\((Sin[thetacm])\)^3)\)\ //. {z\ -> 2, en\ -> \ 5000, beta\ -> \ 0.0004, thetacm\ -> \ relangle*Pi/180})\), {relangle, 0, 90}, \ \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ AxesLabel\ -> \ {"\", \ \ "\"}, PlotRange\ -> \ {{\(-10\), 90}, {0, 0.002}}, \n\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> \ Identity];\)\)\)], "Input"], Cell[BoxData[ \(\(angleqdiff5000\ = \ Part[Transpose[dataqdiff5000], 1];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(61\)], "Output"] }, Open ]], Cell[BoxData[ \(\(listqdiff5000\ = \ \ Part[Transpose[dataqdiff5000], 2];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Length[%]\)], "Input"], Cell[BoxData[ \(61\)], "Output"] }, Open ]], Cell[BoxData[ \(\(integrandqv5000\ = \ Table[\((Part[listqdiff5000, j]*\((Sin[ Part[angleqdiff5000, j]* Pi/180])\)^3)\), \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {j, 1, Length[angleqdiff5000]}];\)\)], "Input"], Cell[BoxData[ \(\(tableintegrandqv5000\ = \ Transpose[{angleqdiff5000, integrandqv5000}];\)\)], "Input"], Cell[BoxData[ \(\(\ \ \ \ \ \ \ \ \ \ \ \(integrand5000plot\ = \ LogLogListPlot[tableintegrandqv5000, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{0.1, 10}, {1*^-5, 0.001}}, \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity]; \)\n\)\)], "Input", Evaluatable->False], Cell[BoxData[ \(\(\(\ \ \ \ \ \ \ \ \ \ \ \)\(integrand5000plot\ = \ ListPlot[ tableintegrandqv5000, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotStyle\ -> {PointSize[0.025], Hue[0.55]}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ PlotRange\ -> \ {{\(-10\), 90}, {0, 0.002}}, \n\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ DisplayFunction\ -> Identity];\)\(\n\)\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Show[integrand5000plot, coulomb5000plot, \ DisplayFunction\ -> \ $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.1 0.01 2.00577e-018 309.017 [ [.3 -0.0125 -6 -9 ] [.3 -0.0125 6 0 ] [.5 -0.0125 -6 -9 ] [.5 -0.0125 6 0 ] [.7 -0.0125 -6 -9 ] [.7 -0.0125 6 0 ] [.9 -0.0125 -6 -9 ] [.9 -0.0125 6 0 ] [.0875 .15451 -36 -4.5 ] [.0875 .15451 0 4.5 ] [.0875 .30902 -30 -4.5 ] [.0875 .30902 0 4.5 ] [.0875 .46353 -36 -4.5 ] [.0875 .46353 0 4.5 ] [.0875 .61803 -30 -4.5 ] [.0875 .61803 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .3 0 m .3 .00625 L s [(20)] .3 -0.0125 0 1 Mshowa .5 0 m .5 .00625 L s [(40)] .5 -0.0125 0 1 Mshowa .7 0 m .7 .00625 L s [(60)] .7 -0.0125 0 1 Mshowa .9 0 m .9 .00625 L s [(80)] .9 -0.0125 0 1 Mshowa .125 Mabswid .15 0 m .15 .00375 L s .2 0 m .2 .00375 L s .25 0 m .25 .00375 L s .35 0 m .35 .00375 L s .4 0 m .4 .00375 L s .45 0 m .45 .00375 L s .55 0 m .55 .00375 L s .6 0 m .6 .00375 L s .65 0 m .65 .00375 L s .75 0 m .75 .00375 L s .8 0 m .8 .00375 L s .85 0 m .85 .00375 L s .05 0 m .05 .00375 L s .95 0 m .95 .00375 L s .25 Mabswid 0 0 m 1 0 L s .1 .15451 m .10625 .15451 L s [(0.0005)] .0875 .15451 1 0 Mshowa .1 .30902 m .10625 .30902 L s [(0.001)] .0875 .30902 1 0 Mshowa .1 .46353 m .10625 .46353 L s [(0.0015)] .0875 .46353 1 0 Mshowa .1 .61803 m .10625 .61803 L s [(0.002)] .0875 .61803 1 0 Mshowa .125 Mabswid .1 .0309 m .10375 .0309 L s .1 .0618 m .10375 .0618 L s .1 .09271 m .10375 .09271 L s .1 .12361 m .10375 .12361 L s .1 .18541 m .10375 .18541 L s .1 .21631 m .10375 .21631 L s .1 .24721 m .10375 .24721 L s .1 .27812 m .10375 .27812 L s .1 .33992 m .10375 .33992 L s .1 .37082 m .10375 .37082 L s .1 .40172 m .10375 .40172 L s .1 .43262 m .10375 .43262 L s .1 .49443 m .10375 .49443 L s .1 .52533 m .10375 .52533 L s .1 .55623 m .10375 .55623 L s .1 .58713 m .10375 .58713 L s .25 Mabswid .1 0 m .1 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath 0 .7 1 r .025 w .10036 .08029 Mdot .10052 .15477 Mdot .10066 .15587 Mdot .1008 .17959 Mdot .10096 .14971 Mdot .1011 .19134 Mdot .10124 .21614 Mdot .1014 .27612 Mdot .10154 .28502 Mdot .1017 .32286 Mdot .10184 .34541 Mdot .10198 .36346 Mdot .10214 .38159 Mdot .10228 .39626 Mdot .10242 .46743 Mdot .10258 .46695 Mdot .1026 .42447 Mdot .10272 .49096 Mdot .10286 .51309 Mdot .10302 .57017 Mdot .10316 .61561 Mdot .10346 .60225 Mdot .1036 .61129 Mdot .109 .56585 Mdot .11 .52567 Mdot .111 .45313 Mdot .112 .42719 Mdot .113 .38524 Mdot .114 .39888 Mdot .1152 .38069 Mdot .1166 .33432 Mdot .1204 .33389 Mdot .124 .31265 Mdot .1276 .31233 Mdot .1314 .31246 Mdot .135 .31829 Mdot .1388 .32572 Mdot .1424 .31556 Mdot .1462 .31838 Mdot .1498 .29592 Mdot .1536 .3171 Mdot .1572 .29866 Mdot .161 .29428 Mdot .1646 .31076 Mdot .1684 .27063 Mdot .172 .28194 Mdot .1756 .26623 Mdot 0 g .5 Mabswid .1 0 m .11777 .18609 L .12 .21919 L .12243 .25015 L .12456 .27275 L .12683 .29199 L .12808 .30068 L .12946 .30875 L .13064 .31446 L .1319 .31953 L .13314 .32354 L .13428 .32641 L .13542 .32862 L .13651 .33014 L .13769 .33122 L .13835 .33158 L .13895 .33177 L .14012 .3318 L .14121 .33145 L .14184 .33111 L .1425 .33063 L .14368 .32953 L .14635 .32603 L .15072 .31809 L .15551 .30745 L .17633 .25704 L .19631 .21608 L .21453 .18681 L .2513 .14473 L .26933 .12972 L .28895 .11624 L .32749 .09583 L .3646 .08128 L .40027 .07037 L .43826 .06105 L .47482 .0537 L .5137 .04715 L .55114 .04181 L .58716 .03735 L .62549 .03321 L .66239 .02969 L .7016 .02638 L .73938 .02354 L .77573 .02108 L .8144 .01872 L .85164 .01667 L .88744 .01488 L .92556 .01315 L .96225 .01163 L 1 .01022 L Mistroke Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgoooo000/oooo1P2co`03oooo0000oooo01Woool00`000?ooooooo`3?oooo000/oooo1P2c o`03oooo0000oooo01Woool00`000?ooooooo`3?oooo000/oooo1P2co`03oooo0000oooo01Soool0 0`000?ooooooo`3@oooo000/oooo1P2co`03oooo0000oooo01Ooool00`000?ooooooo`3Aoooo000/ oooo1P2co`03oooo0000oooo01Koool00`000?ooooooo`3Boooo000/oooo1P2co`03oooo0000oooo 01Koool00`000?ooooooo`3Boooo000/oooo1P2co`03oooo0000oooo01Goool00`000?ooooooo`3C oooo000]oooo102co`;oool00`000?ooooooo`0Coooo00<0003oooooool0e?ooo`00;_ooo`030000 oooooooo00;oool00`000?ooooooo`0Coooo00<0003oooooool0e?ooo`00;Oooo`@0/ol2oooo00<0 003oooooool04_ooo`030000oooooooo0=Goool002coool60;?o00?oool0003oool04oooo`030000 oooooooo0=Goool002coool60;?o00?oool0003oool04_ooo`030000oooooooo0=Koool002coool6 0;?o00?oool0003oool04Oooo`030000oooooooo0=Ooool002coool60;?o00?oool0003oool04Ooo o`030000oooooooo0=Ooool002goool40;?o0_ooo`030000oooooooo00ooool00`000?ooooooo`3H oooo000^oooo00<0003oooooool00_ooo`030000oooooooo00ooool00`000?ooooooo`3Hoooo000^ oooo00<0003oooooool00oooo`030000oooooooo00goool00`000?ooooooo`3Ioooo000^oooo00<0 003oooooool00oooo`030000oooooooo00goool00`000?ooooooo`3Ioooo000^oooo00<0003ooooo ool00oooo`030000oooooooo00coool00`000?ooooooo`3Joooo000^oooo0P0000Coool00`000?oo ooooo`0 oooo00@0003oooooool00008oooo00@0003oooooool00002oooo00@0003oooooool00004oooo00<0 003oooooool00oooo`<000000oooo`2co`2co`020;?o00<00000/ol0/ol00`2co`0300000;?o0;?o 00<0/ooNoooo000>oooo00@0003oooooool00008oooo00@0003oooooool00002oooo00@0003ooooo ool00004oooo00<0003oooooool00_ooo`@0/ol00oooo`2co`2co`020;?o00<00000/ol0/ol00P2c o`0300000;?o0;?o00@0/ooNoooo000>oooo00@0003oooooool00008oooo00@0003oooooool00002 oooo00@0003oooooool00004oooo00@0003oooooooooool:0;?o00<00000/ol0/ol00P2co`030000 0;?o0;?o00@0/ooNoooo000?oooo0P0000[oool200001?ooo`800004oooo0P0000?oool:0;?o00D0 0000/ol0/ol0/ol000001P2comooool002coool;0;?o00@00000/ol0/ol000040;?ohOooo`00;?oo o`/0/ol00`00002co`0000040;?oh_ooo`00;Oooo`/0/ol2000000<0/oooooooool0hoooo`00;Ooo o`d0/ooVoooo000/oooo2`2conWoool002coool:0;?oj_ooo`00;?ooo`H0/oo^oooo000/oooo1P2c onkoool002goool40;?okoooo`00;?ooo`H0/oo^oooo000/oooo1P2conkoool002coool60;?ok_oo o`00;?ooo`P0/oo/oooo000/oooo2@2con_oool002coool90;?ojoooo`00;?ooo`T0/oo[oooo000/ oooo2@2con_oool002coool90;?ojoooo`00;?ooo`T0/oo[oooo000/oooo2@2con_oool002coool9 0;?ojoooo`00;Oooo`P0/oo[oooo000^oooo00<0003oool0/ol00`2concoool002koool00`000?oo ooooo`3_oooo000]oooo1P2congoool002coool80;?ok?ooo`00;?ooo`P0/oo/oooo000/oooo202c oncoool002coool80;?ok?ooo`00;Oooo`H0/oo]oooo000^oooo00<0003oooooool0koooo`00;_oo o`0300000;?o0;?o0080/oo]oooo000^oooo1P2concoool002koool60;?ok?ooo`002Oooo`800004 oooo0P0000Coool200001?ooo`800003oooo1@0000;oool200001Oooo`H0/oo/oooo0008oooo00@0 003oooooool00008oooo00@0003oooooool00002oooo00@0003oooooool00004oooo00D0003ooooo ooooool000000_ooo`040000oooooooooooo1`2concoool000Soool010000?ooooooo`0000Soool0 10000?ooooooo`0000;oool010000?ooooooo`0000Coool00`000?ooooooo`04oooo00<0003ooooo ool01`2congoool000Soool010000?ooooooo`0000Soool010000?ooooooo`0000;oool010000?oo ooooo`0000Coool00`000?ooooooo`02oooo0P0000?oool60;?ok_ooo`002?ooo`040000oooooooo 00002?ooo`040000oooooooo00000_ooo`040000oooooooo00001?ooo`030000oooooooo00;oool0 0`000?ooooooo`02oooo1P2conkoool000Woool200002_ooo`800004oooo0P0000Coool200001?oo o`<00002oooo1P2conkoool002goool40;?okoooo`00;Oooo`@0/oo_oooo000/oooo1P2conkoool0 02coool60;?ok_ooo`00;?ooo`H0/oo^oooo000/oooo1P2conkoool002goool40;?okoooo`00;?oo o`H0/oo^oooo000/oooo1P2conkoool002coool60;?ok_ooo`00;?ooo`L0/oo]oooo000]oooo1`2c oncoool002koool60;?ok?ooo`00;_ooo`H0/oo/oooo000^oooo1P2concoool002koool00`00002c o`2co`020;?okOooo`00;_ooo`030000oooooooo0>ooool002koool00`000?ooooooo`3_oooo000^ oooo00<0003oooooool0koooo`00;_ooo`030000oooooooo0>ooool002koool00`00002co`2co`02 0;?okOooo`00;Oooo`L0/oo/oooo000/oooo202concoool002coool80;?ok?ooo`00;?ooo`P0/oo/ oooo000/oooo1`2congoool002goool40;?okoooo`00;_ooo`030000oooooooo0>ooool002koool2 0000l?ooo`00;_ooo`030000oooooooo0>ooool002goool40;?okoooo`00;?ooo`H0/oo^oooo000/ oooo1P2conkoool002coool60;?ok_ooo`003oooo`800004oooo0P0000Coool200001?ooo`800003 oooo100000;oool60;?ok_ooo`003_ooo`040000oooooooo00002?ooo`040000oooooooo00000_oo o`040000oooooooo00000_ooo`030000oooooooo00?oool60;?ok_ooo`003_ooo`040000oooooooo 00002?ooo`040000oooooooo00000_ooo`040000oooooooo00000oooo`030000oooooooo00;oool6 0;?ok_ooo`003_ooo`040000oooooooo00002?ooo`040000oooooooo00000_ooo`040000oooooooo 00001?ooo`040000oooooooooooo1P2conkoool000koool010000?ooooooo`0000Soool010000?oo ooooo`0000;oool010000?ooooooo`0000;oool010000?ooooooo`0000?oool40;?okoooo`003ooo o`80000:oooo0P0000Coool200001?ooo`80003goooo003ooooo8Oooo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-17.9312, -0.000152823, \ 0.385797, 1.24846*^-05}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqt5000\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 2, en\ -> \ 5000, beta\ -> \ 0.0004, thetacm\ -> \ relangle* Pi/180})\)\ *\ \n\t\t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ Sin[relangle*3.1416/180]*Pi/180, {relangle, 7.56, 90}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.0003790660705133348`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqm5000\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 2, en\ -> \ 5000, beta\ -> \ 0.0004, thetacm\ -> \ relangle*Pi/180})\)\ \ *\ \ \((1 - Cos[relangle*3.1416/180])\)*Sin[relangle*3.1416/180]* Pi/180, \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {relangle, 7.56, 90}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.000016618129950105238`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqv5000\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 2, en\ -> \ 5000, beta\ -> \ 0.0004, thetacm\ -> \ relangle* Pi/180})\)*\n\t\t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \((Sin[relangle*3.1416/180])\)^3* Pi/180, \ {relangle, 7.56, 90}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.000029661951070467243`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(delqv5000\ = \ 2*Pi*NIntegrate[\((icoulombcm\ //. {z\ -> 2, en\ -> \ 5000, beta\ -> \ 0.0004, thetacm\ -> \ relangle* Pi/180})\)*\n\t\t\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \((Sin[relangle*3.1416/180])\)^3*Pi/180, \ {relangle, 0, 90}]*3.1416/180\)], "Input"], Cell[BoxData[ \(0.00004096095299074583`\)], "Output"] }, Open ]], Cell["Summary of cross sections from differential cross sections", "Subsection"], Cell[CellGroupData[{ Cell[BoxData[ \(qtrice\ = \ {{500, \((qt500 + delqt500)\)}, {1500, \((qt1500 + delqt1500)\)}, {5000, \((qt5000 + delqt5000)\)}}\)], "Input"], Cell[BoxData[ \({{500, 14.493110704789743`}, {1500, 11.957617274302008`}, {5000, 7.727983793993276`}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qmrice\ = \ {{500, \((qm500 + delqm500)\)}, {1500, \((qm1500 + delqm1500)\)}, {5000, \((qm5000 + delqm5000)\)}}\)], "Input"], Cell[BoxData[ \({{500, 0.0720081412850742`}, {1500, 0.010441258811749576`}, {5000, 0.000509362941601795`}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(qvrice\ = \ {{500, \((qv500 + delqv500)\)}, {1500, \((qv1500 + delqv1500)\)}, {5000, \((qv5000 + delqv5000)\)}}\)], "Input"], Cell[BoxData[ \({{500, 0.11733024641145282`}, {1500, 0.021256812077749813`}, {5000, 0.0010252789130211588`}}\)], "Output"] }, Open ]], Cell[BoxData[ \(Quit[]\)], "Input"], Cell["Not converted to He-He beyond here", "Subtitle", Background->RGBColor[1, 0, 0]], Cell["\<\ Examples of relative differential scattering cross sections for the \ parameters of our fit.\ \>", "Subsubsection", Evaluatable->False], Cell["f[x_] := a*(1 -(1 - 2*beta)*Cos[x])^-2", "Input", PageWidth->Infinity, FontFamily->"Times New Roman", FontSize->12, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["\<\ Plot[{(f[theta]/f[0] /. beta -> empiricalbeta /. enrel -> 0.01 (*eV*)), (f[theta]/f[0] /. beta -> empiricalbeta /. enrel -> 1 (*eV*)), (f[theta]/f[0] /. beta -> empiricalbeta /. enrel -> 10 (*eV*)), (f[theta]/f[0] /. beta -> empiricalbeta /. enrel -> 100 (*eV*)), (f[theta]/f[0] /. beta -> empiricalbeta /. enrel -> 1000 (*eV*))}, {theta,0,Pi}, PlotRange -> {0,1.2}, PlotStyle -> {Dashing[{ }], (*0.01 eV - solid*) Dashing[{0.003,0.003}], (*1 eV - dots*) Dashing[{0.01,0.01}], (*10 eV - short \ dashes*) Dashing[{0.05,0.01}], (*100 eV - long \ dashes*) Dashing[{0.05,0.01,0.01,0.01}]}, (*1000 \ eV - long & short dashes*) AxesLabel -> {\"theta (rd)\",\"i(theta)/i(0)\"}];\ \>", "Input", PageWidth->Infinity, FontFamily->"Times New Roman", FontSize->12, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["Export data to file.", "Section", Evaluatable->False], Cell["Parameters and and resuktant cross sections. ", "Subsubsection"], Cell[BoxData[ \(\(fittedqdqvqt\ = \ Transpose[{energies, empiricalbetalist, empiricalalist, empiricalblist, \n\t\t\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ qvfitlist, qdfitlist, qtfitlist}]; \)\)], "Input"], Cell[BoxData[ \(\(fittedresult\ = \ Map[ToString[SetPrecision[#, 4.2], FormatType -> FortranForm]&, fittedqdqvqt, {2}]; \)\)], "Input"], Cell["\<\ fitteddataout1 = TableForm[fittedresult,TableHeadings -> {None, {\"%rel. en.\",\"beta\",\"a\",\"b\",\"Qv\", \ \"Qd\",\"Qt\"}}, TableSpacing -> {0,1}]\ \>", "Input", PageWidth->Infinity, FontFamily->"Arial", FontSize->11, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["\<\ Block[ {stream1,theoryfile,today}, SetOptions[$Output,PageWidth->120]; outputfile := \"c:\\\\proplot\\\\arctprop\\\\ararxsec.dat\"; (*DeleteFile[outputfile];*) stream1 = OpenWrite[outputfile,PageWidth -> Infinity]; WriteString[stream1, \"%\", StringForm[\"Diffusion, viscosity, and total cross \ sections\"], \"\\n\", ToString[fitteddataout1]] ; Close[stream1]; SetOptions[$Output,PageWidth->53]; ]\ \>", "Input", PageWidth->Infinity, FontFamily->"Arial", FontSize->11, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["!!c:\\proplot\\arctprop\\ararxsec.dat", "Input", PageWidth->Infinity], Cell["Diffferential cross sections ", "Subsubsection"], Cell[BoxData[ \(\(ithetadatafunction\ = \ Interpolation[ithetadataout]; \)\)], "Input"], Cell[BoxData[ \(\(ithetadatatable\ = \ Table[ithetadatafunction[2*j], {j, 2, 44/2}]; \)\)], "Input"], Cell[BoxData[ \(\(thetatable\ = \ Table[2*j, {j, 2, 44/2}]; \)\)], "Input"], Cell[BoxData[ \(\(ithetafittable\ = \ \(i\ \ //. \ \ {\ beta\ -> \ empiricalbeta, \ a\ -> \ empiricala, \ b\ -> \ empiricalb\ }\)\ //. \ \ \n \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ { enrel\ \ -> \ 0.03, \ theta\ -> \ thetatable*Pi/90}; \)\)], "Input"], Cell[BoxData[ \(\(fittedithetaresult\ = \ Transpose[{thetatable, ithetadatatable, ithetafittable}]; \)\)], "Input"], Cell["\<\ fittedithetadataout = TableForm[fittedithetaresult,TableHeadings -> {None, {\"%theta\",\"i expt.\",\"i fit\"}}, TableSpacing -> {0,1}]\ \>", "Input", PageWidth->Infinity, FontFamily->"Arial", FontSize->11, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["\<\ Block[ {stream1,theoryfile,today}, SetOptions[$Output,PageWidth->120]; outputfile := \"c:\\\\proplot\\\\arctprop\\\\tparson.dat\"; (*DeleteFile[outputfile];*) stream1 = OpenWrite[outputfile,PageWidth -> Infinity]; WriteString[stream1, \"%\", StringForm[\"Diffferential scattering I(theta)\"], \"\\n\", StringForm[\"Parson,Sioska, and Lee, J. Chem. Phys. 56, \ 1511 (1972)\"], \"\\n\", ToString[fittedithetadataout]]; Close[stream1]; SetOptions[$Output,PageWidth->53]; ]\ \>", "Input", PageWidth->Infinity, Evaluatable->False, FontFamily->"Arial", FontSize->11, FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["!!c:\\proplot\\arctprop\\tparson.dat", "Input", PageWidth->Infinity, Evaluatable->False], Cell["Discussion", "Subsubsection", Evaluatable->False], Cell["\<\ 1. The objective of this notebook is to use measured and/or calculated cross \ sections to infer the effectve width of the forward peaked differential cross \ section for the elastic scattering of He by He. 2. At high energies the width is small and we have shown that the momentum \ tramsfer and viscosity cross sections weight the differential cross section \ the same way, i.e., theta^2, except that the ratio of the cross sections at \ high energies approaches 2:1. Therefore we learn very little from their \ ratio. 3. We are trying to fit three observables, i.e., Qv(enrel), Qm(enrel), and \ Qt(enrel) with a formula having three parameters a(enrel), b(enrel), and \ beta(enrel). 4. At low energies one one cannot asssume that the differential cross section \ approaches isotropic and beta approaches 0.5. In fact our analysis give a \ higher degree of anisotropy, i.e., larger beta, than at high energies. 5. We have emphasized the fit to the experimental viscosity coefficient data \ at low energies and to the theoretical viscosity cross section at high \ energies because we plan to use the results for the modeling of the transport \ of He atoms through He where the viscosity angular weighting is the important \ consideration. 6. We get a reasonable fit between the measured differential cross sections \ from Rice and the formula for Couomb scattering when we use the Newman et al \ (1985) values of the effective z, e.g., z = 2 at 5000 eV.\ \>", "SmallText", Evaluatable->False], Cell["Summary", "Subsubsection", Evaluatable->False, FontFamily->"Helv", FontWeight->"Bold", FontColor->GrayLevel[0], Background->GrayLevel[1]], Cell["\<\ We have determined total, diffusion, and viscosity cross sections from \ differential scattering cross sections for He-He at 500, 1500, and 5000 eV.\ \>", "SmallText", Evaluatable->False], Cell["runtime = SessionTime[] - startclock", "Input", PageWidth->Infinity] }, FrontEndVersion->"4.0 for Microsoft Windows", ScreenRectangle->{{0, 1024}, {0, 695}}, WindowToolbars->{}, CellGrouping->Manual, WindowSize->{819, 630}, WindowMargins->{{2, Automatic}, {Automatic, 5}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PrivateNotebookOptions->{"ColorPalette"->{RGBColor, -1}}, ShowCellLabel->True, ShowCellTags->False, RenderingOptions->{"ObjectDithering"->True, "RasterDithering"->False}, StyleDefinitions -> "AvpStyle.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[1717, 49, 170, 5, 64, "Subtitle", Evaluatable->False], Cell[1890, 56, 423, 7, 74, "Text", Evaluatable->False], Cell[CellGroupData[{ Cell[2338, 67, 42, 0, 34, "Text"], Cell[2383, 69, 41, 0, 30, "Input"], Cell[2427, 71, 39, 0, 30, "Input"], Cell[2469, 73, 44, 0, 30, "Input"], Cell[2516, 75, 36, 0, 30, "Input"], Cell[2555, 77, 37, 0, 30, "Input"], Cell[2595, 79, 38, 0, 30, "Input"], Cell[2636, 81, 38, 0, 30, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[2711, 86, 215, 4, 50, "Input"], Cell[2929, 92, 221, 4, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3187, 101, 81, 1, 47, "Subsection", Evaluatable->False], Cell[3271, 104, 93, 1, 43, "Subsubsection", Evaluatable->False], Cell[3367, 107, 121, 3, 31, "SmallText", Evaluatable->False], Cell[3491, 112, 115, 2, 30, "Input", Evaluatable->False], Cell[3609, 116, 183, 4, 31, "SmallText", Evaluatable->False], Cell[3795, 122, 105, 3, 34, "Text"], Cell[3903, 127, 2012, 41, 508, "Input"], Cell[5918, 170, 394, 11, 79, "Input", Evaluatable->False], Cell[6315, 183, 114, 2, 30, "Input", Evaluatable->False], Cell[6432, 187, 1940, 40, 492, "Input"], Cell[8375, 229, 447, 13, 79, "Input", Evaluatable->False] }, Open ]], Cell[8837, 245, 45, 0, 47, "Subsection"], Cell[8885, 247, 73, 0, 43, "Subsubsection"], Cell[8961, 249, 102, 2, 31, "SmallText"], Cell[9066, 253, 118, 2, 30, "Input"], Cell[9187, 257, 142, 3, 31, "SmallText"], Cell[CellGroupData[{ Cell[9354, 264, 138, 3, 30, "Input"], Cell[9495, 269, 66, 1, 40, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9598, 275, 180, 4, 30, "Input"], Cell[9781, 281, 199, 3, 42, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10017, 289, 159, 4, 30, "Input"], Cell[10179, 295, 164, 2, 29, "Output"] }, Open ]], Cell[10358, 300, 161, 3, 31, "SmallText"], Cell[CellGroupData[{ Cell[10544, 307, 295, 6, 70, "Input"], Cell[10842, 315, 123, 2, 44, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11002, 322, 90, 1, 30, "Input"], Cell[11095, 325, 99, 2, 44, "Output"] }, Open ]], Cell[11209, 330, 148, 3, 31, "SmallText"], Cell[11360, 335, 260, 4, 30, "Input"], Cell[11623, 341, 120, 3, 31, "SmallText"], Cell[11746, 346, 466, 9, 90, "Input"], Cell[12215, 357, 74, 0, 43, "Subsubsection"], Cell[12292, 359, 130, 3, 30, "Input"], Cell[12425, 364, 214, 5, 50, "Input"], Cell[12642, 371, 451, 8, 90, "Input"], Cell[CellGroupData[{ Cell[13118, 383, 200, 3, 50, "Input"], Cell[13321, 388, 27931, 1041, 186, 17190, 904, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[41289, 1434, 432, 7, 90, "Input"], Cell[41724, 1443, 13427, 371, 186, 3091, 239, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[55188, 1819, 78, 1, 30, "Input"], Cell[55269, 1822, 125, 3, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[55431, 1830, 130, 3, 30, "Input"], Cell[55564, 1835, 53, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[55654, 1841, 178, 4, 50, "Input"], Cell[55835, 1847, 237, 4, 44, "Message"], Cell[56075, 1853, 53, 1, 29, "Output"] }, Open ]], Cell[56143, 1857, 96, 2, 31, "SmallText"], Cell[CellGroupData[{ Cell[56264, 1863, 178, 4, 50, "Input"], Cell[56445, 1869, 237, 4, 44, "Message"], Cell[56685, 1875, 53, 1, 29, "Output"] }, Open ]], Cell[56753, 1879, 421, 7, 69, "SmallText"], Cell[CellGroupData[{ Cell[57199, 1890, 316, 6, 70, "Input"], Cell[57518, 1898, 237, 4, 44, "Message"], Cell[57758, 1904, 55, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[57850, 1910, 327, 8, 50, "Input"], Cell[58180, 1920, 237, 4, 44, "Message"], Cell[58420, 1926, 55, 1, 29, "Output"] }, Open ]], Cell[58490, 1930, 212, 4, 50, "SmallText"], Cell[58705, 1936, 90, 1, 30, "Input"], Cell[CellGroupData[{ Cell[58820, 1941, 42, 1, 30, "Input"], Cell[58865, 1944, 36, 1, 29, "Output"] }, Open ]], Cell[58916, 1948, 91, 1, 30, "Input"], Cell[CellGroupData[{ Cell[59032, 1953, 42, 1, 30, "Input"], Cell[59077, 1956, 36, 1, 29, "Output"] }, Open ]], Cell[59128, 1960, 334, 7, 50, "Input"], Cell[59465, 1969, 115, 2, 30, "Input"], Cell[59583, 1973, 580, 12, 90, "Input", Evaluatable->False], Cell[60166, 1987, 558, 9, 110, "Input"], Cell[60727, 1998, 107, 3, 31, "SmallText"], Cell[60837, 2003, 899, 17, 110, "Input", Evaluatable->False], Cell[61739, 2022, 341, 7, 102, "Input"], Cell[CellGroupData[{ Cell[62105, 2033, 131, 2, 30, "Input"], Cell[62239, 2037, 17631, 487, 186, 4826, 324, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[79907, 2529, 435, 9, 70, "Input"], Cell[80345, 2540, 53, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[80435, 2546, 237, 5, 66, "Input"], Cell[80675, 2553, 54, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[80766, 2559, 410, 8, 70, "Input"], Cell[81179, 2569, 53, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[81269, 2575, 429, 9, 70, "Input"], Cell[81701, 2586, 54, 1, 29, "Output"] }, Open ]], Cell[81770, 2590, 103, 2, 55, "Subsubtitle"], Cell[81876, 2594, 46, 0, 47, "Subsection"], Cell[CellGroupData[{ Cell[81947, 2598, 115, 2, 30, "Input"], Cell[82065, 2602, 3131, 45, 333, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[85233, 2652, 196, 4, 50, "Input"], Cell[85432, 2658, 2409, 35, 238, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[87878, 2698, 42, 1, 30, "Input"], Cell[87923, 2701, 36, 1, 29, "Output"] }, Open ]], Cell[87974, 2705, 453, 8, 90, "Input"], Cell[88430, 2715, 598, 10, 110, "Input"], Cell[CellGroupData[{ Cell[89053, 2729, 124, 2, 30, "Input"], Cell[89180, 2733, 27534, 1019, 186, 16805, 882, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[116751, 3757, 433, 7, 90, "Input"], Cell[117187, 3766, 13458, 372, 186, 3093, 239, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[130682, 4143, 80, 1, 30, "Input"], Cell[130765, 4146, 124, 3, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[130926, 4154, 136, 3, 30, "Input"], Cell[131065, 4159, 53, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[131155, 4165, 179, 4, 50, "Input"], Cell[131337, 4171, 237, 4, 44, "Message"], Cell[131577, 4177, 53, 1, 29, "Output"] }, Open ]], Cell[131645, 4181, 421, 7, 69, "SmallText"], Cell[CellGroupData[{ Cell[132091, 4192, 317, 6, 70, "Input"], Cell[132411, 4200, 237, 4, 44, "Message"], Cell[132651, 4206, 56, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[132744, 4212, 328, 8, 50, "Input"], Cell[133075, 4222, 237, 4, 44, "Message"], Cell[133315, 4228, 54, 1, 29, "Output"] }, Open ]], Cell[133384, 4232, 635, 11, 90, "Input", Evaluatable->False], Cell[134022, 4245, 775, 12, 110, "Input"], Cell[134800, 4259, 92, 1, 30, "Input"], Cell[CellGroupData[{ Cell[134917, 4264, 42, 1, 30, "Input"], Cell[134962, 4267, 36, 1, 29, "Output"] }, Open ]], Cell[135013, 4271, 93, 1, 30, "Input"], Cell[CellGroupData[{ Cell[135131, 4276, 42, 1, 30, "Input"], Cell[135176, 4279, 36, 1, 29, "Output"] }, Open ]], Cell[135227, 4283, 338, 7, 50, "Input"], Cell[135568, 4292, 118, 2, 30, "Input"], Cell[135689, 4296, 493, 10, 110, "Input", Evaluatable->False], Cell[136185, 4308, 522, 9, 110, "Input"], Cell[CellGroupData[{ Cell[136732, 4321, 133, 2, 30, "Input"], Cell[136868, 4325, 15940, 463, 186, 4774, 321, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[152845, 4793, 405, 8, 70, "Input"], Cell[153253, 4803, 53, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[153343, 4809, 450, 8, 90, "Input"], Cell[153796, 4819, 55, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[153888, 4825, 444, 8, 90, "Input"], Cell[154335, 4835, 54, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[154426, 4841, 441, 8, 90, "Input"], Cell[154870, 4851, 55, 1, 29, "Output"] }, Open ]], Cell[154940, 4855, 46, 0, 47, "Subsection"], Cell[CellGroupData[{ Cell[155011, 4859, 115, 2, 30, "Input"], Cell[155129, 4863, 3205, 47, 352, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[158371, 4915, 196, 4, 50, "Input"], Cell[158570, 4921, 2117, 31, 219, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[160724, 4957, 42, 1, 30, "Input"], Cell[160769, 4960, 36, 1, 29, "Output"] }, Open ]], Cell[160820, 4964, 453, 8, 90, "Input"], Cell[161276, 4974, 597, 10, 110, "Input"], Cell[CellGroupData[{ Cell[161898, 4988, 124, 2, 30, "Input"], Cell[162025, 4992, 27016, 998, 186, 16469, 863, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[189078, 5995, 433, 7, 90, "Input"], Cell[189514, 6004, 13409, 371, 186, 3089, 239, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[202960, 6380, 80, 1, 30, "Input"], Cell[203043, 6383, 124, 3, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[203204, 6391, 135, 3, 30, "Input"], Cell[203342, 6396, 52, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[203431, 6402, 179, 4, 50, "Input"], Cell[203613, 6408, 237, 4, 44, "Message"], Cell[203853, 6414, 52, 1, 29, "Output"] }, Open ]], Cell[203920, 6418, 421, 7, 69, "SmallText"], Cell[CellGroupData[{ Cell[204366, 6429, 317, 6, 70, "Input"], Cell[204686, 6437, 237, 4, 44, "Message"], Cell[204926, 6443, 56, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[205019, 6449, 328, 8, 50, "Input"], Cell[205350, 6459, 237, 4, 44, "Message"], Cell[205590, 6465, 55, 1, 29, "Output"] }, Open ]], Cell[205660, 6469, 213, 4, 50, "SmallText"], Cell[205876, 6475, 632, 11, 110, "Input", Evaluatable->False], Cell[206511, 6488, 746, 12, 90, "Input"], Cell[207260, 6502, 92, 1, 30, "Input"], Cell[CellGroupData[{ Cell[207377, 6507, 42, 1, 30, "Input"], Cell[207422, 6510, 36, 1, 29, "Output"] }, Open ]], Cell[207473, 6514, 93, 1, 30, "Input"], Cell[CellGroupData[{ Cell[207591, 6519, 42, 1, 30, "Input"], Cell[207636, 6522, 36, 1, 29, "Output"] }, Open ]], Cell[207687, 6526, 338, 7, 50, "Input"], Cell[208028, 6535, 118, 2, 30, "Input"], Cell[208149, 6539, 494, 10, 110, "Input", Evaluatable->False], Cell[208646, 6551, 465, 8, 110, "Input"], Cell[CellGroupData[{ Cell[209136, 6563, 130, 2, 30, "Input"], Cell[209269, 6567, 14076, 419, 186, 3967, 290, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[223382, 6991, 478, 9, 90, "Input"], Cell[223863, 7002, 56, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[223956, 7008, 436, 7, 90, "Input"], Cell[224395, 7017, 58, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[224490, 7023, 426, 8, 90, "Input"], Cell[224919, 7033, 58, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[225014, 7039, 411, 8, 70, "Input"], Cell[225428, 7049, 57, 1, 29, "Output"] }, Open ]], Cell[225500, 7053, 83, 1, 47, "Subsection"], Cell[CellGroupData[{ Cell[225608, 7058, 160, 2, 30, "Input"], Cell[225771, 7062, 128, 2, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[225936, 7069, 160, 2, 30, "Input"], Cell[226099, 7073, 133, 2, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[226269, 7080, 160, 2, 30, "Input"], Cell[226432, 7084, 135, 2, 29, "Output"] }, Open ]], Cell[226582, 7089, 39, 1, 30, "Input"], Cell[226624, 7092, 87, 1, 80, "Subtitle"], Cell[226714, 7095, 147, 4, 43, "Subsubsection", Evaluatable->False], Cell[226864, 7101, 205, 6, 45, "Input"], Cell[227072, 7109, 1093, 23, 249, "Input"], Cell[228168, 7134, 61, 1, 53, "Section", Evaluatable->False], Cell[228232, 7137, 70, 0, 43, "Subsubsection"], Cell[228305, 7139, 308, 5, 50, "Input"], Cell[228616, 7146, 156, 3, 30, "Input"], Cell[228775, 7151, 351, 11, 76, "Input"], Cell[229129, 7164, 623, 21, 236, "Input"], Cell[229755, 7187, 77, 1, 30, "Input"], Cell[229835, 7190, 54, 0, 43, "Subsubsection"], Cell[229892, 7192, 92, 1, 30, "Input"], Cell[229987, 7195, 110, 2, 30, "Input"], Cell[230100, 7199, 80, 1, 30, "Input"], Cell[230183, 7202, 422, 8, 50, "Input"], Cell[230608, 7212, 127, 2, 30, "Input"], Cell[230738, 7216, 336, 10, 76, "Input"], Cell[231077, 7228, 757, 24, 268, "Input", Evaluatable->False], Cell[231837, 7254, 98, 2, 30, "Input", Evaluatable->False], Cell[231938, 7258, 57, 1, 43, "Subsubsection", Evaluatable->False], Cell[231998, 7261, 1511, 24, 278, "SmallText", Evaluatable->False], Cell[233512, 7287, 153, 5, 59, "Subsubsection", Evaluatable->False], Cell[233668, 7294, 199, 4, 50, "SmallText", Evaluatable->False], Cell[233870, 7300, 76, 1, 30, "Input"] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)