{VERSION 3 0 "DEC ALPHA UNIX" "3.0" }
{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 polyno
m with respect to \",vars,\" but dismantle received \",Phi, vars)\n  f
i;\n  Psi:=collect(Phi,vars);\n  if not degree(Psi,vars)<=2  then\n   \+
 ERROR(\"Phi is expected to have degree at most 2, but dismantle recei
ved \",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:
=collect(subs(seq(v=t*v,v=V),Psi),t);\n  for d from 0 to 2 do\n    Psi
.d:=coeff(Psi,t,d)\n  od;\n  A:=array(1..n,1..n,symmetric);\n  for i t
o 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$#\"
\"(\"\"),$*$-%%sqrtG6#\"\"$\"\"\"#!\"\"F*7$F+#\"\"&F*-F$6#7$7#\"\"#7#,
$F,F;F)7$%\"xG%\"yG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 71 "dism
antle(z-(7/8*x^2-sqrt(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%#!\"(\"\"),$*$-%%s
qrtG6#\"\"$\"\"\"#\"\"\"F*\"\"!7%F+#!\"&F*F47%F4F4F4-F$6#7%7#!\"#7#,$F
,F=7#F3F)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 "1 0 0" 0 }
{VIEWOPTS 1 1 0 1 1 1803 }