Hello guys, can you help me with my problem in this thread?

http://community.wolfram.com/groups/-/m/t/455321

in that discussion, the formula is just an example of my work using ParametricNDSolveValue, the problem surfaced when FindRoot have an answer of complex value. is there anyway to force FindRoot to present only Real numbers?

However, here you can download the obsolete add-on package which did the job I actually need. Works for version 6

http://library.wolfram.com/infocenter/MathSource/6771/

At a certain point it says

"The RealOnly package is obsolete. Use instead functionality provided by \
Reduce[equations, variables, Reals] and by Assuming, Refine, and similar \
functions."

I tried several combinations of Reduce or Assuming without success :-(

Attachments:

Luther, no, probably not an error on your part. Some integrals do not converge in the limit x -> infinity. For instance, solutions for an oscillating system will cycle around a stable state (and have "imaginary" parts). Solutions for an oscillating system, driven at resonant frequency (without proper damping) will diverge! In engineering, this can have catastrophic results.

Have an example, or a description of the system you are trying to model?

My take away is that I CAN ignore the imaginary component, even if it is very large.

Not always, especially if it is large. You can safely ignore imaginary components whose magnitude is small in comparison to your expected/desired numerical precision. On the other hand, a large imaginary component in a result you'd expect to be real is likely a symptom that you've wandered off onto the wrong branch of a multivalued function. The result may be mathematically correct, but in a sense that doesn't apply to your physical problem.

I'm a physicist. I'm used to mathematics needing guidance to properly represent physical reality. It starts with very simple high school problems. Consider dropping a ball: how long will it take to hit the floor. Well, it's the solution to an algebra problem:

In[24]:= Solve[ 1/2 a t^2 == h, t]

Out[24]= {{t -> -((Sqrt[2] Sqrt[h])/Sqrt[a])}, {t -> (
   Sqrt[2] Sqrt[h])/Sqrt[a]}}

Oops, there are two solutions! However, the ball can't hit the floor before you drop it, so we throw out the negative one. This is completely normal.

Now, for algebra problems there are "root isolation" methods that find all solutions rigorously, and Mathematica uses them. So, you throw out the solutions you don't want and you have the solutions for your problem. However, once you go beyond algebra, you generally can't have confidence that your method, even if it has found a solution, has found every solution to your mathematical problem, nor can you be certain that the solution it has found is the solution that you want for your physical problem. You're in the realm of "experimental verification".

Comparing results from Integrate and NIntegrate is a good sort of mathematical experiment, as they have different sorts of problems. You've already encountered Integrate's branch cut trouble. NIntegrate has the problem that if its sampling of points in the function in question accidentally misses the exciting parts, it may yield an entirely bogus result.

Mathematica will not do your thinking for you.

The basic mathematical issue is that the grand principles of algebra and analysis operate in the complex domain. Often, you have to reason through the complex domain to get a real answer. Any attempt to restrict Mathematica to operate in the real domain will generally impact its ability to break a problem down into subproblems. It is, however, often possible to ask for a real result:

In[20]:= Solve[x^3 == -10, Element[x, Reals]]

Out[20]= {{x -> Root[10 + #1^3 &, 1]}}

In[21]:= % // N

Out[21]= {{x -> -2.15443}}

The Root object is yet another way to express "the real cube root of -10".

Yes, that is the answer. Thanks!

I understand that this is a little off-topic, but for some root operations you can use

In[1]:= ?Surd                                                                                                      
Surd[x, n] gives the real-valued n'th root of x.

In[2]:= ?CubeRoot                                                                       
CubeRoot[x] gives the real-valued cube root of x.

However,

psurd2 = Plot[y = Surd[x, 3]; y, {x, -10, 10}]

Thank you, John! That solves a mystery.

Another way to look at it: for negative real numbers, there are several solutions coming from the complex roots. For instance, in the example above...

in:= solutions = Solve[x^3 == -10, x] // N
out:= {{x -> 1.07722 + 1.8658 I}, {x -> -2.15443}, {x -> 
   1.07722 - 1.8658 I}}

These solutions all have the same magnitude r[x]:

in:= Clear[r]
r[x_] := Sqrt[Re[x]^2 + Im[x]^2]
r@Values@solutions
out:= {{2.15443}, {2.15443}, {2.15443}}

If you don't care to differentiate from among these solutions (though in my opinion, you'd be missing out on the best part!!) you could keep a single, real-valued magnitude:

in:= -1*Flatten@Union@%
out:= {-2.15443}

This yields a similar result as as plotting psurd[] above...

Plot[Piecewise[{{-1*Sqrt[Re[x]^2 + Im[x]^2]^(1/3), x < 0}, {x^(1/3), 
    x > 0}}], {x, -10, 10}]

If we have to cite Euler every time we use his discoveries we will all spend half our time writing citations. ;-)

Happy Pi Day!

Thanks.

What I would like to obtain, if it exists, is a general command that sets Mathematica to work only on the real branches of solutions. In the old Derive program there is such a command : "Branch ? Real"

I know that I could arrange every plot or calculation for every single case, for example in the plot I mentioned I could do

Plot[Sign[x]*Abs[x]^(1/3), {x, -10, 10}]

and I would obtain my result.

But I'd like something general, so to set it up and go ahead without particular arrangements of the functions involved.

Yes, that was my first try to solve the problem......

Ah . . . It chooses one of the complex roots:

In[16]:= x /. Solve[x^3 == -10, x] // N

Out[16]= {1.07722 + 1.8658 I, -2.15443, 1.07722 - 1.8658 I}