我需要更多的信息关于这个
sys1:= {diff (Eim (t), t) = 0.5 *(α* Er (t) + Eim (t)) * (g * (N (t) - N0) - 1 / taup),
diff (Er (t), t) = 0.5 * (Er (t) -α* Eim (t)) * (g * (N (t) - N0) - 1 / taup), diff (N (t), t) =
stroom1 /电荷- N (t) / tauc (Er (t) ^ 2 + Eim (t) ^ 2) * g * (N (t) - N0)};