Simple Integration
1 次查看(过去 30 天)
显示 更早的评论
Hi there,
I am new to MATLAB and have a simple question regarding the integration.
I have two expressions:
syms c L;
first = int( exp(-(L-3)^2/2, c, Inf);
second = int (first * exp(-(c-5)^2/.5), -Inf, Inf);
It gives me the integral of the first expression however on the second it says "Explicit integral could not be found". Then I also tried with quad but no avail.
what should I do to evaluate the second integral?
Many Thanks.
0 个评论
采纳的回答
Teja Muppirala
2011-5-13
QUADGK works better for integrals that go to infinity:
syms c L;
first = int( exp(-(L-3)^2/2), c, Inf);
second = matlabFunction (first * exp(-(c-5)^2/.5));
quadgk(second,-inf,inf)
更多回答(1 个)
Ivan van der Kroon
2011-5-13
Bibha,
It get this answer for your third line (note also that on line two a right bracket is missing).
second =
infinity
/
| 2
| exp(-2.000000000 (c - 5) )
|
/
-infinity
1/2 1/2 1/2 1/2 1/2 1/2
(- 1/2 pi 2 erf(1/2 2 c - 3/2 2 ) + 1/2 pi 2 ) dc
I really do not think this integral can be evaluated explicitly. It is an error-function multiplied with a exponential function.
Ciao, Ivan
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!