MathGroup Archive 2006

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

Search the Archive

Re: solving with inverse functions/inexact coefficients

  • To: mathgroup at smc.vnet.net
  • Subject: [mg67403] Re: [mg67397] solving with inverse functions/inexact coefficients
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Thu, 22 Jun 2006 06:21:18 -0400 (EDT)
  • Reply-to: hanlonr at cox.net
  • Sender: owner-wri-mathgroup at wolfram.com

eqns={F2 == Exp[a* NU]*F1,F1 == F2/2* Exp[a* NU], F2/2+NU == 15};

soln=Reduce[eqns,{F1,F2,NU},Reals]//Simplify

(F1 == 0 && F2 == 0 && NU == 15) || (F2 + 2*NU == 30 && (Sqrt[2]*(30*a - Log[2]))/a == 
    2*F1 && ((Sqrt[2]*Sqrt[F1^2] == F2 && (30*a > Log[2] || a < 0)) || 
    (F2 + Sqrt[2]*Sqrt[F1^2] == 0 && 0 < a && 30*a < Log[2])))

soln2=Solve[soln[[2]]/.a->76/100,{F1,F2,NU}][[1]]

{NU -> (25*Log[2])/38, F2 -> (5/19)*(114 - 5*Log[2]), 
  F1 -> (5/38)*(114*Sqrt[2] - 5*Sqrt[2]*Log[2])}

%//N

{NU -> 0.45601788194733245, F2 -> 29.087964236105336, F1 -> 20.56829676226186}


Bob Hanlon

---- Dave Spence <dspence at ics.mq.edu.au> wrote: 
> I'm trying to solve a set of simultaneous equations describing a laser cavity. A much simplified set of equations in Real unknowns F1, F2 and NU looks like this:
> 
> F2  == Exp[0.76 NU] F1, F1 == 0.5 F2 Exp[0.76 NU], 0.5 F2 + NU == 15
> 
> They can be solved by hand very easily with two solutions, {F1=F2=0 && NU = 15} and {NU=0.46, F1=14.5,  F2=29.1}
> 
> Now Solve cannot solve these equations unless I change both instances of 0.76 to integers. Reduce seems to hang (or at least take too long to return and answer). Reduce can solve the problem quickly if I change both instances of 0.76 to 0.75! Reduce gives the message "Reduce::ratnz: Reduce was unable to solve the system with inexact coefficients. The answer was obtained by solving a corresponding exact system and numericizing the result"
> 
> I guess this is just a difficulty with the way mathematica attempts to solve problem. There are references in the manuals to difficulties solving multiple simultaneous equations that all invoke inverse functions. Are there any options I can set to enable mathematica to quickly solve these equations?
> 
> I note that while these equations can be manipulated by hand into a form that can be solved by mathematica, by full set of equations are far more complex and can't easily be edited by hand.
> 


  • Prev by Date: Re: Re: Package writing
  • Next by Date: Style sheet for 11 x 17
  • Previous by thread: solving with inverse functions/inexact coefficients
  • Next by thread: How to simplify exponents?