Question

我有一个整体，MATLAB无法得到我的问题的答案。我有两个问题。

有没有办法在MATLAB中计算积分？
这是可以使用其他软件，如Mathematica或Maple 计算这个积分？

b=5;
t0=6;
syms s t

Z=8/5*(4/5+1/5*(25*(911860908983391/2251799813685248*s-911860908983391/2251799813685248)^2*exp(-5505815247359351/72057594037927936*s)^2-20*(911860908983391/2251799813685248*s-911860908983391/2251799813685248)*exp(-5505815247359351/72057594037927936*s)+16)^(1/2))^2*(((911860908983391/2251799813685248*s-911860908983391/2251799813685248)*exp(-5505815247359351/72057594037927936*s)+4/5)^2-(911860908983391/2251799813685248*s-911860908983391/2251799813685248)^2*exp(-5505815247359351/72057594037927936*s)^2+4/5*(911860908983391/2251799813685248*s-911860908983391/2251799813685248)*exp(-5505815247359351/72057594037927936*s)-16/25)^(1/2)/((911860908983391/2251799813685248*s-911860908983391/2251799813685248)*exp(-5505815247359351/72057594037927936*s)+1/5*(25*(911860908983391/2251799813685248*s-911860908983391/2251799813685248)^2*exp(-5505815247359351/72057594037927936*s)^2-20*(911860908983391/2251799813685248*s-911860908983391/2251799813685248)*exp(-5505815247359351/72057594037927936*s)+16)^(1/2))^2

V=int(Z,s,t0-b,t)

Answer 1

假设我已正确放置指数函数括号，Mma 10.4只能找到数值近似的答案，例如

z = 8/5*(4/5 + 
      1/5*(25*(911860908983391/2251799813685248*s - 
              911860908983391/2251799813685248)^2*
           Exp[-5505815247359351/72057594037927936*s]^2 - 
          20*(911860908983391/2251799813685248*s - 
             911860908983391/2251799813685248)*
           Exp[-5505815247359351/72057594037927936*s] + 16)^(1/
          2))^2*(((911860908983391/2251799813685248*s - 
             911860908983391/2251799813685248)*
           Exp[-5505815247359351/72057594037927936*s] + 
          4/5)^2 - (911860908983391/2251799813685248*s - 
           911860908983391/2251799813685248)^2*
        Exp[-5505815247359351/72057594037927936*s]^2 + 
       4/5*(911860908983391/2251799813685248*s - 
          911860908983391/2251799813685248)*
        Exp[-5505815247359351/72057594037927936*s] - 16/25)^(1/
       2)/((911860908983391/2251799813685248*s - 
          911860908983391/2251799813685248)*
        Exp[-5505815247359351/72057594037927936*s] + 
       1/5*(25*(911860908983391/2251799813685248*s - 
               911860908983391/2251799813685248)^2*
            Exp[-5505815247359351/72057594037927936*s]^2 - 
           20*(911860908983391/2251799813685248*s - 
              911860908983391/2251799813685248)*
            Exp[-5505815247359351/72057594037927936*s] + 16)^(1/2))^2;

NIntegrate[z, {s, 1, 2}]

2.67636

MATLAB中一个令人不安的缺陷

1 个答案: