Help required for logarithmic equation

Status
Not open for further replies.

Delta11

Junior Member level 3
Joined
May 14, 2009
Messages
25
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,454
Need help solving a logarithmic equation:
(1/x) . exp(-x^2) = 10^-6.
Thanks.
 

If you just need a numerical result, you can use the free math tool Scilab. Code you have to write is:

deff('[y]=funct(x)','y=exp(-x^2)/x-1e-6');

[xres]=fsolve(1,funct)

from which you will obtain x=3.5426862
 

    Delta11

    Points: 2
    Helpful Answer Positive Rating
You may solve it by iteratives methods.

(1/x) . exp(-x^2) = 10^-6 --> x_n = SQRT(LN(10^6/x_n_1))

Starting with x = 1 converge in few steps to 3.542686172...
 

    Delta11

    Points: 2
    Helpful Answer Positive Rating
Thanx for the help guys. Is there any other approach to solving this problem other than iterative method? Numerical iteration obviously fits but i am looking for a more conventional solution. It is a second order logarithmic equation after all!!!
 

It cannot be expressed by elementary functions.

You need special functions like the Lambert W-function (see Google).

If exp(-x^2)/x = a then x = sqrt( W(2/a^2) / 2 )

with W() the Lambert W-function

But... to obtain W() you need numerical methods or need software like Mathematica with the ProductLog[z] function...
 

use first few terms(depending on the accuracy) of the taylor's series of the exponential and divide by x. Now it's just solving a polynomial.
 

Status
Not open for further replies.
Cookies are required to use this site. You must accept them to continue using the site. Learn more…