MathGroup Archive 2011

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

Search the Archive

Re: Roots of a Jacobi polynomial

  • To: mathgroup at smc.vnet.net
  • Subject: [mg120644] Re: Roots of a Jacobi polynomial
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Mon, 1 Aug 2011 06:57:55 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Reply-to: hanlonr at cox.net

Restart Mathematica

eqlist = Rationalize[{
    101.74 == 
     a + ((0.055/(479248/c))*(352343 + ((479248/c)*a) - (158194*b/c))),
    47205 + ((158194/c)*a) - ((61842/c)*b) == 68043.98 + (299.73*b),
    460584 + (479248*(a^2)/(c^2)) - (158194*b*
        a/(c^2)) + (61842*(b^2)/(c^2)) == 21986 + (3576*c)}];

NSolve[eqlist, {a, b, c}, Reals][[1]]

{a -> 92.1341, b -> 103.075, c -> 158.517}

And @@ (eqlist /. %)

True


Bob Hanlon

---- Richard Kendall-Smith <richpks at gmail.com> wrote: 

=============
Hi,

I'm trying to solve a set of non-linear equations by doing the following:

eqlist = {101.74 == a + ((0.055/(479248/c))*(352343 + ((479248/c)*a) -
(158194*b/c))),
 47205 + ((158194/c)*a) - ((61842/c)*b) == 68043.98 + (299.73*b),
 460584 + (479248*(a^2)/(c^2)) - (158194*b*a/(c^2)) + (61842*(b^2)/(c^2))
== 21986 + (3576*c)}

NSolve[eqlist, {a, b, c}, Reals]

But it interprets c as Musical Notation. I have tried altering the variables
to x y and z but then it comes up with a Visual Form error. I'm new to
Mathematica but I don't see why this is happening because aren't we
specifying a, b, c as variables to be solved?

Thanks,

Richard





  • Prev by Date: Re: Multiple Styles
  • Next by Date: Re: Compile and Total
  • Previous by thread: Re: Multiple Styles
  • Next by thread: Re: Compile and Total