MATLAB Answers

0

Integral returning infinity for a definite integral with a finite value

Asked by Jack Joynson on 8 Dec 2018 at 20:19
Latest activity Answered by Star Strider
on 8 Dec 2018 at 21:18
Hi,
I am trying to calculate an integral in a loop however it is just returning infinity for every value; The integral in question is finite - I have tested it in wolfram alpha.
Here is my code:
ytheory = 1:10;
s = 14;
func = @(L) s.^(-0.3) .* L.^(-0.7) .*exp(-L./s);
for(i = 1:10)
ytheory(i) = integral(func,0,i);
end
Is this an issue with Matlab or am I doing something wrong here?
Cheers

  0 Comments

Sign in to comment.

1 Answer

Answer by Star Strider
on 8 Dec 2018 at 21:18

I can’t reproduce that in R2018b.
Your code works as expected:
s = 14;
func = @(L) s.^(-0.3) .* L.^(-0.7) .*exp(-L./s);
for i = 1:10
ytheory(i) = integral(func,0,i);
end
producing:
ytheory =
1.4858 1.8004 2.0019 2.1496 2.2646 2.3576 2.4346 2.4994 2.5547 2.6023

  0 Comments

Sign in to comment.