{VERSION 2 3 "DEC ALPHA UNIX" "2.3" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 
0 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }
{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 
0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Output" 0 11 1 
{CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 
0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 852 "dismantle:=proc(Phi
,vars::\{set,list\})\n  local Psi,n,V,d,A,i,j,b;\n  if not type(Phi,po
lynom(anything,vars)) then\n    ERROR(`Phi is expected to be a polynom
 with respect to `,vars,` but dismantle received `,Phi, vars)\n  fi;\n
  Psi:=collect(Phi,vars);\n  if not degree(Psi,vars)<=2  then\n    ERR
OR(`Phi is expected to have degree at most 2, but dismantle received `
,Phi, vars)\n  fi;\n  n:=nops(vars);\n  if type(vars,set) then \n    V
:=[op(vars)]\n  else \n    V:=[op(\{op(vars)\})]\n  fi;\n  Psi:=collec
t(subs(seq(v=t*v,v=V),Psi),t);\n  for d from 0 to 2 do\n    Psi.d:=coe
ff(Psi,t,d)\n  od;\n  A:=array(1..n,1..n,symmetric);\n  for i to n do
\n    A[i,i]:=coeff(Psi2,V[i],2);\n    for j from i+1 to n do \n      \+
A[i,j]:=coeff(coeff(Psi2,V[i],1),V[j],1)/2\n    od\n  od;\n  b:=array(
1..n,1..1);\n  for i to n do\n    b[i,1]:=coeff(Psi1,V[i],1)\n  od;\n \+
 eval(A),eval(b),Psi0,V\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 65 "dismantle(7/8*x^2-sqrt(3)/4*x*y+5/8*y^2+2*x+2*sqrt(3)*y+7,\{x,y
\});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6&-%'MATRIXG6#7$7$#\"\"(\"\"),$*
$\"\"$#\"\"\"\"\"##!\"\"F*7$F+#\"\"&F*-F$6#7$7#F07#,$F,F0F)7$%\"xG%\"y
G" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 71 "dismantle(z-(7/8*x^2-s
qrt(3)/4*x*y+5/8*y^2+2*x+2*sqrt(3)*y+7),\{x,y,z\});" }}{PARA 11 "" 1 "
" {XPPMATH 20 "6&-%'MATRIXG6#7%7%#!\"(\"\"),$*$\"\"$#\"\"\"\"\"##F/F*
\"\"!7%F+#!\"&F*F27%F2F2F2-F$6#7%7#!\"#7#,$F,F;7#F/F)7%%\"xG%\"yG%\"zG
" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "dismantle(x^2-2*x*y*cos
(theta)+y^2+x*cos(theta)+y*sin(theta),[x,y]);" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6&-%'MATRIXG6#7$7$\"\"\",$-%$cosG6#%&thetaG!\"\"7$F)F(-F$
6#7$7#F*7#-%$sinGF,\"\"!7$%\"xG%\"yG" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 0 "" }}}}{MARK "4 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }