除非你真的理解了局部结构的用法，否则不要乱用。
z1 = With[{eps = $MachineEpsilon},
NDSolve[{x y[x] + 2 y'[x] + x y''[x] == 0, y[eps] == 1,
y'[eps] == 0}, y, {x, eps, 10}]];
Plot[y[x] /. z1, {x, 0, 10}, Method -> Evaluated, PlotStyle -> Black]
内容全部来源于网络收集，如有侵权，请联系网站删除！