{VERSION 5 0 "SUN SPARC SOLARIS" "5.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {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 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 2 6 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 2 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "W arning" 2 7 1 {CSTYLE "" -1 -1 "" 0 1 0 0 255 1 0 0 0 0 0 0 1 0 0 0 } 0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Error" 7 8 1 {CSTYLE "" -1 -1 "" 0 1 255 0 255 1 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 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 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 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 43 "Fitting a circle to data w ith Least-Squares" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "read( \"/home/mat331/Worksheets/lsq_data.txt\");" }}{PARA 6 "" 1 "" {TEXT -1 82 "defined line_pts(), bad_line_pts(), quadratic_pts(), cubic_pts( ), and circle_pts()" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "data := circle_pts();" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%dataG777$$\"*yK #[$)!\"*$\"+E$e>N)F)7$$\"+nE0&Q)F)$\"+M!e>n%F)7$$\"+E!*Qj()F)$\"+n:WmC F)7$$\"+6%*4NmF)$\"+iWWOsF)7$$\"+H(=`$>F)$\"+Vps%[)F)7$$!+9!y3=#F)$!+E x(yG\"F)7$$\"+jL84bF)$!+dfu)z#F)7$$\"+4Sm7lF)$!+k)QUP#F)7$$\"+dJU&G)F) $\"*WR#QaF)7$$\"+#pIgX)F)$\"+7t!fA$F)7$$\"+tpcF')F)$\"+*p**3@#F)7$$!+G *3x>#F)$!+=cn)G\"F)7$$\"+2=0VvF)$\"+PVl,jF)7$$!+%RclS$F)$\"+L:n\\LF)7$ $\"*$RB-')F)$!+3qv0LF)7$$\"+(3;#ouF)$\"+kP\")>jF)7$$\"+HOP*\\'F)$\"+X& ***frF)7$$!+7;nmLF)$\"*m\"fLdF)7$$!*Zws;'F)$\"+)>/+`(F)7$$\"+#z(Ry[F)$ \"+1u0%Q)F)7$$!+b#*)y&GF)$!*7Kju'F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "nops(data);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"#@" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 44 "plot(data,style=point, sc aling=constrained);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6'-%'CURVESG6$777$$\"3-+++!yK#[$)!#=$\"34+++E$e>N)!#<7$$ \"31+++nE0&Q)F-$\"3>+++M!e>n%F-7$$\"3Z+++E!*Qj()F-$\"3&******pcTkY#F-7 $$\"3n*****4T*4NmF-$\"3p*****>YWkB(F-7$$\"3,+++H(=`$>F-$\"33+++Vps%[)F -7$$!38+++9!y3=#F-$!30+++Ex(yG\"F-7$$\"37+++jL84bF-$!3%******p&fu)z#F- 7$$\"3')******3Sm7lF-$!30+++k)QUP#F-7$$\"3J+++dJU&G)F-$\"3/+++S%R#QaF* 7$$\"3c+++#pIgX)F-$\"3#******>J2fA$F-7$$\"3U+++tpcF')F-$\"30+++*p**3@# F-7$$!3%******z#*3x>#F-$!32+++=cn)G\"F-7$$\"3K+++2=0VvF-$\"3W+++PVl,jF -7$$!3))*****RRclS$F-$\"3(******H`r'\\LF-7$$\"3e+++IRB-')F*$!30+++3qv0 LF-7$$\"31+++(3;#ouF-$\"31+++kP\")>jF-7$$\"3#*******GOP*\\'F-$\"3!)*** **\\a***frF-7$$!3>+++7;nmLF-$\"3J+++g;fLdF*7$$!3X+++qkFnhF*$\"3K+++)>/ +`(F-7$$\"3M+++#z(Ry[F-$\"3:+++1u0%Q)F-7$$!3!******\\D*)y&GF-$!3P+++?@ LYnF*-%'COLOURG6&%$RGBG$\"#5!\"\"$\"\"!FdrFcr-%&STYLEG6#%&POINTG-%(SCA LINGG6#%,CONSTRAINEDG-%+AXESLABELSG6$Q!6\"F`s-%%VIEWG6$%(DEFAULTGFes" 1 5 0 1 10 0 2 9 1 4 1 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "cdist:=(a,b,r,p) -> ( (p[1] \+ -a)^2 + (p[2]-b)^2 - r^2 )^2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%&cd istGf*6&%\"aG%\"bG%\"rG%\"pG6\"6$%)operatorG%&arrowGF+*$),(*$),&&9'6# \"\"\"F89$!\"\"\"\"#F8F8*$),&&F66#F;F89%F:F;F8F8*$)9&F;F8F:F;F8F+F+F+ " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "E:= (a,b,r,data) -> sum ( cdist(a,b,r,data[i]), i=1...nops(data) );" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"EGf*6&%\"aG%\"bG%\"rG%%dataG6\"6$%)operatorG%&arrow GF+-%$sumG6$-%&cdistG6&9$9%9&&9'6#%\"iG/F;;\"\"\"-%%nopsG6#F9F+F+F+" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 122 "soln:=solve(\{diff(E(a,b, r,data), a)=0,\n diff(E(a,b,r,data), b)=0,\n diff(E(a,b,r, data), r)=0\},\n \{a,b,r\});" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 1 " " }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%solnG6)<%/%\"rG$\"\"!F*/%\" bG^$$\"+#ol0." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "soln[7];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<%/% \"aG$\"+.(Qfe#!\"*/%\"bG$\"+k=W0DF(/%\"rG$\"+j:#*GhF(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 125 "plots[display](\n plot( subs(soln[ 7], [a+ r*cos(t), b+ r*sin(t), t=0...2*Pi]),\n plot(data,style=POINT, \+ scaling=constrained) );" }}{PARA 8 "" 1 "" {TEXT -1 22 "Error, `;` une xpected\n" }}}}{MARK "10 1 0" 22 }{VIEWOPTS 1 1 0 2 1 1805 1 1 1 1 } {PAGENUMBERS 0 1 2 33 1 1 }