(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 40790, 1165]*) (*NotebookOutlinePosition[ 41436, 1187]*) (* CellTagsIndexPosition[ 41392, 1183]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ \(\(\( (*\ Preliminary\ Definitions*) \)\(\[IndentingNewLine]\)\(<< Calculus`VectorAnalysis`\[IndentingNewLine] SetCoordinates[Cylindrical[R, theta, Z]]\[IndentingNewLine] ClearAll[m, kz, lambda]\[IndentingNewLine] B = B0 {0, BesselJ[1, lambda\ R], BesselJ[0, lambda\ R]}\[IndentingNewLine] Curl[B] - lambda\ B // FullSimplify\[IndentingNewLine] Div[B]\[IndentingNewLine] k = {0, m/R, \ kz}\[IndentingNewLine] k . B\[IndentingNewLine] F[R_] = R\ \(\((k . B)\)^2/k . k\)\ /B0^2 // Simplify\[IndentingNewLine] G[R_] = \((\((m^2 - 1)\)/R^2\ + \ kz^2)\) F[R] + \ \((2/R)\) \(kz^2/\((k . k)\)^2\)/ B0^2\ \((kz^2\ B[\([3]\)]^2\ - \((m/R)\)^2\ B[\([2]\)]^2)\)\ // Simplify\[IndentingNewLine]\[IndentingNewLine]\[IndentingNewLine]\ \[IndentingNewLine]\[IndentingNewLine]\[IndentingNewLine]\[IndentingNewLine] \)\)\)], "Input"], Cell[BoxData[ \(Cylindrical[R, theta, Z]\)], "Output"], Cell[BoxData[ \({0, B0\ BesselJ[1, lambda\ R], B0\ BesselJ[0, lambda\ R]}\)], "Output"], Cell[BoxData[ \({0, 0, 0}\)], "Output"], Cell[BoxData[ \(0\)], "Output"], Cell[BoxData[ \({0, m\/R, kz}\)], "Output"], Cell[BoxData[ \(B0\ kz\ BesselJ[0, lambda\ R] + \(B0\ m\ BesselJ[1, lambda\ R]\)\/R\)], "Output"], Cell[BoxData[ \(\(R\ \((kz\ R\ BesselJ[0, lambda\ R] + m\ BesselJ[1, lambda\ R])\)\^2\)\ \/\(m\^2 + kz\^2\ R\^2\)\)], "Output"], Cell[BoxData[ \(\(\(1\/\((m\^2 + kz\^2\ R\^2)\)\^2\)\((R\ \((2\ kz\^4\ R\^2\ BesselJ[0, \ lambda\ R]\^2 - 2\ kz\^2\ m\^2\ BesselJ[1, lambda\ R]\^2 + \((kz\^2 + \(\(-1\) + \ m\^2\)\/R\^2)\)\ \((m\^2 + kz\^2\ R\^2)\)\ \((kz\ R\ BesselJ[0, lambda\ R] + m\ \ BesselJ[1, lambda\ R])\)\^2)\))\)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\( (*\ Determine\ MHD\ Stability\ Boundaries\ \ *) \)\(\[IndentingNewLine]\)\(\(ClearAll[kappa1];\)\[IndentingNewLine] \(kappa1[x_] = \(\(-BesselJ[1, x]\)/x\)/ BesselJ[0, x];\)\[IndentingNewLine] \(atable = Table[{kappa1[5.2 - \((i - 1)\)/100], 5.20 - \((i - 1)\)/100}, {i, 301}];\)\[IndentingNewLine] \(a1 = Interpolation[atable];\)\[IndentingNewLine] \(ClearAll[lambda];\)\[IndentingNewLine] btable = Table[{0, 0}, {i, 161}]; Do[kzr = \(-0.7\) + \((i - 1)\)/100; btable[\([i]\)] = {kzr, \(\(FindRoot[ BesselJ[1, Sqrt[lambda^2 - kzr^2\ lambda^2]] + kzr\ lambda\ \ Sqrt[ lambda^2 - kzr^2\ lambda^2]/\((lambda - kzr\ lambda)\) BesselJ[0, Sqrt[lambda^2 - kzr^2\ lambda^2]] \[Equal] 0, {\ lambda, 3 + \((kzr - 0.2)\)^2\ 8\ + If[kzr < 0, \((kzr + .2)\)^3\ 20, 0]}]\)[\([1]\)]\)[\([2]\)]}, {i, 161}];\[IndentingNewLine] \(a2 = Interpolation[btable];\)\[IndentingNewLine] kappaminus = \(\(FindRoot[ a1[x] - a2[x] \[Equal] 0, {x, \(- .3\)}]\)[\([1]\)]\)[\([2]\)]\[IndentingNewLine] kappaplus = \(\(FindRoot[ a1[x] - a2[x] \[Equal] 0, {x, .3}]\)[\([1]\)]\)[\([2]\)]\)\)\)], "Input"], Cell[BoxData[ \(\(-0.23650606600691793`\)\)], "Output"], Cell[BoxData[ \(0.27236682023963227`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\( (*Main\ program\ for\ calculating\ radial\ displacement\ ksi \((x)\ \)\ inside\ and\ outside\ singular\ surface \((s)\) . \ Set\ parameters\ m, \ kz, \ and\ lambda\ here\ *) \)\(\[IndentingNewLine]\)\(\(m = 1; lambda = 3.6; kz = 0.4 lambda; ClearAll[ksi];\)\(\[IndentingNewLine]\) \(MHDStable = If[kappaminus < kz/lambda && kappaplus > kz/lambda && lambda > a1[kz/lambda], 0, 1];\)\(\[IndentingNewLine]\) \(Print[ If[MHDStable \[Equal] 0, "\", "\"]]\)\(\[IndentingNewLine]\) \(rsep1 = \(\(FindRoot[ F[R] \[Equal] 0, {R, 0.8}]\)[\([1]\)]\)[\([2]\)]\)\(\[IndentingNewLine]\)\( (*\ Note\ that\ for\ kz < \(-0.5\)\ lambda\ a\ second\ Singular\ surface\ \ appears*) \)\(\[IndentingNewLine]\) \(rsep2a = \(\(FindRoot[ F[R] \[Equal] 0, {R, 0.3}]\)[\([1]\)]\)[\([2]\)];\)\(\[IndentingNewLine]\) \(rsep2 = If[rsep2a < 0.1, 0.001, If[rsep2a > 1, 0.001, rsep2a]];\)\(\[IndentingNewLine]\) \(rsep = If[rsep1 < 1 && rsep1 > 0.1, rsep1, 1];\)\(\[IndentingNewLine]\) \(rsep = If[rsep == 1 && rsep2 > 0.1, rsep2, rsep1];\)\(\[IndentingNewLine]\) \(rsep1 = rsep;\)\(\[IndentingNewLine]\) \(rsep2 = If[Abs[rsep2 - rsep1] < 0.1, 0, rsep2]\)\(\[IndentingNewLine]\) \(Print[ If[rsep1 < 1 && rsep > .1, StringJoin["\", ToString[If[rsep2 > 0.1, rsep2 "\< and \>", "\< \>"]], ToString[ rsep1]], "\"]]\)\(\ \[IndentingNewLine]\) \(epsilon = 0.0005;\)\(\[IndentingNewLine]\)\( (*Calculate\ inner\ ksi\ as\ ksix\ \ and\ outer\ ksi\ as\ ksix1\ *) \)\(\[IndentingNewLine]\) \(rmin = If[rsep2 > 0.1 && rsep1 == 1, 0.01, rsep2 + epsilon]\)\(\[IndentingNewLine]\) \(rmax = If[rsep2 > 0.1 && rsep1 == 1, 1, rsep1 - epsilon]\)\(\[IndentingNewLine]\) \(sol = NDSolve[{F[R] \(ksi'\)[R] \[Equal] z[R], \(z'\)[R] - G[R] ksi[R] \[Equal] 0, ksi[rmin] \[Equal] rmin, z[rmin] \[Equal] 0}, {ksi, z}, {R, rmin, rmax}]\)\(\[IndentingNewLine]\) \(ksix[ x_] = \(Evaluate[ ksi[x] /. sol]\)[\([1]\)];\)\(\[IndentingNewLine]\) \(sol1 = NDSolve[{F[R] \(ksi'\)[R] \[Equal] z[R], \(z'\)[R] - G[R] ksi[R] \[Equal] 0, ksi[1] \[Equal] 0, z[1] \[Equal] 1}, {ksi, z}, {R, If[rsep < 1, rsep + epsilon, 0.001], 1}]\)\(\[IndentingNewLine]\) \(ksix1[ x_] = \(Evaluate[ ksi[x] /. sol1]\)[\([1]\)];\)\(\[IndentingNewLine]\)\( (*Normalization\ \ of\ Inner\ and\ Outer\ solutions, \ so\ that\ delta\ B\ radial\ is\ continuopus\ across\ Singular\ Surface\ *) \)\(\[IndentingNewLine]\) \(y0 = 1/\((\((k . B)\)/B0 /. R \[Rule] \ rsep1 - epsilon)\)\ ksix[ If[rsep1 > 1, 1, rsep1 - epsilon]];\)\(\[IndentingNewLine]\) \(y1 = 1/\((\((k . B)\)/B0 /. R \[Rule] \ rsep1 + epsilon)\)\ ksix1[ If[rsep1 > 1, 0.01, rsep1 + epsilon]];\)\(\[IndentingNewLine]\)\( (*\ Optional\ delta\ W\ calculation\ for\ MHD\ Stability\ *) \)\(\ \[IndentingNewLine]\) \( (*deltaW1 = \ \ \(If[rsep1 > 1, NIntegrate[\ F[x] \((\(ksix1'\)[x])\)^2 + \ \ G[x]\ ksix1[x]^2, {x, 0.02, 1}]/ y1^2, \ \ NIntegrate[\ F[x] \((\(ksix1'\)[x])\)^2 + \ \ G[x]\ ksix1[x]^2, {x, rsep + \ 20 epsilon, 1}]/ y1^2]\[IndentingNewLine]deltaW2 = \(If[rsep1 > 1, 0, \ \ NIntegrate[\ F[x] \((\(ksix'\)[x])\)^2 + \ \ G[x]\ ksix[x]^2, {x, 0.1, rsep - 20 epsilon}]/ y0^2]\[IndentingNewLine]\[IndentingNewLine]deltaW = \((Pi/ 2)\) \((deltaW1 + deltaW2)\)\)\)*) \)\)\)\)], "Input"], Cell[BoxData[ \("MHD Stable"\)], "Print"], Cell[BoxData[ \(0.8418181581881599`\)], "Output"], Cell[BoxData[ \(0.001`\)], "Output"], Cell[BoxData[ \("Singular Surface at R= 0.841818"\)], "Print"], Cell[BoxData[ \(0.0015`\)], "Output"], Cell[BoxData[ \(0.8413181581881599`\)], "Output"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"ksi", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.0015`, 0.8413181581881599`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"z", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.0015`, 0.8413181581881599`}}, "<>"]\), False, Editable->False]}]}], "}"}], "}"}]], "Output"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"ksi", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.8423181581881598`, 1.`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"z", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.8423181581881598`, 1.`}}, "<>"]\), False, Editable->False]}]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[BoxData[""], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\( (*Plot\ bradial\ and\ do\ Delta\ Prime\ calculation . \ Do\ only\ if\ there\ is\ a\ singular\ surface\ inside\ plasma\ \ \((rsep1 < 1)\)\ and\ the\ plasma\ is\ MHD\ Stable\ *) \)\(\ \[IndentingNewLine]\)\(\(Rmin = If[rsep1 > 1, 0.5, rsep1 - 0.05];\)\[IndentingNewLine] Rmin = rsep - 0.05\[IndentingNewLine] \(y1 = 1/\((\((k . B)\)/B0 /. R \[Rule] \ rsep1 + epsilon)\)\ ksix1[ If[rsep1 > 1, 0.01, rsep1 + epsilon]];\)\[IndentingNewLine] \(gg1 = Plot[\ \ \((k . B)\)/B0\ ksix1[R]/y1, {R, If[rsep1 > 1, 0.01, rsep1 + epsilon], 1}, PlotRange \[Rule] All, DisplayFunction \[Rule] Identity];\)\[IndentingNewLine] \(y0 = 1/\((\((k . B)\)/B0 /. R \[Rule] \ rsep1 - epsilon)\)\ ksix[ If[rsep1 > 1, 1, rsep1 - epsilon]];\)\[IndentingNewLine] \(gg2 = Plot[\ \ \((k . B)\)/B0\ ksix[R]/y0, {R, Rmin, If[rsep1 > 1, 1, rsep1 - epsilon]}, PlotRange \[Rule] All, DisplayFunction \[Rule] Identity];\)\[IndentingNewLine] \(Brsep = \((\((k . B)\)/B0\ ksix[R]/y0 /. R \[Rule] rsep1 - epsilon)\);\)\[IndentingNewLine] \(DeltaPrime = \(\(\(-\((\((\((\((k . B)\)/B0\ ksix[R]/y0 /. R \[Rule] rsep1 - epsilon)\) - \((\((k . B)\)/B0\ ksix[R]/ y0 /. R \[Rule] rsep1 - 1.1 epsilon)\))\) - \ \[IndentingNewLine]\((\((\((k . B)\)/B0\ ksix1[R]/y1 /. R \[Rule] rsep1 + 1.1\ epsilon)\) - \((\((k . B)\)/ B0\ ksix1[R]/y1 /. R \[Rule] rsep1 + epsilon)\))\))\)\)/ epsilon\)/ .1\)/Brsep;\)\[IndentingNewLine] Show[gg1, gg2, Frame \[Rule] True, DisplayFunction \[Rule] $DisplayFunction, PlotLabel \[Rule] StringJoin["\<\[CapitalDelta]'=\>", ToString[DeltaPrime]]]\[IndentingNewLine] Print[ If[DeltaPrime > 0, "\", "\"]]\[IndentingNewLine] \(lambdamin = a2[kz/lambda];\)\[IndentingNewLine] Print[ If[lambda > lambdamin, "\", "\"]]\)\)\)], "Input"], Cell[BoxData[ \(0.7918181581881598`\)], "Output"], 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 -3.59856 4.57476 0.0147151 50345 [ [.06124 -0.0125 -9 -9 ] [.06124 -0.0125 9 0 ] [.28998 -0.0125 -12 -9 ] [.28998 -0.0125 12 0 ] [.51871 -0.0125 -9 -9 ] [.51871 -0.0125 9 0 ] [.74745 -0.0125 -12 -9 ] [.74745 -0.0125 12 0 ] [.97619 -0.0125 -3 -9 ] [.97619 -0.0125 3 0 ] [ 0 0 -0.125 0 ] [-0.0125 .01472 -6 -4.5 ] [-0.0125 .01472 0 4.5 ] [-0.0125 .11541 -33 -6.25 ] [-0.0125 .11541 0 6.25 ] [-0.0125 .2161 -33 -6.25 ] [-0.0125 .2161 0 6.25 ] [-0.0125 .31679 -33 -6.25 ] [-0.0125 .31679 0 6.25 ] [-0.0125 .41748 -33 -6.25 ] [-0.0125 .41748 0 6.25 ] [-0.0125 .51816 -42 -4.5 ] [-0.0125 .51816 0 4.5 ] [ 0 0 -0.125 0 ] [ 0 .61803 .125 0 ] [ 1 0 .125 0 ] [.5 .63053 -32 0 ] [.5 .63053 32 12.5625 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .06124 0 m .06124 .00625 L s [(0.8)] .06124 -0.0125 0 1 Mshowa .28998 0 m .28998 .00625 L s [(0.85)] .28998 -0.0125 0 1 Mshowa .51871 0 m .51871 .00625 L s [(0.9)] .51871 -0.0125 0 1 Mshowa .74745 0 m .74745 .00625 L s [(0.95)] .74745 -0.0125 0 1 Mshowa .97619 0 m .97619 .00625 L s [(1)] .97619 -0.0125 0 1 Mshowa .125 Mabswid .10699 0 m .10699 .00375 L s .15273 0 m .15273 .00375 L s .19848 0 m .19848 .00375 L s .24423 0 m .24423 .00375 L s .33572 0 m .33572 .00375 L s .38147 0 m .38147 .00375 L s .42722 0 m .42722 .00375 L s .47297 0 m .47297 .00375 L s .56446 0 m .56446 .00375 L s .61021 0 m .61021 .00375 L s .65596 0 m .65596 .00375 L s .70171 0 m .70171 .00375 L s .7932 0 m .7932 .00375 L s .83895 0 m .83895 .00375 L s .8847 0 m .8847 .00375 L s .93044 0 m .93044 .00375 L s .01549 0 m .01549 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .01472 m .00625 .01472 L s [(0)] -0.0125 .01472 1 0 Mshowa 0 .11541 m .00625 .11541 L s gsave -0.0125 .11541 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (2) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore 0 .2161 m .00625 .2161 L s gsave -0.0125 .2161 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (4) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore 0 .31679 m .00625 .31679 L s gsave -0.0125 .31679 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (6) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore 0 .41748 m .00625 .41748 L s gsave -0.0125 .41748 -33 -4.25 Mabsadd m 1 1 Mabs scale /Courier findfont 10 scalefont setfont (8) show /Math1Mono findfont 10 scalefont setfont (\\264) show /Courier findfont 10 scalefont setfont (10) show 0 5.25 rmoveto /Courier findfont 7.5 scalefont setfont (-6) show grestore 0 .51816 m .00625 .51816 L s [(0.00001)] -0.0125 .51816 1 0 Mshowa .125 Mabswid 0 .03989 m .00375 .03989 L s 0 .06506 m .00375 .06506 L s 0 .09023 m .00375 .09023 L s 0 .14058 m .00375 .14058 L s 0 .16575 m .00375 .16575 L s 0 .19092 m .00375 .19092 L s 0 .24127 m .00375 .24127 L s 0 .26644 m .00375 .26644 L s 0 .29161 m .00375 .29161 L s 0 .34196 m .00375 .34196 L s 0 .36713 m .00375 .36713 L s 0 .3923 m .00375 .3923 L s 0 .44265 m .00375 .44265 L s 0 .46782 m .00375 .46782 L s 0 .49299 m .00375 .49299 L s 0 .54334 m .00375 .54334 L s 0 .56851 m .00375 .56851 L s 0 .59368 m .00375 .59368 L s .25 Mabswid 0 0 m 0 .61803 L s .06124 .61178 m .06124 .61803 L s .28998 .61178 m .28998 .61803 L s .51871 .61178 m .51871 .61803 L s .74745 .61178 m .74745 .61803 L s .97619 .61178 m .97619 .61803 L s .125 Mabswid .10699 .61428 m .10699 .61803 L s .15273 .61428 m .15273 .61803 L s .19848 .61428 m .19848 .61803 L s .24423 .61428 m .24423 .61803 L s .33572 .61428 m .33572 .61803 L s .38147 .61428 m .38147 .61803 L s .42722 .61428 m .42722 .61803 L s .47297 .61428 m .47297 .61803 L s .56446 .61428 m .56446 .61803 L s .61021 .61428 m .61021 .61803 L s .65596 .61428 m .65596 .61803 L s .70171 .61428 m .70171 .61803 L s .7932 .61428 m .7932 .61803 L s .83895 .61428 m .83895 .61803 L s .8847 .61428 m .8847 .61803 L s .93044 .61428 m .93044 .61803 L s .01549 .61428 m .01549 .61803 L s .25 Mabswid 0 .61803 m 1 .61803 L s .99375 .01472 m 1 .01472 L s .99375 .11541 m 1 .11541 L s .99375 .2161 m 1 .2161 L s .99375 .31679 m 1 .31679 L s .99375 .41748 m 1 .41748 L s .99375 .51816 m 1 .51816 L s .125 Mabswid .99625 .03989 m 1 .03989 L s .99625 .06506 m 1 .06506 L s .99625 .09023 m 1 .09023 L s .99625 .14058 m 1 .14058 L s .99625 .16575 m 1 .16575 L s .99625 .19092 m 1 .19092 L s .99625 .24127 m 1 .24127 L s .99625 .26644 m 1 .26644 L s .99625 .29161 m 1 .29161 L s .99625 .34196 m 1 .34196 L s .99625 .36713 m 1 .36713 L s .99625 .3923 m 1 .3923 L s .99625 .44265 m 1 .44265 L s .99625 .46782 m 1 .46782 L s .99625 .49299 m 1 .49299 L s .99625 .54334 m 1 .54334 L s .99625 .56851 m 1 .56851 L s .99625 .59368 m 1 .59368 L s .25 Mabswid 1 0 m 1 .61803 L s gsave .5 .63053 -93 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20.5625 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def /MISOfy { /newfontname exch def /oldfontname exch def oldfontname findfont dup length dict begin {1 index /FID ne {def} {pop pop} ifelse} forall /Encoding WindowsANSIEncoding def currentdict end newfontname exch definefont pop } def 63.000 12.813 moveto %%IncludeResource: font Courier %%IncludeFont: Courier %%BeginResource: font Courier-MISO %%BeginFont: Courier-MISO /Courier /Courier-MISO MISOfy %%EndFont %%EndResource %%IncludeResource: font Courier-MISO %%IncludeFont: Courier-MISO /Courier-MISO findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 63.000 12.813 moveto %%IncludeResource: font Mathematica1Mono %%IncludeFont: Mathematica1Mono /Mathematica1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (D) show 69.000 12.813 moveto %%IncludeResource: font Courier-MISO %%IncludeFont: Courier-MISO /Courier-MISO findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (') show 75.000 12.813 moveto %%IncludeResource: font Mathematica1Mono %%IncludeFont: Mathematica1Mono /Mathematica1Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (=) show 81.000 12.813 moveto %%IncludeResource: font Courier-MISO %%IncludeFont: Courier-MISO /Courier-MISO findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (37.3528) show 123.000 12.813 moveto %%IncludeResource: font Courier-MISO %%IncludeFont: Courier-MISO /Courier-MISO findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore .5 Mabswid .25483 .19416 m .2841 .18652 L .31601 .1782 L .34598 .17039 L .37481 .1629 L .40549 .15495 L .43503 .14732 L .46642 .13923 L .49666 .13148 L .52576 .12404 L .55671 .11617 L .58652 .10863 L .61517 .10142 L .64569 .09378 L .67505 .08648 L .70628 .07876 L .73636 .07138 L .76528 .06434 L .79607 .0569 L .82571 .0498 L .8572 .04232 L .88755 .03517 L .91675 .02837 L .94781 .0212 L .97619 .01472 L s .02381 .60332 m .033 .58612 L .04301 .56741 L .05242 .5499 L .06147 .53311 L .0711 .51529 L .08038 .49819 L .09023 .48007 L .09972 .46267 L .10886 .44599 L .11858 .4283 L .12793 .41132 L .13693 .39505 L .14651 .37778 L .15573 .36121 L .16553 .34366 L .17497 .32682 L .18405 .31067 L .19372 .29355 L .20302 .27712 L .21291 .25973 L .22243 .24304 L .2316 .22703 L .24135 .21007 L .25026 .19463 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgT0oooo0P0000030?ooo`8000000@3oool0 02D0oooo00@000000?ooo`3oool00000103oool00`000000oooo0?ooo`3T0?ooo`@000001P3oool1 00000040oooo000U0?ooo`040000003oool0oooo000000@0oooo00<000000?ooo`3oool0g`3oool5 000000X0oooo0@0000010?ooo`009P3oool2000000D0oooo00<000000?ooo`3oool0f@3oool60000 00l0oooo0@0000010?ooo`00;@3oool00`000000oooo0?ooo`3E0?ooo`@000005@3oool100000040 oooo000]0?ooo`030000003oool0oooo0=80oooo0`00000I0?ooo`4000000@3oool002d0oooo0P00 003?0?ooo`@000006`3oool200000040oooo000]0?ooo`030000003oool0oooo0`3oool400000:l0oooo0P00 00010?ooo`00;@3oool00`000000oooo0?ooo`0i0?ooo`030000003oool0oooo0;80oooo0@000001 0?ooo`00;@3oool00`000000oooo0?ooo`0h0?ooo`030000003oool0oooo0;<0oooo0@0000010?oo o`00303oool3000000030?ooo`000000oooo0080oooo1P0000020?ooo`8000003`3oool00`000000 oooo0?ooo`0h0?ooo`030000003oool0oooo0;<0oooo0@0000010?ooo`003@3oool01`000000oooo 0?ooo`3oool000000?ooo`0000000`3oool01@000000oooo0?ooo`3oool000000080oooo00<00000 0?ooo`3oool0303oool00`000000oooo0?ooo`0g0?ooo`030000003oool0oooo0;@0oooo0@000001 0?ooo`002P3oool5000000<0oooo00<000000?ooo`3oool00P3oool01@000000oooo0?ooo`3oool0 00000080oooo00<000000?ooo`3oool0303oool00`000000oooo0?ooo`0f0?ooo`030000003oool0 oooo0;D0oooo0@0000010?ooo`002P3oool010000000oooo0?ooo`0000030?ooo`030000003oool0 000000<0oooo00D000000?ooo`3oool0oooo000000020?ooo`030000003oool0oooo00`0oooo0P00 000g0?ooo`030000003oool0oooo0;@0oooo0P0000010?ooo`002`3oool00`000000oooo00000002 0?ooo`050000003oool0oooo0?ooo`0000000P3oool01@000000oooo0?ooo`3oool000000080oooo 00<000000?ooo`3oool0303oool00`000000oooo0?ooo`0e0?ooo`030000003oool0oooo0;H0oooo 0@0000010?ooo`00303oool2000000P0oooo0P0000040?ooo`800000203oool3000000@0oooo00<0 00000?ooo`3oool0=@3oool00`000000oooo0?ooo`2f0?ooo`4000000@3oool002H0oooo00<00000 0?ooo`000000103oool00`000000oooo0?ooo`0d0?ooo`030000003oool0oooo0;L0oooo0@000001 0?ooo`008@3oool300000080oooo0`0000040?ooo`030000003oool0oooo03<0oooo00<000000?oo o`3oool0^03oool100000040oooo000V0?ooo`030000003oool0oooo00@0oooo00<000000?ooo`3o ool0<`3oool00`000000oooo0?ooo`2h0?ooo`4000000@3oool002L0oooo0P0000040?ooo`800000 <`3oool00`000000oooo0?ooo`2h0?ooo`8000000@3oool002d0oooo00<000000?ooo`3oool0<@3o ool00`000000oooo0?ooo`2j0?ooo`4000000@3oool002d0oooo00<000000?ooo`3oool0<@3oool0 0`000000oooo0?ooo`2j0?ooo`4000000@3oool002d0oooo00<000000?ooo`3oool0<03oool00`00 0000oooo0?ooo`2k0?ooo`4000000@3oool002d0oooo00<000000?ooo`3oool0<03oool00`000000 oooo0?ooo`2k0?ooo`4000000@3oool002d0oooo00<000000?ooo`3oool0<03oool00`000000oooo 0?ooo`2k0?ooo`4000000@3oool002d0oooo0P00000`0?ooo`030000003oool0oooo0;/0oooo0P00 00010?ooo`00;@3oool00`000000oooo0?ooo`0_0?ooo`030000003oool0oooo0;`0oooo0@000001 0?ooo`00;@3oool00`000000oooo0?ooo`0^0?ooo`030000003oool0oooo0;d0oooo0@0000010?oo o`00;@3oool00`000000oooo0?ooo`0^0?ooo`030000003oool0oooo0;d0oooo0@0000010?ooo`00 ;@3oool00`000000oooo0?ooo`0]0?ooo`030000003oool0oooo0;h0oooo0@0000010?ooo`00;@3o ool00`000000oooo0?ooo`0/0?ooo`030000003oool0oooo0;l0oooo0@0000010?ooo`00;@3oool2 000002d0oooo00<000000?ooo`3oool0_P3oool200000040oooo000]0?ooo`030000003oool0oooo 02/0oooo00<000000?ooo`3oool0`03oool100000040oooo000]0?ooo`030000003oool0oooo02X0 oooo00<000000?ooo`3oool0`@3oool100000040oooo000:0?ooo`<000000`3oool010000000oooo 0?ooo`3oool600000080oooo0P00000?0?ooo`030000003oool0oooo02X0oooo00<000000?ooo`3o ool0`@3oool100000040oooo000:0?ooo`040000003oool0oooo000000<0oooo00<000000?ooo`00 00000`3oool01@000000oooo0?ooo`3oool000000080oooo00<000000?ooo`3oool0303oool00`00 0000oooo0?ooo`0Y0?ooo`030000003oool0oooo0<80oooo0@0000010?ooo`002P3oool010000000 oooo0?ooo`0000040?ooo`030000003oool0oooo0080oooo00D000000?ooo`3oool0oooo00000002 0?ooo`030000003oool0oooo00`0oooo00<000000?ooo`3oool0:@3oool00`000000oooo0?ooo`32 0?ooo`4000000@3oool000X0oooo0`0000040?ooo`030000003oool0000000<0oooo00D000000?oo o`3oool0oooo000000020?ooo`030000003oool0oooo00`0oooo0P00000Y0?ooo`030000003oool0 oooo0<80oooo0P0000010?ooo`002`3oool00`000000oooo0?ooo`020?ooo`050000003oool0oooo 0?ooo`0000000P3oool01@000000oooo0?ooo`3oool000000080oooo00<000000?ooo`3oool0303o ool00`000000oooo0?ooo`0X0?ooo`030000003oool0oooo0<<0oooo0@0000010?ooo`002`3oool3 000000P0oooo0P0000040?ooo`800000203oool3000000@0oooo00<000000?ooo`3oool09`3oool0 0`000000oooo0?ooo`340?ooo`4000000@3oool002H0oooo00<000000?ooo`000000103oool00`00 0000oooo0?ooo`0W0?ooo`030000003oool0oooo0<@0oooo0@0000010?ooo`008@3oool300000080 oooo0`0000040?ooo`030000003oool0oooo02H0oooo00<000000?ooo`3oool0a@3oool100000040 oooo000V0?ooo`030000003oool0oooo00@0oooo00<000000?ooo`3oool09P3oool00`000000oooo 0?ooo`350?ooo`4000000@3oool002L0oooo0P0000040?ooo`8000009P3oool00`000000oooo0?oo o`350?ooo`8000000@3oool002d0oooo00<000000?ooo`3oool09@3oool00`000000oooo0?ooo`36 0?ooo`4000000@3oool002d0oooo00<000000?ooo`3oool0903oool00`000000oooo0?ooo`370?oo o`4000000@3oool002d0oooo00<000000?ooo`3oool08`3oool00`000000oooo0?ooo`380?ooo`40 00000@3oool002d0oooo00<000000?ooo`3oool08`3oool00`000000oooo0?ooo`380?ooo`400000 0@3oool002d0oooo00<000000?ooo`3oool08P3oool00`000000oooo0?ooo`390?ooo`4000000@3o ool002d0oooo0P00000R0?ooo`030000003oool0oooo00?ooo`030000003oool0oooo0=d0oooo0@0000010?oo o`0000D0oooo0000003oool0oooo000000080?ooo`040000003oool0oooo00000080oooo00@00000 0?ooo`3oool000000P3oool010000000oooo0?ooo`0000020?ooo`040000003oool0oooo000000@0 oooo00<000000?ooo`3oool00`3oool00`000000oooo0?ooo`0>0?ooo`030000003oool0oooo0=d0 oooo0@0000010?ooo`0000D0oooo0000003oool0oooo000000080?ooo`040000003oool0oooo0000 0080oooo00@000000?ooo`3oool000000P3oool010000000oooo0?ooo`0000020?ooo`040000003o ool0oooo000000@0oooo00<000000?ooo`3oool00`3oool2000000h0oooo00<000000?ooo`3oool0 g@3oool200000040oooo00001@3oool000000?ooo`3oool0000000P0oooo00@000000?ooo`3oool0 00000P3oool010000000oooo0?ooo`0000020?ooo`040000003oool0oooo00000080oooo00@00000 0?ooo`3oool00000103oool00`000000oooo0?ooo`030?ooo`030000003oool0oooo00`0oooo00<0 00000?ooo`3oool0g`3oool100000040oooo00001@3oool000000?ooo`3oool0000000P0oooo00@0 00000?ooo`3oool000000P3oool010000000oooo0?ooo`0000020?ooo`040000003oool0oooo0000 0080oooo00@000000?ooo`3oool00000103oool00`000000oooo0?ooo`030?ooo`030000003oool0 oooo00`0oooo00<000000?ooo`3oool0g`3oool100000040oooo00020?ooo`8000002P3oool20000 00@0oooo0P0000040?ooo`800000103oool2000000@0oooo0P0000050?ooo`030000003oool0oooo 00/0oooo00<000000?ooo`3oool0h03oool100000040oooo000]0?ooo`030000003oool0oooo00X0 oooo00<000000?ooo`3oool0h@3oool100000040oooo000]0?ooo`030000003oool0oooo00X0oooo 00<000000?ooo`3oool0h@3oool100000040oooo000]0?ooo`8000002P3oool00`000000oooo0?oo o`3Q0?ooo`8000000@3oool002d0oooo00<000000?ooo`3oool02@3oool00`000000oooo0?ooo`3R 0?ooo`4000000@3oool002d0oooo00<000000?ooo`3oool02@3oool00`000000oooo0?ooo`3R0?oo o`4000000@3oool002d0oooo00<000000?ooo`3oool0203oool00`000000oooo0?ooo`3S0?ooo`40 00000@3oool002d0oooo00<000000?ooo`3oool0203oool00`000000oooo0?ooo`3S0?ooo`400000 0@3oool002d0oooo00<000000?ooo`3oool01`3oool00`000000oooo0?ooo`3T0?ooo`4000000@3o ool002d0oooo0P0000080?ooo`030000003oool0oooo0><0oooo0P0000010?ooo`00;@3oool00`00 0000oooo0?ooo`060?ooo`030000003oool0oooo0>D0oooo0@0000010?ooo`00;@3oool00`000000 oooo0?ooo`060?ooo`030000003oool0oooo0>D0oooo0@0000010?ooo`00;@3oool00`000000oooo 0?ooo`050?ooo`030000003oool0oooo0>H0oooo0@0000010?ooo`00;@3oool00`000000oooo0?oo o`050?ooo`030000003oool0oooo0>H0oooo0@0000010?ooo`00;@3oool00`000000oooo0?ooo`04 0?ooo`030000003oool0oooo0>L0oooo0@0000010?ooo`00;@3oool2000000@0oooo00<000000?oo o`3oool0i`3oool200000040oooo000]0?ooo`030000003oool0oooo00<0oooo00<000000?ooo`3o ool0j03oool100000040oooo000]0?ooo`030000003oool0oooo0080oooo00<000000?ooo`3oool0 j@3oool100000040oooo000]0?ooo`030000003oool0oooo0>h0oooo0@0000010?ooo`00;@3oool0 0`000000oooo0?ooo`0;0?ooo`030000003oool0oooo03D0oooo00<000000?ooo`3oool0=03oool0 0`000000oooo0?ooo`0d0?ooo`030000003oool0oooo03D0oooo00<000000?ooo`3oool00P3oool1 00000040oooo000]0?ooo`040000003oool0oooo000000X0oooo00<000000?ooo`3oool0203oool0 0`000000oooo0?ooo`090?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool0203oool0 0`000000oooo0?ooo`080?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool0203oool0 0`000000oooo0?ooo`080?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool0203oool0 0`000000oooo0?ooo`080?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool0203oool0 0`000000oooo0?ooo`080?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool0203oool0 0`000000oooo0?ooo`080?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool0203oool0 0`000000oooo0?ooo`090?ooo`030000003oool0oooo0080oooo0@0000010?ooo`00;@3oooob0000 0040oooo003o0?ooob40oooo003o0?ooob40oooo003o0?ooob40oooo003o0?ooob40oooo003o0?oo ob40oooo003o0?ooob40oooo003o0?ooob40oooo00270?ooo`D00000203oool4000000<0oooo0P00 00050?ooo`030000003oool0oooo0080oooo0P0000040?ooo`800000103oool2000000<0oooo1000 00030?ooo`800000G`3oool008P0oooo00<000000?ooo`0000003`3oool010000000oooo0?ooo`00 00040?ooo`030000003oool0oooo00L0oooo00@000000?ooo`3oool000000P3oool010000000oooo 0?ooo`0000020?ooo`030000003oool0oooo00<0oooo00@000000?ooo`3oool00000GP3oool008P0 oooo00<000000?ooo`0000002@3oool4000000@0oooo00<000000?ooo`3oool00`3oool00`000000 oooo0?ooo`090?ooo`030000003oool0oooo00@0oooo00D000000?ooo`3oool0oooo000000040?oo o`040000003oool0oooo000005h0oooo00290?ooo`030000003oool0oooo00@0oooo00<000000?oo o`3oool02P3oool00`000000oooo0?ooo`030?ooo`030000003oool0oooo00T0oooo00@000000?oo o`3oool0oooo0P0000050?ooo`030000003oool0oooo0080oooo0P00001O0?ooo`00R@3oool00`00 0000oooo0?ooo`040?ooo`030000003oool0oooo00L0oooo00@000000?ooo`3oool000000P3oool0 10000000oooo0?ooo`0000080?ooo`040000003oool0oooo000000<0oooo00<000000?ooo`3oool0 0P3oool010000000oooo0?ooo`0000020?ooo`040000003oool0oooo000005h0oooo00290?ooo`03 0000003oool0oooo00@0oooo00<000000?ooo`3oool0203oool2000000<0oooo100000090?ooo`80 0000103oool3000000<0oooo0P0000040?ooo`800000G`3oool00?l0oooo8@3oool00?l0oooo8@3o ool00001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {0.745781, -1.28048*^-06, \ 0.000904662, 8.2205*^-08}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"], Cell[BoxData[ \("Computational Resistive MHD Unstable"\)], "Print"], Cell[BoxData[ \("Theoretical Resistive MHD Unstable"\)], "Print"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(lambdamin\)], "Input"], Cell[BoxData[ \(3.112411350941672`\)], "Output"] }, Open ]] }, FrontEndVersion->"5.2 for Microsoft Windows", ScreenRectangle->{{0, 1024}, {0, 685}}, WindowSize->{499, 527}, WindowMargins->{{Automatic, 160}, {Automatic, 51}} ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1776, 53, 933, 17, 430, "Input"], Cell[2712, 72, 58, 1, 29, "Output"], Cell[2773, 75, 91, 1, 29, "Output"], Cell[2867, 78, 43, 1, 29, "Output"], Cell[2913, 81, 35, 1, 29, "Output"], Cell[2951, 84, 47, 1, 40, "Output"], Cell[3001, 87, 112, 2, 42, "Output"], Cell[3116, 91, 131, 2, 46, "Output"], Cell[3250, 95, 340, 6, 138, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3627, 106, 1442, 30, 490, "Input"], Cell[5072, 138, 59, 1, 29, "Output"], Cell[5134, 141, 54, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5225, 147, 4022, 84, 1170, "Input"], Cell[9250, 233, 45, 1, 25, "Print"], Cell[9298, 236, 53, 1, 29, "Output"], Cell[9354, 239, 40, 1, 29, "Output"], Cell[9397, 242, 66, 1, 25, "Print"], Cell[9466, 245, 41, 1, 29, "Output"], Cell[9510, 248, 53, 1, 29, "Output"], Cell[9566, 251, 515, 15, 48, "Output"], Cell[10084, 268, 507, 15, 48, "Output"] }, Open ]], Cell[10606, 286, 26, 0, 30, "Input"], Cell[CellGroupData[{ Cell[10657, 290, 2442, 49, 690, "Input"], Cell[13102, 341, 53, 1, 29, "Output"], Cell[13158, 344, 27203, 798, 186, 8800, 566, "GraphicsData", "PostScript", \ "Graphics"], Cell[40364, 1144, 130, 3, 29, "Output"], Cell[40497, 1149, 71, 1, 25, "Print"], Cell[40571, 1152, 69, 1, 25, "Print"] }, Open ]], Cell[CellGroupData[{ Cell[40677, 1158, 42, 1, 30, "Input"], Cell[40722, 1161, 52, 1, 29, "Output"] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)