Re: RE: [Q] Equation solving?

*To*: mathgroup at smc.vnet.net*Subject*: [mg23335] Re: [mg23318] RE: [mg23282] [Q] Equation solving?*From*: Andrzej Kozlowski <andrzej at tuins.ac.jp>*Date*: Thu, 4 May 2000 02:59:05 -0400 (EDT)*Sender*: owner-wri-mathgroup at wolfram.com

David, In[1]:= eqns = {Q1 == (2*z + 1 - 1)^2 + 4*1*1, k1 == 1/Sqrt[Q1]*(v*z - 3/2*(2*z + 1 + 1 - Sqrt[Q1])), Q2 == (2*z + 1 - 1)^2 + 4*1*1, k2 == 1/Sqrt[Q2]*((6.1 - v)*z - 2/2*(2*z + 1 + 1 - Sqrt[Q2])), 3/2*(2*z + 1 + 1) + (k1 - 3/2)*Sqrt[Q1] == 6.1*z - (2/2*(2*z + 1 + 1) + (k2 - 2/2)*Sqrt[Q2])}; In[2]:= eqns2 = Eliminate[eqns, {Q1, Q2}]; In[3]:= ans = Solve[eqns2, {k1, k2, z}] Out[3]= {{z -> 0., k1 -> 0., k2 -> 0.}, {z -> 0., k1 -> 0., k2 -> 2.}, {z -> 0., k1 -> 3., k2 -> 0.}, {z -> 0., k1 -> 3., k2 -> 2.}} However, this is not the end of the story! Now let's try to verify the answer. In[4]:= eqns /. ans // FullSimplify Out[4]= 2. 2. {{Q1 == 4., 1. == --------, Q2 == 4., 1. == --------, Sqrt[Q1] Sqrt[Q2] 5. == 1.5 Sqrt[Q1] + 1. Sqrt[Q2]}, 2. 2. {Q1 == 4., 1. == --------, Q2 == 4., 1. + -------- == 0, Sqrt[Q1] Sqrt[Q2] 5. + 1. Sqrt[Q2] == 1.5 Sqrt[Q1]}, 2. 2. {Q1 == 4., 1. + -------- == 0, Q2 == 4., 1. == --------, Sqrt[Q1] Sqrt[Q2] 5. + 1.5 Sqrt[Q1] == 1. Sqrt[Q2]}, 2. {Q1 == 4., 1. + -------- == 0, Q2 == 4., Sqrt[Q1] 2. 1. + -------- == 0, 5. + 1.5 Sqrt[Q1] + 1. Sqrt[Q2] == 0}} Sqrt[Q2] Just look at the last equation: 5. + 1.5 Sqrt[Q1] + 1. Sqrt[Q2] == 0 . Clearly this is not possible. In general you can't rely on any solutions Mathematica gives you to equations of this sort. This is not a bug: spurious solutions naturally arise as part of the slving process and since (in general) it is impossible for Mathematica to verify their correctness (in equations with parameters) it does not even try to do so. Andrzej Kozlowski on 5/2/00 7:00 AM, David Park at djmp at earthlink.net wrote: > > >> From: James [mailto:research at proton.csl.uiuc.edu] To: mathgroup at smc.vnet.net > To: mathgroup at smc.vnet.net >> Hello! >> >> I'm a newbie to Mathematica >> and have a question about what's wrong in the following >> expression. >> (The parameters are generated by the program, >> so there're some terms that can be cancelled out.) >> >> Solve[{ Q1==(2*z+1-1)^2+4*1*1 , >> k1==(1/Sqrt[Q1])*(v*z-(3/2)*(2*z+1+1-Sqrt[Q1])) , >> Q2==(2*z+1-1)^2+4*1*1 , >> k2==(1/Sqrt[Q2])*((6.10-v)* z-(2/2)*(2*z+1+1-Sqrt[Q2])) , >> (3/2)*(2*z+1+1)+(k1-(3/2))*Sqrt[Q1] == >> 6.10*z-((2/2)*(2*z+1+1)+(k2-(2/2))*Sqrt[Q2])}, {z,v}] >> >> I want to solve the above equations according to 'z', >> and I tried several ways in expressing {z,v} or z, {v}, etc >> to get an answer(s), but all in vain. >> >> Is there anyone who can tell me how to express the equations? >> Sorry about the question, if this is trivial. >> Thank you. >> >> >> James. >> > > This seems to work: > > eqns = {Q1 == (2*z + 1 - 1)^2 + 4*1*1, > k1 == 1/Sqrt[Q1]*(v*z - 3/2*(2*z + 1 + 1 - > Sqrt[Q1])), Q2 == (2*z + 1 - 1)^2 + 4*1*1, > k2 == 1/Sqrt[Q2]*((6.1 - v)*z - > 2/2*(2*z + 1 + 1 - Sqrt[Q2])), > 3/2*(2*z + 1 + 1) + (k1 - 3/2)*Sqrt[Q1] == > 6.1*z - (2/2*(2*z + 1 + 1) + (k2 - 2/2)*Sqrt[Q2])}; > > Eliminate[eqns, {Q1, Q2}] > Solve[%, {k1, k2, z}] > > -3.*k1 + k1^2 == 0. && -2.*k2 + k2^2 == 0. && z == 0. > > {{z -> 0., k1 -> 0., k2 -> 0.}, {z -> 0., k1 -> 0., > k2 -> 2.}, {z -> 0., k1 -> 3., k2 -> 0.}, > {z -> 0., k1 -> 3., k2 -> 2.}} > > Or... > > Eliminate[eqns, {k1, k2}] > Solve[%, {Q1, Q2, z}] > > Q1 == 4. && Q2 == 4. && z == 0. && Q1 != 0. && Q2 != 0. > > {{Q1 -> 4., Q2 -> 4., z -> 0.}} > > It appears that v drops out of the equations. > > David Park > djmp at earthlink.net > http://home.earthlink.net/~djmp/ >