MathGroup Archive 2006

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

Search the Archive

Re: Re: Strange empty set of solutions

  • To: mathgroup at smc.vnet.net
  • Subject: [mg71871] Re: [mg71814] Re: Strange empty set of solutions
  • From: Andrzej Kozlowski <akoz at mimuw.edu.pl>
  • Date: Fri, 1 Dec 2006 06:22:29 -0500 (EST)
  • References: <ekh7pg$sgs$1@smc.vnet.net> <ekjfrf$d7h$1@smc.vnet.net> <200611301105.GAA08424@smc.vnet.net>

On 30 Nov 2006, at 20:05, José Carlos Santos wrote:

> On 29-11-2006 8:19, Jens-Peer Kuska wrote:
>
>> and we can't read you mind and the memory of your computer
>> so we must imagine a matrix M and write down
>
> Do you want an example? Here it is:
>
> M = {{4/5(c - 1), -2/Sqrt[5]s, 2/5(1 - c)},
>       {2/5(1 - c), -1 + s/Sqrt[5], 4/5 + c/5},
>       {-2s/Sqrt[5], -c, -1 + s/Sqrt[5]}} // N
>
> with c = Cos[Sqrt[5]] and s = Sin[Sqrt[5]].
>
> If I type
>
> Solve[{M.{x, y, z} == {0, 0, 0}}, {x, y, z}]
>
> I get
>
> {{x -> 0. - 0.0438861 z, y -> 0. + 1. z}}
>
> but if I type
>
> Solve[{M.{x, y, z} == {0, 0, 0}, x^2 + y^2 + z^2 == 1}, {x, y, z}]
>
> then I get the empty set. Why is that?
>
> Best regards,
>
> Jose Carlos Santos
>

First,  is a bad idea to mix arbitrary precision numbers with  
algebraic functions such as Solve. Use the definition:

M = {{4/5(c - 1), -2/Sqrt[5]s, 2/5(1 - c)},
       {2/5(1 - c), -1 + s/Sqrt[5], 4/5 + c/5},
       {-2s/Sqrt[5], -c, -1 + s/Sqrt[5]}} ;


Second, you have three equations and 5 "unknowns". which you can  
arbitrarily divide into "variables" (in your case x,y,z} and  
parameters (s,c}.  This means that you will only get solutions  
provided certain relations hold between the parameters. But Solve  
gives only "generic" solutions for the variables, where generic means  
solutions valid for any values of the parameters. But your equation  
obviously has no such solutions.

There are two things you can do. One is to use Reduce:

Reduce[{M.{x, y, z} == {0, 0, 0}, x^2 + y^2 + z^2 == 1}, {x, y, z}]

will show you also the relations that have to hold between the  
parameters. The other possibility is to treat one of the "parameters"  
as a variable:

Solve[{M.{x, y, z} == {0, 0, 0}, x^2 + y^2 + z^2 == 1}, {x, y, z, c}]

will give you a (pretty awful) list of solutions for the "variables"  
x, y, z and c in terms of the "parameter" s.

Andrzej Kozlowski
Tokyo, Japan


  • Prev by Date: Re: Re: Strange empty set of solutions
  • Next by Date: Re: please explain why ToExpression[SubscriptBox["x","1"],StandardForm] ...
  • Previous by thread: Re: Re: Strange empty set of solutions
  • Next by thread: Re: Strange empty set of solutions