求解求解，我用mma编程解微分方程，求解出来后，我手动改变步长，超过0.001就不运算了，为什么啊

发代码

Clear["Global`*"]\[CapitalLambda] = \!\(TraditionalForm\`2.4*10^\((\(-3\))\)\);M = 1.0*10^23;R1 = 0.01;\[Rho]2 = 5.0*10^1;x1 = 0 10^-6;rou[x_] := If[x <= R1, 1.2*10^19, \[Rho]2]eq = \[Phi]''[x] + 2/x \[Phi]'[ x] == ((5.07*10^6)^2)*(-\[CapitalLambda]^5/\[Phi][x] ^2 + rou[x] /M)esol = Block[{\[Epsilon] = $MachineEpsilon}, NDSolve[{eq, \[Phi]'[\[Epsilon]] == 0, \[Phi][0.1] == Sqrt[( M \[CapitalLambda]^5)/\[Rho]2 ] }, \[Phi], {x, \[Epsilon], 0.1}, MaxStepSize -> 0.0001]] // Quiet
LogLogPlot[Evaluate[\[Phi][x] /. First[esol]], {x, x1, 0.1},
PlotRange -> All]

抱歉，下午一直在帮别人算CFD，算到现在才有时间帮你看看，首先，我不明白你为什么要用，帮你改了一下，最大步长可以变到0.0001

当精度小于0.0001时，就会出现一些数字太小无法计算的问题

\[CapitalLambda] = \!\(TraditionalForm\`2.4*10^\((\(-3\))\)\); M =
1.0*10^23; R1 = 0.01; \[Rho]2 = 5.0*10^1; x1 = 0 10^-6;
rou[x_] := If[x <= R1, 1.2*10^19, \[Rho]2];
eq = \[Phi]''[x] +
2/x \[Phi]'[
x] == ((5.07*10^6)^2)*(-\[CapitalLambda]^5/\[Phi][x] ^2 +
rou[x] /M);
s = NDSolve[{eq, \[Phi]'[$MachineEpsilon] == 0, \[Phi][0.1] ==
Sqrt[( M \[CapitalLambda]^5)/\[Rho]2 ] }, \[Phi][
x], {x, $MachineEpsilon, 0.1}, MaxStepSize -> 0.0001];
Plot[Evaluate[\[Phi][x] /. s], {x, 0, 30}, PlotRange -> All]

1. 谁告诉你们MaxStepSize代表精度了。
2. MaxStepSize没事别动，或者说具体点，这选项对非线性常微分方程边值问题的求解没帮助，瞎动反而容易导致问题。
3. 不要在没搞清楚代码正误的情况下往里加Quiet。
4. 四楼的画图范围错了。