Reduction of a quadratic form : gauss

gauss takes two arguments : a symbolic expression representing a quadratic form q and a vector of variable names.
gauss returns q written as sum or difference of squares using Gauss algorithm.
Input :

gauss(2*x*y,[x,y])

Output :

(y+x)^2/2+(-(y-x)^2)/2