如何利用c++编程求解下列微分方程：
（1）dN1/dt=(k1-N1/s1)(1-N1/k2-a*N2/k1)N1(2)dN2/dt=(k2-N2/s1)(1-N2/k1-b*N1/k2)N2其中N1(0)=50;N2(0)=50;k1,k2,s1,s2均已知；
http://www.veryhuo.com/a/view/34327.html