Question

用mathmetica求解常微分方程组的问题！！ 在线等啊 急 急 急

r=1,d=0.5,a=0.1,b=0.02,x0=25,y0=2;
x'[t] == r*x[t] - a*x[t]*y[t]
y'[t] == -d*y[t] + b*x[t]*y[t]求解x[0]=x0 y[0]=yo时的数值解的图形以及相轨线y[x] 用mathmetica求解 求源代码 最好是调试过的 先谢了！！！在线等啊 可以追加分数

Answer 1

r = 1; d = 0.5; a = 0.1; b = 0.02; x0 = 25; y0 = 2;
sol = Flatten[NDSolve[{x'[t] == r*x[t] - a*x[t]*y[t],
y'[t] == -d*y[t] + b*x[t]*y[t], x[0] == x0, y[0] == y0}, {x,
y}, {t, -10, 10}]]
X[t_] = Evaluate[x[t] /. sol];
Y[t_] = Evaluate[y[t] /. sol];
Plot[{X[t], Y[t]}, {t, -10, 10}, PlotRange -> All]
ParametricPlot[{X[t], Y[t]}, {t, -10, 10}, PlotRange -> All]
最下边的命令是相图，倒数第二个是函数图象，两个。