Wolfram mathematica Mathematica中光学布洛赫方程的代码中有什么错误/缺失？

为了工作，我必须解2能级系统的光学布洛赫方程，我似乎真的被Mathematica中的代码卡住了：

O=1;
g=1;
d=0;
sol3=NDSolve
[
{
x'[t]==g y[t] + I/2 (O* b[t] - O a[t]),
y'[t]==-g y[t]+ I/2 (O a[t]-O* b[t]),
a'[t]==-(g/2+I d) a[t] + I/2 =O* (y[t]-x[t]),
b'[t]==-(g/2-I d) b[t] + I/2 O* (x[t]-y[t]),
x[0]==1,
y[0]==0,
b[0]==0,
a[0]==0

},
{x,y},{t,0,100}
]

我得到的错误是：Syntax:：tsntxi：“整个DE系统”不完整；需要更多的投入

如果您能指出我的错误，我将不胜感激

谢谢大家：）

a[t]

和

b[t]

非常复杂，因此它们不会出现在绘图中。可以分别绘制实部和虚部

---

a[t]

和

b[t]

非常复杂，因此它们不会出现在绘图中。您可以分别绘制实部和虚部。

uhm 1错误可能与使用大写字母有关。Mathematica不知道I是什么，如果你是指虚单位，那么正确的实现方法是：ii：第三个方程在rhsGuys上有一个“=”，非常感谢。很抱歉，这是我的错误，因为我在这里编写代码而不是直接复制它（由于某些原因，它的效果不好）。为了澄清，O是一个ω，代码中不存在第三个等式中的“=”。我已经把I改为：ii:但它仍然会产生同样的错误。谢谢你们。我设法让它使用“导数[1]”而不是“x”。而@Bill an“Omega”是O的希腊字母。uhm 1错误可能与大写字母的使用有关。Mathematica不知道I是什么，如果你是指虚单位，那么正确的实现方法是：ii：第三个方程在rhsGuys上有一个“=”，非常感谢。很抱歉，这是我的错误，因为我在这里编写代码而不是直接复制它（由于某些原因，它的效果不好）。为了澄清，O是一个ω，代码中不存在第三个等式中的“=”。我已经把I改为：ii:但它仍然会产生同样的错误。谢谢你们。我设法让它使用“导数[1]”而不是“x”。而@billan“Omega”是希腊字母O。

w = 1;
g = 1;
d = 0;
swl3 = NDSolve[{
        x'[t] == g y[t] + I/2 (w*b[t] - w a[t]),
        y'[t] == -g y[t] + I/2 (w a[t] - w*b[t]), 
        a'[t] == -(g/2 + I d) a[t] + I/2 w*(y[t] - x[t]), 
        b'[t] == -(g/2 - I d) b[t] + I/2 w*(x[t] - y[t]),
        x[0] == 1, y[0] == 0, b[0] == 0, a[0] == 0},
     {x, y, a, b},
     {t, 0, 100}]

Plot[Evaluate[{x[t], y[t], a[t], b[t]} /. swl3], {t, 0, 100}]