How to treat with logarithm numerically when it becomes negative?
I have a term log(f(x)) in my equations which I am trying to solve them numerically. But the solver does not solve, perhaps when the argument of the logarithm becomes negative. What should I write in the program when f(x)<0 in the loop in order not to loose the iteration?
Thanks
I think you have to provide more details. In general, log of a negative number does make sense, but it is multi-valued.
you mean detail in the program?
Have you tried: find_root(log(x*x-x),-2,0) -0.6180339887498988 find_root(log(x*x-x),-2,2) 1.6180339887498947