MathGroup Archive 2002

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

Search the Archive

Re: Solving problem

  • To: mathgroup at smc.vnet.net
  • Subject: [mg32387] Re: Solving problem
  • From: Erk Jensen <Erk.Jensen at cern.ch>
  • Date: Fri, 18 Jan 2002 00:17:32 -0500 (EST)
  • References: <a25uuu$dk4$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Nicolas Le Roux wrote:
> 
> I have to solve a system of 4 equations with 4 variables. The problem is
> that one of a variable is a born(?) of an integrale and so Mathematica
> could not know whether or not the integrale is convergent. If you could
> help me to solve this system, it would really help me.
> You'll find the worksheet at www.via.ecp.fr/~genji/test.nb
> 
> Thank you
> 
> --
> Genji
> "La femme ne voit jamais ce qu'on fait pour elle, elle ne voit que ce
> qu'on ne fait pas."
> Georges Courteline

I'm not a guru, but why do you use "set delayed"? If you evaluate the integral
BEFORE you try and solve the equations, you don't get the problem you encounter.
However, I don't believe that you can solve the equations all analytically. What
I did is

\!\(eqns = 
    Simplify[{\@\(h\^2 - e\^2\) - a*e\^2 - b\ e - c == 0, 
        2\ a\ e + b + \(2\ e\)\/\@\(h\^2 - e\^2\) == 0, 
        a\ d\^2 + b\ d + c == 0, 
        L == \[Integral]\_e\%d\(\@\( 1 + b\^2 + 4\ a\ x\ b\  + \ 
                      4\ a\^2\ x\^2\)\) \[DifferentialD]x + 
            h\ \((\[Pi]\/2 - ArcCos[e\/h])\)}]\)

and then solved the 1st 3 equations for {a,b,c}, like

Solve[eqns[[{1, 2, 3}]], {a, b, c}] // Simplify

and insert the result into the forth like

eq4=eqns[[4]] /. %

Now you can plot the RHS - L versus e. This doesn't solve your real problem
however, the function doesn't seem to be well-behaved there ...

Good luck,
     -erk-


  • Prev by Date: Re: Taylor Series Expansions
  • Next by Date: Re: Wrestling with Mathematica on Trig Simplification
  • Previous by thread: Solving problem
  • Next by thread: Comparison of Mathematica on Various Computers