dN1/dt=-Wa*N1+We*N3+A*N3                   (1)
dN3/dt=Wa*N1-We*N3-A*N3                    (2)
de/dz=(Ge*N3-Ga*N1-0.1)*e                  (3)
dea/dz=(Ge*N3-Ga*N1-0.1)*e+V*H*Ge*N3       (4)
P=e+ea                                    （5）
N1+N3-N=0                                 （6）
Wa=P*Ga/(H*Va*S)                           (7)
We=P*Ge/(H*Ve*S)                           (8)
 V=Va-Ve                 Va=3*10^8/360*10^(-9)     Ve=3*10^8/450*10^(-9)
 Ge=Ga=10^(-20)
 A=10000
S=1
Z=1
H=6.624*10^(-34)

大神们 求帮助啊！！！！先谢过啦

采用ode求解微分方程

再具体点啊？老大。