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-