Efficient Numerical Method for Inverse Integration with Known Constants a, b, c

  • Context: Graduate 
  • Thread starter Thread starter jimbo_durham
  • Start date Start date
  • Tags Tags
    Integral Inverse
Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
3 replies · 6K views
jimbo_durham
Messages
13
Reaction score
0
I have a complicated integral which i need to compute numerically. I can do this in C++ using a version of Simpson's rule. I also need to compute the inverse of this integral (presumably this is what it is called) ie i have

d=integral(f(x)dx)

and i need to be able to compute x given a value of d.

The function f(x) is too complicated to simply rearrange, so i need an iterative way of guessing x, running it in the program, seeing what value of d is given, and then improving my guess of x until i reach my starting d value.

This seems on the surface to be a simple iteration problem, however i cannot find an efficent way of doing this.

Can anyone tell me of a nice way of solving this problem (and even if it is called an 'inverse integration' or not?) and if there is a name given to the numerical method of solving it.

It has been mentioned that there is a 'press' method however i can find no mention of this. Prehaps it involves finding the roots of an equation which takes the value of the integral evaluated at some guessed x, and the roots give the point where the guessed x corresponds to the known d, thus solving the problem? does this ring any bells?

sorry for the long post, thanks in advance for your help
 
Last edited:
Physics news on Phys.org
Well, one way would be:
Let:
[tex]y(x)=y(0)+\int_{x_{0}}^{x}f(t)dt[/tex]
Or, that is:
[tex]y=y_{0}+\int_{x(0)}^{x(y)}f(t)dt[/tex]
Differentiating the latter expression with respect to y yields:
[tex]1=f(x(y))\frac{dx}{dy}[/tex]
That is, you may solve the following differential equation numerically:
[tex]\frac{dx}{dy}=\frac{1}{f(x)}, x(y_{0})=x(y_{0})[/tex]
 
i am not sure how to apply that, i think it is worth me giving you the integral;



[tex] \begin{equation}<br /> <br /> d_{M}=c_{1}\cdot sinh\left[c_{2}\int^{b}_{a}\left[\left(1+z\right)^{2}\cdot\left(1+c_{3}\cdot z\right)-z\cdot\left(2+z\right)\cdot c_{4}\right]^{-\frac{1}{2}}dz\right]<br /> <br /> \end{equation}[/tex]

i am trying to find z given a value of [tex]d_{M}[/tex].
 
note in formula, [tex]a, b, c[/tex] are known constants