Question

我试图以数字方式求解微分方程但我需要为我的情节改变y0并查看常数x的结果。我可以像我预期的那样正常地解决我的等式：但是当我尝试真正的目的时，我无法得到结果，你可以看到

`\[Sigma] = 1;
n = 23.04;
Rop = y[x];
R = 0.5;
sz = R/(Rop + R);
F = -n*\[Sigma]*y[x]*(1 - 2*sz);
s = NDSolve[{y'[x] == F, y[0] == 0.8}, y, {x, 0, 0.07}]
Plot[Evaluate[y[x] /. s], {x, 0, 0.07}, PlotRange -> All,]`


`[Sigma] = 1;
 n = 23.04;
 Rop = y[x];
 R = 0.5;
 sz = R/(Rop + R);
 F = -n*\[Sigma]*y[x]*(1 - 2*sz);
 y0 = 0.8;
 \!\(\*
 ButtonBox["Array",
 BaseStyle->"Link",
 ButtonData->"paclet:ref/Array"]\)[s, 140]
 i = 1;
 For[i < 140,
  s = NDSolve[{y'[x] == F, y[0] == y0}, y, {x, 0, 0.07}]
      Plot[Evaluate[y[] /. s], x = 0.07, {y0, 0.8, 2.2}] // print
      y0 == y0 + i*0.01];`

Answer 1

各种拼写错误或误解

\[Sigma] = 1;
n = 23.04;
Rop = y[x];
R = 0.5;
sz = R/(Rop + R);
F = -n*\[Sigma]*y[x]*(1 - 2*sz);
y0 = 0.8;
For[i = 1, i < 140, i++, 
  s = NDSolve[{y'[x] == F, y[0] == y0}, y, {x, 0, 0.07}]; 
  Plot[Evaluate[y[x] /. s], {x, 0, 0.07}] // Print;
  y0 = y0 + i*0.01
];

仔细检查并将其与原始角色进行比较。在弄清楚每个更改的原因之后，您可以尝试决定是否将按钮放回去。

ND的图解用另一个参数求解微分方程

1 个答案: