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

10 views (last 30 days)

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

1 Answer

Link

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

You are now following this question

Jack Joynson (view profile)