{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 } {CSTYLE "" -1 256 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 257 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 0 0 0 1 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 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 "Warni ng" 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 "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 }{PSTYLE "" 0 256 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 "" {MPLTEXT 1 0 103 "Seed:=readlib(rando mize)():\nwith(stats):\nmyrand:=(a,b)->evalf(a+(b-a)*rand()/10^12):\np oints_wanted:=30:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 258 "line_ pts:= proc()\n local fuzz, m, b, x,i;\n m:= myrand(-5,5);\n b:= my rand(-20,20);\n fuzz := [ random[normald[0,0.5]](points_wanted)];\n \+ x := [ random[uniform[-10,10]](points_wanted)];\n RETURN( [ [x[ i], m*x[i]+b+fuzz[i]] $i=1..points_wanted]);\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "lineardata:=line_pts();" }}{PARA 12 "" 1 " " {XPPMATH 20 "6#>%+lineardataG7@7$$!+#eid\"Q!\"*$\"+[?s))>!#57$$\"*&o [^O!\")$!+68;LAF07$$\"*At%)4\"F0$!+%*4rN9F07$$\"*Ja\\)>F0$!+C&[hz\"F07 $$\"*-3qx(F0$!+ibm*\\$F07$$!+yXDr'*F)$\"+,Yc\\ 9F07$$\"*ug+K$F0$!+&y7KC#F07$$!++c*yy'F)$\"+EyP\\**F)7$$!+-#)=1tF)$\"+ A;2t5F07$$!+UF6\"o*F)$\"+A=8+=F07$$!+E5n*z&F)$\"+&z$o.hF)7$$\"*K)ek\") F0$!+RT9YOF07$$\"*-\"GW!)F0$!+:C)fn$F07$$\"(R)[jF0$!+W]$)=6F07$$!+3c:: HF)$!+YgR5LF)7$$!+ccx!Q'F)$\"+?')R%3(F)7$$!+GF')\\7F)$!+HR6D$)F)7$$\"* Z\"eHfF0$!+?))[QIF07$$!+1LZhAF)$!+hLl:YF)7$$\"*B`^:\"F0$!+Z:;h:F07$$\" *9%3v%*F0$!+DX!=6%F07$$!+K;-$o*F)$\"+Ab\"Gw\"F07$$\"*Mq0C#F0$!+9CD&*=F 07$$!+e;#Q2)F)$\"+-0GB8F07$$\"*$)o*y=F0$!+X+$3u\"F07$$\"**[prPF0$!++&* e1CF07$$!+/$*>x^F)$\"+?[A`TF)7$$\"*2DB[*F0$!+#zQ!)3%F07$$!*1F3b)F)$!+s 8=V\"*F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "with(plots):" } }{PARA 7 "" 1 "" {TEXT -1 50 "Warning, the name changecoords has been \+ redefined\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "with(CurveFi tting);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7*%(BSplineG%-BSplineCurveG %,InteractiveG%-LeastSquaresG%8PolynomialInterpolationG%6RationalInter polationG%'SplineG%4ThieleInterpolationG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "polyfit:=PolynomialInterpolation(lineardata,x);" }} {PARA 12 "" 1 "" {XPPMATH 20 "6#>%(polyfitG,hn*&$\"+2/)zD\"!\"(\"\"\"% \"xGF*!\"\"*&$\"*jxGN(!#6F*)F+\"#:F*F**&$\"+HhL1BF)F*)F+\"\"'F*F**&$\" +$\\TG'))!\")F*)F+\"\"#F*F**&$\"+w!=V5$F)F*)F+\"\"$F*F**&$\"+*Q9dj#F)F *)F+\"\"%F*F,*&$\"+l`Y`CF)F*)F+\"\"&F*F,*&$\"+!)zL^zF;F*)F+\"\"(F*F**& $\"+lR%yu)F;F*)F+\"\")F*F,*&$\"+^fK=8F;F*)F+\"\"*F*F,*&$\"+v!Q$p\")!#7F* )F+\"#;F*F,*&$\"+!ze;-%!#8F*)F+\"#F*F**&$\"+0'4,&e!#:F*)F+\"#?F*F,*&$\"+z8n^P!#;F*) F+\"#@F*F,*&$\"+h,(oH)!#F*)F+\"#CF*F,*&$\"+[ZBhb!#?F*)F+\"#DF*F,*&$\"+.ZLaO!#@F*)F+ \"#EF*F**&$\"+A0^HH!#AF*)F+\"#FF*F**&$\"+qR1Ty!#CF*)F+\"#GF*F,*&$\"+3p Ebl!#DF*)F+\"#HF*F,$\"+([)zLOF^oF," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 126 "display(plot(lineardata, style=point, color=black, s ymbol=BOX,axes=boxed), plot(polyfit,x=-10..10, y=-200..200, thickness= 2));" }}{PARA 13 "" 1 "" {GLPLOT2D 427 427 427 {PLOTDATA 2 "6'-%'CURVE SG6&7@7$$!32+++#eid\"Q!#<$\"31+++[?s))>!#=7$$\"3&*******\\o[^OF*$!3#** ****4JhJB#!#;7$$\"3-+++?KZ)4\"F*$!33+++%*4rN9F37$$\"3++++5V&\\)>F*$!3= +++C&[hz\"F37$$\"3C+++?!3qx(F*$!3')*****>cl'*\\$F37$$!3?*****zda7n*F*$ \"3:+++,Yc\\++++<^#3)F-$!3/+++XE!)>9F37$$\"3=+++S21?LF*$!31+ ++&y7KC#F37$$!3p*******f&*yy'F*$\"3?*****f#yP\\**F*7$$!3T+++-#)=1tF*$ \"3(******>irI2\"F37$$!3#)*****>u76o*F*$\"3<+++A=8+=F37$$!3S+++E5n*z&F *$\"3C+++&z$o.hF*7$$\"3++++?$)ek\")F*$!3u******QT9YOF37$$\"3n******>5G W!)F*$!3++++:C)fn$F37$$\"3d++++!R)[j!#>$!32+++W]$)=6F37$$!31+++3c::HF* $!3:+++YgR5LF*7$$!3&******flv2Q'F*$\"3(*******>')R%3(F*7$$!34+++GF')\\ 7F*$!3Y+++HR6D$)F*7$$\"3Y+++q9eHfF*$!3*)******>))[QIF37$$!32+++1LZhAF* $!3o*****4O`ch%F*7$$\"3.+++IK:b6F*$!3.+++Z:;h:F37$$\"3))******RT3v%*F* $!3;+++DX!=6%F37$$!3(3++?j@Io*F*$\"3++++Ab\"Gw\"F37$$\"35+++S.dSAF*$!3 9+++9CD&*=F37$$!3&******zl@Q2)F*$\"32+++-0GB8F37$$\"3.+++I)o*y=F*$!3() *****\\/I3u\"F37$$\"3-+++!*[prPF*$!37++++&*e1CF37$$!3;+++/$*>x^F*$\"3: +++?[A`TF*7$$\"3\"*******p]K#[*F*$!3p*****>zQ!)3%F37$$!3D+++gq#3b)F-$! 3!3++?P\"=V\"*F*-%'COLOURG6&%$RGBG\"\"!F_uF_u-%&STYLEG6#%&POINTG-%'SYM BOLG6#%$BOXG-F$6%7bo7$$!#5F_u$\"3++++X)eLM\"!\"*7$$!3!pmmm\"p0k&*F*$!3 ++voa?fB5!#67$$!3uKL$3s%HaF*$\"34U$) Gg>3\"3$F\\y7$$!3]******\\$*4)*\\F*$!3&R4(z#QyY]$F\\y7$$!3o******\\_& \\c%F*$!3)p+s[K4L3)F\\y7$$!3$)******\\1aZTF*$!3SsCfTd*)fSF\\y7$$!3Imm; /#)[oPF*$\"3s;&4VKtr!QFgy7$$!3%HLLL=exJ$F*$\"3q'pd$Q7Q\\6F\\y7$$!3lKLL L2$f$HF*$\"3UC=h;wCUUF37$$!3%)****\\PYx\"\\#F*$!3D'GW#4SS\"F37$$\"3% \\LL$e\"*[H7F*$!3?D!4xJxmt\"F37$$\"3=++++dxd;F*$!3]nsz***\\F#>F37$$\"3 e+++D0xw?F*$!3]$>0;40h*=F37$$\"34,+]i&p@[#F*$!3%3(G=#*o(42\"F*7$$\"3++ ++vgHKHF*$\"3[l`Ld#QX9(F37$$\"3ElmmmZvOLF*$!3*>_%)*\\=YEGF37$$\"3%4+++ v+'oPF*$!3[N@mhw'QA#F37$$\"3UKL$eR<*fTF*$!3[qd-(3:)p[F\\y7$$\"3K-++])H xe%F*$!3^8=jW]]^gFax7$$\"3!fmm\"H!o-*\\F*$!3xR_;(\\#G*e#F[x7$$\"3X,+]7 k.6aF*$!3kfQ9P1^veF[x7$$\"3#emmmT9C#eF*$!3%G9$4&)*ebQ$F[x7$$\"32****\\ i!*3`iF*$\"3T02&4TE<7#Ffv7$$\"3;NLLL*zym'F*$\"3Q\"Qv%*y?+#oFfv7$$\"3'e LL$3N1#4(F*$\"3c/\"os%oaF#*Ffv7$$\"3,pm;HYt7vF*$\"3)*y!4]gwA!RFfv7$$\" 37-+++xG**yF*$!3wM\")3f;)3t%F[x7$$\"3gpmmT6KU$)F*$!3.2(*>UCG:hFfv7$$\" 3qNLLLbdQ()F*$!34YAbetJL5F`v7$$\"33om\"zW?)\\*)F*$!3W[@qpASW?F`v7$$\"3 Z++]i`1h\"*F*$!3ilZ8`7>#Q#F`v7$$\"3\"4+]7`f@E*F*$!3voaN*R@O&=F`v7$$\"3 M,+++PDj$*F*$!3+v=nj!H%4')F]v7$$\"3X-]P%y+QT*F*$!3+v=#\\>-RW$F]v7$$\"3 y,+voyMk%*F*$!3,vVBS_1$)>Ffv7$$\"36,]7`\\*[^*F*$!3+]7.#[=mo\"F]v7$$\"3 @-+]P?Wl&*F*$!3+v$f$o8$QA\"F`v7$$\"3],D\"Gyh(>'*F*$!3]7G8E.-1TF`v7$$\" 3w+]7G:3u'*F*$!3+]7`R9*4%**F`v7$$\"3#=]PMF,%G(*F*$!3+DJ&pgB-/#!\")7$$ \"35,+v=5s#y*F*$!3](oaI&Q:&y$F_hl7$$\"3k^iS\"*3))4)*F*$!3]7`Wugw<]F_hl 7$$\"3Q+D1k2/P)*F*$!3]Pf3x2iclF_hl7$$\"3\"4v=nj+U')*F*$!3]PM-f)\\]7.#Q?&=**F*$!3voagRh6m8Fiil7$$\"3cu$f$=.5K**F*$!3](=UL -&GJ:Fiil7$$\"3s+voa-oX**F*$!3voa&R\"o*Hr\"Fiil7$$\"3(oi:5>g#f**F*$!3D \"y+&3uk7>Fiil7$$\"3D^PMF,%G(**F*$!3Dc^e*QM<8#Fiil7$$\"3iv=nj+U')**F*$ !3vV)R+aj=P#Fiil7$$\"#5F_u$!3++D\"*HdvMEFiil-F\\u6&F^u$F_\\m!\"\"$F_uF _uFf\\m-%*THICKNESSG6#\"\"#-%*AXESSTYLEGFfu-%+AXESLABELSG6%Q!6\"F`]m-% %FONTG6#%(DEFAULTG-%%VIEWG6$;F\\vF^\\m;$!$+#F_u$\"$+#F_u" 1 2 0 1 10 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "easydata:=[[1,1],[2,2.5],[3, 3]];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)easydataG7%7$\"\"\"F'7$\"\" #$\"#D!\"\"7$\"\"$F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "pol yfit2:=PolynomialInterpolation(easydata,x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)polyfit2G,(*&$\"+++++]!#5\"\"\")%\"xG\"\"#F*!\"\"*&$ \"+++++I!\"*F*F,F*F*$\"+++++:F2F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 121 "display(plot(easydata, style=point, color=black, sym bol=BOX,axes=boxed), plot(polyfit2,x=-1..4, y=-20..20, thickness=1)); " }}{PARA 13 "" 1 "" {GLPLOT2D 426 426 426 {PLOTDATA 2 "6'-%'CURVESG6& 7%7$$\"\"\"\"\"!F(7$$\"\"#F*$\"3++++++++D!#<7$$\"\"$F*F2-%'COLOURG6&%$ RGBGF*F*F*-%&STYLEG6#%&POINTG-%'SYMBOLG6#%$BOXG-F$6%7S7$$!\"\"F*$!\"&F *7$$!3Fnmm\"HU,\"*)!#=$!3SN()yle**pXF07$$!3)HL$3FH'='zFK$!3H?LZ%=:b?%F 07$$!35mm;/OU&*oFK$!3u$\\aYThj!QF07$$!3Plm;H_\">#eFK$!358!GM0[gT$F07$$ !3SLL3_!4Nv%FK$!3B\">mr&>.RIF07$$!3YmmTg(fHw$FK$!3Kvv\"*eso*p#F07$$!3] ++vVLIPFFK$!3'\\_.6:b'eBF07$$!3\"em;/,oln\"FK$!3MiFf/[-'!#>$!3BS9E(f(o(o\"F07$$\"3B%ommTIOo%F^p$!3p&*R>2zeg8F07$$\"3eM L$3_>jU\"FK$!3(3bF13wA3\"F07$$\"3E,++D;v/DFK$!31/q9gTV*z(FK7$$\"3x+++v =h(e$FK$!3=[?GBBr!)[FK7$$\"3U+++v$[6j%FK$!37bKcQ;$*y@FK7$$\"3DMLe*[z(y bFK$\"3;v@!3R%*>!=F^p7$$\"3knmmTXg0nFK$\"3-o*Q:,d&oGFK7$$\"3OommmJ+](=xpe =\"F0$\"3'zc,xvlWN\"F07$$\"37nm\"H28IH\"F0$\"3$**H$Q:y4V:F07$$\"3um;zp SS\"R\"F0$\"3D$femc4iq\"F07$$\"3GLL3_?`(\\\"F0$\"3'4aTF.&Hr=F07$$\"3fL $e*)>pxg\"F0$\"3iRaJ)pY3.#F07$$\"33+]Pf4t.))G#F07$$\"3/+++DRW9>F0$\"3-(H:S+%y5CF07$$\"39++ DJE>>?F0$\"3o]Gq_%3!>DF07$$\"3F+]i!RU07#F0$\"3T0ilcrF8EF07$$\"3+++v=S2 LAF0$\"3Fi=kpA\"fq#F07$$\"3Jmmm\"p)=MBF0$\"3$ysZ3lZ$yFF07$$\"3B++](=]@ W#F0$\"3%z&oM$z,W%GF07$$\"35L$e*[$z*RDF0$\"3#)*f2.]!>%*GF07$$\"3e++]iC $pk#F0$\"31V=)plrw$HF07$$\"3Zm;H2qcZFF0$\"3Gdh4#zQ\"oHF07$$\"3O+]7.\"f F&GF0$\"3&GsV\"f+;*)HF07$$\"3Xmm;/OgbHF0$\"3y/&=+[9!**HF07$$\"3w**\\il AFjIF0$\"3+T8?5$)*z*HF07$$\"3yLLL$)*pp;$F0$\"3KKGL701')HF07$$\"3(RL$3x e,tKF0$\"3)z?Ia;JF'HF07$$\"3Cn;HdO=yLF0$\"3_[(y1c)[GHF07$$\"3`+++D>#[Z $F0$\"3AZ'p(p?F()GF07$$\"3RnmT&G!e&e$F0$\"3amcmkyaGGF07$$\"3\"RLLL)Qk% o$F0$\"3wEr]w8jlFF07$$\"37+]iSjE!z$F0$\"3C'ePabRxo#F07$$\"3a+]P40O\"*Q F0$\"3&*oL;@#QFg#F07$$\"\"%F*F.-F56&F7$\"#5FE$F*F*Ff\\l-%*THICKNESSG6# F)-%*AXESSTYLEGF>-%+AXESLABELSG6%Q!6\"F_]l-%%FONTG6#%(DEFAULTG-%%VIEWG 6$;FDF`\\l;$!#?F*$\"#?F*" 1 2 0 1 10 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}{EXCHG {PARA 256 "" 0 "" {TEXT 256 7 "SPLINES" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "data:=[[0,1],[2,4],[5,-9],[10,1],[11,6]];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%dataG7'7$\"\"!\"\"\"7$\"\"#\"\"%7$\"\"&! \"*7$\"#5F(7$\"#6\"\"'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 58 "s pline1:=Spline(data,x,degree=1);\nspline2:=Spline(data,x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(spline1G-%*PIECEWISEG6&7$,&\"\"\"F**(\"\" $F*\"\"#!\"\"%\"xGF*F*2F/F-7$,&#\"#QF,F**(\"#8F*F,F.F/F*F.2F/\"\"&7$,& \"#>F.*&F-F*F/F*F*2F/\"#57$,&\"#\\F.*&F8F*F/F*F*%*otherwiseG" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(spline2G-%*PIECEWISEG6&7$,(\"\"\"F**&#\"% ')p\"%VBF*%\"xGF*F**&#\"%Vp\"&W(=F**$)F/\"\"$F*F*!\"\"2F/\"\"#7$,*#\"& Ev$\"%HqF7*&#\"&4&>\"%i:F*F/F*F**&#\"&bX%\"%s$*F**$)F/F9F*F*F7*&#\"&j= \"\"&;\"GF*F4F*F*2F/\"\"&7$,*#\"&hm%\"$\"yF**&#\"%`<\"#mF*F/F*F7*&#\"% Z&*\"%CJF*FGF*F**&#\"%FY\"&go%F*F4F*F72F/\"#57$,*#\"$\\&\"#rF7*&#\"&.' H\"%'o%F*F/F*F7*&#\"$$H\"$%GF*FGF*F**&#FgoFFF*F4F*F7%*otherwiseG" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 115 "display(plot(data,style=poi nt,symbol=BOX,axes=boxed,color=black),plot([spline1,spline2],x=0..11,c olor=[red,blue]));" }}{PARA 13 "" 1 "" {GLPLOT2D 426 426 426 {PLOTDATA 2 "6(-%'CURVESG6&7'7$$\"\"!F)$\"\"\"F)7$$\"\"#F)$\"\"%F)7$$ \"\"&F)$!\"*F)7$$\"#5F)F*7$$\"#6F)$\"\"'F)-%'COLOURG6&%$RGBGF)F)F)-%&S TYLEG6#%&POINTG-%'SYMBOLG6#%$BOXG-F$6$7fnF'7$$\"3ALLLepo(R#!#=$\"3%)** *\\P/`'f8!#<7$$\"3`mmTg:!R[%FP$\"3.+D1M_es;FS7$$\"3yKL$32o+$oFP$\"3#** *\\i5-^C?FS7$$\"3'RLLe\\'y\">*FP$\"34+]PuzwyBFS7$$\"3^m;a3!GU:\"FS$\"3 m*\\7G,U8t#FS7$$\"3>L$3F&)[@P\"FS$\"3o*\\i!zKAeIFS7$$\"3!***\\PkKz(f\" FS$\"33+Dc'*)*o'R$FS7$$\"3cm;/^TZ9FS$\"3$)*\\P4f#>@RFS7$$\"3*)*\\(o*RQl(>FS$\"3&*\\7`*f2['RFS7$$\"3 BLL3s]h0?FS$\"3%>As4-oc(RFS7$$\"3ym\"zWu\"pM?FS$\"3LVpD2\"p'\\QFS7$$\" 3!***\\(oToP1#FS$\"3]m;a$>qOs$FS7$$\"3HL3-VTS$=#FS$\"3I)QwNr[_?$FS7$$ \"3qmm;p)RII#FS$\"3&=66OBFoo#FS7$$\"3ZLLe%H!z8DFS$\"3K*))Q,RvNx\"FS7$$ \"3B++]d`/^FFS$\"3H+++v6qauFP7$$\"3f****\\7YF*)HFS$!3rTmmTlcoGFP7$$\"3 #*****\\UE&)=KFS$!3rlm;%y%p\"G\"FS7$$\"3'HL3x[JtU$FS$!3&\\x-M6.^=#FS7$ $\"3cmm;**HBvOFS$!3-))))Q'*HMfKFS7$$\"3Tmmm'4Q_)QFS$!3%Gbbb3l$pTFS7$$ \"3y**\\P\\R_HTFS$!3\\J$eR6PzA&FS7$$\"3vlmm@$edM%FS$!3?;AAFF&\\;'FS7$$ \"3,+]P*p,Ie%FS$!3iJ$eRO2I>(FS7$$\"3N+]7)\\8*3[FS$!3A,](=\\e><)FS7$$\" 3gL3FH6xE\\FS$!3aVp]$*[n#o)FS7$$\"3'om;/wGY/&FS$!3Gmm;zCu5*)FS7$$\"3!p ;zp&)eG:&FS$!3=m;/'G#G%p)FS7$$\"3%pmTN&*)3h_FS$!36mm\"H4AyZ)FS7$$\"3yK Le90d%\\&FS$!3UML$3(*e3,)FS7$$\"3mK$3xB#4PdFS$!3pMLeCb\"e_(FS7$$\"3)** *\\i5\"3#[fFS$!3-++vyPe.rFS7$$\"3TLL3P!>i<'FS$!3;LL$e#>cZmFS7$$\"3&*)* ***\\jwloFS$!3G,+D\"[8'p_FS7$$\"3Z***\\7%Gw7rFS$!31,+]R*)FS$!3!G+]7.-;7\"FS7$$\"3;MLLjR Ln\"*FS$!3k;LLL2K`mFP7$$\"38LLeH\\j+%*FS$!3MPLL39I()>FP7$$\"3rm;/YS+K' *FS$\"3BML$3#43SEFP7$$\"3\"3++]B3Y%)*FS$\"3M;+++Z;#*oFP7$$\"3nm;ziw#)3 5!#;$\"3'pL$eRJQT9FS7$$\"3ILLLa;iI5F`^l$\"3&pmmmr#3JDFS7$$\"3%**\\P\\f eQ0\"F`^l$\"33&*\\(ouHHp$FS7$$\"32+D17$*4w5F`^l$\"3?0]7.c'\\![FSF9-F?6 &FA$\"*++++\"!\")F(F(-F$6$7ZF'7$FN$\"3?K8GS**z4(QFS7$$\" 3')**\\iI%)=j7FS$\"35*oFLI!y>SFS7$Fco$\"3)pUeIj:V8%FS7$$\"3Wm;ae5(\\[ \"FS$\"3cc^nq9s9UFS7$Fho$\"3$*oj^:#=JD%FS7$F]p$\"3ckH**)zV_C%FS7$Fbp$ \"3A:j]xL\\&=%FS7$F\\q$\"3e+M!R3)zqSFS7$F`r$\"3%ot4%G*fx*QFS7$Fjr$\"3- &yGm)4r=FFS7$Fds$\"3[:h#=FS7$Fis$\"3\")G[g\" =f%eyFP7$F^t$!3Z&Gnd4eC@$FP7$Fct$!36nM$)3VZ*Q\"FS7$Fht$!3Bp#*[)>e[q#FS 7$F]u$!3'3HUsq/0$QFS7$Fbu$!3ZlQwE26?^FS7$Fgu$!3g&)R(RY,d@'FS7$F\\v$!3q e'QA]VkL(FS7$Fav$!3QX(Rp&yf%H)FS7$F[w$!3zk;#\\Sx)[\"*FS7$Few$!3e#o@yFo Oy*FS7$Fjw$!3ms<^')G7J5F`^l7$F_x$!3M%Q**oGU%p5F`^l7$Fdx$!3k'>vi%zx*3\" F`^l7$$\"3qmT&Qd8A1'FS$!3'4:-TcFf4\"F`^l7$Fix$!3SPlq?7!))4\"F`^l7$$\"3 =m;/Oy*RH'FS$!3O&=>A!GU)4\"F`^l7$F^y$!3]+y#oYQZ4\"F`^l7$$\"3d)*\\(=@+q _'FS$!3Y*yk:%z-)3\"F`^l7$Fcy$!3sh?zdmLy5F`^l7$Fhy$!3,W\"*G$ob90\"F`^l7 $F]z$!3/'>=oGF(45F`^l7$Fbz$!3@cK%zaMGi*FS7$Fgz$!3Q+,4A#[*>!*FS7$F\\[l$ !3C5Fbr#HUR)FS7$Fa[l$!3K'=hhaf4j(FS7$Ff[l$!3K\\_%3eaR%oFS7$F[\\l$!3aAn \\!Q$4dfFS7$F`\\l$!3OI;1pF^L]FS7$Fe\\l$!3/?>!H4mU,%FS7$Fj\\l$!3&)4DvwS B*)HFS7$F_]l$!3pDZ1hcR/>FS7$Fd]l$!37=RUIP[$*zFP7$Fi]l$\"31QOJ%4DRN#FP7 $F^^l$\"30JvCsGdO9FS7$Fd^l$\"3t*3q.*G$)>DFS7$Fi^l$\"3bD?U,cd\"o$FS7$F^ _l$\"3AU^x\"H?zz%FS7$F:$\"3?9++++++gFS-F?6&FAF(F(Fd_l-%*AXESSTYLEGFH-% +AXESLABELSG6%Q!6\"F`[m-%%FONTG6#%(DEFAULTG-%%VIEWG6$;F(F:Fe[m" 1 2 0 1 10 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" "Curve 3" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT 257 27 "LEAST SQUARES, LINE AR CASE." }{MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "lpts:=line_pts();" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%lptsG7@7$$ \"*CJJ\"y!\")$!+]Y&Hk&!\"*7$$!+@V4#Q)F,$\"+CYzCZF)7$$\"*R.\"QEF)$\"+F% pj?\"F)7$$\"*nT9%>F)$\"+wso(R\"F)7$$\"*(yQ*=&F)$\"+tQ9$)HF,7$$\"*_(yyM F)$\"+Zx&\\%zF,7$$!+yQd&p(F,$\"+&e.5\\%F)7$$\"*`%pW\\F)$\"+vPS[KF,7$$! +@et`!*F,$\"+(Hgy'\\F)7$$!+!*puM^F,$\"+08f&e$F)7$$\").4b:F)$\"+Z\"eb)= F)7$$!+m`=QZF,$\"+rXk4NF)7$$\"*&*fBe*F)$!+m4UH7F)7$$\"*S<#F)7$$\"*i\\zp\"F)$\"+RZ 7s8F)7$$!++'z(=XF,$\"+%H#GmNF)7$$!+3-h:GF,$\"+9E2)*GF)7$$\"*17R\\*F)$! +vOv67F)7$$\"*/JCX\"F)$\"+$\\xUd\"F)7$$!+yeRT')F,$\"+QehE[F)" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "m:='m';b:='b';" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"mGF$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\" bGF$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 59 "E(m,b):=sum((m*lpts [i][1]+b-lpts[i][2])^2,i=1..nops(lpts));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>-%\"EG6$%\"mG%\"bG,hn*$),(*&$\"+yeRT')!\"*\"\"\"F'F1!\"\"F(F1$ \"+QehE[!\")F2\"\"#F1F1*$),(*&$\"*/JCX\"F5F1F'F1F1F(F1$\"+$\\xUd\"F5F2 F6F1F1*$),(*&$\"*17R\\*F5F1F'F1F1F(F1$\"+vOv67F5F1F6F1F1*$),(*&$\"+3-h :GF0F1F'F1F2F(F1$\"+9E2)*GF5F2F6F1F1*$),(*&$\"++'z(=XF0F1F'F1F2F(F1$\" +%H#GmNF5F2F6F1F1*$),(*&$\"*i\\zp\"F5F1F'F1F1F(F1$\"+RZ7s8F5F2F6F1F1*$ ),(*&$\"*=v\\`)F0F1F'F1F2F(F1$\"+?y$><#F5F2F6F1F1*$),(*&$\"+-/`%o&F0F1 F'F1F2F(F1$\"+F5F1F'F1F1F(F1$ \"+wso(R\"F5F2F6F1F1*$),(*&$\"*R.\"QEF5F1F'F1F1F(F1$\"+F%pj?\"F5F2F6F1 F1*$),(*&$\"+@V4#Q)F0F1F'F1F2F(F1$\"+CYzCZF5F2F6F1F1*$),(*&$\"*CJJ\"yF 5F1F'F1F1F(F1$\"+]Y&Hk&F0F1F6F1F1" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "solve(\{diff(E(m,b),m)=0,diff(E(m,b),b)=0 \},\{m,b\}) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$/%\"mG$\"+M3Ze;!\"*/%\"bG$!+0B] <9!\")" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "assign(%);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "f:=x->m*x+b;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fGf*6#%\"xG6\"6$%)operatorG%&arrowGF(,&*&%\"mG \"\"\"9$F/F/%\"bGF/F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 89 "display(plot(f(x), x=-10..10, color=blue,axes=boxed), plot(lpts,st yle=point,symbol=BOX));" }}{PARA 13 "" 1 "" {GLPLOT2D 567 567 567 {PLOTDATA 2 "6'-%'CURVESG6$7S7$$!#5\"\"!$!3*)******QJ(f2$!#;7$$!3\"pmm m\"p0k&*!#<$!3Vj25]Kn.IF-7$$!3uKL$3+(GF-7$$!3/nmm\"4m(G$)F1$!3#RtB'pQ!))z#F-7$$!3OLL$3i. 9!zF1$!3(yS\"G]q#zs#F-7$$!3fmm;/R=0vF1$!3-nm)e;:Am#F-7$$!3k++]P8#\\4(F 1$!3K-ox:V<%f#F-7$$!3Kmm;/siqmF1$!31Q(f7P1Q_#F-7$$!3Q****\\(y$pZiF1$!3 GuJj(4kOX#F-7$$!3jKLL$yaE\"eF1$!3+#*oiZT^\"Q#F-7$$!3s%HaF1$!3[A $*QKW'zJ#F-7$$!3]******\\$*4)*\\F1$!3sN4u/DUYAF-7$$!3p******\\_&\\c%F1 $!3soeu@F-7$$!3$)******\\1aZTF1$!3JSa3D)f`5#F-7$$!3Imm;/#)[oPF1$!3 IX$)[#3&\\U?F-7$$!3%HLLL=exJ$F1$!3)eJKJ#Gun>F-7$$!3lKLLL2$f$HF1$!3Qcx= ayT/>F-7$$!3%)****\\PYx\"\\#F1$!3?n%>6'evI=F-7$$!3gLLLL7i)4#F1$!3:W'4d ^_bw\"F-7$$!3n)***\\P'psm\"F1$!3#*\\2i6T,%p\"F-7$$!3?****\\74_c7F1$!3b J#pNj#*ei\"F-7$$!3M:LL$3x%z#)!#=$!3)ozMw,:[b\"F-7$$!3()HLL3s$QM%Fer$!3 =[m9F-7$$\"3fVLLezw5VFer$!3 :C'*ov%4gM\"F-7$$\"3..++v$Q#\\\")Fer$!3L*[dLc\\BG\"F-7$$\"3%\\LL$e\"*[ H7F1$!3c&\\=S6&f87F-7$$\"3>++++dxd;F1$!3f1OAS]cU6F-7$$\"3f+++D0xw?F1$! 33jw`qf2t5F-7$$\"35,+]i&p@[#F1$!3$37bDsTe+\"F-7$$\"3++++vgHKHF1$!3kB)f H`&*=J*F17$$\"3ElmmmZvOLF1$!3IhGF*e76k)F17$$\"3&4+++v+'oPF1$!3.tW8@'3 \\#zF17$$\"3UKL$eR<*fTF1$!3E&=hrP@fF(F17$$\"3L-++])Hxe%F1$!3'4!Q]**oSm lF17$$\"3!fmm\"H!o-*\\F1$!3aBWyp!4))*eF17$$\"3X,+]7k.6aF1$!3`an:yp(4?& F17$$\"3#emmmT9C#eF1$!3x5s\\c&=(=XF17$$\"33****\\i!*3`iF1$!3iZ$Ru@dW!Q F17$$\"38)F1$!3&** ****4_xbx#F-7$$\"3<+++I\"Q%pnF1$!3$)*****H\\>-n#F17$$!3X+++uN')\\$*F1$ !3*)*****fzE&>HF-7$$\"3/+++5BEh!*F1$\"3G+++sGu\"o(Fer7$$\"3,+++?HT9KF1 $!3K*****\\nx%Q&)F17$$!36*****>#*))3O*F1$!3++++BP%o-$F-7$$!3<+++#elRB# F1$!3-+++35b==F-7$$\"3))******pW_u`F1$!3K+++Y(*f0v)F1$!37+++lkhW-!475F1$!3,+++TcPI:F-7$$!3z* ****4C_9b*F1$!3;+++CF-7$$\"3(* ******>8Lo5F1$!3********f4U97F-7$$!3$******ReCXI\"F1$!3/+++?OHw:F-7$$ \"3#*******464[8F1$!3(*******)=,M7\"F-7$$\"3)*******HN,glF1$!3\"****** >P>58$F17$$\"36+++?YSjEF1$!3E+++Y'**G))*F17$$\"3%********p6j'>Fer$!3&* *****faqJS\"F-7$$!3*)*****>Y09(\\F1$!31+++8B:n@F--F\\[l6&F^[l$Fhz!\"\" F_[lF_[l-%&STYLEG6#%&POINTG-%'SYMBOLG6#%$BOXG-%*AXESSTYLEGFfel-%+AXESL ABELSG6%Q\"x6\"Q!F^fl-%%FONTG6#%(DEFAULTG-%%VIEWG6$;F(FgzFcfl" 1 2 0 1 10 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "LeastSquares(lpts,x);" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#,&$\"+xqEe>!\")\"\"\"*&$\"3NKUpPhA3L! # " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "" 0 "" {TEXT 258 10 "Cubic Case" }{MPLTEXT 1 0 0 "" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 381 "cubic_pts := proc()\n loca l fuzz, a, b, c, d,s,x,i;\n a:= myrand(-3,3);\n b:= myrand(-3,3); \n c:= myrand(-3,3);\n d:= myrand(-10,10);\n if (myrand(-10,10) \+ > 0) then s:=1 else s:=-1 fi;\n fuzz := [ random[normald[0,0.5]](poi nts_wanted)];\n x := [ random[uniform[-3,3]](points_wanted)];\n \+ RETURN( [ [x[i], s*(a-x[i])*(b-x[i])*(c-x[i])+d+fuzz[i]] $i=1..points _wanted]);\nend:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "cpts:=cubic_pts();" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%cptsG7@7$$\"+HUES9!\"*$\"+7KrM(*!#57$$\"+&y2$>C F)$\"+tZVlOF)7$$\"+6tC=F,7$$!*%*HGi$F)$\"+5Qn8!\")7$$!+/`6_=F)$!+DVhJwF)7$ $\"+K\"e::\"F)$\"+!*3o()>F)7$$!+JB%*G=F)$!+G)46F(F)7$$!+Mui)>\"F)$\"+v vB=JF,7$$\"+PTlX=F)$\"+X>Xw:F)7$$\"+ffu\"Q\"F)$\"+9@)*\\AF)7$$!)bO$G\" F)$\"+uKk#4*>F)$!+7b^l5FU7$ $\"+xeXw7F)$\"+U0xE:F)7$$\"+Nkl&R\"F)$\"+P&[!p:F)7$$\"*91\"**)*F)$\"+I #*)pt#F)7$$!+'>Tg:#F)$!+1\"Q;Z\"FU7$$!+:d^(f\"F)$!+R(H'oRF)7$$!)LIc'*F )$\"+QVY;PF)7$$\"+YI0v?F)$\"+hm&e5#F)7$$\"*<2IZ)F)$\"+:mi4HF)7$$\"+`([ B#GF)$\"+9zv^hF)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "cub:=(x ,a,b,c,d)->a*x^3+b*x^2+c*x+d;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$cu bGf*6'%\"xG%\"aG%\"bG%\"cG%\"dG6\"6$%)operatorG%&arrowGF,,**&9%\"\"\") 9$\"\"$F3F3*&9&F3)F5\"\"#F3F3*&9'F3F5F3F39(F3F,F,F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 74 "E:=(a,b,c,d)->sum((cub(cpts[i][1],a,b,c,d )-cpts[i][2])^2,i=1..nops(cpts));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %\"EGf*6&%\"aG%\"bG%\"cG%\"dG6\"6$%)operatorG%&arrowGF+-%$sumG6$*$),&- %$cubG6'&&%%cptsG6#%\"iG6#\"\"\"9$9%9&9'F>&F96#\"\"#!\"\"FEF>/F<;F>-%% nopsG6#F:F+F+F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "a:='a';b :='b';c:='c';d:='d';" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"aGF$" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bGF$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cGF$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"dGF$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 103 "solve(\{diff(E(a,b,c,d),a)=0,diff( E(a,b,c,d),b)=0,diff(E(a,b,c,d),c)=0,diff(E(a,b,c,d),d)=0\},\{a,b,c,d \});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<&/%\"aG$\"+klGN(*!#5/%\"dG$\" +V8b\"G%!\"*/%\"bG$!+Y4z*=#F-/%\"cG$!+ev#4#oF(" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 10 "assign(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 15 "cub(x,a,b,c,d);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,* *&$\"+klGN(*!#5\"\"\")%\"xG\"\"$F(F(*&$\"+Y4z*=#!\"*F()F*\"\"#F(!\"\"* &$\"+ev#4#oF'F(F*F(F2$\"+V8b\"G%F/F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 97 "display(plot(cub(x,a,b,c,d), x=-5..5, color=blue,axes =boxed), plot(cpts,style=point,symbol=BOX));" }}{PARA 13 "" 1 "" {GLPLOT2D 567 567 567 {PLOTDATA 2 "6'-%'CURVESG6$7S7$$!\"&\"\"!$!3.++y 0%Quo\"!#:7$$!3XLLLe%G?y%!#<$!3+bHGv/#**[\"F-7$$!3OmmT&esBf%F1$!3YvO>7 EdI8F-7$$!3@LL$3s%3zVF1$!3:'=h2^eZ;\"F-7$$!3^LL$e/$QkTF1$!3%p-*3'\\3;, \"F-7$$!3ommT5=q]RF1$!3N>/7C%\\Ks)!#;7$$!3HLL3_>f_PF1$!3W>\"GOPISa(FH7 $$!3K++vo1YZNF1$!3'f!f0+0tJkFH7$$!3;LL3-OJNLF1$!3S!y@t.ICR&FH7$$!3o*** \\P*o%Q7$F1$!3X5KGiaMjWFH7$$!3Kmmm\"RFj!HF1$!3)>*RO7+=8OFH7$$!33LL$e4O Zr#F1$!3rVvt3yC[HFH7$$!3u*****\\n\\!*\\#F1$!3es'*)R-q$)G#FH7$$!3%)**** *\\ixCG#F1$!3iI;5`6g9JVg$F17$$!3mdmmT&Q(RTFdq$\"3nVM1,%z&>TF17$$!3$\\ mmTg=><#Fdq$\"3K([hHiCkJ%F17$$!3uDMLLe*e$\\!#?$\"3b#*fp9Y'[G%F17$$\"3z rm;zRQb@Fdq$\"3-%fv>k^D/%F17$$\"3^,+](=>Y2%Fdq$\"3j\"=@4?Afq$F17$$\"3r umm\"zXu9'Fdq$\"3%Qe=^ah3E$F17$$\"3#4+++]y))G)Fdq$\"3li!y$oB3mFF17$$\" 3H++]i_QQ5F1$\"3ed<`\")=9-BF17$$\"3a++D\"y%3T7F1$\"3%yQf3&*GJ#>F17$$\" 3+++]P![hY\"F1$\"3f*=GbBRDk\"F17$$\"3iKLL$Qx$o;F1$\"3y_Mu+rGp:F17$$\"3 Y+++v.I%)=F1$\"3T\")>W*Q4Xt\"F17$$\"3?mm\"zpe*z?F1$\"3C*y$39\"o%\\@F17 $$\"3;,++D\\'QH#F1$\"3w9RS$\\M]%HF17$$\"3%HL$e9S8&\\#F1$\"3Ylu?rmSpSF1 7$$\"3s++D1#=bq#F1$\"3Za7hsb$po&F17$$\"3\"HLL$3s?6HF1$\"3.h_I\">4ov(F1 7$$\"3a***\\7`Wl7$F1$\"3e)GU+v%p\\5FH7$$\"3enmmm*RRL$F1$\"3\"y@(H/MRu8 FH7$$\"3$zmmTvJga$F1$\"3ag9hNKjtCF1$\"3')*****HxMam$F17$$\"30+++6tC=xLh)=F17$$!3\"******4Z1Ug\"F1$!3))*****zRsu8%F17$$!3!******R)=aY6F1$ \"3++++c72&*>Fdq7$$!3C+++S*HGi$Fdq$\"3R+++5Qn8FH7$$!3,+++/`6_=F1$!39+++DVhJwF17$$\"3!******>8e::\" F1$\"3.+++!*3o()>F17$$!3#******4LU*G=F1$!3;+++G)46F(F17$$!3/+++Mui)>\" F1$\"3?+++vvB=JFdq7$$\"3\"******p8ac%=F1$\"32+++X>Xw:F17$$\"32+++ffu\" Q\"F1$\"3\")*****R6#)*\\AF17$$!3%********\\lLG\"!#>$\"3#******RFVw\"QF 17$$\"3-+++&))=qP#F1$\"3-+++#H&[vLF17$$\"3.+++rtN8CF1$\"3K+++vS!oN%F17 $$!37++++([oC&Fe`l$\"32+++&[M%eVF17$$\"3M******\\@OCyFdq$\"3))*****HE* oOHF17$$!3-+++%*>#4*>F1$!3*******>^:b1\"FH7$$\"32+++xeXw7F1$\"3&****** >aqn_\"F17$$\"3-+++Nkl&R\"F1$\"3-+++P&[!p:F17$$\"3A+++Sh5**)*Fdq$\"3$) ******H#*)pt#F17$$!3))*****f>Tg:#F1$!3)******f5Q;Z\"FH7$$!35+++:d^(f\" F1$!3#)******Q(H'oRF17$$!3k++++LIc'*Fe`l$\"3?+++QVY;PF17$$\"3%)*****f/ `]2#F1$\"3!)*****4mce5#F17$$\"3!)******pr+t%)Fdq$\"33+++:mi4HF17$$\"3* ******Hv[B#GF1$\"30+++9zv^hF1-F][l6&F_[l$\"#5!\"\"F`[lF`[l-%&STYLEG6#% &POINTG-%'SYMBOLG6#%$BOXG-%*AXESSTYLEGFiel-%+AXESLABELSG6%Q\"x6\"Q!Faf l-%%FONTG6#%(DEFAULTG-%%VIEWG6$;F(FhzFffl" 1 2 0 1 10 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "a:='a';b:='b';c:='c';d:='d';" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%\"aGF$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bGF$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cGF$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"dGF$ " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 370 "circle_pts := proc()\n local dat, i, rad, rads, angles, a, b;\n rad := myrand(3,10 );\n a := myrand(-10,10);\n b:= myrand(-10,10);\n rads := [ random[normald[rad,rad/50]](points_wanted)];\n angles := [ \+ random[uniform[0,2*Pi]](points_wanted) ];\n RETURN ( [ [a+rads[i] *cos(angles[i]), b+rads[i]*sin(angles[i])] $ i=1..points_wanted ]): \n end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "cpt:=circle_pts( );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$cptG7@7$$\"+6<^$***!\"*$\"*Qe 7j(F)7$$\"+DU(F) $\"**)H$faF)7$$\"+f+DD7!\")$\"+\"p8&QRF)7$$\"+MejSgF)$\"+#=>%p6F)7$$\" +?'f%zbF)$\"+:l$3')F)$\"+u]xix F)7$$\"+H\"\\RA\"F>$\"+([4Zz$F)7$$\"+)yy:<\"F>$\"+Dc(>$>F)7$$\"+#G(>c] F)$\"+MtozIF)7$$\"+JwXu6F>$\"+af'f6'F)7$$\"+g&30?\"F>$\"+ldc3CF)7$$\"+ 'p\\%[nF)$\"*_vre)F)7$$\"+c`[]bF)$\"+TSTHjF)7$$\"+\"4Co?\"F>$\"+IZdSaF )7$$\"+)y`\"*p(F)$\"+@3zgwF)7$$\"+Q)4O/\"F>$\"+4v9OtF)7$$\"+ATwupF)$\" *Io'3xF)7$$\"+UfK27F>$\"+V*=x3$F)7$$\"+/F:]\"*F)$\"*N5q=$F)7$$\"+&p.U7 \"F>$\"+v2^\"\\\"F)7$$\"+0!GKX'F)$\"+L[EprF)7$$\"+KHM$\"+o*[eX%F)7$$\"+Ai<'*pF)$\"**H%>(fF)" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "epsilon:=(pt,a,b,c)->(pt[1]- a)^2+(pt[2]-b)^2-r^2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(epsilonGf* 6&%#ptG%\"aG%\"bG%\"cG6\"6$%)operatorG%&arrowGF+,(*$),&&9$6#\"\"\"F69% !\"\"\"\"#F6F6*$),&&F46#F9F69&F8F9F6F6*$)%\"rGF9F6F8F+F+F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "E:=(a,b,r,cpt)->sum(epsilon(cpt[i], a,b,r)^2,i=1..nops(cpt));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"EGf*6 &%\"aG%\"bG%\"rG%$cptG6\"6$%)operatorG%&arrowGF+-%$sumG6$*$)-%(epsilon G6&&9'6#%\"iG9$9%9&\"\"#\"\"\"/F:;F?-%%nopsG6#F8F+F+F+" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 101 "solutions:=\{solve(\{diff(E(a,b,r, cpt),a)=0, diff(E(a,b,r,cpt),b)=0, diff(E(a,b,r,cpt),r)=0\},\{a,b,r\}) \};" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%*solutionsG<)<%/%\"rG$\"\"!F* /%\"aG$\"+_&)4x')!\"*/%\"bG$\"+T5$3%RF/<%F'/F,^$$\"+F8SP%)F/$!+\"*=ieb F//F1^$$\"+..O_gF/$!+T^r=\"*!#5<%/F,^$F7$\"+\"*=iebF/F'/F1^$F=$\"+T^r= \"*FA<%/F,^$$\"+Op&\\`*F/$!+$3\\5a\"F/F'/F1^$$\"+kXaYRF/$\"+%fi!o^F/<% /F,^$FN$\"+$3\\5a\"F//F1^$FT$!+%fi!o^F/F'<%/F,$\"+F@Ct&)F//F1$\"+kAnoS F//F($!+Ld9CPF/<%F\\oF_o/F($\"+Ld9CPF/" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "finalsol:=op(select(s->if(subs(s,r)>0)then true else \+ false fi,solution));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)finalsolG<% /%\"aG$\"+F@Ct&)!\"*/%\"bG$\"+kAnoSF*/%\"rG$\"+Ld9CPF*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 121 "display(plot(subs(finalsol,[a+r*co s(t),b+r*sin(t),t=0..2*Pi]), color=blue,axes=boxed), plot(cpt,style=po int,symbol=BOX));" }}{PARA 13 "" 1 "" {GLPLOT2D 600 600 600 {PLOTDATA 2 "6'-%'CURVESG6$7S7$$\"31+++'yQ(H7!#;$\"3#******REs'oS!#<7$$\"3d*QXVe ^iA\"F*$\"3;7Izb<7xXF-7$$\"3R&=#y+2f<7F*$\"3YWJ$4*f57]F-7$$\"3A=`T^]v, 7F*$\"3b&>WG(f+&[&F-7$$\"3;Pe4%=x&z6F*$\"3@llOp3PNfF-7$$\"3A#)y]x+p^6F *$\"3.#>9vjH*\\jF-7$$\"3oUzQ'f*3@6F*$\"3[efkY)[xp'F-7$$\"3h%Q(QQ*4^3\" F*$\"350w1sF'\\,(F-7$$\"3w+Y!ykKR/\"F*$\"3qvnz:xb\"H(F-7$$\"35)\\;F!y \"f***F-$\"3su!*F-$\"3X$>)zy>(*exF-7$$\"3Nczf8%=5d)F-$\"3'ea-JLC5!R[`F-$\"3W#[,H/c8$fF-7$$\"3Y'R!oLp]N^F- $\"3[6**p'p')3]&F-7$$\"3\"e\\87!zOu\\F-$\"3D]HG&Qkk-&F-7$$\"3:O&['p(>P )[F-$\"372JwhVJvXF-7$$\"3taPo.V6\\[F-$\"3,DYI#*>A!3%F-7$$\"3:ex=5a>$)[ F-$\"3Dy[([[jec$F-7$$\"3e.l$R;z/(\\F-$\"3mM*\\0;8c7$F-7$$\"32Xsc#GoM7& F-$\"32Di>=LqlEF-7$$\"3=W)*))*z[GM&F-$\"3d\"Q\">C_h:AF-7$$\"3,1fa)o5Sh &F-$\"3o<#\\,U_w!=F-7$$\"3o>mhZP;DfF-$\"3=DZ%)oe4]9F-7$$\"3_$\\fDRg1K' F-$\"3cHR#fk7I5\"F-7$$\"3v-ECHqj9nF-$\"3Ae[&y=yYT)!#=7$$\"36\")3L@L?or F-$\"3o^=jH#4u>'Fhu7$$\"3Cl0mZqt,wF-$\"30!yTkibZt%Fhu7$$\"3!y35v$QC#4) F-$\"3hI*=TL#>dPFhu7$$\"3\"[PZb)f&=c)F-$\"3HsxSp$RaW$Fhu7$$\"3AM*HI.4G 0*F-$\"39%4lB\"3LbPFhu7$$\"3c)=C&fnxC&*F-$\"3]D2hT^Q\"o%Fhu7$$\"3$oB-Y $\\<+5F*$\"3U]=&3K\\RH'Fhu7$$\"3s.Rh4U(F-$\"3z*******))H$faFhu7$$\"3/+++f+DD7F*$\"3#) *****4p8&QRF-7$$\"3q*****R$ejSgF-$\"3%******>=>%p6F-7$$\"3t******>'f%z bF-$\"39+++: l$3')F-$\"3L+++u]xixF-7$$\"3-+++H\"\\RA\"F*$\"3)******p[4Zz$F-7$$\"31+ ++)yy:<\"F*$\"3(******\\iv>$>F-7$$\"3H+++#G(>c]F-$\"3%******RL(ozIF-7$ $\"3%******4jdW<\"F*$\"3)******R&f'f6'F-7$$\"3#*******f&30?\"F*$\"33++ +ldc3CF-7$$\"3x*****fp\\%[nF-$\"3%*******>b<(e)Fhu7$$\"3E+++c`[]bF-$\" 3r*****4/9%HjF-7$$\"3,+++\"4Co?\"F*$\"3!)******HZdSaF-7$$\"3E+++)y`\"* p(F-$\"3(******4#3zgwF-7$$\"32+++Q)4O/\"F*$\"3!*******3v9OtF-7$$\"3%)* ****>7kZ(pF-$\"3.++++$o'3xFhu7$$\"3%******>%fK27F*$\"3/+++V*=x3$F-7$$ \"3\"******Rq_,:*F-$\"3@+++].,(=$Fhu7$$\"3#******\\p.U7\"F*$\"31+++v2^ \"\\\"F-7$$\"3r*****\\+GKX'F-$\"3s*****H$[EprF-7$$\"3.+++KHMAwh*pF-$\"3$)*******)H%>(fFhu -Fhz6&Fjz$\"#5!\"\"F[[lF[[l-%&STYLEG6#%&POINTG-%'SYMBOLG6#%$BOXG-%*AXE SSTYLEGFdel-%+AXESLABELSG6%Q!6\"F[fl-%%FONTG6#%(DEFAULTG-%%VIEWG6$F`fl F`fl" 1 2 0 1 10 0 2 9 1 2 2 1.000000 45.000000 45.000000 0 0 "Curve 1 " "Curve 2" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK " 54 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }