Question

如何用matlab求解二阶常微分方程组

y 的各返回值各表示什么

Answer 1

mathematica的解法要不要……

sol = NDSolve[{D[x[t], t, t] == x[t]/(x[t] + y[t])^2, 

   D[y[t], t, t] == y[t]/(x[t] + y[t])^2, x[0] == 1, x'[0] == 1, 

   y[0] == 2, y'[0] == 1}, {x[t], y[t]}, {t, 0, 100}];

Plot[{x[t], y[t]} /. Flatten@sol, {t, 0, 100}]

解的图形：