{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 }{PSTYLE "Maple Plot" 0 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 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 47 "A:=array(1..2,1..2,[[7/8,-sqrt(3)/4],[0,5/8]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"AG-%'MATRIXG6#7$7$#\"\"(\"\"),$*$\"\"$#\"\"\" \"\"##!\"\"\"\"%7$\"\"!#\"\"&F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "b:=array(1..2,1..1,[[2],[2*sqrt(3)]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bG-%'MATRIXG6#7$7#\"\"#7#,$*$\"\"$#\"\"\"F*F*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "c:=7;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cG\"\"(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "X:=array(1..2,1..1,[[x],[y]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%\"XG-%'MATRIXG6#7$7#%\"xG7#%\"yG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 90 "Phi:=collect(evalm(transpose(X) &* A &* X + transpose (b) &* X +c)[1,1],\{x,y\},distributed);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$PhiG,.*$%\"xG\"\"##\"\"(\"\")*(\"\"$#\"\"\"F(F'F/%\"yGF/#!\" \"\"\"%*$F0F(#\"\"&F+F'F(*&F-F.F0F/F(F*F/" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 21 "dismantle(Phi,[x,y]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6&-%'MATRIXG6#7$7$#\"\"(\"\"),$*$\"\"$#\"\"\"\"\"##!\"\"F*7$F+#\" \"&F*-F$6#7$7#F07#,$F,F0F)7$%\"xG%\"yG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "A1:=matrix([[7/8, -1/8*sqrt(3)], [-1/8*sqrt(3), 5/8]] );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#A1G-%'MATRIXG6#7$7$#\"\"(\"\" ),$*$\"\"$#\"\"\"\"\"##!\"\"F,7$F-#\"\"&F," }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 26 "evalm((transpose(A)+A)/2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7$7$#\"\"(\"\"),$*$\"\"$#\"\"\"\"\"##!\"\"F *7$F+#\"\"&F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "linalg[jor dan](A1,P);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7$7$\"\"\" \"\"!7$F)#F(\"\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "eval(P) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7$7$#\"\"$\"\"%#\"\" \"F*7$,$*$F)#F,\"\"##!\"\"F*,$F/F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "tPP:=evalm(transpose(P) &* P);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$tPPG-%'MATRIXG6#7$7$#\"\"$\"\"%\"\"!7$F-#\"\"\"F," } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "Q:=array(1..2,1..2,diagona l):\nQ[1,1]:=1/sqrt(tPP[1,1]):\nQ[2,2]:=1/sqrt(tPP[2,2]):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "P1:=evalm(P&*Q);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%#P1G-%'MATRIXG6#7$7$,$*$\"\"$#\"\"\"\"\"#F-F-7$#!\" \"F/F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 27 "evalm(transpose(P 1) &* P1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%'MATRIXG6#7$7$\"\"\"\" \"!7$F)F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "newX:=array(1. .2,1..1,[[newx],[newy]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%newXG- %'MATRIXG6#7$7#%%newxG7#%%newyG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "XX:=evalm(P1&*newX);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#XXG -%'MATRIXG6#7$7#,&*&\"\"$#\"\"\"\"\"#%%newxGF.F-%%newyGF-7#,&F0#!\"\"F /*&F,F-F1F.F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "tosubs1:= \{x=XX[1,1],y=XX[2,1]\}:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "newPhi:=collect(subs(tosubs1,Phi),\{newx,newy\},distributed);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%'newPhiG,**$%%newxG\"\"#\"\"\"\"\"(F )*$%%newyGF(#F)F(F,\"\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "\{diff(newPhi,newx),diff(newPhi,newy)\};" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$,$%%newxG\"\"#,&%%newyG\"\"\"\"\"%F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "center:=solve(\{diff(Phi,x),diff(Ph i,y)\},\{x,y\});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'centerG<$/%\"yG ,$*$\"\"$#\"\"\"\"\"#!\"#/%\"xGF." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "newcenter:=solve(\{diff(newPhi,newx),diff(newPhi,newy )\},\{newx,newy\});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*newcenterG<$ /%%newxG\"\"!/%%newyG!\"%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "subs(newcenter,tosubs1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$/%\" yG,$*$\"\"$#\"\"\"\"\"#!\"#/%\"xGF," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 84 "tosubs2:=subs(newxx=newx,newyy=newy,subs(newcenter,\{ newxx=xi+newx,newyy=eta+newy\}));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %(tosubs2G<$/%%newxG%#xiG/%%newyG,&%$etaG\"\"\"!\"%F-" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "newnewPhi:=collect(subs(tosubs2,new Phi),\{xi,eta\},distributed);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*ne wnewPhiG,(*$%$etaG\"\"##\"\"\"F(*$%#xiGF(F*!\"\"F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "tosubs:=subs(tosubs2,tosubs1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'tosubsG<$/%\"xG,(*&\"\"$#\"\"\"\"\"#%#xiGF,F +%$etaGF+!\"#F,/%\"yG,&F.#!\"\"F-*&F*F+,&F/F,!\"%F,F,F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 69 "arc:=subs(subs(xi=cos(t),eta=sqrt(2 )*sin(t),tosubs),[x,y,t=0..2*Pi]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%$arcG7%,(*&\"\"$#\"\"\"\"\"#-%$cosG6#%\"tGF*F)*&F+F)-%$sinGF.F*F)!\" #F*,&F,#!\"\"F+*&F(F),&F0F*!\"%F*F*F)/F/;\"\"!,$%#PiGF+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "plot(arc,scaling=constrained,labels =[x,y]);" }}{PARA 13 "" 1 "" {INLPLOT "6%-%'CURVESG6$7S7$$!3Shb@'fuR8 \"!#<$!3Rax8:;5kRF*7$$!3q!G'R1TaX5F*$!3:wZi))z?#z$F*7$$!35x8?He$4$)*!# =$!3!zAxRXGvj$F*7$$!3*ee-*fSG,$*F5$!3*oQ5D&pugMF*7$$!3KRX7FR!>'*)F5$!3 #pW#>x9'GG$F*7$$!3%*RIgvPFB))F5$!3!\\Y=MD'34JF*7$$!3L()Q*)H.Yu))F5$!3 \"p$4%=X;O&HF*7$$!3V#RNp#f'*)4!GF*7$$!3%)[zQsRAT&*F5$! 31$\\=*zFuaEF*7$$!3ztg+Q))o:5F*$!3Cz7/,EDBDF*7$$!3z+l'e$z5(4\"F*$!3n\"F*$!3+Ta9*GE79#F*7$$!3,,#[X;Ug#=F*$!3*G.O5vMV:#F*7 $$!3&p]k/wa%f>F*$!3hZvq.\"QF* 7$$!3?-+cXW)H6$F*$!3'**Gja.M(pRF*7$$!3>NB2:`\"p3$F*$!3C/UMyo(p8%F*7$$! 3?rT,.%e\\/$F*$!3E*[.t(\\)eF%F*7$$!38UNVYSe\")HF*$!3=xb6f()p4WF*7$$!3 \"R*p,$=S&3HF*$!3B=Eeci+;XF*7$$!32IG_fr.8GF*$!3kOX#)G*4Sh%F*7$$!399W-[ 7v4FF*$!3:BFfq?J(o%F*7$$!3WNY_N%*p*e#F*$!31cvQ\\c.VZF*7$$!3U]\\R&RjBY# F*$!3m_Emwn%fx%F*7$$!3Nmw;'oG3K#F*$!37*H2XF *7$$!3:2O.yC.u9F*$!38HSgM$*4)R%F*7$$!3S)**)GLP4[8F*$!3kDjE8Y!fE%F*7$$! 3BkSG/uEQ7F*$!3sY&pZ;Lh7%F*7$$!3BTaj&fuR8\"F*$!3\"o&H89;5kRF*-%'COLOUR G6&%$RGBG$\"#5!\"\"\"\"!F_[l-%(SCALINGG6#%,CONSTRAINEDG-%+AXESLABELSG6 $%\"xG%\"yG" 2 258 258 258 2 0 1 0 2 9 0 4 1 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 -24000 0 0 0 0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "7 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 }