[求助]求解常微分方程组
 

有一组常微分方程组,如下:
dn1/dt=-28.0807*n1-1.86e-17*n3*n1+3650*n3+n2*1e6/a;
dn2/dt=350*n3+2*1.86e-17*n3*n1-n2*1e6/a;
dn3/dt=28.0807*n1-1.86e-17*n3*n1-n3*1e6/250;
n1+n2+n3=1e20
n1(0)=1e20,n2(0)=n3(0)=0

我想问一下在matlab里面能通过这四个方程,解出一个确定的a值吗?

我看的文献中是说以a作为变量用四阶Runge-Kutta方法可以解出a值,但是我使用ode45时.不知道如何将n1+n2+n3=1e20这个条件加进去,望达人解答,小弟先多谢了!:lol: