[求助] 请大家帮忙看看我这个Mathematica程序有什么问题，为什么给不出一个完整的波形？

请大家帮忙看看我这个程序有什么问题，为什么给不出一个完整的波形？ 我这个用的是夹逼原理，就是通过调节参数EE来最终达到一个完整的正弦或者余弦波形，为什么我调节EE，波形都不发生变化呢？ 是不是这种类型的方程，我的程序本身有问题。 第一步：这一步表示先定义一个方程； Equation[\[Sigma]0_, \[Delta]0_, R_, a_, M_, fv_, k_, EE_] :=   G[r] (-((-1 + k) k)/        r^2 - (M + EE - \[Delta]0/(1 + Exp[(r - R)/a])) (M -          EE + \[Sigma]0/(1 + Exp[(r - R)/a])) -       2 k (fv/(         4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] +           D[\[Delta]0/(1 + Exp[(r - R)/a]), r]))/r -       D[fv/(4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] +           D[\[Delta]0/(1 + Exp[(r - R)/a]), r]),        r] - (fv/(          4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] +            D[\[Delta]0/(1 + Exp[(r - R)/a]), r]))^2 + (k/r +          fv/(4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] +             D[\[Delta]0/(1 + Exp[(r - R)/a]), r]))        D[\[Sigma]0/(1 + Exp[(r - R)/a]), r]/(        M - EE + \[Sigma]0/(1 + Exp[(r - R)/a]))) -    D[\[Sigma]0/(1 + Exp[(r - R)/a]), r]/(     M - EE + \[Sigma]0/(1 + Exp[(r - R)/a])) D[G[r], r] +    D[G[r], {r, 2}]; 第二步：用NDSolve数值求解该方程。 sol = NDSolve[{Equation[-66.596, 706.4, 7.030, 0.584, 939, 0,      2, -10] == 0, G[0.001] == 0, G'[0.001] == 1},   G[r], {r, 0.001, 2}] 第三步：Plot图象 g1 = Plot[(G[r])/1 /. sol, {r, 0.001, 2},   PlotStyle -> {Thin, RGBColor[1, 0, 0]},   FrameLabel -> {r, "\!\(\*FormBox[ RowBox[{\" \", RowBox[{\"G\", \"(\", \"r\", \"\"}]}], TraditionalForm]\)"}, Frame -> True, PlotRange -> Automatic] 输出的图象不随参数EE改变，让人失望。

回复此楼