求解Mathematica中的积分方程

时间:2012-12-06 22:14:22

标签: wolfram-mathematica integration

亲爱的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非常新,我不太确定我做的是正确的事情......是否有更好的方法可以解决这种方程?

非常感谢你!

0 个答案:

没有答案