MATLAB中一个麻烦的积分

Question

我有一个积分，MATLAB无法获得我的问题的答案。我有两个问题。

1. 有没有一种方法可以在MATLAB中计算积分？

2. 这是吗可以使用其他软件，如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