mathematica中对定义的函数作运算

Question

举个简单的例子：我要解一个微分方程组{dx1/dt=x1+H, dx2/dt=x1*x2-H} 其中H=x1^2+x1*x2，而x1，x2都是t的函数。我要先定义H，然后在用DSolve解这个方程组的时候会引用H。求问怎样实现？

Answer 1

H = x1[t]^2 + x2[t]^2;
equ = {D[x1[t], t] == x1[t] + H, D[x2[t], t] == x1[t]*x2[t] - H};
sol = DSolve[equ, {x1, x2}, t]
但这样求不出来 这可能这方程组没有解析解

你可以尝试数值解法 带入x1 ,x2的初始值，在一定范围的t内求数值解 就是
H = x1[t]^2 + x2[t]^2;
equ = {D[x1[t], t] == x1[t] + H, D[x2[t], t] == x1[t]*x2[t] - H, x1[0] == -0.1, x2[0] == 0.1};
sol = NDSolve[equ, {x1, x2}, {t, 0, 5}]
Plot[{sol[[1, 1, 2]][t], sol[[1, 2, 2]][t]}, {t, 0,5}, PlotRange -> All]
这里假设 x1[0] == -0.1, x2[0] == 0.1，并画出了t在[0,5]之间的值
你改变下初始条件 你可以看到这方程增长很快 经常趋向无穷大

望采纳！！

Answer 2

这样？：

h := x[t]^2 + x[t] y[t]
DSolve[{x'[t] == x[t] + h, y'[t] == x[t] y[t] - h}, {x[t], y[t]}, t]

不过你这个方程好像不能DSolve……