Solving an Improper Integral

49 次查看（过去 30 天）
Ahmad El Kari
Ahmad El Kari 2020-4-8
回答： John D'Errico 2021-3-21
Hello,
I've tried to solve the following integral as such:
syms x;
f = 1/(sqrt(exp(x)-x))
A = int(f, 1, inf)
and the answer was the function itself:
int(1/(exp(x) - x)^(1/2), x, 1, Inf)
Did I make an error because why is the answer not coming out?
Thank you.
  2 个评论
darova
darova 2020-4-8
Maybe it's too complicated for symbolic calculations. What about numerical solution?
syms x;
f = 1/(sqrt(exp(x)-x))
F = matlabFunction(f);
integral(F,1,inf)
2/sqrt(exp(1))
output
ans =
1.3725
ans =
1.2149
doesn't look like truth

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

Deepak Meena
Deepak Meena 2021-3-21
Hi,
You did not make any mistake, To evaluate the A use vpa fucntion :
>> vpa(A)
ans =
1.3725437567734794327982409932998
John D'Errico
John D'Errico 2021-3-21
Ran in:
Think of it like this. Do ALL problems you might pose to a computer have a solution? If so, then I might try
solve('Peace in the middle East')
To no surprise, MATLAB would have a fit if I try it, generating some randomly useless error message.
But many (even most) problems you might pose in mathematics actually have no analytical solution. It is frighteningly easy to pose such a problem. In fact, it looks like you did exactly that.
One solution, and a reason why there are numerical analysis courses taught at many schools, is to learn to use numerical methods to solve intractable problems, where no simple algebraic solution is available. Often numerical analysis is valuable even for problems where a solution exists in theory, yet it is numerically difficult to work with. But this is a probem of the first sort, where no analytical solution apparently exists. At least neither of MATLAB or Wolfram Alpha are able to offer one.
In MATLAB, you can then often tell MATLAB to look for a numerical solution. You can use vpa for that, or you can use vpaintegral.
syms x;
f = 1/(sqrt(exp(x)-x));
int(f,1,inf)
ans = 
But this will work
vpa(int(f,1,inf))
ans = 
1.3725437567734794327982409932998
Or you might do it as
vpaintegral(f,1,inf)
ans = 
1.37254

类别

在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by