MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Mathematica issue

  • To: mathgroup at smc.vnet.net
  • Subject: [mg127861] Re: Mathematica issue
  • From: "Dr. Wolfgang Hintze" <weh at snafu.de>
  • Date: Tue, 28 Aug 2012 04:50:22 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-newout@smc.vnet.net
  • Delivered-to: mathgroup-newsend@smc.vnet.net
  • References: <k1epvg$rnq$1@smc.vnet.net>

On 27 Aug., 05:36, Matthias Bode <lvs... at hotmail.com> wrote:
> Hola:
> As neither Solve[] nor Reduce[] find an analytical solution there isn't any - most probably.
> One might, however, do - inter alia - this:
> Clear[a0, a1, b0, b1, c0, d0, eq01, sol01]eq01 = Expand[a0 + b0 +   x^a1*(c0 - d0*x^b1) == 0]; sol01 = FindInstance[eq=
01, {a0, b0, c0, d0,        a1, b1, x}, Reals]a0 + b0 + x^a1*(c0 - d0*x^b1) == 0 /. sol01N[sol01]
> Out[3]= {{a0 -> (46314691670528 - 1684956133536*           134^(7/10) - 15*134^(4/5))/         57893364588160, b0 -> -(4/5),      c0 -> 39/10, d0 -> -(3/2),      a1 -> -(3/10), b1 -> -(49/10), x -> 134}}
> Out[4]= {True}
> Out[5]= {{a0 -> -0.09729070208467958, b0 -> -0.8,      c0 -> 3.9,d0 -> -1.5, a1 -> -0.3,      b1 -> -4.9, x -> 134.}}
> Best regards,
> MATTHIAS BODE
> S 17.35775=B0, W 066.14577=B0
> 2'740 m
> AMSL.
>
>
>
>
>
> > From: n... at 12000.org
> > Subject: Re: Mathematica issue
> > To: mathgr... at smc.vnet.net
> > Date: Sun, 26 Aug 2012 05:45:30 -0400
>
> > On 8/26/2012 3:18 AM, Nasser M. Abbasi wrote:
>
> > >> I need analytical expression for the x in terms of A,B,C,D and a,b,
>
> > > try this:
>
> > > parms = {A0 -> 1, B0 -> 2, a -> 3, C0 -> 4, D0 -> 5, b -> 6}
> > > eq    = A0 + B0 + x^a (C0 - D0 x^b) == 0
> > > sol   = Solve[eq /. parms, x]
>
> > Opps, just noticed you want symbolic solution.
> > Mathematica 8.04 does not do it. I doubt this can be solved symbolicall=
y.
> > Need to use some numbers for the parameters.
>
> > Clear[A0, B0, a , C0, D0 , b]
> > eq = A0 + B0 + x^a (C0 - D0 x^b) == 0
> > Solve[eq, x]
>
> > Solve::nsmet: This system cannot be solved with the methods available t=
o =
> Solve. >>
>
> > --Nasser

If you restrict yourself to integer exponents, the solubility of the
equation is governed by the Abel=96Ruffini theorem:

In algebra, the Abel=96Ruffini theorem (also known as Abel's
impossibility theorem) states that there is no general algebraic
solution=97that is, solution in radicals=97 to polynomial equations of
degree five or higher.

You might want to watch this in Mathematica usind this command

Table[{n, eq = x^n + a*x + b == 0; Simplify[Solve[eq, x]]}, {n, 1, 5}]

More generally you are facing a transcendental quation which has no
closed solution in most cases but must be soved numerically.
Here's an example:

In[57]:= eq1 = x^(-E) + x^Pi == 3

Out[57]= x^(-E) + x^Pi == 3

Watch the behaviour of the function

In[66]:= Plot[{3, x^(-E) + x^Pi}, {x, 0, 2}]
[... the graph of the function]

and then find the solutions (roots) numerically

In[61]:= FindRoot[eq1, {x, 1}]

Out[61]= {x -> 1.34813204709293}

In[62]:= FindRoot[eq1, {x, 0.5}]

Out[62]= {x -> 0.6958252736289695}

Regards,
Wolfgang



  • Prev by Date: Re: Ambiguity of "Plot"
  • Next by Date: Re: ContourPlot non rectangular evaluation?
  • Previous by thread: Re: Mathematica issue
  • Next by thread: Re: Mathematica issue