| 1 |
(* Content-type: application/vnd.wolfram.mathematica *) |
| 2 |
|
| 3 |
(*** Wolfram Notebook File ***) |
| 4 |
(* http://www.wolfram.com/nb *) |
| 5 |
|
| 6 |
(* CreatedBy='Mathematica 9.0' *) |
| 7 |
|
| 8 |
(*CacheID: 234*) |
| 9 |
(* Internal cache information: |
| 10 |
NotebookFileLineBreakTest |
| 11 |
NotebookFileLineBreakTest |
| 12 |
NotebookDataPosition[ 157, 7] |
| 13 |
NotebookDataLength[ 18049, 524] |
| 14 |
NotebookOptionsPosition[ 17236, 490] |
| 15 |
NotebookOutlinePosition[ 17591, 506] |
| 16 |
CellTagsIndexPosition[ 17548, 503] |
| 17 |
WindowFrame->Normal*) |
| 18 |
|
| 19 |
(* Beginning of Notebook Content *) |
| 20 |
Notebook[{ |
| 21 |
Cell[BoxData[{ |
| 22 |
RowBox[{ |
| 23 |
RowBox[{"rc", "=."}], ";"}], "\[IndentingNewLine]", |
| 24 |
RowBox[{ |
| 25 |
RowBox[{"\[Kappa]", "=."}], ";"}], "\[IndentingNewLine]", |
| 26 |
RowBox[{ |
| 27 |
RowBox[{ |
| 28 |
RowBox[{"f", "[", |
| 29 |
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
| 30 |
RowBox[{ |
| 31 |
RowBox[{"Erfc", "[", |
| 32 |
RowBox[{"\[Kappa]", " ", "r"}], "]"}], "/", "r"}]}], |
| 33 |
";"}], "\[IndentingNewLine]", |
| 34 |
RowBox[{ |
| 35 |
RowBox[{ |
| 36 |
RowBox[{"g", "[", |
| 37 |
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
| 38 |
RowBox[{ |
| 39 |
RowBox[{ |
| 40 |
RowBox[{"(", |
| 41 |
RowBox[{ |
| 42 |
RowBox[{"-", "2"}], "*", "\[Kappa]"}], ")"}], "/", |
| 43 |
RowBox[{"(", |
| 44 |
RowBox[{ |
| 45 |
RowBox[{"E", "^", |
| 46 |
RowBox[{"(", |
| 47 |
RowBox[{ |
| 48 |
RowBox[{"r", "^", "2"}], "*", |
| 49 |
RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", |
| 50 |
RowBox[{"Sqrt", "[", "Pi", "]"}], "*", "r"}], ")"}]}], "-", |
| 51 |
RowBox[{ |
| 52 |
RowBox[{"Erfc", "[", |
| 53 |
RowBox[{"r", "*", "\[Kappa]"}], "]"}], "/", |
| 54 |
RowBox[{"r", "^", "2"}]}]}]}], " ", ";"}], "\[IndentingNewLine]", |
| 55 |
RowBox[{ |
| 56 |
RowBox[{ |
| 57 |
RowBox[{"h", "[", |
| 58 |
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", |
| 59 |
RowBox[{ |
| 60 |
RowBox[{ |
| 61 |
RowBox[{"(", |
| 62 |
RowBox[{"4", "*", "\[Kappa]"}], ")"}], "/", |
| 63 |
RowBox[{"(", |
| 64 |
RowBox[{ |
| 65 |
RowBox[{"E", "^", |
| 66 |
RowBox[{"(", |
| 67 |
RowBox[{ |
| 68 |
RowBox[{"r", "^", "2"}], "*", |
| 69 |
RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", |
| 70 |
RowBox[{"Sqrt", "[", "Pi", "]"}], "*", |
| 71 |
RowBox[{"r", "^", "2"}]}], ")"}]}], "+", |
| 72 |
RowBox[{ |
| 73 |
RowBox[{"(", |
| 74 |
RowBox[{"4", "*", |
| 75 |
RowBox[{"\[Kappa]", "^", "3"}]}], ")"}], "/", |
| 76 |
RowBox[{"(", |
| 77 |
RowBox[{ |
| 78 |
RowBox[{"E", "^", |
| 79 |
RowBox[{"(", |
| 80 |
RowBox[{ |
| 81 |
RowBox[{"r", "^", "2"}], "*", |
| 82 |
RowBox[{"\[Kappa]", "^", "2"}]}], ")"}]}], "*", |
| 83 |
RowBox[{"Sqrt", "[", "Pi", "]"}]}], ")"}]}], "+", |
| 84 |
RowBox[{ |
| 85 |
RowBox[{"(", |
| 86 |
RowBox[{"2", "*", |
| 87 |
RowBox[{"Erfc", "[", |
| 88 |
RowBox[{"r", "*", "\[Kappa]"}], "]"}]}], ")"}], "/", |
| 89 |
RowBox[{"r", "^", "3"}]}]}]}], ";"}], "\[IndentingNewLine]", |
| 90 |
RowBox[{ |
| 91 |
RowBox[{ |
| 92 |
RowBox[{"f1", "[", |
| 93 |
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
| 94 |
RowBox[{ |
| 95 |
RowBox[{"f", "[", |
| 96 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ", |
| 97 |
RowBox[{"f", "[", |
| 98 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-", |
| 99 |
RowBox[{ |
| 100 |
RowBox[{"(", |
| 101 |
RowBox[{"r", "-", "rc"}], ")"}], |
| 102 |
RowBox[{"g", "[", |
| 103 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], "-", |
| 104 |
RowBox[{ |
| 105 |
RowBox[{ |
| 106 |
RowBox[{"(", |
| 107 |
RowBox[{"r", "-", "rc"}], ")"}], "^", "2"}], " ", |
| 108 |
RowBox[{ |
| 109 |
RowBox[{"h", "[", |
| 110 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "/", "2"}]}]}]}], " ", |
| 111 |
";"}], "\[IndentingNewLine]", |
| 112 |
RowBox[{ |
| 113 |
RowBox[{ |
| 114 |
RowBox[{"g1", "[", |
| 115 |
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", " ", |
| 116 |
RowBox[{ |
| 117 |
RowBox[{"g", "[", |
| 118 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ", |
| 119 |
RowBox[{"g", "[", |
| 120 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}], "-", |
| 121 |
RowBox[{ |
| 122 |
RowBox[{"(", |
| 123 |
RowBox[{"r", "-", "rc"}], ")"}], |
| 124 |
RowBox[{"h", "[", |
| 125 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}]}], |
| 126 |
";"}], "\[IndentingNewLine]", |
| 127 |
RowBox[{ |
| 128 |
RowBox[{ |
| 129 |
RowBox[{"h1", "[", |
| 130 |
RowBox[{"r_", ",", "\[Kappa]_"}], "]"}], ":=", |
| 131 |
RowBox[{ |
| 132 |
RowBox[{"h", "[", |
| 133 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", |
| 134 |
RowBox[{"h", "[", |
| 135 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], " ", |
| 136 |
";"}], "\[IndentingNewLine]", |
| 137 |
RowBox[{ |
| 138 |
RowBox[{ |
| 139 |
RowBox[{"v11", "[", |
| 140 |
RowBox[{"r_", ",", "method_"}], "]"}], ":=", |
| 141 |
RowBox[{"Switch", "[", |
| 142 |
RowBox[{"method", ",", "\[IndentingNewLine]", "BARE", ",", " ", |
| 143 |
RowBox[{"g", "[", |
| 144 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
| 145 |
"TSF", ",", " ", |
| 146 |
RowBox[{"g1", "[", |
| 147 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", "\[IndentingNewLine]", "SP", |
| 148 |
",", " ", |
| 149 |
RowBox[{ |
| 150 |
RowBox[{"g", "[", |
| 151 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", " ", |
| 152 |
RowBox[{"g", "[", |
| 153 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}], ",", "\[IndentingNewLine]", |
| 154 |
"GSF", ",", " ", |
| 155 |
RowBox[{ |
| 156 |
RowBox[{"g", "[", |
| 157 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], " ", "-", " ", |
| 158 |
RowBox[{"g", "[", |
| 159 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}], " ", "-", " ", |
| 160 |
RowBox[{ |
| 161 |
RowBox[{"(", |
| 162 |
RowBox[{"r", "-", "rc"}], ")"}], |
| 163 |
RowBox[{"h", "[", |
| 164 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}]}], "\[IndentingNewLine]", |
| 165 |
"]"}]}], ";"}], "\[IndentingNewLine]", |
| 166 |
RowBox[{ |
| 167 |
RowBox[{ |
| 168 |
RowBox[{"v11prime", "[", |
| 169 |
RowBox[{"r_", ",", " ", "method_"}], "]"}], ":=", " ", |
| 170 |
RowBox[{"Switch", "[", |
| 171 |
RowBox[{"method", ",", " ", "\[IndentingNewLine]", "BARE", ",", |
| 172 |
RowBox[{"h", "[", |
| 173 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
| 174 |
"TSF", ",", " ", |
| 175 |
RowBox[{"h1", "[", |
| 176 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
| 177 |
"SP", ",", " ", |
| 178 |
RowBox[{"h", "[", |
| 179 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], ",", " ", "\[IndentingNewLine]", |
| 180 |
"GSF", ",", " ", |
| 181 |
RowBox[{ |
| 182 |
RowBox[{"h", "[", |
| 183 |
RowBox[{"r", ",", "\[Kappa]"}], "]"}], "-", |
| 184 |
RowBox[{"h", "[", |
| 185 |
RowBox[{"rc", ",", "\[Kappa]"}], "]"}]}]}], "\[IndentingNewLine]", |
| 186 |
"]"}]}], ";"}], "\[IndentingNewLine]", |
| 187 |
RowBox[{ |
| 188 |
RowBox[{ |
| 189 |
RowBox[{"R", "[", "\[Alpha]_", "]"}], " ", ":=", |
| 190 |
RowBox[{"r", " ", |
| 191 |
RowBox[{"Switch", "[", |
| 192 |
RowBox[{"\[Alpha]", ",", " ", "\[IndentingNewLine]", "1", ",", |
| 193 |
RowBox[{ |
| 194 |
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
| 195 |
RowBox[{"Cos", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]", |
| 196 |
"2", ",", |
| 197 |
RowBox[{ |
| 198 |
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
| 199 |
RowBox[{"Sin", "[", "\[Phi]", "]"}]}], " ", ",", "\[IndentingNewLine]", |
| 200 |
"3", ",", |
| 201 |
RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}]}], |
| 202 |
";"}], "\[IndentingNewLine]", |
| 203 |
RowBox[{ |
| 204 |
RowBox[{ |
| 205 |
RowBox[{"T", "[", |
| 206 |
RowBox[{"\[Alpha]_", ",", "\[Beta]_", ",", " ", "method_"}], "]"}], " ", ":=", |
| 207 |
" ", |
| 208 |
RowBox[{ |
| 209 |
RowBox[{ |
| 210 |
RowBox[{"KroneckerDelta", "[", |
| 211 |
RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], "*", " ", |
| 212 |
RowBox[{"(", |
| 213 |
RowBox[{ |
| 214 |
RowBox[{"v11", "[", |
| 215 |
RowBox[{"r", ",", "method"}], "]"}], "/", "r"}], ")"}]}], "+", " ", |
| 216 |
RowBox[{ |
| 217 |
RowBox[{"R", "[", "\[Alpha]", "]"}], |
| 218 |
RowBox[{"R", "[", "\[Beta]", "]"}], |
| 219 |
RowBox[{ |
| 220 |
RowBox[{"(", |
| 221 |
RowBox[{ |
| 222 |
RowBox[{"v11prime", "[", |
| 223 |
RowBox[{"r", ",", "method"}], "]"}], "-", |
| 224 |
RowBox[{ |
| 225 |
RowBox[{"v11", "[", |
| 226 |
RowBox[{"r", ",", "method"}], "]"}], "/", "r"}]}], ")"}], "/", |
| 227 |
RowBox[{"r", "^", "2"}]}]}]}]}], ";"}], "\[IndentingNewLine]", |
| 228 |
RowBox[{ |
| 229 |
RowBox[{"TBare", " ", "=", " ", |
| 230 |
RowBox[{"MatrixForm", "[", |
| 231 |
RowBox[{"Table", "[", |
| 232 |
RowBox[{ |
| 233 |
RowBox[{"Table", "[", |
| 234 |
RowBox[{ |
| 235 |
RowBox[{"Integrate", "[", |
| 236 |
RowBox[{ |
| 237 |
RowBox[{ |
| 238 |
RowBox[{"r", "^", "2"}], |
| 239 |
RowBox[{"Integrate", "[", |
| 240 |
RowBox[{ |
| 241 |
RowBox[{"Integrate", "[", |
| 242 |
RowBox[{ |
| 243 |
RowBox[{ |
| 244 |
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
| 245 |
RowBox[{"T", "[", |
| 246 |
RowBox[{"i", ",", "j", ",", "BARE"}], "]"}]}], ",", |
| 247 |
RowBox[{"{", |
| 248 |
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
| 249 |
",", |
| 250 |
RowBox[{"{", |
| 251 |
RowBox[{"\[Phi]", ",", "0", ",", |
| 252 |
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
| 253 |
RowBox[{"{", |
| 254 |
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
| 255 |
RowBox[{"{", |
| 256 |
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
| 257 |
RowBox[{"{", |
| 258 |
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
| 259 |
";"}], "\[IndentingNewLine]", |
| 260 |
RowBox[{ |
| 261 |
RowBox[{"TSP", " ", "=", " ", |
| 262 |
RowBox[{"MatrixForm", "[", |
| 263 |
RowBox[{"Table", "[", |
| 264 |
RowBox[{ |
| 265 |
RowBox[{"Table", "[", |
| 266 |
RowBox[{ |
| 267 |
RowBox[{"Integrate", "[", |
| 268 |
RowBox[{ |
| 269 |
RowBox[{ |
| 270 |
RowBox[{"r", "^", "2"}], |
| 271 |
RowBox[{"Integrate", "[", |
| 272 |
RowBox[{ |
| 273 |
RowBox[{"Integrate", "[", |
| 274 |
RowBox[{ |
| 275 |
RowBox[{ |
| 276 |
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
| 277 |
RowBox[{"T", "[", |
| 278 |
RowBox[{"i", ",", "j", ",", "SP"}], "]"}]}], ",", |
| 279 |
RowBox[{"{", |
| 280 |
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
| 281 |
",", |
| 282 |
RowBox[{"{", |
| 283 |
RowBox[{"\[Phi]", ",", "0", ",", |
| 284 |
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
| 285 |
RowBox[{"{", |
| 286 |
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
| 287 |
RowBox[{"{", |
| 288 |
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
| 289 |
RowBox[{"{", |
| 290 |
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
| 291 |
";"}], "\[IndentingNewLine]", |
| 292 |
RowBox[{ |
| 293 |
RowBox[{"TGSF", " ", "=", " ", |
| 294 |
RowBox[{"MatrixForm", "[", |
| 295 |
RowBox[{"Table", "[", |
| 296 |
RowBox[{ |
| 297 |
RowBox[{"Table", "[", |
| 298 |
RowBox[{ |
| 299 |
RowBox[{"Integrate", "[", |
| 300 |
RowBox[{ |
| 301 |
RowBox[{ |
| 302 |
RowBox[{"r", "^", "2"}], |
| 303 |
RowBox[{"Integrate", "[", |
| 304 |
RowBox[{ |
| 305 |
RowBox[{"Integrate", "[", |
| 306 |
RowBox[{ |
| 307 |
RowBox[{ |
| 308 |
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
| 309 |
RowBox[{"T", "[", |
| 310 |
RowBox[{"i", ",", "j", ",", "GSF"}], "]"}]}], ",", |
| 311 |
RowBox[{"{", |
| 312 |
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
| 313 |
",", |
| 314 |
RowBox[{"{", |
| 315 |
RowBox[{"\[Phi]", ",", "0", ",", |
| 316 |
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
| 317 |
RowBox[{"{", |
| 318 |
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
| 319 |
RowBox[{"{", |
| 320 |
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
| 321 |
RowBox[{"{", |
| 322 |
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
| 323 |
";"}], "\[IndentingNewLine]", |
| 324 |
RowBox[{ |
| 325 |
RowBox[{"TTSF", " ", "=", " ", |
| 326 |
RowBox[{"MatrixForm", "[", |
| 327 |
RowBox[{"Table", "[", |
| 328 |
RowBox[{ |
| 329 |
RowBox[{"Table", "[", |
| 330 |
RowBox[{ |
| 331 |
RowBox[{"Integrate", "[", |
| 332 |
RowBox[{ |
| 333 |
RowBox[{ |
| 334 |
RowBox[{"r", "^", "2"}], |
| 335 |
RowBox[{"Integrate", "[", |
| 336 |
RowBox[{ |
| 337 |
RowBox[{"Integrate", "[", |
| 338 |
RowBox[{ |
| 339 |
RowBox[{ |
| 340 |
RowBox[{"Sin", "[", "\[Theta]", "]"}], |
| 341 |
RowBox[{"T", "[", |
| 342 |
RowBox[{"i", ",", "j", ",", "TSF"}], "]"}]}], ",", |
| 343 |
RowBox[{"{", |
| 344 |
RowBox[{"\[Theta]", ",", "0", ",", "Pi"}], "}"}]}], "]"}], |
| 345 |
",", |
| 346 |
RowBox[{"{", |
| 347 |
RowBox[{"\[Phi]", ",", "0", ",", |
| 348 |
RowBox[{"2", " ", "Pi"}]}], "}"}]}], "]"}]}], ",", |
| 349 |
RowBox[{"{", |
| 350 |
RowBox[{"r", ",", "0", ",", "rc"}], "}"}]}], "]"}], ",", |
| 351 |
RowBox[{"{", |
| 352 |
RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], ",", |
| 353 |
RowBox[{"{", |
| 354 |
RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], |
| 355 |
";"}], "\[IndentingNewLine]", |
| 356 |
RowBox[{ |
| 357 |
RowBox[{"QBARE", " ", "=", |
| 358 |
RowBox[{"3", |
| 359 |
RowBox[{ |
| 360 |
RowBox[{"Part", "[", |
| 361 |
RowBox[{"TBare", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
| 362 |
RowBox[{"(", |
| 363 |
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", |
| 364 |
RowBox[{ |
| 365 |
RowBox[{"QTSP", " ", "=", " ", |
| 366 |
RowBox[{"3", |
| 367 |
RowBox[{ |
| 368 |
RowBox[{"Part", "[", |
| 369 |
RowBox[{"TSP", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
| 370 |
RowBox[{"(", |
| 371 |
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", |
| 372 |
RowBox[{ |
| 373 |
RowBox[{"QGSF", " ", "=", " ", |
| 374 |
RowBox[{"3", |
| 375 |
RowBox[{ |
| 376 |
RowBox[{"Part", "[", |
| 377 |
RowBox[{"TGSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
| 378 |
RowBox[{"(", |
| 379 |
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", |
| 380 |
RowBox[{ |
| 381 |
RowBox[{"QTSF", " ", "=", " ", |
| 382 |
RowBox[{"3", |
| 383 |
RowBox[{ |
| 384 |
RowBox[{"Part", "[", |
| 385 |
RowBox[{"TTSF", ",", "1", ",", "1", ",", "1"}], "]"}], "/", |
| 386 |
RowBox[{"(", |
| 387 |
RowBox[{"4", " ", "Pi"}], ")"}]}]}]}], ";"}]}], "Input", |
| 388 |
CellChangeTimes->{{3.613731558054296*^9, 3.613732200209785*^9}, |
| 389 |
3.6137322781282454`*^9, {3.61373233169027*^9, 3.613732397354877*^9}, { |
| 390 |
3.613732488348338*^9, 3.6137327713906393`*^9}, {3.613732810835518*^9, |
| 391 |
3.613732821957967*^9}, {3.613732862763562*^9, 3.613732913141625*^9}, { |
| 392 |
3.61373295595298*^9, 3.613733106904414*^9}, {3.613733177027177*^9, |
| 393 |
3.613733199740512*^9}, {3.61373326352699*^9, 3.6137332727080936`*^9}, { |
| 394 |
3.613733425731758*^9, 3.6137335719815893`*^9}, {3.613733605778666*^9, |
| 395 |
3.613733648768111*^9}, {3.613733776819376*^9, 3.6137338150569887`*^9}, { |
| 396 |
3.6137338699918118`*^9, 3.61373393057332*^9}, {3.613733985426153*^9, |
| 397 |
3.613733985628289*^9}, {3.613734054134171*^9, 3.6137340588514023`*^9}, { |
| 398 |
3.6137371511302147`*^9, 3.613737188118826*^9}, {3.613815304031159*^9, |
| 399 |
3.613815306899302*^9}, {3.613815415274508*^9, 3.613815450160325*^9}, { |
| 400 |
3.613815535717202*^9, 3.613815604127789*^9}, {3.613815665097225*^9, |
| 401 |
3.613815673383039*^9}, {3.613816023942026*^9, 3.6138160527230864`*^9}, { |
| 402 |
3.613825536655637*^9, 3.613825555245915*^9}, {3.613825660601698*^9, |
| 403 |
3.6138256692518787`*^9}, {3.613825968874196*^9, 3.613826000721004*^9}, { |
| 404 |
3.6138263115515537`*^9, 3.613826312408416*^9}, {3.613826919678602*^9, |
| 405 |
3.613826920319913*^9}, {3.613827178556149*^9, 3.6138272197490597`*^9}, { |
| 406 |
3.6138410090125504`*^9, 3.613841090576209*^9}, {3.613841142597171*^9, |
| 407 |
3.613841191624243*^9}, {3.6138412274048777`*^9, 3.613841229484283*^9}, { |
| 408 |
3.6138413897003736`*^9, 3.613841438224633*^9}, {3.61384150925097*^9, |
| 409 |
3.6138415133903227`*^9}, {3.613841547354535*^9, 3.613841573983226*^9}, { |
| 410 |
3.613898318732862*^9, 3.613898326606717*^9}, {3.614440844373805*^9, |
| 411 |
3.614440869838758*^9}, {3.61444103194166*^9, 3.614441084299322*^9}, { |
| 412 |
3.614441116069621*^9, 3.614441164401059*^9}, {3.6144412214281187`*^9, |
| 413 |
3.614441279916031*^9}, {3.6144413160863132`*^9, 3.614441341635235*^9}, { |
| 414 |
3.614441394376206*^9, 3.614441421059911*^9}, {3.6144416231023483`*^9, |
| 415 |
3.614441625159938*^9}, {3.614442650450296*^9, 3.6144427006914186`*^9}, { |
| 416 |
3.614442735072257*^9, 3.614442776902771*^9}, {3.614442819684997*^9, |
| 417 |
3.614442867826069*^9}, {3.614442914658215*^9, 3.614442924512073*^9}}], |
| 418 |
|
| 419 |
Cell[CellGroupData[{ |
| 420 |
|
| 421 |
Cell[BoxData["QBARE"], "Input", |
| 422 |
CellChangeTimes->{{3.61389845863691*^9, 3.6138984602794447`*^9}}], |
| 423 |
|
| 424 |
Cell[BoxData[ |
| 425 |
FractionBox[ |
| 426 |
RowBox[{"3", " ", |
| 427 |
RowBox[{"(", |
| 428 |
RowBox[{ |
| 429 |
RowBox[{ |
| 430 |
RowBox[{"-", |
| 431 |
FractionBox["8", "3"]}], " ", |
| 432 |
SuperscriptBox["\[ExponentialE]", |
| 433 |
RowBox[{ |
| 434 |
RowBox[{"-", |
| 435 |
SuperscriptBox["rc", "2"]}], " ", |
| 436 |
SuperscriptBox["\[Kappa]", "2"]}]], " ", |
| 437 |
SqrtBox["\[Pi]"], " ", "rc", " ", "\[Kappa]"}], "+", |
| 438 |
RowBox[{ |
| 439 |
FractionBox["4", "3"], " ", "\[Pi]", " ", |
| 440 |
RowBox[{"Erf", "[", |
| 441 |
RowBox[{"rc", " ", "\[Kappa]"}], "]"}]}]}], ")"}]}], |
| 442 |
RowBox[{"4", " ", "\[Pi]"}]]], "Output", |
| 443 |
CellChangeTimes->{ |
| 444 |
3.61389846058971*^9, {3.6144418380848494`*^9, 3.614441838682383*^9}, { |
| 445 |
3.614443124832337*^9, 3.614443142122366*^9}}] |
| 446 |
}, Open ]], |
| 447 |
|
| 448 |
Cell[CellGroupData[{ |
| 449 |
|
| 450 |
Cell[BoxData["QTSP"], "Input", |
| 451 |
CellChangeTimes->{{3.613898461774685*^9, 3.61389846255107*^9}}], |
| 452 |
|
| 453 |
Cell[BoxData["1"], "Output", |
| 454 |
CellChangeTimes->{3.6138984628297167`*^9, 3.614441838710091*^9, |
| 455 |
3.614443128936681*^9}] |
| 456 |
}, Open ]], |
| 457 |
|
| 458 |
Cell[CellGroupData[{ |
| 459 |
|
| 460 |
Cell[BoxData[ |
| 461 |
RowBox[{"QGSF", " "}]], "Input", |
| 462 |
CellChangeTimes->{{3.613898464699469*^9, 3.613898465498917*^9}, { |
| 463 |
3.614441850430952*^9, 3.6144418663262777`*^9}, 3.61444313552145*^9}], |
| 464 |
|
| 465 |
Cell[BoxData["1"], "Output", |
| 466 |
CellChangeTimes->{ |
| 467 |
3.613898465771669*^9, {3.614441838760746*^9, 3.614441866820423*^9}, { |
| 468 |
3.614443131316518*^9, 3.614443135852601*^9}}] |
| 469 |
}, Open ]], |
| 470 |
|
| 471 |
Cell[CellGroupData[{ |
| 472 |
|
| 473 |
Cell[BoxData["QTSF"], "Input", |
| 474 |
CellChangeTimes->{{3.613898466555105*^9, 3.613898468075191*^9}}], |
| 475 |
|
| 476 |
Cell[BoxData["1"], "Output", |
| 477 |
CellChangeTimes->{3.613898468428026*^9, 3.614441838800758*^9, |
| 478 |
3.614443138734713*^9}] |
| 479 |
}, Open ]], |
| 480 |
|
| 481 |
Cell[CellGroupData[{ |
| 482 |
|
| 483 |
Cell[BoxData[ |
| 484 |
RowBox[{"FullSimplify", "[", "QTSP", "]"}]], "Input", |
| 485 |
CellChangeTimes->{{3.613898477301433*^9, 3.613898481802129*^9}}], |
| 486 |
|
| 487 |
Cell[BoxData["1"], "Output", |
| 488 |
CellChangeTimes->{3.613898482310692*^9, 3.614441838849948*^9}] |
| 489 |
}, Open ]] |
| 490 |
}, |
| 491 |
WindowSize->{884, 1337}, |
| 492 |
WindowMargins->{{500, Automatic}, {Automatic, 0}}, |
| 493 |
FrontEndVersion->"10.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (June 27, \ |
| 494 |
2014)", |
| 495 |
StyleDefinitions->"Default.nb" |
| 496 |
] |
| 497 |
(* End of Notebook Content *) |
| 498 |
|
| 499 |
(* Internal cache information *) |
| 500 |
(*CellTagsOutline |
| 501 |
CellTagsIndex->{} |
| 502 |
*) |
| 503 |
(*CellTagsIndex |
| 504 |
CellTagsIndex->{} |
| 505 |
*) |
| 506 |
(*NotebookFileOutline |
| 507 |
Notebook[{ |
| 508 |
Cell[557, 20, 14650, 396, 811, "Input"], |
| 509 |
Cell[CellGroupData[{ |
| 510 |
Cell[15232, 420, 98, 1, 28, "Input"], |
| 511 |
Cell[15333, 423, 719, 21, 58, "Output"] |
| 512 |
}, Open ]], |
| 513 |
Cell[CellGroupData[{ |
| 514 |
Cell[16089, 449, 95, 1, 28, "Input"], |
| 515 |
Cell[16187, 452, 119, 2, 28, "Output"] |
| 516 |
}, Open ]], |
| 517 |
Cell[CellGroupData[{ |
| 518 |
Cell[16343, 459, 186, 3, 28, "Input"], |
| 519 |
Cell[16532, 464, 169, 3, 28, "Output"] |
| 520 |
}, Open ]], |
| 521 |
Cell[CellGroupData[{ |
| 522 |
Cell[16738, 472, 96, 1, 28, "Input"], |
| 523 |
Cell[16837, 475, 117, 2, 28, "Output"] |
| 524 |
}, Open ]], |
| 525 |
Cell[CellGroupData[{ |
| 526 |
Cell[16991, 482, 134, 2, 28, "Input"], |
| 527 |
Cell[17128, 486, 92, 1, 28, "Output"] |
| 528 |
}, Open ]] |
| 529 |
} |
| 530 |
] |
| 531 |
*) |
| 532 |
|
| 533 |
(* End of internal cache information *) |
