level 5
Clear["@"]
Subscript[K, 1] = 1000;
Subscript[K, 2] = 100;
Subscript[\[Gamma], 1] = Subscript[\[Gamma], 2] = 1;
Subscript[\[CapitalDelta], 1] = Subscript[\[CapitalDelta], 2] = 0;
dd11 = {};
y1 = 5;
Do[
e = 10^-5;
x10 = {2 x2, 0.1, y1/2, 0.1, 0.4, 0.4, 0.2, 0.2, 0.1, 0.3, 0.2, 0.2,
0.1};
F1[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := -2 Subscript[\
\[Gamma], 1]*\[Rho]11 + 2 x1r*\[Rho]12i;
F2[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := -(Subscript[\
\[Gamma], 1] + Subscript[\[Gamma], 2])*\[Rho]12r +
x1r*(\[Rho]22 - \[Rho]11) +
Subscript[\[CapitalDelta], 1]*\[Rho]12i + x2*\[Rho]13i;
F3[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := -(Subscript[\
\[Gamma], 1] + Subscript[\[Gamma], 2])*\[Rho]12i +
x1i*(\[Rho]22 - \[Rho]11) - x2*\[Rho]13r -
Subscript[\[CapitalDelta], 1]*\[Rho]12r;
F4[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := (Subscript[\
\[CapitalDelta], 1] + Subscript[\[CapitalDelta], 2]) \[Rho]13i -
Subscript[\[Gamma], 1]*\[Rho]13r - x1r*\[Rho]23i - x1i*\[Rho]23r +
x2*\[Rho]12i;
F5[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := \
-Subscript[\[Gamma],
1]*\[Rho]13i - (Subscript[\[CapitalDelta], 1] +
Subscript[\[CapitalDelta], 2]) \[Rho]13r + x1r*\[Rho]23r -
x1i*\[Rho]23i - x2*\[Rho]12r;
F6[yr_, yi_, x1r_,
x1i_, \[Rho]11_, \[Rho]22_, \[Rho]33_, \[Rho]12r_, \[Rho]12i_, \
\[Rho]13r_, \[Rho]13i_, \[Rho]23r_, \[Rho]23i_] := \
-Subscript[\[Gamma], 2]*\[Rho]23r +
Subscript[\[CapitalDelta], 2]*\[Rho]23i - x1r*\[Rho]13i -
x1i*\[Rho]13r;
2023年10月31日 08点10分
