Finding the inverse of a nasty 1-to-1 function.

  • Thread starter jforeman83
  • Start date
  • #1
jforeman83
4
0
Hey, everybody.

I have a function:
[tex]
\int\limits_{x}^{x+c} exp(-t^2) dt = y
[/tex]

c is a known constant here.

I am beating my head against the wall trying to find a good way to numerically evaluate the inverse here, i.e. I have y and c and I want to know x. I know that erf^-1 is readily available in mathematica and maple and the like but the limits of integration here make this a bit nastier. Any ideas? I don't need a perfect evaluation, just a moderately good approximation will work.
 

Answers and Replies

  • #2
quasar987
Science Advisor
Homework Helper
Gold Member
4,793
21
If you're interested only in values of x smaller than 1, then you can always Taylor expand exp(-t²) and integrate term by term and keep only the first 2 terms. I get y=c-xc+c²/2.
 
  • #3
lzkelley
277
2
i think the rigorous method would be to use a Fourier inversion expansion; can't remember the details of how to do that - but with a smooth gaussian function i think its not that bad.
 
  • #4
jforeman83
4
0
If you're interested only in values of x smaller than 1, then you can always Taylor expand exp(-t²) and integrate term by term and keep only the first 2 terms. I get y=c-xc+c²/2.

I wish I could say with certainty that this was the case, because you're right, that would be a good idea. Unfortunetly, I think I'll need something that covers a bit more values.
 
  • #5
jforeman83
4
0
i think the rigorous method would be to use a Fourier inversion expansion; can't remember the details of how to do that - but with a smooth gaussian function i think its not that bad.

Do you have any recommendations on sources where I might read up on this?
 
  • #6
ObsessiveMathsFreak
406
8
This function does not have a proper inverse as y(x)=y(-x-c)

I believe that the inverse may exist if you restrict x to be greater than zero.
 
  • #7
Gib Z
Homework Helper
3,352
6
I think we can view this problem as a differential equation, and use some method of numerically solving ODE's such as the midpoint method or Runge-Kutta? Not 100% sure how it'll work out though.
 
  • #8
daudaudaudau
302
0
If you only want to evaluate it numerically, then why don't you just use Newton-Rhapson to determine when y(x) = K ?
 
  • #9
Gib Z
Homework Helper
3,352
6
That involves evaluating the integral, which can't be done analytically though I do know there are many tables of data for that particular function (The error function). So Yes, I guess that method can do it numerically, good idea =]
 
  • #10
jforeman83
4
0
If you only want to evaluate it numerically, then why don't you just use Newton-Rhapson to determine when y(x) = K ?

This is a fantastic idea. Thanks!
 

Suggested for: Finding the inverse of a nasty 1-to-1 function.

  • Last Post
Replies
3
Views
452
  • Last Post
Replies
16
Views
746
  • Last Post
Replies
5
Views
524
Replies
6
Views
941
Replies
5
Views
416
  • Last Post
Replies
6
Views
653
  • Last Post
Replies
2
Views
584
Replies
5
Views
482
Replies
44
Views
2K
Top