您的位置: 首页  >  作业答案  >  其他  >  mathematica 解复杂的四元微分方程组NDsolve[{y1[x]*y2'[x] + y2[x]*y1'[x] == 0,3200*y1[x]*y2[x]*y3[x] == 0.075*y3'[x] + 4/3*25*10^-6*y2[x]*y2'[x],y1[x]*y2[x]*y2[x] == -y4[x] + 4/3*25*10^-6*y2'[x],y4[x] == y1[x]*461.76*y3[x],y1[0] == 67.99,y2[0] == 50,y3[0] == 811,y4[0] == 22.99*10^6},{y1[x],y2[x],y3[x],y4[x]},{x,0,1}]

mathematica 解复杂的四元微分方程组NDsolve[{y1[x]*y2'[x] + y2[x]*y1'[x] == 0,3200*y1[x]*y2[x]*y3[x] == 0.075*y3'[x] + 4/3*25*10^-6*y2[x]*y2'[x],y1[x]*y2[x]*y2[x] == -y4[x] + 4/3*25*10^-6*y2'[x],y4[x] == y1[x]*461.76*y3[x],y1[0] == 67.99,y2[0] == 50,y3[0] == 811,y4[0] == 22.99*10^6},{y1[x],y2[x],y3[x],y4[x]},{x,0,1}]

分类: 作业答案 • 2022-07-10 15:31:33
问题描述：

mathematica 解复杂的四元微分方程组
NDsolve[{y1[x]*y2'[x] + y2[x]*y1'[x] == 0,
3200*y1[x]*y2[x]*y3[x] == 0.075*y3'[x] + 4/3*25*10^-6*y2[x]*y2'[x],
y1[x]*y2[x]*y2[x] == -y4[x] + 4/3*25*10^-6*y2'[x],
y4[x] == y1[x]*461.76*y3[x],y1[0] == 67.99,
y2[0] == 50,y3[0] == 811,y4[0] == 22.99*10^6},{y1[x],y2[x],
y3[x],y4[x]},{x,0,1}]

答

函数名写错了,应该是 NDSolve

相关推荐