|
|
Re: Trying to close some loose plotting
Posted:
Feb 1, 2012 3:55 AM
|
|
Here's a correction. My last posted version didn't plot on its own.
http://home.comcast.net/~cy56/Mma/RootsOfCubicPic3.png
http://home.comcast.net/~cy56/Mma/RootsOfCubic;6.nb
\[ScriptCapitalA][\[Mu]_, rc_] = Sqrt[
4 (-3 + 2 \[Mu]^2)^3 + (27 rc - 18 \[Mu] + 20 \[Mu]^3)^2]^(1/3);
\[Lambda]1[\[Mu]_, rc_] =
2 \[Mu] - (
2*2^(1/3)*(-3 + 2 \[Mu]^2))/\[ScriptCapitalA][\[Mu], rc]^(1/3) +
2^(2/3)*\[ScriptCapitalA][\[Mu], rc]/6;
\[Lambda]2[\[Mu]_, rc_] =
1/12 (4*\[Mu] + (
2*2^(1/3)*(1 + I Sqrt[3])*(-3 +
2*\[Mu]^2))/\[ScriptCapitalA][\[Mu], rc]^(1/3) +
2^(2/3)*(-1 + Sqrt[3] I)*\[ScriptCapitalA][\[Mu], rc]^(1/3));
\[Lambda]3[\[Mu]_, rc_] =
1/12 (4*\[Mu] + (
2*2^(1/3)*(1 - Sqrt[3] I)*(-3 +
2*\[Mu]^2))/\[ScriptCapitalA][\[Mu], rc]^(1/3) -
2^(2/3)*(1 + Sqrt[3] I)* \[ScriptCapitalA][\[Mu], rc]^(1/3));
thick = AbsoluteThickness[1.25];
rctab = {-2, -1, -0.5`, -0.3`, -0.15`, -0.08`, -0.04`, 0, 0.04`,
0.08`, 0.15`, 0.3`, 0.5`, 1, 2};
Np = Length[rctab];
\[Lambda]1tab[\[Mu]_] =
Table[\[Lambda]1[\[Mu], rc] /. rc -> rctab[[i]], {i, 1, Np}];
\[Lambda]2tab[\[Mu]_] =
Table[\[Lambda]2[\[Mu], rc] /. rc -> rctab[[i]], {i, 1, Np}];
\[Lambda]3tab[\[Mu]_] =
Table[\[Lambda]3[\[Mu], rc] /. rc -> rctab[[i]], {i, 1, Np}];
\[Lambda]123[\[Mu]_, rc_] =
Flatten[{\[Lambda]1tab[\[Mu], rc], \[Lambda]2tab[\[Mu],
rc], \[Lambda]3tab[\[Mu], rc]}];
Off[Plot::"plnr"];
Plot[
Evaluate @ \[Lambda]1tab[\[Mu]], {\[Mu], -1, 1},
PlotStyle -> Table[Hue[0.8 (k - 1)/(Np - 1)], {k, 1, Np}],
PlotRange -> {{-1.1, 1.1}, {-0.1, 12}},
PlotPoints -> 300,
MaxRecursion -> 6,
AspectRatio -> 1/5 12.1/2.2,
Frame -> True,
FrameStyle -> LightGray,
Axes -> False,
ImageSize -> 300
]
|
|
