2

亲爱的 Mathematica 用户,

我是一名物理专业的学生(对 Mathematica 不太熟练),我必须求解一个积分方程。下面是我尝试的草图:

iterstepGap := (valuesgap = Table[
    {t,
     NIntegrate[
     Exp[2 (tprime - t)] gapStart[tprime]/(3- tprime) EllipticK[Exp[tprime - t]],
       {tprime, -10, t}, MaxRecursion -> 150, WorkingPrecision -> 30] +
      NIntegrate[
        Exp[tprime - t] gapStart[tprime]/(3- tprime) EllipticK[Exp[t - tprime]],
       {tprime, t, 0},  MaxRecursion -> 150, WorkingPrecision -> 30]
     }, {t, -10, 0, 1}];
  approxgap[t_] = InterpolatingPolynomial[valuesgap, t]
)

gapStart[t_] = 0.002 Exp[-t] ((Exp[t] + t)/(t - 1) + 1);

问题是该函数的振荡非常糟糕。正如我所说,我对mathematica 很陌生,我不太确定我在做正确的事情......是否有更好的方法来解决这种方程?

非常感谢您!

4

0 回答 0