Value of x that maximises the integral in a numerical integrator

Click For Summary
SUMMARY

The discussion centers on the challenge of determining the value of x that maximizes the product of two functions, f(x) and g(x), within the context of a definite integral, specifically $$\int_a^b dx f(x) g(x)$$. Participants emphasize that maximizing an integral over a range does not equate to finding a single x value, as the integral represents an accumulation of values. The GNU Scientific Library (GSL) provides methods for finding maxima of functions, which may be more appropriate than attempting to derive a maximum from the integral itself. Clarification from the original poster (OP) regarding their specific goals is deemed necessary for further discussion.

PREREQUISITES
  • Understanding of definite integrals and their properties
  • Familiarity with the GNU Scientific Library (GSL) and its optimization functions
  • Knowledge of calculus concepts, particularly antiderivatives and maxima
  • Basic programming skills in C++ for implementing numerical methods
NEXT STEPS
  • Explore GSL optimization methods for finding function maxima
  • Study the properties of definite integrals and their applications
  • Learn about antiderivatives and their significance in calculus
  • Investigate numerical methods for adaptive integration in C++
USEFUL FOR

Mathematicians, data scientists, and software developers interested in numerical integration and optimization techniques, particularly those using the GNU Scientific Library in C++.

CAF123
Gold Member
Messages
2,918
Reaction score
87
Consider an integral of form $$\int_a^b dx f(x) g(x).$$ Is it possible to tell a numerical integrator to spit out the value of ##x \in [a,b]## that maximises the value of ##f(x)g(x)##? I'm mostly interested in incorporating this into some code I have for adaptive integrator gsl_qags in C++.

Thanks!
 
Technology news on Phys.org
This link directs me simply to google -> in any case I suppose the gsl_qags will use some stepwise increment over the integration domain. Is it maybe simpler to print out the values sampled in the integration domain together with the value of the integrand? Then I can see the maximum by inspection. Thanks.
 
@DrClaude Thanks, sorry I see I was not careful with my expression, I meant the value of the integral in my above post.
 
Sorry, but I don't understand. Let's start from the beginning: what function do you want to maximise, and with respect to what parameter/variable?
 
  • Like
Likes pbuk
CAF123 said:
I meant the value of the integral in my above post.

That doesn't make sense. You originally asked:

CAF123 said:
Is it possible to tell a numerical integrator to spit out the value of ##x \in [a,b]## that maximises the value of ##f(x)g(x)##?

If you meant maximizes the integral, then you would be asking for a particular value of ##x## that maximizes an integral over a range of values of ##x##. That doesn't make sense.
 
PeterDonis said:
That doesn't make sense. You originally asked:
If you meant maximizes the integral, then you would be asking for a particular value of ##x## that maximizes an integral over a range of values of ##x##. That doesn't make sense.

I suppose it might make sense if he/she wants---for some unknown reason--- to spit up the integral as
$$\int_a^b f(x) g(x) \, dx = \int_a^m f(x) g(x) \, dx + \int_m^b f(x) g(x) \, dx,$$
where ##m = \text{arg max}_{x \in [a,b]} \, f(x) g(x).##

Even then the problem may be ill-posed, because for some nasty functions ##f,g## the integrand could have hundreds of values of ##m##--that is, hundreds of points that maximize the integrand.
 
Ray Vickson said:
I suppose

Rather than suppose, I would like the OP to give more information about what the actual problem is.
 
  • #10
If looking for a maximum x for f(x)g(x), and assuming that f(x)g(x) has a maximum, then the maximum value will occur at one of the roots of the derivative of f(x)g(x).

If looking for the maximum of the indefinite integral of f(x)g(x), and assuming that there is a maximum, then the maximum value will occur at one of the roots of f(x)g(x).

If looking for a maximum integral range, and if the entire range of the integral results in "positive" areas, then the maximum is the entire range. Otherwise you're looking for a range with a mix of "positive" and "negative" areas that has the maximum "positive" area, which seems complicated.
 
  • #11
rcgldr said:
If looking for the maximum of the indefinite integral of f(x)g(x), and assuming that there is a maximum, then the maximum value will occur at one of the roots of f(x)g(x).

This doesn't make sense. You don't maximize an integral over a range of values of ##x## by picking one value of ##x##. See my response to the OP in post #7.

And in fact, you can't maximize an indefinite integral at all since it's indefinite. You might mean maximize the function of ##x## that is the antiderivative of the function ##f(x) g(x)##, but that's not quite the same thing.
 
  • #12
PeterDonis said:
you can't maximize an indefinite integral at all since it's indefinite. You might mean maximize the function of ##x## that is the antiderivative of the function ##f(x) g(x)##, but that's not quite the same thing.
I meant antiderivative. "In calculus, an antiderivative, primitive function, primitive integral or indefinite integral ... "

https://en.wikipedia.org/wiki/Antiderivative
 
  • #13
rcgldr said:
I meant antiderivative. "In calculus, an antiderivative, primitive function, primitive integral or indefinite integral ... "

https://en.wikipedia.org/wiki/Antiderivative

The OP started this all by stating a definite integral ##\int_a^b f(x) g(x) \, dx.##
 
  • #14
I am closing this thread because there is no point in speculating about what the OP meant. The OP needs to tell us, and that hasn't happened.

@CAF123 if you want to clarify what you meant, please PM me; if I get sufficient clarification I can reopen the thread for discussion.
 

Similar threads

  • · Replies 6 ·
Replies
6
Views
7K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 15 ·
Replies
15
Views
3K
  • · Replies 31 ·
2
Replies
31
Views
4K
  • · Replies 50 ·
2
Replies
50
Views
6K
  • · Replies 21 ·
Replies
21
Views
4K
  • · Replies 19 ·
Replies
19
Views
3K
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 6 ·
Replies
6
Views
3K
  • · Replies 4 ·
Replies
4
Views
6K