Problem 58259. Easy Sequences 115: Integral involving square root, floor, and round functions
Given a postive real number n, we are asked to evaluate the following integral:
We may rewrite the above function in Matlab as:
>> S = @(n) round(integral(@(x) sqrt(floor(x))-floor(sqrt(x)),0,n));
Therefore, for , we have:
>> s = S(10*pi)
s =
12
Be careful though, in using the Matlab integral function, as it is only an approximation. For example if :
>> s = S(100000)
Warning: Reached the limit on the maximum number of intervals in use. Approximate bound on error is
8.0e+01. The integral may not exist, or it may be difficult to approximate numerically to the
requested accuracy.
> In integralCalc/iterateScalarValued (line 372)
In integralCalc/vadapt (line 132)
In integralCalc (line 75)
In integral (line 87)
In @(n)round(integral(@(x)sqrt(floor(x))-floor(sqrt(x)),0,n))
s =
49841
The correct answer is . The integral function is off by only 2 units here, but the discrepancy could be even greater for other values of n. integral function can also be quite slow. The challenge is to find an efficient and more accurate algorithm to evaluate the integral.
Solution Stats
Problem Comments
-
1 Comment
Ramon Villamangca
on 18 May 2023
Hi David,
You right I should have use floor (34) instead of round (35). Test suites were adjusted. Thanks.
Solution Comments
Show commentsProblem Recent Solvers1
Suggested Problems
-
1561 Solvers
-
Lychrel Number Test (Inspired by Project Euler Problem 55)
97 Solvers
-
Geometry: Find Circle given 3 Non-Colinear Points
61 Solvers
-
97 Solvers
-
172 Solvers
More from this Author116
Problem Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!