MathGroup Archive 2001

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

Search the Archive

Re: Sum of Squares

  • To: mathgroup at smc.vnet.net
  • Subject: [mg26615] Re: [mg26587] Sum of Squares
  • From: BobHanlon at aol.com
  • Date: Thu, 11 Jan 2001 10:39:22 -0500 (EST)
  • Sender: owner-wri-mathgroup at wolfram.com

Here is an approach

var = {x, y, z};

expr1 = 5 x^2 + 8 x y + 5 y^2 + 2 x z - 2 y z + 2 z^2;

expr2 = Plus @@ ((Partition[coef = Table[Unique[s], {2*Length[var]}], 
                Length[var]].var)^2);

cl1 = Flatten[CoefficientList[expr1, var]];

cl2 = Flatten[CoefficientList[expr2, var]];

soln = Solve[Complement[Thread[cl1 == cl2],{True}], coef];

There are multiple solutions

ans = expr2 /. soln /. coef[[3]] -> 1

{(-2*x - y - z)^2 + (-x - 2*y + z)^2, 
  (-x - 2*y + z)^2 + (2*x + y + z)^2, 
  (x + 2*y - z)^2 + (2*x + y + z)^2, 
  (-x - 2*y + z)^2 + (2*x + y + z)^2}

And @@ ((expr1 == Expand[#])& /@ ans)

True


Bob Hanlon

In a message dated 2001/1/9 2:28:14 AM, wself at msubillings.edu writes:

>Is there, or could there be, a "sum of squares" function for
>multivariate polynomials, which would rewrite an expression as a sum of
>squares, when possible?
>
>It might work like this:
>
>SumOfSquares[5 x^2 + 8 x y + 5 y^2 + 2 x z - 2 y z + 2 z^2]
>
>--->
>
>(x + 2y - z)^2 + (2x + y + z)^2
>


  • Prev by Date: Re: Random Sampling Without Replacement?
  • Next by Date: Re: Expansion Coefficients in Multivariable Series Expansions?
  • Previous by thread: Re: Sum of Squares
  • Next by thread: Re: Sum of Squares