{VERSION 4 0 "IBM INTEL NT" "4.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 "Helvetica" 1 14 128 0 0 1 0 0 2 0 0 0 0 0 0 1 } {CSTYLE "" -1 257 "" 0 24 0 0 255 1 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 258 "Helvetica" 0 1 128 0 0 1 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 259 "Helvetica" 0 1 128 0 0 1 0 0 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 260 " " 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 261 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 263 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }{PSTYLE "Norm al" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 } 1 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE " " -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Map le Plot" -1 13 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }3 1 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Title" -1 256 1 {CSTYLE " " -1 -1 "Times" 1 18 0 0 0 1 2 1 1 2 2 2 1 1 1 1 }1 1 0 0 12 12 1 0 1 0 2 2 19 1 }} {SECT 0 {EXCHG {PARA 256 "" 0 "" {TEXT 257 27 "Some Interpolation Exam ples" }}{PARA 0 "" 0 "" {TEXT 256 18 "Mth 351 Aug 7 2001" }}{PARA 0 " " 0 "" {TEXT 258 16 "Bent E. Petersen" }}{PARA 0 "" 0 "" {TEXT 259 38 "Filename: 351u2001_interp_examples.mws" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 73 "In this worksheet we illustrate some of Maple's interpolation faci lities." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT 260 57 "Example: Cubic Interpolation Spline for \+ sin(x) and cos(x)" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 55 "XX1:=[0,Pi/4,Pi/2,3*Pi/4,Pi,5*Pi/4,3*Pi/2,7* Pi/4,2*Pi];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$XX1G7+\"\"!,$%#PiG# \"\"\"\"\"%,$F(#F*\"\"#,$F(#\"\"$F+F(,$F(#\"\"&F+,$F(#F1F.,$F(#\"\"(F+ ,$F(F." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "YY1:=map(sin,XX1) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$YY1G7+\"\"!,$*$-%%sqrtG6#\"\"# \"\"\"#F-F,F-F'F&,$F(#!\"\"F,F1F/F&" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "ex1:=spline(XX1,YY1,x,cubic):" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 89 "plot(ex1,x=0..2*Pi,thickness=3,color=blue, tit le=\"cubic spline interpolation of sin(x)\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7en7$$\"\"!F)F(7$$\"3 i]cC&eb&p8!#=$\"3%*>gL0PQi8F-7$$\"3YqR*pd)>hDF-$\"3j1h'*>>%)GDF-7$$\"3 Zs[E^dK,RF-$\"3U^?%G;$o)z$F-7$$\"3ab^NehL]_F-$\"3m?qT#eV)4]F-7$$\"35*Q hW&\\$Hf'F-$\"3=Yu'yU*RDhF-7$$\"3zS11.fpPyF-$\"3CX=M%Ro&fqF-7$$\"3upr' fwtl7*F-$\"3%z(fyyTF2zF-7$$\"3!QRa&4L&f/\"!#<$\"3e$4ec(=6X')F-7$$\"3Yj Xwg<#)y6FQ$\"3q`Bk%)y#4B*F-7$$\"36qGW%H$\\:8FQ$\"3-4&44g8zm*F-7$$\"3SH 22vMov8FQ$\"3='>x2yTX!)*F-7$$\"3$*)e)pbO(eV\"FQ$\"3uz)pRk>-OUv**F-7$$\"3!=3SCmpuq\"FQ$\"33h'zc3zL!**F-7$$\"3 3u?N%*p.t!z*F-7$$\"3LmSEEVgQ=FQ$\"31(e(zd<FQ$\"3)=+uoi)F-7$$\"3 Ua[#o!GC>AFQ$\"39B8>:SzlzF-7$$\"3-YwJf&y(eBFQ$\"3w'**p&HC\"G0(F-7$$\"3 oNrpV8H#[#FQ$\"3x]S*GIha7'F-7$$\"3EzY)QW/yh#FQ$\"3)GB?%*QS!**\\F-7$$\" 35v)>8'\\%ou#FQ$\"3vvPd1@OTQF-7$$\"3w$GCS?&[\")GFQ$\"3]v/G]&otc#F-7$$ \"3Z%)='p(p70IFQ$\"3[)f6WSRvN\"F-7$$\"3sA%)\\K8\\QJFQ$\"3'*))fIK=E%4$! #?7$$\"3c&fJlT>qF$FQ$!3)p]w3'[EZ8F-7$$\"39l2\"4_3wR$FQ$!3F`nr%oOy_#F-7 $$\"3MOnGd![y_$FQ$!3G=F'GenFw$F-7$$\"3#*=4JU#)RiOFQ$!3k(G7ReOJ(\\F-7$$ \"3%G,f%[$HSz$FQ$!3#elUncg42'F-7$$\"3KoE+7#*Q@RFQ$!3`q685@VJqF-7$$\"3y (f0ii+G1%FQ$!3WIjn]'Q*ezF-7$$\"3;fORr]')*=%FQ$!3uXDOh[kc')F-7$$\"3R&p% *z[LbK%FQ$!3+\\+)4K$Q]#*F-7$$\"3M'pExBp%[WFQ$!3Cx;\"fX6bk*F-7$$\"3,z&[ 2')pc^%FQ$!3![bq-(4N,)*F-7$$\"3oh/x$[qGe%FQ$!3dD$RL)>:8**F-7$$\"3qQq'H .,hk%FQ$!3=Y]k*4pq(**F-7$$\"3e;O;#eJ$4ZFQ$!3/jdzJ3&*****F-7$$\"3sU'G^s Dax%FQ$!3W-:%[FQ$!3^F(=lzmO\"**F-7$$\"3;@9vt*Q h!\\FQ$!3(R\\()=KHs!)*F-7$$\"3Ou\"4%z!e2(\\FQ$!31j7q&=B+m*F-7$$\"3*zX# [4&eg5&FQ$!39`K%38MWA*F-7$$\"3n*e![/*ojB&FQ$!3![Ya7Q4%\\')F-7$$\"3#ob8 X5I'p`FQ$!30x#>GNEF\"zF-7$$\"3PA\"GX$yy,bFQ$!3MG3(Rb$yUqF-7$$\"39L1/jq ABcFQ$!3YP9%Rch18'F-7$$\"3mweKB,TidFQ$!3is4\"QHsG(\\F-7$$\"3'oIe;6(*o) eFQ$!3))Hz?R/fbQF-7$$\"3nKcu0ii>gFQ$!3;Q?P3bo+EF-7$$\"3w)*zj%)[mYhFQ$! 3%H/*R'y!3e8F-7$$\"3)****>YH&=$G'FQ$!3Gwak5D5`7!#D-%&TITLEG6#QEcubic~s pline~interpolation~of~sin(x)6\"-%'COLOURG6&%$RGBGF(F($\"*++++\"!\")-% +AXESLABELSG6$Q\"xFd]lQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F($\"+3`=$G' !\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 105 "plot(sin(x) -ex1,x=0..2*Pi,thickness=3,color=blue, title=\"error in cubic spline i nterpolation of sin(x)\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7ev7$$\"\"!F)F(7$$\"3;`#Gi#zxZo!#>$\"3)[z$3 BnnH:!#@7$$\"3i]cC&eb&p8!#=$\"31oS\\KVZ%*GF07$$\"3q5)>63x`'>F4$\"3='*3 !3zH0$QF07$$\"3YqR*pd)>hDF4$\"3U%otENMqW%F07$$\"3\\LyxtWrGFF4$\"3#Q>@^ f!RaXF07$$\"3'fph0PIi*GF4$\"3K2+_Jo,JYF07$$\"3WebMniujIF4$\"37fp&\\y#H wYF07$$\"3Y@%HT;i7B$F4$\"3iD7?dI#)*o%F07$$\"3[%G841y()R$F4$\"3#e([lP( \\9n%F07$$\"3'p9(pdRHmNF4$\"3')\\RmX\\E@YF07$$\"3W45[a)4Qt$F4$\"3&*e1n -LiRXF07$$\"3Zs[E^dK,RF4$\"3#p&*)pIV:FWF07$$\"3+9+\"[&4$ed%F4$\"33NDfs @l&o$F07$$\"3ab^NehL]_F4$\"3=W1)=OZtd#F07$$\"3Ls#3kbN;#fF4$\"3\\YptrHe 68F07$$\"35*QhW&\\$Hf'F4$\"39smX#e5B'=!#A7$$\"3J-7h\">DT!pF4$!3iq0(pV, aw\"Fgp7$$\"3]:5wGaJ:sF4$!3d_-wGs08PFgp7$$\"3c@fLZ0\"4P(F4$!3Yzt0Mz`$) QFgp7$$\"3fF3\"fm0l_(F4$!3eob\"4Y`$=MFgp7$$\"3uMd[%y+@o(F4$!3x&3+yL@vB #Fgp7$$\"3zS11.fpPyF4$!3e)yZD0F0d#!#B7$$\"3'QF(yo`\"*f\")F4$\"3`U7K!** fn]'Fgp7$$\"3#e!R^M[8#[)F4$\"3x,iGAf$eh\"F07$$\"3yP0C+VN/))F4$\"3K&*4' =;Yex#F07$$\"3upr'fwtl7*F4$\"32u-#z%)e[/%F07$$\"3#=Oh$)f8)f%*F4$\"39+E U2!=&*Q&F07$$\"3!Rbb2V`Iz*F4$\"3)Gn.F$=r\"p'F07$$\"3gu\\JE$HE,\"!#<$\" 3-Z&p'>(o_)yF07$$\"3!QRa&4L&f/\"Fht$\"3/!*[N]>D9*)F07$$\"31OpNA/$*Q*f`(p9Fht$\"3&)Q%zSc&*G/#F07$$\"3/]+3VNj.:Fht$\"3LY%z R?V&=**Fgp7$$\"3q!yqn[8v`\"Fht$\"3o*)ythO_lEFgp7$$\"396:YIMRr:Fht$\"3` GzM&=+\\K*!#E7$$\"3\"Q:c%)[7ag\"Fht$\"3)GuICRmL(GFgp7$$\"3Z'z]kaJ%R;Fh t$\"3'[W;%ffnJ5F07$$\"39RaW/1Xt;Fht$\"3')yNcED\"*)4#F07$$\"3!=3SCmpuq \"Fht$\"3)3kbQ*o?oLF07$$\"3%z2'RGLDSo%F07$$\"33u?N%*p. t\\u,3O)F07$$\"3')4G5Xd9)*=Fht$\"3w#[jHt.7*)*F07$$\"3S`:%R ;(od>Fht$\"3on^#)e_nh5Fju7$$\"3)>qL/H(Qv>Fht$\"3'y!3#f*)yg1\"Fju7$$\"3 M]e#pT(3$*>Fht$\"3%[!))oF-Xi5Fju7$$\"3:**zTVvy5?Fht$\"3M1V\"pvS30\"Fju 7$$\"3^Z,\"*pw[G?Fht$\"3;)z(e7:RJ5Fju7$$\"3oWW*G#z)Q1#Fht$\"31x%>Zm5)* p*F07$$\"3#=uye<)G*4#Fht$\"3%o#\\CRO46))F07$$\"34q_hLoFH@Fht$\"33ZD1!> FZ(yF07$$\"3!zz^8\\l#f@Fht$\"3r]Jx\")[L0oF07$$\"3;E$)3\\TD*=#Fht$\"3Gj +))**)zLk&F07$$\"3Ua[#o!GC>AFht$\"3Jh_j[<0OWF07$$\"3W_![\\uETD#Fht$\"3 g.eg:\"\\i/$F07$$\"3W]72$o5!*G#Fht$\"3-L3>b_nj@Y*QK#Fht$ \"3))f)3]'G5apFgp7$$\"3-YwJf&y(eBFht$!3(4l@\"y#fL.%Ffr7$$\"3;#3ltlBh=Fgp7$$\"3C 24z$*y/]DFht$\"3q.t>3r:B8F07$$\"3EzY)QW/yh#Fht$\"3^?V?X$o+g#F07$$\"3=x Ag-ZK#o#Fht$\"3YaIw7lyhOF07$$\"35v)>8'\\%ou#Fht$\"3[ye@qt'3R%F07$$\"3S Ezl\"*\\njFFht$\"34FD\"Q-00y#Fht$\"3(4[X%eze-YF07$$ \"37GSL_]L(z#Fht$\"3AO*>[K!RhYF07$$\"3Vz?n#3lT\"GFht$\"3qOI=3.:)o%F07$ $\"3uI,,8^*4$GFht$\"3qWoK2nq#o%F07$$\"3g\"=[L9Dy%GFht$\"3k>/#pLc^k%F07 $$\"3YKiot^lkGFht$\"31n4p17%ed%F07$$\"3w$GCS?&[\")GFht$\"3G(Q47jN`Z%F0 7$$\"3*Q3$\\!41L%HFht$\"3G4!f)f(=N&QF07$$\"3Z%)='p(p70IFht$\"3#>$GYRSm &)GF07$$\"3g`,ta\"4=2$Fht$\"3E*4'Rx:vd:F07$$\"3sA%)\\K8\\QJFht$\"3Wl&4 (\\vGaqFfr7$$\"3#*3]^u`v2KFht$!3m&)QM^=vz9F07$$\"3c&fJlT>qF$Fht$!3E'*y jxh#p'GF07$$\"3e!=@(oRJPLFht$!3W$*Rt5Kq>QF07$$\"39l2\"4_3wR$Fht$!3;#\\ Q^ctiW%F07$$\"3\\hx&HY))QT$Fht$!3C#3g$fr?^XF07$$\"3%yv/]So,V$Fht$!3#=% *[WlRri%F07$$\"3=a<0Z$[kW$Fht$!3!=HX+#>\\tYF07$$\"3_]()4*GGFY$Fht$!31F HEBZ))*o%F07$$\"3)ouX6B3!zMFht$!31$G#=k'\\hn%F07$$\"3?VF>t\")G&\\$Fht$ !3EtPWb+MKYF07$$\"3+S(R_6o:^$Fht$!39d[+V \\4#*RFht$!37RAyTC'Gp\"F07$$\"3)\\'[lsxWFSFht$!3ihnCqU)Q)HF07$$\"3y(f0 ii+G1%Fht$!3o0Sr+[))*Q%F07$$\"378E]Pnc%4%Fht$!3=t+2;_locF07$$\"3ZG'*z[ GLETFht$!3Y4)R@C!R&*oF07$$\"3#Qk'4g*)4eTFht$!3K?kD$H2W,)F07$$\"3;fORr] ')*=%Fht$!3[GLaAd]y*)F07$$\"3))=Rav@yBUFht$!3%R>bN!z`$z*F07$$\"3sxTpz# *pdUFht$!3X5.-HIBN5Fju7$$\"392$p<$yluUFht$!3'\\joJcjF0\"Fju7$$\"3cOW%Q Q;;H%Fht$!3'zBU&Re2j5Fju7$$\"3(fc>f$\\d3VFht$!3->8k!>ag1\"Fju7$$\"3R&p %*z[LbK%Fht$!3E\\MMS%e;1\"Fju7$$\"3'epgGO,qQ%Fht$!39hbXtWx_)*F07$$\"3M 'pExBp%[WFht$!3L/nT[`@P#)F07$$\"3oPwB\\&p?[%Fht$!39%QE9OMO2(F07$$\"3,z &[2')pc^%Fht$!33es98::!y&F07$$\"3O?&fA](o$ed[9YFht$!3*H]=\"=%3^$>F07$$\"3 qQq'H.,hk%Fht$!3;F774L%[o*Fgp7$$\"33G`c2jrxYFht$!3%4;\\n'3L#)GFgp7$$\" 3e;O;#eJ$4ZFht$!3'>`2JB@)HC!#C7$$\"3lHhk`'yBu%Fht$!3sDzn->t#=#Fgp7$$\" 3sU'G^sDax%Fht$!3g_h.xkmN))Fgp7$$\"3zb6h'zs%3[Fht$!3]*[FKv$[t=F07$$\"3 &)oO4o)>:%[Fht$!3SlNUU=1vIF07$$\"3WXD#4UHQ([Fht$!3ar%\\72#QdVF07$$\"3; @9vt*Qh!\\Fht$!3)\\4QB+.=m&F07$$\"3w(H!eE&[%Q\\Fht$!35A#f0oKb\"pF07$$ \"3Ou\"4%z!e2(\\Fht$!3m_%HG\"3Zc!)F07$$\"3t:eW%H3%Q]Fht$!39IHD_[Cl)*F0 7$$\"3*zX#[4&eg5&Fht$!3=O*HNI5V1\"Fju7$$\"3yBs&)3tMA^Fht$!3I/X`#[Ee1\" Fju7$$\"3Y!*>B3hjQ^Fht$!3#R$pn:Wag5Fju7$$\"39dng2\\#\\:&Fht$!3?1)*f#oA &[5Fju7$$\"3#Q_\")pq87<&Fht$!3dQ'[&eR))H5Fju7$$\"3>d5t08z._Fht$!3yYjsV !Qct*F07$$\"3n*e![/*ojB&Fht$!3u*)oU#4&GQ*)F07$$\"3-J))[/Uop_Fht$!3c!*p veZ]9zF07$$\"3Dtq\\/&**HI&Fht$!3Qo'ex%G\"\\s'F07$$\"3[:`]/[JO`Fht$!34! pm7WI^U&F07$$\"3#ob8X5I'p`Fht$!3WUA%)\\$Fgp7$$\"3O9.AXvK ZbFht$\"3MSduZx>#*QFgp7$$\"3wxVy[u]ibFht$\"3[oNFmcv6_yaEF07$$\"3w\"4# \\%[#)[.fFht$!3q1Y)o,&z+XF07$$\"3(y8!=&Qz+#fFht$!3qqKw_@S$f%F07$$\"3Q` 5%>_qm$fFht$!3M/EdMyNbYF07$$\"3yp>qe;E`fFht$!3WgH^cW@'o%F07$$\"3;')GY& z_)pfFht$!37o#Qjt\"y&o%F07$$\"3o,QAKRW')fFht$!37]#Ht&[6aYF07$$\"3<\\UD]\"f%F07$$\"3nKcu0ii>gFht$!3;:k3k]X)\\%F07$$\"3G:=>X b9$3'Fht$!3y'GNx%\\GwQF07$$\"3w)*zj%)[mYhFht$!3))[N&GU\\m)GF07$$\"3P** *G'*3D\\@'Fht$!3j>ifJH*\\_\"F07$$\"3)****>YH&=$G'Fht$!3!=vlr$)Rh&G!#G- %&TITLEG6#QNerror~in~cubic~spline~interpolation~of~sin(x)6\"-%'COLOURG 6&%$RGBGF(F($\"*++++\"!\")-%+AXESLABELSG6$Q\"xFjenQ!6\"-%*THICKNESSG6# \"\"$-%%VIEWG6$;F($\"+3`=$G'!\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 18 "YY2:=map(cos,XX1);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$YY2G7+\"\"\",$*$-%%sqrtG6#\"\"#F&#F&F,\"\"!,$F(#!\"\"F,F1F/F.F'F &" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "ex2:=spline(XX1,YY2,x, cubic):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 88 "plot(ex2,x=0..2* Pi,thickness=3,color=blue,title=\"cubic spline interpolation of cos(x) \");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVE SG6#7W7$$\"\"!F)$\"\"\"F)7$$\"3i]cC&eb&p8!#=$\"3]\"f4,%Rkn'*F/7$$\"3Yq R*pd)>hDF/$\"3?ztpxoL_$*F/7$$\"3Zs[E^dK,RF/$\"37<4%fFr)R*)F/7$$\"3ab^N ehL]_F/$\"3=V([bpQ@V)F/7$$\"35*QhW&\\$Hf'F/$\"30dX$[i7H!yF/7$$\"3zS11. fpPyF/$\"3M1ue!)H]\"3(F/7$$\"3upr'fwtl7*F/$\"3?[wjRlIwhF/7$$\"3!QRa&4L &f/\"!#<$\"3;Z(yFx$R#4&F/7$$\"3YjXwg<#)y6FS$\"3mW08-T6&*QF/7$$\"36qGW% H$\\:8FS$\"33c%z&o>cwDF/7$$\"3$*)e)pbO(eV\"FS$\"3;'\\9KC>.P\"F/7$$\"39 6:YIMRr:FS$!39'[bB$y2fg!#@7$$\"3!=3SCmpuq\"FS$!3Gu*=9!3Kw8F/7$$\"3LmSE EVgQ=FS$!3%4K%f*eUSm#F/7$$\"3S`:%R;(od>FS$!37rwM7Cv*y$F/7$$\"3#=uye<)G *4#FS$!3u(=G%*4%Ra]F/7$$\"3Ua[#o!GC>AFS$!3W$Rg()>%[YgF/7$$\"3-YwJf&y(e BFS$!3CmyR[T;*3(F/7$$\"3oNrpV8H#[#FS$!3)o^]\\H?k*yF/7$$\"3EzY)QW/yh#FS $!3KN8hedSX')F/7$$\"35v)>8'\\%ou#FS$!3o#zV)pb3;#*F/7$$\"3w$GCS?&[\")GF S$!3Ed^tM_0`'*F/7$$\"3*Q3$\\!41L%HFS$!3G57e]'=mz*F/7$$\"3Z%)='p(p70IFS $!3MVUpT_$G!**F/7$$\"3g`,ta\"4=2$FS$!3[6!R-NdV(**F/7$$\"3sA%)\\K8\\QJF S$!3kq&=&*))[*****F/7$$\"3#*3]^u`v2KFS$!3!)e!)*p$y$p(**F/7$$\"3c&fJlT> qF$FS$!3')[DAi_H/**F/7$$\"3e!=@(oRJPLFS$!3y&RFq'ex,)*F/7$$\"39l2\"4_3w R$FS$!3)3chXd!pj'*F/7$$\"3MOnGd![y_$FS$!35UH]Pb`[#*F/7$$\"3#*=4JU#)RiO FS$!31=UV%z.-m)F/7$$\"3%G,f%[$HSz$FS$!3s:%[FS$\"3tp+iwQn68F/7$$\"3Ou\"4%z!e2(\\FS$ \"3Y6o#)3V&og#F/7$$\"3*zX#[4&eg5&FS$\"3SVp0*[%)4\"RF/7$$\"3n*e![/*ojB& FS$\"3Oh3H.r&\\3&F/7$$\"3#ob8X5I'p`FS$\"3_K0%RCs%phF/7$$\"3PA\"GX$yy,b FS$\"3?z$z;(Ql'4(F/7$$\"39L1/jqABcFS$\"3CD@rjpV*z(F/7$$\"3mweKB,TidFS$ \"3)*pohPg!*\\%)F/7$$\"3'oIe;6(*o)eFS$\"3iO^m;l'*=*)F/7$$\"3nKcu0ii>gF S$\"3)*RJ*H%fIJ$*F/7$$\"3w)*zj%)[mYhFS$\"3sW)pFGN(o'*F/7$$\"3)****>YH& =$G'FS$\"21dL+(********FS-%&TITLEG6#QEcubic~spline~interpolation~of~co s(x)6\"-%'COLOURG6&%$RGBGF(F($\"*++++\"!\")-%+AXESLABELSG6$Q\"xFb\\lQ! 6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F($\"+3`=$G'!\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 105 "plot(cos(x)-ex2,x=0..2*Pi,thicknes s=3,color=blue, title=\"error in cubic spline interpolation of cos(x) \");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVE SG6#7es7$$\"\"!F)F(7$$\"3Hjqb\"[W>r\"!#>$\"3MGdED33RR!#?7$$\"3eET6j*)) QU$F-$\"3+?K0`@p\"f(F07$$\"3_*=rYWLe8&F-$\"3YF'Q&f=X'4\"F-7$$\"3;`#Gi# zxZoF-$\"3slQHgIV19F-7$$\"3!oJ&y2Csf&)F-$\"3I`C@S%)z*o\"F-7$$\"3\"zBM* )omr-\"!#=$\"39M]q!>Us%>F-7$$\"3EW**3P6O)>\"FH$\"3sUajgqYz@F-7$$\"3i]c C&eb&p8FH$\"3C%y'\\UV=(Q#F-7$$\"3_IF=LjYn;FH$\"3)QSHn[kBp#F-7$$\"3q5)> 63x`'>FH$\"3Uirt1\"o$HHF-7$$\"3e!*o0HyGjAFH$\"3[+KX0v7-JF-7$$\"3YqR*pd )>hDFH$\"3\"G.u]A_Y@$F-7$$\"3\\LyxtWrGFFH$\"3)=`>dL(*HD$F-7$$\"3'fph0P Ii*GFH$\"3#[rTF'\\KuKF-7$$\"3WebMniujIFH$\"3eYnP$zz$zKF-7$$\"3Y@%HT;i7 B$FH$\"3'>0n519*oKF-7$$\"3'p9(pdRHmNFH$\"3\"G&3%zLnW?$F-7$$\"3Zs[E^dK, RFH$\"3OQozg?:(3$F-7$$\"3CVu.`$y&QUFH$\"3G_)GA$4AF-7$$\"3]9F-7$$\"3Ls#3kbN;#fFH$ \"3B\"GL4#Q!Ri\"F-7$$\"3;I[Vb_GdiFH$\"3p'*R1/nC=8F-7$$\"35*QhW&\\$Hf'F H$\"3:7eSmTR85F-7$$\"3J-7h\">DT!pFH$\"35hI0?ebxtF07$$\"3]:5wGaJ:sFH$\" 3=D=ojHA[ZF07$$\"3fF3\"fm0l_(FH$\"3t>wXW!o[I#F07$$\"3zS11.fpPyFH$\"3-W eaM)*>r5!#@7$$\"3#e!R^M[8#[)FH$!3AyC]dC99NF07$$\"3upr'fwtl7*FH$!3sL*oU zx])fF07$$\"3!Rbb2V`Iz*FH$!3QIr-mac5vF07$$\"3!QRa&4L&f/\"!#<$!3mood>(> A;)F07$$\"3/lc&f'=ci5F`v$!37-yDL%=!3#)F07$$\"31OpNA/a4+E7 #)F07$$\"3G2#e(y*yd4\"F`v$!3?#*pa6JNx\")F07$$\"3_y%f^`(Q76F`v$!3$z?5][ qc5)F07$$\"3+@?'zk/c9\"F`v$!3yzF_lhHhyF07$$\"3YjXwg<#)y6F`v$!3#3XogUVu \\(F07$$\"3!or.w_drC\"F`v$!3e$GZ.Ny$[kF07$$\"36qGW%H$\\:8F`v$!3a$Hl__m v6&F07$$\"3$*)e)pbO(eV\"F`v$!3e%3V+.R#=DF07$$\"396:YIMRr:F`v$\"3#o@(>- @:/))!#B7$$\"3Z'z]kaJ%R;F`v$\"37'R83*)3?W)F]u7$$\"3!=3SCmpuq\"F`v$\"3# *>`))>y%QQ\"F07$$\"3%z2'RGLDSF`v$\"3Q>Z2*)>Nk;F07$$\"3#=uye<)G*4#F`v$\"3B1E3Tk)y?\"F07$$\"3Ua[#o !GC>AF`v$\"3E<6aq+w%*pF]u7$$\"3-YwJf&y(eBF`v$!3CzT$GSH;_\"!#A7$$\"3oNr pV8H#[#F`v$!3!R_c$[e9FF`v$!3\"=M?*f0[6:F07$$\"35v)>8' \\%ou#F`v$!3SHgtz!o][\"F07$$\"3Vz?n#3lT\"GF`v$!3^d[va/RG8F07$$\"3w$GCS ?&[\")GF`v$!3a+^w))e#o0\"F07$$\"3Z%)='p(p70IF`v$!35LJ$o,pZ>%F]u7$$\"3g `,ta\"4=2$F`v$!3dJh#\\cXQI\"F]u7$$\"3sA%)\\K8\\QJF`v$!3@6$GjTd&>I!#C7$ $\"3#er1NNBJ<$F`v$!3?'o/w')[3#HF]]l7$$\"3#*3]^u`v2KF`v$!3QQN,M!=E=\"F] u7$$\"3--L_&R(QUKF`v$!3s**)*QV;r9DF]u7$$\"3c&fJlT>qF$F`v$!3mp#\\A#fwUT F]u7$$\"39l2\"4_3wR$F`v$!3&zU_85\"\\P5F07$$\"3_]()4*GGFY$F`v$!3*)p^S7* RtI\"F07$$\"3MOnGd![y_$F`v$!3kZ$eJQqEZ\"F07$$\"3?#yUNg&[hNF`v$!3!RMm3@ 3&3:F07$$\"3iF))z\\J7&f$F`v$!3M\"4+Y+O\"3:F07$$\"31t[0'pg(GOF`v$!3A#Gk ,:N6Z\"F07$$\"3#*=4JU#)RiOF`v$!3=/10ZwB)R\"F07$$\"36m\\Q&z8#GPF`v$!3C \\i!HI'*)f6F07$$\"3%G,f%[$HSz$F`v$!3m//<21)R=)F]u7$$\"3KoE+7#*Q@RF`v$! 3I(>P\\pY3J$F]]l7$$\"3y(f0ii+G1%F`v$\"3i:,0Y*4E%pF]u7$$\"3;fORr]')*=%F `v$\"3\\nL7&z(*4B\"F07$$\"3R&p%*z[LbK%F`v$\"3i:v#*zNVk;F07$$\"3'epgGO, qQ%F`v$\"3k<&=hl$yp:%[F`v$!3a/!GH_s]R#F07$$\"3Ou\"4%z!e2(\\F`v$!3)>uQha18=&F0 7$$\"3t:eW%H3%Q]F`v$!3;?WMxF]!\\'F07$$\"3*zX#[4&eg5&F`v$!3]Vw\\\"R?'=v F07$$\"3Y!*>B3hjQ^F`v$!3s7b$[L%*4(yF07$$\"3#Q_\")pq87<&F`v$!3f@BFld*y5 )F07$$\"3]!Hcj]-v=&F`v$!3gktw!Rjw<)F07$$\"3>d5t08z._F`v$!3fDD)opW?@)F0 7$$\"3(Q#e50,3?_F`v$!3[b.[$\\,)3#)F07$$\"3n*e![/*ojB&F`v$!3GZa3uBml\") F07$$\"3Dtq\\/&**HI&F`v$!3rk68ri\\CvF07$$\"3#ob8X5I'p`F`v$!3-WIx(\\$>7 gF07$$\"3fR3_p*3dV&F`v$!3%oXAQed&zMF07$$\"3PA\"GX$yy,bF`v$\"3_`qbvAVYE F]u7$$\"31]ilTw9KbF`v$\"3$o(y\\$*)))fU#F07$$\"3wxVy[u]ibF`v$\"3dd#p(yH 9=[F07$$\"3W0D\"fDnGf&F`v$\"3=)Q)>R2w&Q(F07$$\"39L1/jqABcF`v$\"3tID&fc 8u+\"F-7$$\"3CW>6GG-ecF`v$\"3\"QGEL![TB8F-7$$\"3WaK=$f=Gp&F`v$\"3$Hu$Q F<:S;F-7$$\"3blXDeVhFdF`v$\"3?0/s%G9(\\>F-7$$\"3mweKB,TidF`v$\"3bU^QvT JWAF-7$$\"3w$)*3/(=`$z&F`v$\"3O(RKDmk))[#F-7$$\"3w\"4#\\vXOvd&eF`v$\"3yl*z,4_@!HF-7$$\"3'oIe;6(*o)eF`v$\"3 3*=<;)QMgIF-7$$\"3(y8!=&Qz+#fF`v$\"3s=?@Wzr&=$F-7$$\"3yp>qe;E`fF`v$\"3 -O8\\ySTgKF-7$$\"3;')GY&z_)pfF`v$\"3aia'G!R(oF$F-7$$\"3o,QAKRW')fF`v$ \"3?[?S2xTyKF-7$$\"3<gF`v$\" 3([>X$G*GQB$F-7$$\"3)Rsoa(eQ^gF`v$\"3]&[OX*f=FJF-7$$\"3G:=>Xb9$3'F`v$ \"3#GHSM/KH&HF-7$$\"3Y2\\\"\\@0\\6'F`v$\"3:?V(o'\\>1FF-7$$\"3w)*zj%)[m YhF`v$\"3O4l\"QN-AQ#F-7$$\"3O\\d)e$*HP;'F`v$\"3G8!eqTrX<#F-7$$\"31*\\L r)\\z!='F`v$\"3Q2<#\\csD%>F-7$$\"3w[7QQ+'y>'F`v$\"3E.X%zV,bo\"F-7$$\"3 P***G'*3D\\@'F`v$\"3`'=FD6iES\"F-7$$\"3'*\\n(39!*>B'F`v$\"35I]Z9rO$4\" F-7$$\"3o*\\C@>b!\\iF`v$\"353')4,VPpvF07$$\"3Q\\APV-7miF`v$\"3Z-,-\\^- FRF07$$\"3)****>YH&=$G'F`v$\"31WBxAGk'*H!#E-%&TITLEG6#QNerror~in~cubic ~spline~interpolation~of~cos(x)6\"-%'COLOURG6&%$RGBGF(F($\"*++++\"!\") -%+AXESLABELSG6$Q\"xFifmQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F($\"+3`=$ G'!\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 227 "Note because the second derivative of sin(x) vanish es at 0 and at Pi we expect the natural cubic spline for sin(x) \+ to have somewhat smaller error than the natural cubic spline for cos( x). That is what we observe above." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 261 55 "Example: Interpolation Polynom ial for sin(x) and cos(x)" }{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "ex1p:=interp(XX1,YY1 ,x):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 89 "plot(ex1p,x=0..2*Pi ,thickness=3,color=blue, title=\"interpolation polynomial for sin(x)\" );" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG 6#7en7$$\"\"!F)F(7$$\"3i]cC&eb&p8!#=$\"3I*)yI-alv8F-7$$\"3YqR*pd)>hDF- $\"3E4,;2ZIXDF-7$$\"3Zs[E^dK,RF-$\"3P/P)3#f$H\"QF-7$$\"3ab^NehL]_F-$\" 3(==IaK>&=]F-7$$\"35*QhW&\\$Hf'F-$\"3z$Hlv$G6GhF-7$$\"3zS11.fpPyF-$\"3 ub0=H!p&fqF-7$$\"3upr'fwtl7*F-$\"31zJR*3y(4zF-7$$\"3!QRa&4L&f/\"!#<$\" 3KgR(*[d'=l)F-7$$\"3YjXwg<#)y6FQ$\"3k\\l2zlaR#*F-7$$\"36qGW%H$\\:8FQ$ \"3W%paQUhWn*F-7$$\"3SH22vMov8FQ$\"3clcQC+>4)*F-7$$\"3$*)e)pbO(eV\"FQ$ \"3/%zQ1H(R3**F-7$$\"3/]+3VNj.:FQ$\"33140#*H6x**F-7$$\"396:YIMRr:FQ$\" 369on/,++5FQ7$$\"3Z'z]kaJ%R;FQ$\"3!*fVctAvw**F-7$$\"3!=3SCmpuq\"FQ$\"3 m5D')HgF2**F-7$$\"33u?N%*p.tFQ$\"3=PHpkgmh#*F-7$$\"3#=uye<)G*4#FQ$\"3l> A]]fHO')F-7$$\"3Ua[#o!GC>AFQ$\"39Rj?l:dqzF-7$$\"3-YwJf&y(eBFQ$\"3!y_.' pew_qF-7$$\"3oNrpV8H#[#FQ$\"3QFQ]6TPDhF-7$$\"3EzY)QW/yh#FQ$\"3h')4q+y< ,]F-7$$\"35v)>8'\\%ou#FQ$\"3!y;$[B.CXQF-7$$\"3w$GCS?&[\")GFQ$\"3K'))Gw 4:9d#F-7$$\"3Z%)='p(p70IFQ$\"3!407`Rr,O\"F-7$$\"3sA%)\\K8\\QJFQ$\"3D$> @*3Rr+J!#?7$$\"3c&fJlT>qF$FQ$!3YI\\Ww)z)\\8F-7$$\"39l2\"4_3wR$FQ$!35!) =3c'e=`#F-7$$\"3MOnGd![y_$FQ$!3eKq(*f@rmPF-7$$\"3#*=4JU#)RiOFQ$!318fo7 bKv\\F-7$$\"3%G,f%[$HSz$FQ$!3K\\]z5$e42'F-7$$\"3KoE+7#*Q@RFQ$!3M\\u'*y ULJqF-7$$\"3y(f0ii+G1%FQ$!3]@!3&=smjzF-7$$\"3;fORr]')*=%FQ$!3='3u.kXim )F-7$$\"3R&p%*z[LbK%FQ$!3Qfc>S%*zh#*F-7$$\"3M'pExBp%[WFQ$!3)R$[Dap_a'* F-7$$\"3,z&[2')pc^%FQ$!3Cn$f!Gy!y!)*F-7$$\"3oh/x$[qGe%FQ$!36/5H'\\\\n \"**F-7$$\"3qQq'H.,hk%FQ$!3-(=`FUD$y**F-7$$\"3e;O;#eJ$4ZFQ$!3+t642x'** ***F-7$$\"3sU'G^sDax%FQ$!3=gIHt#>)z**F-7$$\"3&)oO4o)>:%[FQ$!3RCsH4]0;* *F-7$$\"3;@9vt*Qh!\\FQ$!32,!e&R5$=\")*F-7$$\"3Ou\"4%z!e2(\\FQ$!3-\\w$p Lcmm*F-7$$\"3*zX#[4&eg5&FQ$!3]fZ2>I0L#*F-7$$\"3n*e![/*ojB&FQ$!3Mt*=u3' =c')F-7$$\"3#ob8X5I'p`FQ$!3$onT_0f_\"zF-7$$\"3PA\"GX$yy,bFQ$!3lQQJ`ryU qF-7$$\"39L1/jqABcFQ$!3%Qxzm0]L8'F-7$$\"3mweKB,TidFQ$!3-)ygFQ$!3ll3[0.;YH&=$G'FQ$!3;EkV*>#Qq7!#D-%&TI TLEG6#QDinterpolation~polynomial~for~sin(x)6\"-%'COLOURG6&%$RGBGF(F($ \"*++++\"!\")-%+AXESLABELSG6$Q\"xFd]lQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG 6$;F($\"+3`=$G'!\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 105 "plot(sin(x)-ex1p,x=0..2*Pi,thickness=3,color=blue, t itle=\"error in interpolation polynomial for sin(x)\");" }}{PARA 13 " " 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7fq7$$\"\"!F)F (7$$\"3Hjqb\"[W>r\"!#>$!3)>NO6'))Q+?!#@7$$\"3eET6j*))QU$F-$!3WpzD(4S,x $F07$$\"3_*=rYWLe8&F-$!3`e26rQ-D`F07$$\"3;`#Gi#zxZoF-$!3G^zg'Gf+o'F07$ $\"3\"zBM*)omr-\"!#=$!3_U(***H%eq%))F07$$\"3i]cC&eb&p8FC$!3h3v$REAx.\" !#?7$$\"3_IF=LjYn;FC$!3J!ofES=l7\"FK7$$\"3q5)>63x`'>FC$!3E31p^g0!=\"FK 7$$\"3x]$)3bCL9@FC$!3Y-(z^>S_>\"FK7$$\"3e!*o0HyGjAFC$!3`RY)G:7N?\"FK7$ $\"3QIa-.KC7CFC$!3wj(3m*HW07FK7$$\"3YqR*pd)>hDFC$!3(RrKTGv:?\"FK7$$\"3 Y@%HT;i7B$FC$!3EbG\\+$\\c7\"FK7$$\"3Zs[E^dK,RFC$!3NAb&4r.c#)*F07$$\"3+ 9+\"[&4$ed%FC$!3&Hc;dKgQ-)F07$$\"3ab^NehL]_FC$!3/bNG^cR)4'F07$$\"3Ls#3 kbN;#fFC$!3l[u:&=Z/B%F07$$\"35*QhW&\\$Hf'FC$!3bi$y\"F07$$\"3!Rbb2V`Iz*FC$\"3US$G2KGW'>F07$$\"3gu\\JE$H E,\"!#<$\"3'eIDV&Qo)3#F07$$\"3!QRa&4L&f/\"Fes$\"3ej/]M'y.;#F07$$\"3/lc &f'=ci5Fes$\"3)3.7)ph$z<#F07$$\"31OpNA//K%HG !#C7$$\"3Z'z]kaJ%R;Fes$!3/7y#)HT&*pHFbr7$$\"3!=3SCmpuq\"Fes$!3/[q6!*GP (G&Fbr7$$\"33u?N%*p.t=QS;3J9)Fbr7$$\"3S`:%R;(od>Fes$!3#=aO[SJ ^*zFbr7$$\"3^Z,\"*pw[G?Fes$!3-_;3/p6'G(Fbr7$$\"3#=uye<)G*4#Fes$!3ymos? )[&3hFbr7$$\"3Ua[#o!GC>AFes$!3Mr.MkE)\\T$Fbr7$$\"3-YwJf&y(eBFes$\"3\\w JswNPEiF_w7$$\"3oNrpV8H#[#Fes$\"3&p7Z'**e9LFFbr7$$\"3EzY)QW/yh#Fes$\"3 HC#)>S1dEYFbr7$$\"35v)>8'\\%ou#Fes$\"32'H%f9Y]E^Fbr7$$\"3w$GCS?&[\")GF es$\"3g,LWKG)yG%Fbr7$$\"3Z%)='p(p70IFes$\"3%Qk*)4RJn`#Fbr7$$\"3sA%)\\K 8\\QJFes$\"3**H*zDh2@-'F_w7$$\"3c&fJlT>qF$Fes$!3-/2e0h5>DFbr7$$\"39l2 \"4_3wR$Fes$!3-5Sf+?-VUFbr7$$\"3MOnGd![y_$Fes$!3Hd1o/v18^Fbr7$$\"3#*=4 JU#)RiOFes$!3!yW67xfDS\" )Fbr7$$\"3M'pExBp%[WFes$\"3=NE*3(*HGy(Fbr7$$\"3,z&[2')pc^%Fes$\"3S+st- G1nnFbr7$$\"3oh/x$[qGe%Fes$\"3](=8iFe^2&Fbr7$$\"3e;O;#eJ$4ZFes$\"3WB8G ltaW9F\\r7$$\"3&)oO4o)>:%[Fes$!3w@deH1\\ooFbr7$$\"3Ou\"4%z!e2(\\Fes$!3 Ag\"fHf>LU\"F07$$\"3t:eW%H3%Q]Fes$!35p:e*[2Ew\"F07$$\"3*zX#[4&eg5&Fes$ !3B3hjQ^Fes$!3K!4U:7#)36#F07$$\"3#Q_\")pq87<&Fes $!3g*zMPGld;#F07$$\"3]!Hcj]-v=&Fes$!3Kf'pY]L)z@F07$$\"3>d5t08z._Fes$!3 !)3@j/-@%=#F07$$\"3(Q#e50,3?_Fes$!3Q`Jh,lJy@F07$$\"3n*e![/*ojB&Fes$!39 !fGP\"*y:;#F07$$\"3-J))[/Uop_Fes$!3y%3e&[!474#F07$$\"3Dtq\\/&**HI&Fes$ !3Q:p=JaUo>F07$$\"3[:`]/[JO`Fes$!3i%=v\"[c!))y\"F07$$\"3#ob8X5I'p`Fes$ !3/[GzlkI[:F07$$\"3fR3_p*3dV&Fes$!3+i_ZIMRy()Fbr7$$\"3PA\"GX$yy,bFes$ \"3?r)H4Nb4]'F\\r7$$\"3wxVy[u]ibFes$\"3))HwH2/Dw6F07$$\"39L1/jqABcFes$ \"3*3mBkO59^#F07$$\"3WaK=$f=Gp&Fes$\"3!o(HiOWryUF07$$\"3mweKB,TidFes$ \"32z,#\\z=+A'F07$$\"3w\"4#\\'e&*Rw3s'*F07$$\"3yp>qe;E`fFes$\"3H!pI_$>g86FK7$$\"3nKcu0ii >gFes$\"3)*='4d?Mw>\"FK7$$\"3)y<21/1b.'Fes$\"3l@WgL![W?\"FK7$$\"3)Rsoa (eQ^gFes$\"3J2^6O(G\\?\"FK7$$\"3e'\\**GzD$fgFes$\"3WhN(pL%f-7FK7$$\"31 q-L5dEngFes$\"3k$GftRK%)>\"FK7$$\"3oU5wFc?vgFes$\"3Y.0X/nN#>\"FK7$$\"3 G:=>Xb9$3'Fes$\"3c7)z;vyU=\"FK7$$\"3Y2\\\"\\@0\\6'Fes$\"38FS%*3&)3I6FK 7$$\"3w)*zj%)[mYhFes$\"3NEzjDj8O5FK7$$\"31*\\Lr)\\z!='Fes$\"3gB_BL'p&H ))F07$$\"3P***G'*3D\\@'Fes$\"3Qra&[VKSm'F07$$\"3'*\\n(39!*>B'Fes$\"34f 5O)fq6J&F07$$\"3o*\\C@>b!\\iFes$\"3ErcsN**efPF07$$\"3Q\\APV-7miFes$\"3 )fTgf<3W*>F07$$\"3)****>YH&=$G'Fes$\"3iE7!)G;+P9!#F-%&TITLEG6#QMerror~ in~interpolation~polynomial~for~sin(x)6\"-%'COLOURG6&%$RGBGF(F($\"*+++ +\"!\")-%+AXESLABELSG6$Q\"xF^]mQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F($ \"+3`=$G'!\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 " ex2p:=interp(XX1,YY2,x):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 89 "plot(ex2p,x=0..2*Pi,thickness=3,color=blue, title=\"interpolation pol ynomial for cos(x)\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7en7$$\"\"!F)$\"\"\"F)7$$\"3;`#Gi#zxZo!#>$ \"3a:\"*[U\"3(y**!#=7$$\"3i]cC&eb&p8F2$\"3Ewc6Y\"='4**F27$$\"3q5)>63x` '>F2$\"3\\+nM\"Q76\")*F27$$\"3YqR*pd)>hDF2$\"3#[QF27$$\"36qGW%H$\\:8Fhn$\"3@d(eBb?^_#F27$$\"3$*)e) pbO(eV\"Fhn$\"3:CiPw)4]M\"F27$$\"396:YIMRr:Fhn$!3c&*zj]tdqf!#@7$$\"3!= 3SCmpuq\"Fhn$!3#Rp;H//CO\"F27$$\"3LmSEEVgQ=Fhn$!3Kt%Q*op2YEF27$$\"3S`: %R;(od>Fhn$!3?w\\_#G6Ix$F27$$\"3#=uye<)G*4#Fhn$!3C[\"z@h\\A/&F27$$\"3U a[#o!GC>AFhn$!3DU)pr)pXRgF27$$\"3-YwJf&y(eBFhn$!33'o0L\"oJ*3(F27$$\"3o NrpV8H#[#Fhn$!3KI`LDFA/zF27$$\"3EzY)QW/yh#Fhn$!3CLbva;Lf')F27$$\"35v)> 8'\\%ou#Fhn$!3MJ:lvq&4B*F27$$\"3w$GCS?&[\")GFhn$!3Q^*G#p\\jj'*F27$$\"3 *Q3$\\!41L%HFhn$!354QT=E1/)*F27$$\"3Z%)='p(p70IFhn$!38U!*HsN.2**F27$$ \"3g`,ta\"4=2$Fhn$!3-xDrd@mv**F27$$\"3sA%)\\K8\\QJFhn$!3iXU+6>&*****F2 7$$\"3#*3]^u`v2KFhn$!3sb)=2K@\"y**F27$$\"3c&fJlT>qF$Fhn$!3!etX!Q:W3**F 27$$\"3e!=@(oRJPLFhn$!3&>:1TY$34)*F27$$\"39l2\"4_3wR$Fhn$!3E'='[im2u'* F27$$\"3MOnGd![y_$Fhn$!3QkgrkDGj#*F27$$\"3#*=4JU#)RiOFhn$!3SGOf86@u')F 27$$\"3%G,f%[$HSz$Fhn$!3]6\"HS-(3YzF27$$\"3KoE+7#*Q@RFhn$!3m&3;y2n06(F 27$$\"3y(f0ii+G1%Fhn$!3![7-]Na&[gF27$$\"3;fORr]')*=%Fhn$!3MCW+i7i!*\\F 27$$\"3R&p%*z[LbK%Fhn$!39c!Q:&RosPF27$$\"3M'pExBp%[WFhn$!3-9!p:tg&3EF2 7$$\"3oh/x$[qGe%Fhn$!3e\\:%[Fhn$\"3YgFhn$\"3 #)[E=!o$Ge'*F27$$\"3G:=>Xb9$3'Fhn$\"3PJyG^@A/)*F27$$\"3w)*zj%)[mYhFhn$ \"3e_\"\\wz1-\"**F27$$\"3P***G'*3D\\@'Fhn$\"3IL5w!y^)y**F27$$\"3)****> YH&=$G'FhnF*-%&TITLEG6#QDinterpolation~polynomial~for~cos(x)6\"-%'COLO URG6&%$RGBGF(F($\"*++++\"!\")-%+AXESLABELSG6$Q\"xFe]lQ!6\"-%*THICKNESS G6#\"\"$-%%VIEWG6$;F($\"+3`=$G'!\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 105 "plot(cos(x)-ex2p,x=0..2*Pi,thickness=3,color=blue, title=\"error in interpolation polynomial for cos(x)\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7es7$$\"\"!F) F(7$$\"3W;`y2Csf&)!#?$!3-OIT1j-uL!#A7$$\"3Hjqb\"[W>r\"!#>$!3-)p+Y\"e1M lF07$$\"3w%fNBs;zc#F4$!3g\\CV)zL%)[*F07$$\"3eET6j*))QU$F4$!3]'RYm:@XA \"!#@7$$\"3SeE*Q?h)zUF4$!3%p%GQS^A\"[\"FA7$$\"3_*=rYWLe8&F4$!3(yzYO9=( >/T2%>FA7$$\"3;`#Gi#zxZoF4$!3xA'Q[C8]9#F A7$$\"3!oJ&y2Csf&)F4$!3g!)>f['yg]#FA7$$\"3\"zBM*)omr-\"!#=$!3-9G&>vd#3 GFA7$$\"3EW**3P6O)>\"Fhn$!3I*z\")[Tzl0$FA7$$\"3i]cC&eb&p8Fhn$!3Kl&3kvr dD$FA7$$\"3_IF=LjYn;Fhn$!3Q$y;_Z<$)\\$FA7$$\"3q5)>63x`'>Fhn$!3HDYsjq%o i$FA7$$\"3t!3/\"oZ&)R?Fhn$!3n@#oVraMk$FA7$$\"3x]$)3bCL9@Fhn$!3_2[%3!*> Wl$FA7$$\"3#3is?95))=#Fhn$!3+g#)yp=,gOFA7$$\"3e!*o0HyGjAFhn$!3%4'*zc(= \\gOFA7$$\"3Mg6/;bwPBFhn$!3Ma:k\\L6cOFA7$$\"3QIa-.KC7CFhn$!3y#oJg1Ark$ FA7$$\"3U+(4+*3s'[#Fhn$!3='GhDFhn$!31(*yQ>qC;OF A7$$\"3'fph0PIi*GFhn$!3))3W__;2#\\$FA7$$\"3Y@%HT;i7B$Fhn$!3,GPl[*HqI$F A7$$\"3'p9(pdRHmNFhn$!3%3:RWu,s2$FA7$$\"3Zs[E^dK,RFhn$!3pU\"QJsAj\"GFA 7$$\"3CVu.`$y&QUFhn$!3Ftz84C7MDFA7$$\"3+9+\"[&4$ed%Fhn$!3>5&y5)RFA7$$\"3ab^NehL]_Fhn$!3!Q`dZfN.m\"FA 7$$\"3]9#\\AN@,_'F07 $$\"3]:5wGaJ:sFhn$!37O,#oItQ\"HF07$$\"3zS11.fpPyFhn$!3'Hp1oZ@'[k!#C7$$ \"3#e!R^M[8#[)Fhn$\"3Tp%F07$$\"3/lc&f'=ci5Fjx$\"3)>nXOmX?p%F07$$\"31OpNA/YF07$$\"3_y%f^`(Q76Fjx$ \"3_Ge`a!e9b%F07$$\"3+@?'zk/c9\"Fjx$\"3[kC-f@oeVF07$$\"3YjXwg<#)y6Fjx$ \"3Cq\\+%epI5%F07$$\"3!or.w_drC\"Fjx$\"3y!Q5_hvYV$F07$$\"36qGW%H$\\:8F jx$\"3A\"z./o4vl#F07$$\"3$*)e)pbO(eV\"Fjx$\"3I)f3(y`wp7F07$$\"396:YIMR r:Fjx$!31HLv3Ig*e%!#D7$$\"3Z'z]kaJ%R;Fjx$!3UA4Zx$)>2Y!#B7$$\"3!=3SCmpu q\"Fjx$!3nL))ei..GyF_]l7$$\"33u?N%*p.tFjx$!3%H\\_*zJ*4w*F_]l7$$\"3#=uye<)G*4#Fjx$!3UsC\"H%RGilF_]l7$$\" 3Ua[#o!GC>AFjx$!3]m=b:?:XKF_]l7$$\"3-YwJf&y(eBFjx$\"3`X^-eJo>]Fi\\l7$$ \"3oNrpV8H#[#Fjx$\"3K+K8ph@b=F_]l7$$\"3EzY)QW/yh#Fjx$\"3u=I_!3$G%\\#F_ ]l7$$\"35v)>8'\\%ou#Fjx$\"33M$p9%4]#3#F_]l7$$\"3w$GCS?&[\")GFjx$\"3%G? IpRgv9\"F_]l7$$\"3Z%)='p(p70IFjx$\"3+(yU@D\"ehNF]w7$$\"3sA%)\\K8\\QJFj x$\"3(ygexwG\"G>!#F7$$\"3c&fJlT>qF$Fjx$\"3?\\!)=jP!*4NF]w7$$\"39l2\"4_ 3wR$Fjx$\"3zs^%H&*)o<6F_]l7$$\"3MOnGd![y_$Fjx$\"3WY$)p^oLK?F_]l7$$\"3# *=4JU#)RiOFjx$\"3q;)G2@pU\\#F_]l7$$\"3%G,f%[$HSz$Fjx$\"3z%[_zxq.#>F_]l 7$$\"3KoE+7#*Q@RFjx$\"3u_#*pYe-#3\"F]w7$$\"3y(f0ii+G1%Fjx$!3j5Wie;S9KF _]l7$$\"3;fORr]')*=%Fjx$!3*)3QM4fNBnF_]l7$$\"3R&p%*z[LbK%Fjx$!3KDzYE%[ :w*F_]l7$$\"3'epgGO,qQ%Fjx$!3g_'y_31d/\"F07$$\"3M'pExBp%[WFjx$!3Cc:\\I q3;\")elf*F_]l7$$\"3oh/x$[qGe%Fjx$!3E3 jck\"y:b(F_]l7$$\"3qQq'H.,hk%Fjx$!3#pX:Niu@Z%F_]l7$$\"3e;O;#eJ$4ZFjx$! 3%pR)o^'Q%RBF]w7$$\"3&)oO4o)>:%[Fjx$\"3L:qF>E.27F07$$\"3Ou\"4%z!e2(\\F jx$\"3I\"[<*[cY$p#F07$$\"3t:eW%H3%Q]Fjx$\"3EL`<\"[I.Y$F07$$\"3*zX#[4&e g5&Fjx$\"3()**z!4%HUB3hjQ^Fjx$\"3+E[K*oYeO%F07$$\"3#Q_\" )pq87<&Fjx$\"3Q$Q_e0:Mb%F07$$\"3]!Hcj]-v=&Fjx$\"3A1_8^'z-i%F07$$\"3>d5 t08z._Fjx$\"3EV&yLXHpm%F07$$\"3(Q#e50,3?_Fjx$\"3G8l%zo>;p%F07$$\"3n*e! [/*ojB&Fjx$\"3+ah#*zWe#p%F07$$\"3!3r%[al-`_Fjx$\"3'ex&3q8:nYF07$$\"3-J ))[/Uop_Fjx$\"3\\CyVRm38YF07$$\"39_H\\a=M'G&Fjx$\"3%)o86GG-ecFjx$!3r$f'\\lq?4))F07$ $\"3WaK=$f=Gp&Fjx$!3'o\"3'RZ-_8\"FA7$$\"3blXDeVhFdFjx$!3nxQ%y_[\"39FA7 $$\"3mweKB,TidFjx$!3kCuYEeE'p\"FA7$$\"3w$)*3/(=`$z&Fjx$!3ag6'za;I'>FA7 $$\"3w\"4#\\vXOvd&eFjx$!3,m%\\!ygn.DFA 7$$\"3'oIe;6(*o)eFjx$!3!Qz5Wwlfw#FA7$$\"3(y8!=&Qz+#fFjx$!3+l@MRI!)GIFA 7$$\"3yp>qe;E`fFjx$!315*)QnPajKFA7$$\"3o,QAKRW')fFjx$!3[d9)>,nsX$FA7$$ \"3nKcu0ii>gFjx$!3mNEXuQ%[f$FA7$$\"3)y<21/1b.'Fjx$!3)4*z[&>#yNOFA7$$\" 3)Rsoa(eQ^gFjx$!3]d1Kg`tdOFA7$$\"3e'\\**GzD$fgFjx$!3E`XeRH)3m$FA7$$\"3 1q-L5dEngFjx$!3CA=wM\\TeOFA7$$\"3oU5wFc?vgFjx$!3Eo**y_$>+l$FA7$$\"3G:= >Xb9$3'Fjx$!38etXpaPNOFA7$$\"3Y2\\\"\\@0\\6'Fjx$!3\"e&)=qHVv]$FA7$$\"3 w)*zj%)[mYhFjx$!3#z8'>A\"z7D$FA7$$\"3O\\d)e$*HP;'Fjx$!3EV[tOwg^IFA7$$ \"31*\\Lr)\\z!='Fjx$!3\\Hrp8e,.GFA7$$\"3w[7QQ+'y>'Fjx$!3@WiAsn\"3]#FA7 $$\"3P***G'*3D\\@'Fjx$!3T7&eypA+9#FA7$$\"3AuGD:wXBiFjx$!3O\"HY]35g$>FA 7$$\"3'*\\n(39!*>B'Fjx$!3mb!\\iFjx$!3#>ljLF07$$ \"3)****>YH&=$G'FjxF(-%&TITLEG6#QMerror~in~interpolation~polynomial~fo r~cos(x)6\"-%'COLOURG6&%$RGBGF(F($\"*++++\"!\")-%+AXESLABELSG6$Q\"xFhf mQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F($\"+3`=$G'!\"*%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 376 "I n the above examples we see the interpolation polynomial has smaller e rror than the interpolation spline. The next example will demonstrate \+ this is not always the case. One problem with interpolation polynomial s is that they are of high degree if there are many nodes. Evaluating \+ a polynomial of high degree runs a serious risk of loss of significanc e errors due to roundoff." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 262 61 "Example. A function pathological for equispaced interpolation" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 68 "We use a do-loop to build a list of 21 equispaced \+ nodes in [-1,1]" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "a:=-1: XX2:=[]:for k from 1 to 21 do XX2:=[op (XX2),a];a:=a+1/10;od:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 31 "Let's see what XX2 looks like" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "XX 2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#77!\"\"#!\"*\"#5#!\"%\"\"&#!\"(F '#!\"$F*#F$\"\"##!\"#F*#F.F'#F$F*#F$F'\"\"!#\"\"\"F'#F8F*#\"\"$F'#F0F* #F8F0#F;F*#\"\"(F'#\"\"%F*#\"\"*F'F8" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 128 "Now here's the standard exampl e of a function whose interpolation polynomials with equispaced nodes \+ yield very bad approximation" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "f:=x->1/(1+20*x^2);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fGR6#%\"xG6\"6$%)operatorG%&arrowG F(*&\"\"\"F-,&F-F-*&\"#?F-)9$\"\"#F-F-!\"\"F(F(F(" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 16 "YY3:=map(f,XX2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$YY3G77#\"\"\"\"#@#\"\"&\"#')#F*\"#p#F*\"#a#F*\"#T#F' \"\"'#F*F(#F*\"#9#F*\"\"*#F*F3F'F9F7F5F4F2F0F.F,F)F&" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 59 "Let's look at t he natural cubic spline interpolation first." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "ex3:=spline(XX2,Y Y3,x,cubic):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 83 "plot(ex3,x= -1..1,thickness=3,color=blue,title=\"cubic interpolation spline for f \");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVE SG6#7ao7$$!\"\"\"\"!$\"35cZ!>w/>w%!#>7$$!3ommm;p0k&*!#=$\"3-jRk0]\"\\> &F-7$$!3wKL$3dP`\")F-7$$!3u****\\P8#\\4(F1$\"3;)=:\"*Q \"eN!*F-7$$!3+nm;/siqmF1$\"3#*eZg`#[+,\"F17$$!3[++](y$pZiF1$\"3zjZx=*H a8\"F17$$!33LLL$yaE\"eF1$\"3Ufv\\,Q0*G\"F17$$!3hmmm\">s%HaF1$\"3)>DeRY '4]9F17$$!3Q+++]$*4)*\\F1$\"3yh?pB9sn;F17$$!39+++]_&\\c%F1$\"3UT3!*p>[ M>F17$$!31+++]1aZTF1$\"3[Qq:hex^AF17$$!3umm;/#)[oPF1$\"3>R%zPeBUg#F17$ $!3hLLL$=exJ$F1$\"3R'\\UsJ(>CJF17$$!3*RLLLtIf$HF1$\"3=Hli>owqOF17$$!3] ++]PYx\"\\#F1$\"3E*zjd9LhX%F17$$!3QnmTNz>&H#F1$\"3GU()[*H')R'[F17$$!3E MLLL7i)4#F1$\"3`wV7A&3MJ&F17$$!3#pm;aVXH)=F1$\"3'e4C(fp]deF17$$!3c**** \\P'psm\"F1$\"3we5&\\a'4YkF17$$!3s*****\\F&*=Y\"F1$\"3&R`/s.Y(HqF17$$! 3')****\\74_c7F1$\"3%3!\\,SJ?=wF17$$!3ZmmT5VBU5F1$\"3mT1S'y*Q=#)F17$$! 3)3LLL3x%z#)F-$\"3-@z;)y\"Q$y)F17$$!3gKL$e9d;J'F-$\"3fq([T92PC*F17$$!3 KMLL3s$QM%F-$\"3<**R`lTW='*F17$$!3'ym;aQdDG$F-$\"3+J+bnR#\\x*F17$$!3T, +]ivF@AF-$\"3AQ9-*GWO*)*F17$$!3=o;/^wj!p\"F-$\"3XR*3e3Ou$**F17$$!3'\\L $eRx**f6F-$\"3u?v'>u(4q**F17$$!3S<+D\"GyNH'!#?$\"3-&f**R$e1\"***F17$$! 3]^omm;zr)*!#@$\"3IwvSQpx****F17$$\"3o'H$3x\"yY_%F\\v$\"3))Ql,^%f`***F 17$$\"3&yK$3_Nl.5F-$\"3'**)e@5g^x**F17$$\"3/E$ekGR[b\"F-$\"3E/h[Uo(o%* *F17$$\"3@CL$3-Dg5#F-$\"3#\\:=NQtS!**F17$$\"3e?Le*['R3KF-$\"37h;jM$*\\ %y*F17$$\"3%pJL$ezw5VF-$\"3l%>%3\\`%Qi*F17$$\"3L`mmmJ+IiF-$\"3GmfiSr:h #*F17$$\"3s*)***\\PQ#\\\")F-$\"3-MN(fc?f\"))F17$$\"3ilm\"z\\1A-\"F1$\" 3/`![ATpIF)F17$$\"3GKLLe\"*[H7F1$\"3#)3_Na')3&p(F17$$\"3ylm;HCjV9F1$\" 31MO\\Mq8#3(F17$$\"3I*******pvxl\"F1$\"3%))=Wx')>FZ'F17$$\"3g)***\\7JF n=F1$\"3!4dTD^4!**eF17$$\"3#z****\\_qn2#F1$\"3wQbdkA2m`F17$$\"3=)**\\P /q%zAF1$\"3OgVbQ8P)*[F17$$\"3U)***\\i&p@[#F1$\"3)\\9mq[W^Z%F17$$\"3B)* ***\\2'HKHF1$\"3%4`Yw-+ln$F17$$\"3ElmmmZvOLF1$\"3)HB\"oj'>(*4$F17$$\"3 i******\\2goPF1$\"33UWiU%3Tg#F17$$\"3UKL$eR<*fTF1$\"3VlUa+>F17$$\"3ckm;H!o-*\\F1$\"39Py'z#p2s;F17$$ \"3y)***\\7k.6aF1$\"3s3eQDrbe9F17$$\"3#emmmT9C#eF1$\"3W]2_uPH&G\"F17$$ \"33****\\i!*3`iF1$\"3[2,`T9pL6F17$$\"3%QLLL$*zym'F1$\"3^=4*G/'z55F17$ $\"3wKLL3N1#4(F1$\"3NT^Dc@?U!*F-7$$\"3Nmm;HYt7vF1$\"3i3U`'RL$Q\")F-7$$ \"3Y*******p(G**yF1$\"3l#)[Nxv*)=uF-7$$\"3]mmmT6KU$)F1$\"3:Kaw+,Z+nF-7 $$\"3fKLLLbdQ()F1$\"3-`fW7$=@9'F-7$$\"3[++]i`1h\"*F1$\"3ha#4g4tki&F-7$ $\"3W++]P?Wl&*F1$\"3?8/G!=,N>&F-7$$\"\"\"F*F+-%&TITLEG6#QBcubic~interp olation~spline~for~~f6\"-%'COLOURG6&%$RGBG$F*F*Ffal$\"*++++\"!\")-%+AX ESLABELSG6$Q\"xFaalQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F(F[al%(DEFAULT G" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" } }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 97 "plot(f(x)-ex3,x=-1..1,thi ckness=2,color=blue,title=\"error in cubic interpolation spline for f \");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVE SG6#7_w7$$!\"\"\"\"!$\"3s#yDJ7:6b&!#M7$$!3%pmm\"HU,\"*)*!#=$!3;]G'pGHi X'!#A7$$!3!RLL$e%G?y*F1$!34hxB0ZUK5!#@7$$!3#om;HdNvs*F1$!3B?F`C\">B9\" F:7$$!3u****\\(oUIn*F1$!3[U1Rdvo.7F:7$$!3xLL3-)\\&='*F1$!3'HIOw0N=A\"F :7$$!3ommm;p0k&*F1$!38BY.!=A??\"F:7$$!3E++vV5Su$*F1$!3aob<'*4_4\"*F47$ $!3wKL$3oz#z,7#F47$$ !3mmmmT%p\"e()F1$\"3sT.x!p')y3$F47$$!3qlm;/'=3l)F1$\"3AOj4G2\\3NF47$$! 3&emmmwnMa)F1$\"3y*)o5\"GP,Y$F47$$!3+mm;Hp6O%)F1$\"3KM!3O^u&RIF47$$!3: mmm\"4m(G$)F1$\"3M/)=\\&4FcBF47$$!3\"QLL3i.9!zF1$!3]Xg^vg1nNF]o7$$!3\" ommT!R=0vF1$!3rIu6IQOVYF]o7$$!3u****\\P8#\\4(F1$!3!494p)zPR=F]o7$$!3+n m;/siqmF1$\"3'ymoOML_3\"F47$$!3[++](y$pZiF1$\"3GBP!e#o%3Y'F]o7$$!33LLL $yaE\"eF1$\"3XjgCsC\"z$RF]o7$$!3hmmm\">s%HaF1$\"31tCzK&eDX&F]o7$$!3]LL $3x&y8_F1$!3\\#>e7.XK#GF]o7$$!3Q+++]$*4)*\\F1$\"3Bh([t0%4G;!#C7$$!3E++ ++t_\"y%F1$\"3-w'[lHVe;$F47$$!39+++]_&\\c%F1$\"3J'QS9vN=\"fF47$$!3E+++ D%yF^%F1$\"3gZzw(=kj;'F47$$!3S++++;ggWF1$\"383V-!p\"\\=iF47$$!3_+++vZU 3WF1$\"359t<\")ynigF47$$!33+++]zCcVF1$\"3!*39:)>oFq&F47$$!3!)*******H% *=D%F1$\"3qa:DFXMMWF47$$!31+++]1aZTF1$\"3_gr@wTDSEF47$$!3umm;/#)[oPF1$ !36[e)Gc`!3EF47$$!3<++v$>BJa$F1$!3IF4^z:Gl\\F47$$!3hLLL$=exJ$F1$!3K'pP yrHDl'F47$$!3!QL$3_(H+F$F1$!3W7vi_oC7mF47$$!3)RLL3K,BA$F1$!3`\"\\dElQ^ K'F47$$!3gLLe*)GduJF1$!3([Vddh:Vt&F47$$!3!QLL$eW%o7$F1$!3/Ag2&yi'zZF47 $$!3hLL$ef(QJIF1$!3K!*Rk;o-G:F47$$!3*RLLLtIf$HF1$\"3Ol>PRx()))QF47$$!3 CnmT&o_Qr#F1$\"3=c*fFf$)RN#F:7$$!3]++]PYx\"\\#F1$\"31:z73\"=;f%F:7$$!3 An\"z>YIEW#F1$\"3OiBKD&\\q)\\F:7$$!3%RLekG'[$R#F1$\"3%[rZ#)eiyH&F:7$$! 3K9*oBF1$\"3:Ql2+m[7aF:7$$!3n+v$46UVM#F1$\"3kv\")=O\"=X\\&F:7$$!3 .%3xJ-q(>BF1$\"3o%zxK.5&RbF:7$$!3QnmTNz>&H#F1$\"3#F1$\"3G2lwCXFG]F:7$$!3EMLLL7i)4#F1$\"3#*f'3\")pp&eLF:7$$!3InT&QG -Z/#F1$\"37,v2l)z/x\"F:7$$!3g+]PMLy!*>F1$!3Z\"3o>!*4RA%F47$$!3!R$e*[Qk o$>F1$!3\")H4U@YPPKF:7$$!3#pm;aVXH)=F1$!33m:4^P'HY'F:7$$!3&**\\Pf[E!H= F1$!3K'eKP>_%p)*F:7$$!3CL$ek`2^x\"F1$!3jY/Z#H/ZK\"!#?7$$!3bm\"zpe)=@F_^l7$$!3u**\\(o/E fh\"F1$!3sfm(Gq<\\8#F_^l7$$!3k***\\iX#ek:F1$!3G)z+YY&4%H#F_^l7$$!3a** \\il)QK^\"F1$!3k)**)[!)Ho)Q#F_^l7$$!3s*****\\F&*=Y\"F1$!3Nw\\K$[eOT#F_ ^l7$$!3*)**\\P%o^0T\"F1$!3c%zVJ%3cmBF_^l7$$!3!)***\\P43#f8F1$!3w6q'\\_ $eZAF_^l7$$!3p**\\7.X'yI\"F1$!3at:%G()\\'f?F_^l7$$!3')****\\74_c7F1$!3 +L1&G_l%3=F_^l7$$!3_m\"z>E\\H?\"F1$!3i%G/7*H>)[\"F_^l7$$!3=L$e9hx$\\6F 1$!3ih*RF4r77\"F_^l7$$!3];z>'y\"fA6F1$!3k5kJ51\\_#*F:7$$!3#)*\\P4'f!e4 \"F1$!3[4Y'[NM)QsF:7$$!39$3xc8?!p5F1$!3e(>Q$p_*o>&F:7$$!3ZmmT5VBU5F1$! 3sg$4Y8+P:$F:7$$!3!)\\i:&[[a,\"F1$!3*f%*)[nF3Q6F:7$$!3IJ$e*)fEm))*!#>$ \"3u#R()Rmz\")>)F47$$!3ckTNY$o(='*Fhcl$\"39U&HzZ!o-FF:7$$!3#y**\\P454N *Fhcl$\"3v>d?9)4$4XF:7$$!3Nk;a)e$>:))Fhcl$\"3-**pV6#Rj)yF:7$$!3)3LLL3x %z#)Fhcl$\"3'QG())H9N#4\"F_^l7$$!3KJ$e*)4Avy(Fhcl$\"3v$)*ekx$eP8F_^l7$ $!3uJLe9rc&H(Fhcl$\"3)H$y$oL+ia\"F_^l7$$!3eFhcl$\"3\\(\\^zx)o0>F_^l7$$!3CL3FpYstbFhcl$ \"3*3)>+b%p6#>F_^l7$$!3XLL3xruF`Fhcl$\"3L(e/.?$*G#>F_^l7$$!3)QL3F>#zN[ Fhcl$\"3IJsm'f'Q%)=F_^l7$$!3KMLL3s$QM%Fhcl$\"3Cn61j^4!z\"F_^l7$$!34,]( oH(>8QFhcl$\"3S*p$GckiG;F_^l7$$!3'ym;aQdDG$Fhcl$\"3`5F=#He=T\"F_^l7$$! 3kM$eRZ<>v#Fhcl$\"3g>%y^^/8:\"F_^l7$$!3T,+]ivF@AFhcl$\"3m\\*3Y]=&Q')F: 7$$!3=o;/^wj!p\"Fhcl$\"3#*ye*[TGPs&F:7$$!3'\\L$eRx**f6Fhcl$\"3!>]#HRNf iIF:7$$!3S<+D\"GyNH'F_^l$\"3i56[[5h=5F:7$$!3]^omm;zr)*F:$\"3DR#p2(\\3; GF]o7$$\"3o'H$3x\"yY_%F_^l$\"3)>&z!*QI#oZ&F47$$\"3&yK$3_Nl.5Fhcl$\"3]S ]H1p3yBF:7$$\"3/E$ekGR[b\"Fhcl$\"3<(y2$>cI0]F:7$$\"3@CL$3-Dg5#Fhcl$\"3 E#>Mg'*f(**zF:7$$\"3SA$3_v5sl#Fhcl$\"3NVZ.ngR,6F_^l7$$\"3e?Le*['R3KFhc l$\"3??@eo.vx8F_^l7$$\"3w=$eRA#efPFhcl$\"3y\"Gp$HN34;F_^l7$$\"3%pJL$ez w5VFhcl$\"3#3R!QI$)z\"y\"F_^l7$$\"3/^mTgnd!z%Fhcl$\"3k!y'\\Z;.y=F_^l7$ $\"39&)**\\ibQq_Fhcl$\"3]7h?VKH@>F_^l7$$\"3=_;aj**G5bFhcl$\"3-B28\"RSH #>F_^l7$$\"3C>LekV>]dFhcl$\"3K66XS=W6>F_^l7$$\"3G')\\il()4!*fFhcl$\"3U ^p4Ma0()=F_^l7$$\"3L`mmmJ+IiFhcl$\"3G)oU!*GA,&=F_^l7$$\"3U()*\\(o>\")4 nFhcl$\"3N'*33ZJJSE9XXC&zbF:7$$\"3g:Hd&\\P**4\"F1$!3e6K#))RD>b(F:7$$\"3&* )*\\7Gy%e7\"F1$!3]iCFxm5%\\*F:7$$\"3il\"HK\\ow<\"F1$!3$R2e?2R*>8F_^l7$ $\"3GKLLe\"*[H7F1$!3qF)yX)>d`;F_^l7$$\"3kl;/w\\-$G\"F1$!39R>e:ucX>F_^l 7$$\"3.***\\PzglL\"F1$!3W+^(**)R\"H<#F_^l7$$\"3UK$e9h'4!R\"F1$!3fy/,]J nFBF_^l7$$\"3ylm;HCjV9F1$!3>w3<4?@0CF_^l7$$\"3aKLekSq]:F1$!3Za_5zwGEBF _^l7$$\"3I*******pvxl\"F1$!3i)=SJ/1D'>F_^l7$$\"37**\\7`+:5 Y7q\"F_^l7$$\"3'*)**\\iSCDw\"F1$!3j9g_2*[5S\"F_^l7$$\"3y)*\\Pf()*[\"=F 1$!3[-3SF%Gl2\"F_^l7$$\"3g)***\\7JFn=F1$!3sB%)yQq:WuF:7$$\"3W)*\\iluk> >F1$!3:TQSwy2MUF:7$$\"3G)**\\(==-s>F1$!3*>Y>fV$=S8F:7$$\"35)*\\(=<'RC? F1$\"3E_c$4y#F:7$$\"3/)*\\P% G?\"y@F1$\"3_Mp\"e>kC\"[F:7$$\"3=)**\\P/q%zAF1$\"3yBlH5mQ?bF:7$$\"3@[P f$[2[I#F1$\"3M&)3d$HQna&F:7$$\"3C)\\PM#\\9IBF1$\"3c^m\"\\TI`_&F:7$$\"3 G[7GjB[bBF1$\"37D//hDkhaF:7$$\"3I)*\\7.)>3Q#F1$\"3sE,p(3d2O&F:7$$\"3P) \\7Go%\\JCF1$\"31QP3A*ze1&F:7$$\"3U)***\\i&p@[#F1$\"3iDv[0TIuYF:7$$\"3 L)**\\(=GB2FF1$\"38w+JFx;ACF:7$$\"3B)****\\2'HKHF1$\"3qTl&zb(zQTF47$$ \"3)\\m;zu5M.$F1$!3'G$F1$!3k+U,))Q!3l'F47$$\"3ElmmmZvOLF1$!3ndU(43(z6mF47$$\"3WKL Lexn_NF1$!3)zeEE>TG'[F47$$\"3i******\\2goPF1$!3M!y1f4[pg#F47$$\"3UKL$e R<*fTF1$\"3eWh&p'QBpGF47$$\"3C**\\P40(oE%F1$\"3enP6L.wbYF47$$\"3/mm\"H iBQP%F1$\"3]28B]n;YeF47$$\"3X*\\(oz,IFWF1$\"3^+[*\\=7I9'F47$$\"3%GLekt w2[%F1$\"3uz\\__U*GA'F47$$\"3Em\"HKH`U`%F1$\"3OA7w/?[&3'F47$$\"3m***** *\\)Hxe%F1$\"3UD0CS9qSdF47$$\"37KLeR*)**)y%F1$\"3y)Q=T/Q;/$F47$$\"3ckm ;H!o-*\\F1$\"3]twVw1M*o)Fis7$$\"3nJL$3A_1?&F1$!3!QLO_k;]A$F]o7$$\"3y)* **\\7k.6aF1$\"368E#=$e\"[%[F]o7$$\"3#emmmT9C#eF1$\"3xWbW%GsHk$F]o7$$\" 33****\\i!*3`iF1$\"3k_&>%p)4Sl'F]o7$$\"3%QLLL$*zym'F1$\"3K#Qmj`694\"F4 7$$\"3wKLL3N1#4(F1$!3!Hk$30]f*z\"F]o7$$\"3Nmm;HYt7vF1$!3[s#eQR)p!p%F]o 7$$\"3Y*******p(G**yF1$!3-74WsFX>OF]o7$$\"3]mmmT6KU$)F1$\"3M4r$H,UHX#F 47$$\"3-LLeRZQT%)F1$\"3Q&e/Tu(GnIF47$$\"3a****\\P$[/a)F1$\"3Aj&\\v(G\" HX$F47$$\"32mmTN>^R')F1$\"3sX;WbW\"R_$F47$$\"3fKLLLbdQ()F1$\"3![&QKu#H O?$F47$$\"36+]i!*z>W))F1$\"31Mjq*f$zaBF47$$\"3amm\"zW?)\\*)F1$\"3*Rb&R *pB$>\"*F]o7$$\"3'HL3_!HWb!*F1$!3W9$e&z6lx6F47$$\"3[++]i`1h\"*F1$!3CgB &e#=#[w$F47$$\"3Y++++PDj$*F1$!3a:()yt1=h))F47$$\"3W++]P?Wl&*F1$!3*RnVN O`H?\"F:7$$\"3Q+D\"Gyh(>'*F1$!3X(QWhfj=A\"F:7$$\"3K+]7G:3u'*F1$!39FmoY &QH?\"F:7$$\"3G+vVt7SG(*F1$!3&)3Rp>n'49\"F:7$$\"3A++v=5s#y*F1$!35zVJ3F pI5F:7$$\"35+]P40O\"*)*F1$!30[?;Qb/SkF47$$\"\"\"F*F+-%&TITLEG6#QKerror ~in~cubic~interpolation~spline~for~~f6\"-%'COLOURG6&%$RGBG$F*F*F`in$\" *++++\"!\")-%+AXESLABELSG6$Q\"xF[inQ!6\"-%*THICKNESSG6#\"\"#-%%VIEWG6$ ;F(Fehn%(DEFAULTG" 1 2 0 1 10 2 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 45 "Now lets look at the interpolation polyno mial" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "poly:=interp(XX2,YY3,x);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%polyG,8*$)%\"xG\"\"#\"\"\"#!)Sx1]\"((eaD*&#\"0+++MS3 $e\",,+M*4>F*)F(\"\")F*F**&#\".+7#>^WL\",.e:R,\"F*)F(\"\"%F*F**&#\"0++ giAkg$\",FAS_7*F**$)F(\"\"'F*F*!\"\"*&#\"2++++++]$H\",4uY " 0 "" {MPLTEXT 1 0 82 "plot(poly,x=-1..1,thickness=3,color=blue,titl e=\"interpolation polynomial for f\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7ct7$$!\"\"\"\"!$\"3)o&Hv!y/>w% !#>7$$!3]LekynP')**!#=$!3wY'RL!4\"=`%!#<7$$!3-n;HdNvs**F1$!3]Hvt))=)*z ')F47$$!3_+v$fLI\"f**F1$!3I\"3Rcp^BC\"!#;7$$!3/MLe9r]X**F1$!3m(fd^#[xy :F?7$$!3an\"HK*Q)=$**F1$!3)>q)[>pnz=F?7$$!3/,](=ng#=**F1$!3Y-!p`[[t9#F ?7$$!3bM3_]uj/**F1$!3y9:Kcg'RQ#F?7$$!3%pmm\"HU,\"*)*F1$!3[/hu%)of\"f#F ?7$$!3M*\\7y+\"Rx)*F1$!3Wa;#zkFHF?7$ $!3OmT5lX9])*F1$!3Ec04#='pfIF?7$$!3()***\\PM@l$)*F1$!3!fIQv2A+<$F?7$$! 3)omT5!\\F4)*F1$!3(4c#=%*[wJLF?7$$!3!RLL$e%G?y*F1$!32g\"4&R,dCMF?7$$!3 Sn\"zpB0%o(*F1$!31)yPiSEXPiSW$F?7$$!3sK$3-8*G+( *F1$!3I-Xe&=n$)Q$F?7$$!3u****\\(oUIn*F1$!3s(>!GT)f\"*H$F?7$$!3wm;zWizX '*F1$!3e(o?=lMH=$F?7$$!3xLL3-)\\&='*F1$!3s#)*oCF?7$$!33;HKkf$H\\*F1$ !3MOSz\\%pHE#F?7$$!3#HL3-)*G#p%*F1$!3]1QJ=+b0@F?7$$!3w\\P4'*>_X%*F1$!3 M_Ffwm8\\>F?7$$!3gm\"z>,:=U*F1$!3_]'4([:&\\z\"F?7$$!3U$eky-3\")R*F1$!3 S=)HuEJSk\"F?7$$!3E++vV5Su$*F1$!3k!yWV$\\C(\\\"F?7$$!3&R$3_vq)pK*F1$!3 %z(>x]7\")=7F?7$$!3_m;H2Jdz#*F1$!3E)pW6!z_Q'*F47$$!33*\\i!R\"f@B*F1$!3 eCsQ#)QH[tF47$$!3wKL$3XF17$$!39***\\iId9(*)F1$\"3YcI'fi*QOdF17$$!3YL$eRP8['))F1$\"3Sj6yT,C '*=F47$$!3mmmmT%p\"e()F1$\"3k0:'*\\d.dCF47$$!3umm\"H-%\\/()F1$\"3Gy4zc h5FDF47$$!3qlm;/'=3l)F1$\"3/iW?Su'>\\#F47$$!3wlmT&=Vrf)F1$\"3kXJ(z%z7w BF47$$!3&emmmwnMa)F1$\"3tr\"QKUk7?#F47$$!3+mm;Hp6O%)F1$\"3=_Ls\\\"pou \"F47$$!3:mmm\"4m(G$)F1$\"3)fq=MNriC\"F47$$!3)****\\i&[3:\")F1$\"3%H_R 1_Jr(QF17$$!3\"QLL3i.9!zF1$!3!o+j?%fa#3\"F17$$!3e++Dch(=&yF1$!3\"RK+w \"*3?q\"F17$$!3Mnmm\"p[B!yF1$!39$)\\$GYzv9#F17$$!35ML3F7#Gv(F1$!32**yV $*e:OCF17$$!3'3++Dw$H.xF1$!3%RH1-!e#ee#F17$$!3RMLLL)QUg(F1$!3)H')R#4n( Ga#F17$$!3\"ommT!R=0vF1$!3A;JG/ank@F17$$!3tKL$3i_+I(F1$!3!R&4uNRGb')F- 7$$!3u****\\P8#\\4(F1$\"3K(zgk1%>#\\%F-7$$!3+nm;/siqmF1$\"3'\\x>J@?)G< F17$$!3[++](y$pZiF1$\"3uqs%HaF1$\"39HV00$\\PE\"F17$$!3Q+++]$*4)*\\F1$\"3C4$yCLy%o;F1 7$$!39+++]_&\\c%F1$\"3\\8K\"46g3-#F17$$!31+++]1aZTF1$\"3AR2Ablp!G#F17$ $!3umm;/#)[oPF1$\"3Qc-aUo!3d#F17$$!3hLLL$=exJ$F1$\"3w\\V#o=r74$F17$$!3 *RLLLtIf$HF1$\"3e7[C(>pun$F17$$!3]++]PYx\"\\#F1$\"3g7$[&H%ok[%F17$$!3E MLLL7i)4#F1$\"38YMhKkQB`F17$$!3c****\\P'psm\"F1$\"3Zs')))=Q[6kF17$$!3' )****\\74_c7F1$\"3u,'\\r$*)*)*e(F17$$!3)3LLL3x%z#)F-$\"3%\\X'z.f&o_1$F17$$\"3i******\\2goPF1$\"3cZ*z\"fEqqDF17$$\"3UKL$eR<*fT F1$\"3[e(yu>#zsAF17$$\"3m******\\)Hxe%F1$\"32U$4mYLa+#F17$$\"3ckm;H!o- *\\F1$\"3[y)p'p^#fn\"F17$$\"3y)***\\7k.6aF1$\"3/KSI!QN$y7F17$$\"3#emmm T9C#eF1$\"3vzR.eCbJ6F17$$\"33****\\i!*3`iF1$\"3q-#H4p\"*)y9F17$$\"3%QL LL$*zym'F1$\"3g#=?eo(\\I>b&e3l%F-7$$\"3c ***\\(o!*R-tF1$!3]wgb#*zC<))F-7$$\"3Nmm;HYt7vF1$!3b+\"[%)3#>-AF17$$\"3 j**\\(o*GP4wF1$!336#GlyhXb#F17$$\"3\"HL$ek6,1xF1$!3G[EOqL!4e#F17$$\"3a *\\P%)HIVv(F1$!3\\/+#oc3&HCF17$$\"3>m;HK%\\E!yF1$!3EE(y9Zl`9#F17$$\"3% G$e9m&o4&yF1$!3%>TkJ&*)o6#F 47$$\"3!GLek8!)**e)F1$\"3%\\fg7-ldN#F47$$\"32mmTN>^R')F1$\"3/qR_)*[]tC F47$$\"3K**\\PMP/*o)F1$\"3+lHsFXkEDF47$$\"3fKLLLbdQ()F1$\"3#[Za%*Hpk\\ #F47$$\"36+]i!*z>W))F1$\"3)46v>wLj0#F47$$\"3amm\"zW?)\\*)F1$\"3/po*y1$ od\"*F17$$\"3I+Dcw;j-!*F1$\"3a\"QtkB5ry'!#?7$$\"3'HL3_!HWb!*F1$!3ogI7$ Qz\"o6F47$$\"3rmT&Q8a#3\"*F1$!3*)H$*Rvl_OEF47$$\"3[++]i`1h\"*F1$!37\") )f*zh-@WF47$$\"3Z+](oW7;@*F1$!3([s!GdD_RkF47$$\"3Y++DJ&f@E*F1$!31P]Ar( Gkw)F47$$\"3Y+]i:mq7$*F1$!3F?7$$\"3Y++voyMk%*F1$!3[o0]Y-At? F?7$$\"3+,v$4T@'*[*F1$!36bRO;0&4C#F?7$$\"3W+]7`\\*[^*F1$!3pbj:S6W3CF?7 $$\"3*)*\\7`\\o,a*F1$!3kWxx6&>Pd#F?7$$\"3W++]P?Wl&*F1$!3`Q&4vN%\\MFF?7 $$\"3'4Dc,\">g#f*F1$!39Y[r#f?#**GF?7$$\"3Q+D\"Gyh(>'*F1$!31nK8]c!>0$F? 7$$\"3!)\\(oal@pk*F1$!3-)>uh*Q>)=$F?7$$\"3K+]7G:3u'*F1$!3$R]LY@&4.LF?7 $$\"3'3D\"y+9C,(*F1$!3#oy$=M\">4R$F?7$$\"3G+vVt7SG(*F1$!3%41&p+N?XMF?7 $$\"3UCcw47)>u*F1$!3!z0!ze!>vX$F?7$$\"3q\\P4Y6cb(*F1$!3W!Q@uJP'eMF?7$$ \"3'\\(=U#3T\"p(*F1$!3I*y[/4)\\ZMF?7$$\"3A++v=5s#y*F1$!30%4y[NxHU$F?7$ $\"3v]iS\"*3))4)*F1$!3/PGk3\"e*GLF?7$$\"3;+D1k2/P)*F1$!37Xk6fY>mJF?7$$ \"3KC1R+2i])*F1$!3_JqJvoXbIF?7$$\"3e\\(=nj+U')*F1$!3EQNJhz2BHF?7$$\"3% [(o/t0yx)*F1$!3O'Rze,(RnFF?7$$\"35+]P40O\"*)*F1$!3#*35#o1ime#F?7$$\"3P DJqX/%\\!**F1$!3[,pUog-zBF?7$$\"3k]7.#Q?&=**F1$!3p'ol3vUD9#F?7$$\"3!fP f$=.5K**F1$!3o@H-=B;v=F?7$$\"31+voa-oX**F1$!3Mi-M))psu:F?7$$\"3ACc,\"> g#f**F1$!3W#yaQPm*Q7F?7$$\"3[\\PMF,%G(**F1$!33*>%QCp#\\l)F47$$\"3uu=nj +U')**F1$!3')p#RN3ez^%F47$$\"\"\"F*F+-%&TITLEG6#Q@interpolation~polyno mial~for~~f6\"-%'COLOURG6&%$RGBG$F*F*Fb[n$\"*++++\"!\")-%+AXESLABELSG6 $Q\"xF][nQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F(Fgjm%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 98 "Th ere really isn't much need to plot the error here - it is obviously ba d. Nonetheless, here it is" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 96 "plot(poly-f(x),x=-1..1,thickness=3, color=blue,title=\"error in interpolation polynomial for f\");" }} {PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7ct7 $$!\"\"\"\"!$\"3#G\\zE_>[)=!#F7$$!3]LekynP')**!#=$!3G4BYaPbzX!#<7$$!3- n;HdNvs**F1$!3cE?p=!\\ys)F47$$!3_+v$fLI\"f**F1$!3=n)*G%)3:Z7!#;7$$!3/M Le9r]X**F1$!3e#owu`'e$e\"F?7$$!3an\"HK*Q)=$**F1$!3;\\8t.7]%)=F?7$$!3/, ](=ng#=**F1$!3%pFM+R&=_@F?7$$!3bM3_]uj/**F1$!3q%z`,j:))Q#F?7$$!3%pmm\" HU,\"*)*F1$!3%[gOn\\kqtN\\<$F?7$$!3)omT5!\\F4)*F1$!3J[a&z`/nL$F?7$$!3!RLL$ e%G?y*F1$!3g)*4eqf`HMF?7$$!3Sn\"zpB0%o(*F1$!3w(G`x<8MX$F?7$$!3#4+Dc,#y a(*F1$!3=C\"\\@W')QY$F?7$$!3UM3F%ze6u*F1$!376u#))*o3iMF?7$$!3#om;HdNvs *F1$!3WQpq(>\"3\\MF?7$$!3sK$3-8*G+(*F1$!3KNx\"y\"GT$R$F?7$$!3u****\\(o UIn*F1$!3Yf]5%\\KUI$F?7$$!3wm;zWizX'*F1$!3-(Gl*QX.)=$F?7$$!3xLL3-)\\&= '*F1$!317%y\"fX`]IF?7$$!3w+]PfLI\"f*F1$!3'[=s,qTn*GF?7$$!3ommm;p0k&*F1 $!3z!=u^MD5t#F?7$$!3_$3_D$*\\.a*F1$!3;KjFS+5!e#F?7$$!3D*\\P%[Hk;&*F1$! 3\\)eu_\"**>DCF?7$$!33;HKkf$H\\*F1$!3%R%3:)>E#oAF?7$$!3#HL3-)*G#p%*F1$ !3*4g_!H<$36#F?7$$!3w\\P4'*>_X%*F1$!3!fX1Y_VW&>F?7$$!3gm\"z>,:=U*F1$!3 oIi\\5PG+=F?7$$!3U$eky-3\")R*F1$!3AU.6@*)Q\\;F?7$$!3E++vV5Su$*F1$!3)oq \"of#GE]\"F?7$$!3&R$3_vq)pK*F1$!3OvvOlkCC7F?7$$!3_m;H2Jdz#*F1$!3#=TL=L 1Mp*F47$$!33*\\i!R\"f@B*F1$!3!)oH5Ghq.uF47$$!3wKL$3O2UO=F47$$!3mmmmT%p\"e()F1$\"3ktA9$QSeR# F47$$!3umm\"H-%\\/()F1$\"3igiE<1?lCF47$$!3qlm;/'=3l)F1$\"3Y)\\G_uQLF17$$!3\"ommT!R=0vF1$!3#GzJEam*zHF17$$!3tKL$3i_+I(F 1$!3CSMI`rHB7$$!3+nm;/siqmF1 $\"3!=\\1=OMm=(F]^l7$$!3[++](y$pZiF1$\"3C*[([&GP2P$F]^l7$$!33LLL$yaE\" eF1$!3y%[s;9$3&f\"F]^l7$$!3hmmm\">s%HaF1$!3k2D*[%o,k=F]^l7$$!3Q+++]$*4 )*\\F1$\"3]y5!>#p!Gb(!#A7$$!39+++]_&\\c%F1$\"3'G!GsVdpy&)!#?7$$!31+++] 1aZTF1$\"3+-zgkomlGF]`l7$$!3umm;/#)[oPF1$!3vA%*G)e$f:LF]`l7$$!3hLLL$=e xJ$F1$!3n&GOYu(3EKF]`l7$$!3*RLLLtIf$HF1$\"3_4iMW)*[8j!#@7$$!3]++]PYx\" \\#F1$\"3!>aQwcmVd#F]`l7$$!3EMLLL7i)4#F1$\"3(p$4'4zS$>mFbal7$$!3c**** \\P'psm\"F1$!3!4#e9M#\\Ea\"F]`l7$$!3')****\\74_c7F1$!32y#zOwa>-\"F]`l7 $$!3)3LLL3x%z#)F]^l$\"3Bz5m]`E.^Fbal7$$!3KMLL3s$QM%F]^l$\"3'yG2w#\\eA` Fbal7$$!3]^omm;zr)*Fbal$\"3E[9maXo0R!#C7$$\"3%pJL$ezw5VF]^l$\"3%)H=1&[ w1EmMgDF]`l7$$ \"3B)****\\2'HKHF1$\"3u\"*3I4koamFbal7$$\"3ElmmmZvOLF1$!3?jt^$\\#**yLF ]`l7$$\"3i******\\2goPF1$!35%Q!\\j)3XJ$F]`l7$$\"3UKL$eR<*fTF1$\"3SW(zZ NdL6$F]`l7$$\"3m******\\)Hxe%F1$\"3#\\jfH[9Ye)F]`l7$$\"3ckm;H!o-*\\F1$ \"3J1\\OhBbRQFbal7$$\"3y)***\\7k.6aF1$!36\\4SK>q-=F]^l7$$\"3#emmmT9C#e F1$!3_^h4iuxP:F]^l7$$\"33****\\i!*3`iF1$\"3A5S+$4P8X$F]^l7$$\"3%QLLL$* zym'F1$\"3[uiv<]#f>(F]^l7$$\"3wKLL3N1#4(F1$!3;;FcTO;\"R%F]^l7$$\"3c*** \\(o!*R-tF1$!38^R>]'*)*Q$F17$$\"3>m;HK%\\E!yF1$!3&=fsFt0V!HF17$$\"3%G$e9 m&o4&yF1$!3+Wm>')[,iCF17$$\"3Y*******p(G**yF1$!3#ROLC'zia=F17$$\"3)HLL 3U/37)F1$\"3;aAXyfveLF17$$\"3]mmmT6KU$)F1$\"30yna#=u>C\"F47$$\"3-LLeRZ QT%)F1$\"3<5C'=-naq\"F47$$\"3a****\\P$[/a)F1$\"3_tsf.$=f7#F47$$\"3!GLe k8!)**e)F1$\"3Gkt5'Q.BH#F47$$\"32mmTN>^R')F1$\"3!H>OnHB2T#F47$$\"3K** \\PMP/*o)F1$\"3uIW/EB`kCF47$$\"3fKLLLbdQ()F1$\"3*>FP+3;]V#F47$$\"36+]i !*z>W))F1$\"3yU,'R$>D'*>F47$$\"3amm\"zW?)\\*)F1$\"3YlZ#G2M,d)F17$$\"3I +Dcw;j-!*F1$!3(z4d7MT?8&F]^l7$$\"3'HL3_!HWb!*F1$!3=^WQc*\\cA\"F47$$\"3 rmT&Q8a#3\"*F1$!3E]6#eSqLp#F47$$\"3[++]i`1h\"*F1$!3w9>/&o^\\'F47$$\"3Y++DJ&f@E*F1$!3(\\d!3sB]@))F47$$\"3Y+] i:mq7$*F1$!3M'y0A;f[9\"F?7$$\"3Y++++PDj$*F1$!3ZXyz#4e_V\"F?7$$\"3+,v=U s_)Q*F1$!3a(z2FcO&*e\"F?7$$\"3X+]P%y+QT*F1$!3=l2];7')[F?7$$\"3Y++voyMk%*F1$!3?g\")[=r]y?F?7$$\"3+,v$ 4T@'*[*F1$!3QE\"QVu5iC#F?7$$\"3W+]7`\\*[^*F1$!3;x[L@\\n8CF?7$$\"3*)*\\ 7`\\o,a*F1$!3%HIu9/F*yDF?7$$\"3W++]P?Wl&*F1$!3W$=c\"HenRFF?7$$\"3'4Dc, \">g#f*F1$!3#\\lI\\FuV!HF?7$$\"3Q+D\"Gyh(>'*F1$!3-3`(=wJq0$F?7$$\"3!) \\(oal@pk*F1$!3#zg[Tl#H$>$F?7$$\"3K+]7G:3u'*F1$!3w`PKLo;3LF?7$$\"3'3D \"y+9C,(*F1$!3sH&ye#Q'fR$F?7$$\"3G+vVt7SG(*F1$!35ub^v9A]MF?7$$\"3UCcw4 7)>u*F1$!3sM!*H`P_iMF?7$$\"3q\\P4Y6cb(*F1$!370'[LyGOY$F?7$$\"3'\\(=U#3 T\"p(*F1$!3=k1EzjZ_MF?7$$\"3A++v=5s#y*F1$!3(4Gzx^UzU$F?7$$\"3v]iS\"*3) )4)*F1$!3>+yHsr*QL$F?7$$\"3;+D1k2/P)*F1$!3)zVV_#y5rJF?7$$\"3KC1R+2i])* F1$!3QHronrNgIF?7$$\"3e\\(=nj+U')*F1$!3R`ZfHa'z#HF?7$$\"3%[(o/t0yx)*F1 $!33c\"*[4T&Q&z:F?7$$\"3ACc,\">g#f**F 1$!3].FOVawV7F?7$$\"3[\\PMF,%G(**F1$!3M(R\"4jKz-()F47$$\"3uu=nj+U')**F 1$!3#Q\"*\\0a+dc%F47$$\"\"\"F*F+-%&TITLEG6#QIerror~in~interpolation~po lynomial~for~~f6\"-%'COLOURG6&%$RGBG$F*F*Fe[n$\"*++++\"!\")-%+AXESLABE LSG6$Q\"xF`[nQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F(Fjjm%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 91 "We expect the interpolation polynomial with Chebyshev nodes to much bett er. Let's check it." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 81 "n:=20: XX3:=[]: for k from 0 to n do XX3: =[op(XX3),-cos(Pi*(2*k+1)/(2*n+2))]: od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 4 "XX3;" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#77,$-%$cosG6#, $%#PiG#\"\"\"\"#U!\"\",$-F&6#,$F)#F+\"#9F-,$-F&6#,$F)#\"\"&F,F-,$*$-%% sqrtG6#\"\"$F+#F-\"\"#,$-F&6#,$F)#F?F3F-,$-F&6#,$F)#\"#6F,F-,$-F&6#,$F )#\"#8F,F-,$-F&6#,$F)#F9F3F-,$-F&6#,$F)#\"#F,F-\" \"!FinFYFTFNFHFC,$F;#F+FAF5F/F%" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "YY4:=map(f,XX3):" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 292 "Now if we ask Maple to compute the inter polation polynomial for the nodes determined by XX3 and YY4 we wil l have a very long wait. Maple will try to do the calculation exactly. The resulting expression is huge. To keep the time required reasonabl e we will ask Maple to use floating point." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "XX5:=evalf(XX3,16 );" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$XX5G77$!1,=\"=(z.s**!#;$!1O#= =7z#\\(*F($!1V?W'[P(3$*F($!1&QWy.a-m)F($!1)H!oC[J=yF($!1%>4xPF4xPF " 0 "" {MPLTEXT 1 0 19 "YY5:=evalf( YY4,16);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$YY5G77$\"1f!**>6nty%!#< $\"1nq0.!zv*\\F($\"1`OD*)HQbaF($\"1++++++]iF($\"1pxG%)pJhvF($\"1Hm " 0 "" {MPLTEXT 1 0 55 "Digits:=16: polyc:=sort(interp(XX5,YY5,x)); Digits:=10:" }} {PARA 12 "" 1 "" {XPPMATH 20 "6#>%&polycG,L*$)%\"xG\"#?\"\"\"$\"1[M\"Q rR3X%!#7*&$\"&^_%F-F*)F(\"#>F*!\"\"*&$\"1`kK[]%*eB!#6F*)F(\"#=F*F3*&$ \"&'>=F7F*)F(\"#\")fX&F7F*)F(\"#7F*F**&$\"&]R\"F7F*)F(\"#6F*F3*&$\"1@\"*4'ygYw#F7F* )F(\"#5F*F3*&$\"&bU%F-F*)F(\"\"*F*F**&$\"1tI " 0 "" {MPLTEXT 1 0 93 "plot(polyc,x=-1..1,th ickness=3,color=blue,title=\"Chebyshev interpolation polynomial for f \");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVE SG6#7_q7$$!\"\"\"\"!$\"3')oF.)Gdh;&!#>7$$!3-n;HdNvs**!#=$\"31&4c&zU\" \\z%F-7$$!3/MLe9r]X**F1$\"3A+')3-4$[d%F-7$$!3/,](=ng#=**F1$\"3SZU-tq-r WF-7$$!3%pmm\"HU,\"*)*F1$\"33iD!oG;VX%F-7$$!3()***\\PM@l$)*F1$\"3%3VqB 9<(*e%F-7$$!3!RLL$e%G?y*F1$\"3!H-&e:A@O[F-7$$!3u****\\(oUIn*F1$\"3C%[U Co:@L&F-7$$!3ommm;p0k&*F1$\"3do:Czx&Hg&F-7$$!3E++vV5Su$*F1$\"3\\&RNmZD [a&F-7$$!3wKL$3`F-7$$!3'fmT&Q75y!*F1$\"3-`)eEfr2x4 .\"F17$$!3[++](y$pZiF1$\"3elm&[\"Ro&>\"F17$$!33LLL$yaE\"eF1$\"3ypj.?(3 tJ\"F17$$!3hmmm\">s%HaF1$\"3N2)o#yCJ<9F17$$!3Q+++]$*4)*\\F1$\"3YzDP!oh pf\"F17$$!39+++]_&\\c%F1$\"3AeKZc@9(*=F17$$!31+++]1aZTF1$\"3s!>bs@wgG# F17$$!3umm;/#)[oPF1$\"31_lu&4jWo#F17$$!3hLLL$=exJ$F1$\"3UHhKMvx(=$F17$ $!3*RLLLtIf$HF1$\"3-g.z1'>)oOF17$$!3]++]PYx\"\\#F1$\"35z2=cKf\"Q%F17$$ !3QnmTNz>&H#F1$\"3c\"[wu$fsvZF17$$!3EMLLL7i)4#F1$\"3W![\")o#R`D_F17$$! 3#pm;aVXH)=F1$\"3#*eBqdxi\"y&F17$$!3c****\\P'psm\"F1$\"3iM(3Tp)\\$R'F1 7$$!3s*****\\F&*=Y\"F1$\"31`v_hC(3,(F17$$!3')****\\74_c7F1$\"37;(yB0%p QwF17$$!3ZmmT5VBU5F1$\"3'e6f/;:VF)F17$$!3)3LLL3x%z#)F-$\"3QDLOszAb))F1 7$$!3gKL$e9d;J'F-$\"3[99F'***F17$$\"3&yK$3_Nl.5F-$\"3m*G603t;)**F1 7$$\"3/E$ekGR[b\"F-$\"3uJ1P=c4c**F17$$\"3@CL$3-Dg5#F-$\"38!RRoCe'>**F1 7$$\"3e?Le*['R3KF-$\"3&4OH\"*[I\\\")*F17$$\"3%pJL$ezw5VF-$\"3]t*flTB%p '*F17$$\"3L`mmmJ+IiF-$\"3Iu#\\()emqK*F17$$\"3s*)***\\PQ#\\\")F-$\"3!=% o5Iz(z)))F17$$\"3ilm\"z\\1A-\"F1$\"3[r=njm1(F17$$\"3I*******pvxl\"F1$ \"3C*zVbk_9U'F17$$\"3g)***\\7JFn=F1$\"3VHH*R\"eQCeF17$$\"3#z****\\_qn2 #F1$\"3w`IoO(z*y_F17$$\"3=)**\\P/q%zAF1$\"3zS'px@W'4[F17$$\"3U)***\\i& p@[#F1$\"3YI\"eII,'*R%F17$$\"3B)****\\2'HKHF1$\"3!=/@kcYQn$F17$$\"3Elm mmZvOLF1$\"3#3WS#eWklJF17$$\"3i******\\2goPF1$\"3k@X'3WTVo#F17$$\"3UKL $eR<*fTF1$\"3aEw/DEetAF17$$\"3m******\\)Hxe%F1$\"3U373x%f$y=F17$$\"3ck m;H!o-*\\F1$\"3?PJBk#p7g\"F17$$\"3y)***\\7k.6aF1$\"3Br]mw].B9F17$$\"3# emmmT9C#eF1$\"3cygi*GO\\J\"F17$$\"33****\\i!*3`iF1$\"3A^PsX3(Q>\"F17$$ \"3%QLLL$*zym'F1$\"3@kVmsV9K5F17$$\"3wKLL3N1#4(F1$\"3$QAp&ox!Qj)F-7$$ \"3c***\\(o!*R-tF1$\"3h4F(GC.K/)F-7$$\"3Nmm;HYt7vF1$\"3Qo'*4'G!>,xF-7$ $\"3\"HL$ek6,1xF1$\"3If'>oq;/e(F-7$$\"3Y*******p(G**yF1$\"3gUd5XDi`vF- 7$$\"3)HLL3U/37)F1$\"3!HX7%z\"o%puF-7$$\"3]mmmT6KU$)F1$\"3G?GAfBEcrF-7 $$\"3a****\\P$[/a)F1$\"3))H2?]AYNmF-7$$\"3fKLLLbdQ()F1$\"3!\\x^:0kg*fF -7$$\"36+]i!*z>W))F1$\"3?_7Z\\-([o&F-7$$\"3amm\"zW?)\\*)F1$\"3)*e\\3&f $=ZaF-7$$\"3'HL3_!HWb!*F1$\"3lT5UVuG;`F-7$$\"3[++]i`1h\"*F1$\"3$)\\Y4$ =,iI&F-7$$\"3Y++++PDj$*F1$\"3*38)=!Q$>IbF-7$$\"3W++]P?Wl&*F1$\"3t4[w$Q 57g&F-7$$\"3K+]7G:3u'*F1$\"3V$e#>gICG`F-7$$\"3A++v=5s#y*F1$\"3?)Rkp+8G $[F-7$$\"3;+D1k2/P)*F1$\"3#yo,$3+s(e%F-7$$\"35+]P40O\"*)*F1$\"3Eydb4#* 3aWF-7$$\"3k]7.#Q?&=**F1$\"3s8([Gyk:Z%F-7$$\"31+voa-oX**F1$\"3!p#386-$ ed%F-7$$\"3[\\PMF,%G(**F1$\"3#ec=?1Oez%F-7$$\"\"\"F*$\"3A\")H " 0 "" {MPLTEXT 1 0 107 "plot(f(x)-polyc,x=-1..1,thickness=3,color=blue,title =\"error in Chebyshev interpolation polynomial for f\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6#7i_l7$$!\"\" \"\"!$!3(y,t>EDD/%!#?7$$!3]LekynP')**!#=$!3Op*R-6C$\\=F-7$$!3-n;HdNvs* *F1$!3OrTX#e%G,#)!#A7$$!3_+v$fLI\"f**F1$\"3&3.N3i())38F-7$$!3/MLe9r]X* *F1$\"3e4iWqT\")oBF-7$$!3/,](=ng#=**F1$\"3EC7HBfxeOF-7$$!3%pmm\"HU,\"* )*F1$\"3u'*3kSqvzSF-7$$!3'GLekynP')*F1$\"3'[:$e=xytQF-7$$!3()***\\PM@l $)*F1$\"3(y\"QRQQURKF-7$$!3)omT5!\\F4)*F1$\"32lh1mG[PBF-7$$!3!RLL$e%G? y*F1$\"3E#=GQc)='H\"F-7$$!3#4+Dc,#ya(*F1$\"3SWqI%HOx:#!#@7$$!3#om;HdNv s*F1$!3'*R*\\\\&p8u#)Fao7$$!3sK$3-8*G+(*F1$!31v#eQGcsx\"F-7$$!3u****\\ (oUIn*F1$!3MLLox&GYf#F-7$$!3xLL3-)\\&='*F1$!3)4:3wC$eVPF-7$$!3ommm;p0k &*F1$!3E,:=U*F1$!3EDk_W*G]n#F-7$$!3E++vV5Su$*F1$!3 TW9u&R,]h\"F-7$$!3&R$3_vq)pK*F1$!3`Hgp&pD^`%Fao7$$!3_m;H2Jdz#*F1$\"322 $*y41AUrFao7$$!33*\\i!R\"f@B*F1$\"3YB&o3lTm!=F-7$$!3wKL$3WF-7$$!37;HK*G=\")**)F1$\"3x4G.1:(fV%F-7$$!39***\\iId9(*)F1$\"3e CG5]%p)zVF-7$$!3YL$eRP8['))F1$\"3^7jt.!\\X]$F-7$$!3mmmmT%p\"e()F1$\"3W ?:lmB-\\=F-7$$!3umm\"H-%\\/()F1$\"3KxRS*p,h^)Fao7$$!3qlm;/'=3l)F1$!3Ox xE)GF#>=Fao7$$!3wlmT&=Vrf)F1$!3;*H]!*4Y))>\"F-7$$!3&emmmwnMa)F1$!3_sHM c.J^@F-7$$!3+mm;Hp6O%)F1$!3:Q.oXh&\\q$F-7$$!3:mmm\"4m(G$)F1$!30(*[KE#z jg%F-7$$!3-$eRs$f0-$)F1$!3(3QYYn(Q9ZF-7$$!3))*\\7GyX`F)F1$!3i,yzx*fYx% F-7$$!3w;aQGcj[#)F1$!33S[Q!yUty%F-7$$!3iL$eRZD>A)F1$!3Y)H;7rKIv%F-7$$! 3OnT5l^]o\")F1$!3'*Qge.6#za%F-7$$!3)****\\i&[3:\")F1$!3!p63%4.6sTF-7$$ !3Mm;aQUC3!)F1$!3c)H/))*od&)HF-7$$!3\"QLL3i.9!zF1$!37/jPx)=_Q\"F-7$$!3 Mnmm\"p[B!yF1$\"3ey%e:%[(yp#Fao7$$!3'3++Dw$H.xF1$\"3!pOXUV#Q**=F-7$$!3 RMLLL)QUg(F1$\"3!*4\"y2*zSMLF-7$$!3\"ommT!R=0vF1$\"3.'oS+chlV%F-7$$!37 MLL$3,RX(F1$\"3jst*)Q%GC%[F-7$$!3K++]i#=ES(F1$\"3()f_*f'RRB^F-7$$!3ULL 3_o(pP(F1$\"3GJhyeBA:_F-7$$!3ammmTaL^tF1$\"3_BJK*RnSF&F-7$$!3j***\\7.% pDtF1$\"3o:u@_&y(*H&F-7$$!3tKL$3i_+I(F1$\"3M'=B_t3CH&F-7$$!3Cmm;zp[(>( F1$\"3!pmTyg)**R\\F-7$$!3u****\\P8#\\4(F1$\"3zZ!>,F086%F-7$$!3%ommT!y% )))pF1$\"3%>lf\")yo\\%GF-7$$!3#RLL3FuF)oF1$\"3Y>4$[%*4,H\"F-7$$!3-,+]P 2qwnF1$!3;N*zxY2*RSFao7$$!3+nm;/siqmF1$!3mQ7,A+.#3#F-7$$!3#)******\\Q* [c'F1$!352V(\\t0=f$F-7$$!3uLL$e\\g\"fkF1$!3&**Q_yu[:\"[F-7$$!3mnmmTrU` jF1$!3@KGi,N;UcF-7$$!3[++](y$pZiF1$!3:KoOh(Q*=gF-7$$!3\"[eR(\\Q]?iF1$! 3!=E#Q6QIQgF-7$$!3-o\"z>\"RJ$>'F1$!3Xj@CKLhDgF-7$$!3E^(=U(R7mhF1$!3Q\" zFNS!)4)fF-7$$!3YM$ek.M*QhF1$!36;D<32o/fF-7$$!3!4]P4;aX3'F1$!3/(z\"zXg (*ecF-7$$!3MnmT&Gu,.'F1$!3QsC\\%p7QH&F-7$$!3?+]PMXT@fF1$!3B*4y:a_xB%F- 7$$!33LLL$yaE\"eF1$!3f$f49%Rb@GF-7$$!3YmmTN\"for&F1$!3='HP5#[jj8F-7$$! 3%)****\\([j5i&F1$\"3>\\_/7p&)4?Fao7$$!3BLLeRyEDbF1$\"3+U^S0Re!y\"F-7$ $!3hmmm\">s%HaF1$\"3A?%zu7^QG$F-7$$!3]***\\7)*G;K&F1$\"3g9%)*eD7\"yZF- 7$$!3]LL$3x&y8_F1$\"3A.$H$fC#*ffF-7$$!3]nmTgD%f5&F1$\"39&>SbyJ[u'F-7$$ !3Q+++]$*4)*\\F1$\"3#eN`89Ph2(F-7$$!3]++vVy-r\\F1$\"3#o#\\zx#GV3(F-7$$ !3j++]Pj&R%\\F1$\"39*zeRzb?1(F-7$$!3?++DJ[)o\"\\F1$\"3y;$\\Tpa$4qF-7$$ !3K+++DL\")*)[F1$\"3'pU)>bpQEpF-7$$!3,++]7.nN[F1$\"3#\\ezySE4n'F-7$$!3 E++++t_\"y%F1$\"3#*>GBl-X*H'F-7$$!3?+++v7CtYF1$\"3$RKXhYjHB&F-7$$!39++ +]_&\\c%F1$\"3w%>t-q*4$z$F-7$$!3E+++D%yF^%F1$\"3MPXAFq'=*HF-7$$!3S++++ ;ggWF1$\"32X:#=3Bj8#F-7$$!3_+++vZU3WF1$\"3NRT%*p=**Q7F-7$$!33+++]zCcVF 1$\"3-(o2myj,8$Fao7$$!3m*****\\7rSI%F1$!3=!\\Kv=I.G'Fao7$$!3!)*******H %*=D%F1$!3o.`qlWNq:F-7$$!3#******\\Z<(*>%F1$!3**RZ\"Q=7,]#F-7$$!31+++] 1aZTF1$!3vdL0oIj.MF-7$$!3Ym;aQvx_SF1$!34\\F%HONX$\\F-7$$!3SLL3FW,eRF1$ !3giP&f)prgiF-7$$!3M+]i:8DjQF1$!3Mv2.bY)oJ(F-7$$!3umm;/#)[oPF1$!3a,-MK <.]!)F-7$$!3'**\\i:&p97PF1$!3?2f4EwV;$)F-7$$!3=L$e*)p0el$F1$!3]M%or!zZ \\%)F-7$$!32]ils]jFOF1$!3%)*>@Tr3\\Y)F-7$$!3%p;ajWk%*f$F1$!31Hqxx2)fW) F-7$$!3%Q3_+#QHrNF1$!3+OhZzQk#R)F-7$$!3<++v$>BJa$F1$!3$GfAXq()\\I)F-7$ $!3s$pFF1$\"3M --O0TdJNF-7$$!3CnmT&o_Qr#F1$\"3))*4OnKH=c%F-7$$!3U+vVtJLeEF1$\"3=$pHxb 8*HbF-7$$!3:M$e9m8Gg#F1$\"3=*)[wliU?kF-7$$!3Kn\"z%\\THZDF1$\"3tLW8c67> sF-7$$!3]++]PYx\"\\#F1$\"3G'p9\"p2:8zF-7$$!3%RLekG'[$R#F1$\"3Cr;!*)*Hl ^))F-7$$!3QnmTNz>&H#F1$\"3A%R'=!y..Q*F-7$$!3u]ilZeiqAF1$\"3#)4VG2#*yW% *F-7$$!37Me*)fP0YAF1$\"3GA'R)f_m\"[*F-7$$!3Z#F1$\"3gz;aoQBs%*F-7$$!3cnT&)3awZ@F1$\"3=6\"Ro0?AN*F-7$ $!3EMLLL7i)4#F1$\"3@'e+T\\;L7*F-7$$!3g+]PMLy!*>F1$\"3Cm]!GqwqD)F-7$$!3 #pm;aVXH)=F1$\"3(okkq#QiTpF-7$$!3&**\\Pf[E!H=F1$\"3vo*Q$>x#>9'F-7$$!3C L$ek`2^x\"F1$\"3MG_wJ*G^E&F-7$$!3bm\"zpe)=@7KNwhK%F-7$$!3c** **\\P'psm\"F1$\"3\\G)R@rh6M$F-7$$!3u**\\(o/Efh\"F1$\"3On$)QPZ)eP#F-7$$ !3k***\\iX#ek:F1$\"3!e'Qs'fO$*R\"F-7$$!3a**\\il)QK^\"F1$\"3M+,uDv\\oUF ao7$$!3s*****\\F&*=Y\"F1$!3M#4ri:4IE&Fao7$$!3*)**\\P%o^0T\"F1$!35s#3'z Nu!z!fF-7$$!3IJ$e*)fEm))*!#>$!3v#4B&>*>8 8'F-7$$!3#y**\\P454N*Fbjm$!30\"f$p9dQMiF-7$$!33Je9T=0$3*Fbjm$!31LvnY%f 9C'F-7$$!3Nk;a)e$>:))Fbjm$!3ZhTATJg>iF-7$$!3i(\\PfLNta)Fbjm$!3)>'*oF%f cphF-7$$!3)3LLL3x%z#)Fbjm$!3/RIl*)pE#4'F-7$$!3uJLe9rc&H(Fbjm$!3cA(y.&p l$f&F-7$$!3gKL$e9d;J'Fbjm$!3A^8Y3X>:[F-7$$!3XLL3xruF`Fbjm$!35*[(f%Ra6& QF-7$$!3KMLL3s$QM%Fbjm$!3P&Go&R,%4\"GF-7$$!3'ym;aQdDG$Fbjm$!3c[N2$*)f_ t\"F-7$$!3T,+]ivF@AFbjm$!3w8QGN0Y)R)Fao7$$!3'\\L$eRx**f6Fbjm$!3]r7swi) )oBFao7$$!3]^omm;zr)*Fao$!3Pun!o)e_EqM@%Fao7$$\"3@CL$3-Dg5#Fbjm$!3?<5?bL5&e(Fao7$$\"3e?Le*['R3KFbjm $!3+vz=\"3l`m\"F-7$$\"3%pJL$ezw5VFbjm$!3G1w><%3gx#F-7$$\"39&)**\\ibQq_ Fbjm$!3;a\"[8R&\\\"z$F-7$$\"3L`mmmJ+IiFbjm$!3\">+)GBD#3u%F-7$$\"3_@L$3 x?'*=(Fbjm$!31CPPL`x@bF-7$$\"3s*)***\\PQ#\\\")Fbjm$!3w:4))G4:XgF-7$$\" 3.B3_+'F-7$$\" 3m*[i:N[l#*)Fbjm$!3'pJ-;ne@B'F-7$$\"3)HK$3x;l&=*Fbjm$!3jf=)R_.AC'F-7$$ \"3HcTg-]vW%*Fbjm$!3Mbp5@?&\\A'F-7$$\"3h*)\\7G$eQq*Fbjm$!3]e%Gc$4%*zhF -7$$\"3%H#ek`;'H'**Fbjm$!3#yDQ)z4'o5'F-7$$\"3ilm\"z\\1A-\"F1$!3'[FWaHd b+'F-7$$\"3&*)*\\7Gy%e7\"F1$!3RqtJRcO?`F-7$$\"3GKLLe\"*[H7F1$!3W8S.AEM 5UF-7$$\"3kl;/w\\-$G\"F1$!3woO'QI)p)[$F-7$$\"3.***\\PzglL\"F1$!3]mI@hP M!o#F-7$$\"3UK$e9h'4!R\"F1$!3Ui]KcA[)z\"F-7$$\"3ylm;HCjV9F1$!3!Q()z+-% Gy&)Fao7$$\"3;**\\(oCor\\\"F1$\"3c6\\)G')4sD\"Fao7$$\"3aKLekSq]:F1$\"3 /9AV)eFa8\"F-7$$\"3#fm\"H#))RUg\"F1$\"3qZT#\\>cS:#F-7$$\"3I*******pvxl \"F1$\"3#yV)*o<;U;$F-7$$\"37**\\7`+:5fF-7$$\" 3g)***\\7JFn=F1$\"3!eI&)4:azr'F-7$$\"3G)**\\(==-s>F1$\"3m2qI'G%Qe!)F-7 $$\"3#z****\\_qn2#F1$\"35SxYajM()*)F-7$$\"3*z\\(o/aWF@F1$\"3)Hg%*oRH1F *F-7$$\"3/)*\\P%G?\"y@F1$\"3nUjH!yY$R%*F-7$$\"37)\\iS;&zGAF1$\"39M4\\O e*4\\*F-7$$\"3=)**\\P/q%zAF1$\"3(Qm,&R2vC%*F-7$$\"3I)*\\7.)>3Q#F1$\"3w Y]%=YpN%*)F-7$$\"3U)***\\i&p@[#F1$\"3]E/N6)[<-)F-7$$\"3S)\\ilPN%QDF1$ \"3+U.kM2>PtF-7$$\"3Q)*\\i!>,Zf#F1$\"3ly]AM<=VlF-7$$\"3O)\\(o/q'4l#F1$ \"3Yj]>@(\\Gl&F-7$$\"3L)**\\(=GB2FF1$\"3opB+BV$3o%F-7$$\"3I)\\7Gj)\\jF F1$\"3\\[/(e#y-VOF-7$$\"3F)*\\(oWk(>GF1$\"3+>EU_0IcDF-7$$\"3E)\\P4EIg( GF1$\"3z*)e\\E^@Q9F-7$$\"3B)****\\2'HKHF1$\"3#)RJG\")4MnIFao7$$\"3gJ$e 9T`G)HF1$!3q&RRvL_[1(Fao7$$\"3)\\m;zu5M.$F1$!3W&*e#e%3%Hq\"F-7$$\"3Q)* \\P%3oR3$F1$!3)ym]LpL*pEF-7$$\"3uJL$3UDX8$F1$!35(yzc;<_f$F-7$$\"3])** \\P4ScB$F1$!3#=M(G?\\&[F&F-7$$\"3ElmmmZvOLF1$!3m%pInU(femF-7$$\"3&))** *\\iirWMF1$!3`*R%fUS)Rt(F-7$$\"3WKLLexn_NF1$!3([Wraa_&Q$)F-7$$\"3ql;HK \"o'zNF1$!3:$QbI!347%)F-7$$\"3&*)**\\i]emg$F1$!3i:\")G&\\'3a%)F-7$$\"3 @K$3-))[Oj$F1$!3)pQm+v%\\k%)F-7$$\"3-mm;a#R1m$F1$!3Js?r`&)QV%)F-7$$\"3 5LL3-+i9PF1$!3%G&[`g6^2$)F-7$$\"3i******\\2goPF1$!3eYs'\\Ip$\\!)F-7$$ \"3#GLe9\"*Hk'QF1$!3I'pET)eV'G(F-7$$\"3-mm\"H2fU'RF1$!3CQ]'4O!)4='F-7$ $\"3A**\\PM#)3iSF1$!3_t!*\\y6/#z%F-7$$\"3UKL$eR<*fTF1$!35#)om6,S#>$F-7 $$\"3$e;/E&RR8UF1$!3_'HTz\"=ReAF-7$$\"3C**\\P40(oE%F1$!3P-fQ,9c+8F-7$$ \"3kKe9mqM?VF1$!3$ozcy0rpL$Fao7$$\"3/mm\"HiBQP%F1$\"3!4b0KT`HF'Fao7$$ \"3X*\\(oz,IFWF1$\"3I:'Qam:wc\"F-7$$\"3%GLektw2[%F1$\"3GorPd@zsCF-7$$ \"3Em\"HKH`U`%F1$\"3BF.uH=!*GLF-7$$\"3m******\\)Hxe%F1$\"3)H\"Q&)pWyAT F-7$$\"3!fm\"z%Rk$)o%F1$\"36XLB%*z71aF-7$$\"37KLeR*)**)y%F1$\"3/r$pgPB uN'F-7$$\"3^l\"z>@;$R[F1$\"3oj>3mW#=p'F-7$$\"3M)*\\P%[L'*)[F1$\"3c)H#R b^tDpF-7$$\"3/:Hd?@z9\\F1$\"3>Fe,Lv,/qF-7$$\"3=J3xc2&*R\\F1$\"3-xrt]N< cqF-7$$\"3)yuoHR4^'\\F1$\"3EA5'pvh?3(F-7$$\"3ckm;H!o-*\\F1$\"3y-)QvIN; 3(F-7$$\"3c(****\\7ga4&F1$\"3!HH*H3Vf(z'F-7$$\"3nJL$3A_1?&F1$\"3q6\"pc -Z\"ygF-7$$\"3ylmm;V%eI&F1$\"37(*o6bX!H(\\F-7$$\"3y)***\\7k.6aF1$\"3w! \\.zo\\qb$F-7$$\"3al;a84)Q^&F1$\"3S)G*Q(*e,l>F-7$$\"3IKLe9as;cF1$\"3SD .`_^!*GFFao7$$\"31**\\i:*p&>dF1$!3KI>u:Iy19F-7$$\"3#emmmT9C#eF1$!39>9D Q&31'HF-7$$\"38**\\7yI3IfF1$!3sfuWg=WOVF-7$$\"3WKLeR'F1$!3xDV$Gx>6.'F- 7$$\"3>:a8(*=')=R7#F-7$$\"3e LL3Fe#Rx'F1$!3#3Rw4#fk([%Fao7$$\"3ILL$3sr*zoF1$\"3+\"*R!\\BekC\"F-7$$ \"3-LLe9w,')pF1$\"3&=rV;t(f1GF-7$$\"3wKLL3N1#4(F1$\"3eEw#=GW@3%F-7$$\" 3gl;a)GJs>(F1$\"3rw%>7f'\\Q\\F-7$$\"3c***\\(o!*R-tF1$\"3'ymt]m_WH&F-7$ $\"3@#3-Q,\"pGtF1$\"3q=(z\\=$[)H&F-7$$\"3)f;a)eH)\\N(F1$\"3kl:h$R;xE&F -7$$\"3u\\i!R!\\F\"Q(F1$\"3oP\\()=M7-_F-7$$\"3^L$e*[oc2uF1$\"3zR\\Tc4! >5&F-7$$\"3#**\\i!R2:guF1$\"3g/TI*z0&*z%F-7$$\"3Nmm;HYt7vF1$\"3mbO$fSS nO%F-7$$\"3j**\\(o*GP4wF1$\"3)eA\">#GvrE$F-7$$\"3\"HL$ek6,1xF1$\"3!R2) >@wbc=F-7$$\"3>m;HK%\\E!yF1$\"3W%f/G$)*3ZEFao7$$\"3Y*******p(G**yF1$!3 !=$pTI\"p3N\"F-7$$\"3ymmTgg/5!)F1$!3w^i*e6F'4IF-7$$\"3)HLL3U/37)F1$!3: WXlE(***>UF-7$$\"3km;/,O=w\")F1$!3)>e/`_s#)e%F-7$$\"3=***\\7yi:B)F1$!3 GB`%\\(pxqZF-7$$\"3cmTNrBDf#)F1$!3YYPp3M(zy%F-7$$\"3%GLe9'>%pG)F1$!3'o ?wb=fVv%F-7$$\"36*\\i:bJYJ)F1$!37[Hlr^ZpYF-7$$\"3]mmmT6KU$)F1$!3')z\" \\Y;$RLXF-7$$\"3-LLeRZQT%)F1$!3O$3\\)o:`UOF-7$$\"3a****\\P$[/a)F1$!3tB !)ot)*4-AF-7$$\"3!GLek8!)**e)F1$!3)H=(yw)e1L\"F-7$$\"32mmTN>^R')F1$!3I XMSG%4@*RFao7$$\"3K**\\PMP/*o)F1$\"39Qi;;)zTb&Fao7$$\"3fKLLLbdQ()F1$\" 3gMWC,*yD\\\"F-7$$\"3Mm\"z>w'Q\"z)F1$\"3kJq!fK]RU#F-7$$\"36+]i!*z>W))F 1$\"3WTYr?b7LKF-7$$\"3yK3F>#4q*))F1$\"3=1+O=aXoQF-7$$\"3amm\"zW?)\\*)F 1$\"37_@YbN1$G%F-7$$\"3'GeRA1Ei(*)F1$\"3G!y)[_\"3^R%F-7$$\"3I+Dcw;j-!* F1$\"3R)HouZV$QWF-7$$\"3,fRs$[Me,*F1$\"3b))*R#Q!3KV%F-7$$\"3i;a)3HP!H! *F1$\"3W*o551])4WF-7$$\"3Buo/)4SA/*F1$\"3#*oyR)Hg!oVF-7$$\"3'HL3_!HWb! *F1$\"3I)H:t#))p2VF-7$$\"3rmT&Q8a#3\"*F1$\"3cD&*=^FO!)QF-7$$\"3[++]i`1 h\"*F1$\"3%pJ)=742lJF-7$$\"3Z+](oW7;@*F1$\"3.-()y\"\\Y&RAF-7$$\"3Y++DJ &f@E*F1$\"3G.b(\\i0!G6F-7$$\"3Y+]i:mq7$*F1$!3#*[&34&3WX)*F97$$\"3Y++++ PDj$*F1$!3i6eAiFMZ8F-7$$\"3X+]P%y+QT*F1$!3Ox8otDR2DF-7$$\"3Y++voyMk%*F 1$!3#)[!y_+zQX$F-7$$\"3W+]7`\\*[^*F1$!3S7>L!)QWdSF-7$$\"3W++]P?Wl&*F1$ !3'G+1>P(Q(>%F-7$$\"3Q+D\"Gyh(>'*F1$!3vmJDRfbDPF-7$$\"3K+]7G:3u'*F1$!3 u:s?*pVic#F-7$$\"3'3D\"y+9C,(*F1$!3'e!=S@47YZd\\YfBF-7$$\"3;+D1k2/P)*F1$\" 3+CS^)=hWD$F-7$$\"3e\\(=nj+U')*F1$\"3]$G=j$p'3)QF-7$$\"35+]P40O\"*)*F1 $\"3k%[3:(*)yySF-7$$\"3k]7.#Q?&=**F1$\"3[%[#eZ$))4l$F-7$$\"31+voa-oX** F1$\"3o%o<6SEsN#F-7$$\"3ACc,\">g#f**F1$\"39k2%zOWsH\"F-7$$\"3[\\PMF,%G (**F1$!3(y*=>,4>-#*F97$$\"3uu=nj+U')**F1$!3iOsm4#3c&=F-7$$\"\"\"F*$!3S T^P$)HaUSF--%&TITLEG6#QSerror~in~Chebyshev~interpolation~polynomial~fo r~~f6\"-%'COLOURG6&%$RGBG$F*F*Fbdq$\"*++++\"!\")-%+AXESLABELSG6$Q\"xF] dqQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F(Fecq%(DEFAULTG" 1 2 0 1 10 3 2 6 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 149 "Even for this \+ example the cubic spline does a little better than the Chebyshev inter polation polynomial. Of course, that will not always be the case." }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 263 38 "Example. A sample natural cubic spline " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 168 "In the examples above I suppresse d the formulae for the cubic splines because they a very large. To see what a cubic spline \"looks like\" let's look at a simple example." } }{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "XX6:=[0,1,2,3,4,5] : YY6:=[-1,-1,1,0,2,2]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 " ex6:=spline(XX6,YY6,x,cubic);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$ex 6G-%*PIECEWISEG6'7$,(!\"\"\"\"\"*&#\"\"*\"#6F+%\"xGF+F**&#F.F/F+)F0\" \"$F+F+2F0F+7$,*#\"#@F/F+*&#\"$0\"F/F+F0F+F**&#\"#'*F/F+)F0\"\"#F+F+*& #\"#BF/F+*$F3F+F+F*2F0FA7$,*#!$(QF/F+*&#\"$2&F/F+F0F+F+*&#\"$5#F/F+*$F @F+F+F**&#\"#GF/F+F3F+F+2F0F47$,*\"#!*F+*&#\"$q)F/F+F0F+F**&#\"$\\#F/F +F@F+F+*&#FDF/F+FEF+F*2F0\"\"%7$,*#!%e5F/F+*&#\"$m'F/F+F0F+F+*&#\"$N\" F/F+FQF+F**&F2F+F3F+F+%*otherwiseG" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 104 "To see how it behaves we will \+ plot it together with a piecewise linear function joining the data poi nts." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "ex6lin:=spline(XX6,YY6,x,linear);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%'ex6linG-%*PIECEWISEG6'7$!\"\"2%\"xG\"\"\"7$,&!\"$F ,*&\"\"#F,F+F,F,2F+F17$,&\"\"$F,F+F)2F+F57$,&!\"'F,*&F1F,F+F,F,2F+\"\" %7$F1%*otherwiseG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 130 "plot( [ex6,ex6lin],x=0..5,thickness=3,color=[blue,red], title=\"Natural cubi c spline compared with piecewise linear interpolation\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6$7_o7$$\"\"!F) $!\"\"F)7$$\"3GLLL3x&)*3\"!#=$!3G/y&=-6\")3\"!#<7$$\"3umm\"H2P\"Q?F/$! 3!o(y*HfH)f6F27$$\"3MLL$eRwX5$F/$!3tR;=L#G&H7F27$$\"33ML$3x%3yTF/$!3-d I&=zp@G\"F27$$\"3h+]PfyG7ZF/$!3[M?jJn$**H\"F27$$\"3emm\"z%4\\Y_F/$!31= 42D@568F27$$\"32++v$fl2'o[) y!G\"F27$$\"3Qnm\"HKkIz(F/$!3[)Q.5Nz.D\"F27$$\"3>MLe*)>VB$)F/$!3R&RQz% *3#47F27$$\"3Y++DJbw!Q*F/$!3))G%3Qc6@4\"F27$$\"3%ommTIOo/\"F2$!3+l$G^6 #*==*F/7$$\"3YLL3_>jU6F2$!3VHC-*fWtA(F/7$$\"37++]i^Z]7F2$!3s&yF/7$$\"3/++]P[6j9F2$\"3%pvcsPk!ew!#>7 $$\"3UL$e*[z(yb\"F2$\"3j\\%H(4WvPJF/7$$\"3wmm;a/cq;F2$\"3%o.w!**zA0dF/ 7$$\"3%ommmJF2$\"3#*3b4$y!=*f*F/7$$\"3gmmm\"pW`(>F2$\"3*pq?9]h)4**F /7$$\"3?]iSmTI-?F2$\"3];Pto\"31+\"F27$$\"3dLe9TOEH?F2$\"3NEI*=/vZ+\"F2 7$$\"3'pT&)e6Bi0#F2$\"3CL+IVlr.5F27$$\"3K+]i!f#=$3#F2$\"3sjmtW$>t(**F/ 7$$\"3/++D\"=EX8#F2$\"3PTb(*HE(yt*F/7$$\"3?+](=xpe=#F2$\"3EX=+s4G^$*F/ 7$$\"37nm\"H28IH#F2$\"3iT`r.jKh\")F/7$$\"3um;zpSS\"R#F2$\"3%3K&e$=DWu' F/7$$\"3GLL3_?`(\\#F2$\"35iyk>UQS]F/7$$\"3fL$e*)>pxg#F2$\"37gYnwUOoKF/ 7$$\"33+]Pf4t.FF2$\"3y)*GkKtY\")=F/7$$\"3uLLe*Gst!GF2$\"33;\\r8Tn%4(Fb q7$$\"3)om\"H2\"34'GF2$\"3$)eB-UA[3HFbq7$$\"30+++DRW9HF2$\"3B6+IBR!4-$ !#?7$$\"3S+Dc,6jSHF2$!3!eFRi?_jE$Fgw7$$\"3K+]7y#=o'HF2$!34T$ox/O'Q\\Fg w7$$\"3C+voaa+$*HF2$!3uZfjt[n@>IF2$\"3dh!R(\\h-EmFgw7 $$\"3A+v$4^n)pIF2$\"3Qf[i#\\**zp$Fbq7$$\"3F+]i!RU07$F2$\"3A3*[!egGp%)F bq7$$\"3+++v=S2LKF2$\"3ivCzl!)3XCF/7$$\"3Jmmm\"p)=MLF2$\"3#>BZg'4E&R%F /7$$\"3B++](=]@W$F2$\"3U$)p&eGFH'oF/7$$\"35L$e*[$z*RNF2$\"3:hHU/ie8$*F /7$$\"3e++]iC$pk$F2$\"30A=X)=.$37F27$$\"3[m;H2qcZPF2$\"3qK0@v=:k9F27$$ \"3O+]7.\"fF&QF2$\"3w%R\\Fi?Dr\"F27$$\"3Ymm;/OgbRF2$\"36_)>WJ'pA>F27$$ \"3w**\\ilAFjSF2$\"3#)z1H0s\"R4#F27$$\"3yLLL$)*pp;%F2$\"3hBf.%*>g3AF27 $$\"3)QL3-$H**>UF2$\"3yV!QmA2*\\AF27$$\"3)RL$3xe,tUF2$\"37P\\_HrW!G#F2 7$$\"3h+v=n(*fDVF2$\"32\"4$[-H#3I#F27$$\"3Cn;HdO=yVF2$\"3v0B2&RW?J#F27 $$\"3MMe9\"z-lU%F2$\"3A'f(fOt*[J#F27$$\"3a+++D>#[Z%F2$\"3g0sorr<6BF27$ $\"3TM$3_5,-`%F2$\"3yXLVF*7$FHF*7$FRF* 7$FfnF*7$F`oF*7$$\"3wmmTg()4_))F/F*7$FeoF*7$$\"3=nT&)3\\m_'*F/F*7$$\"3 +N$ekGkX#**F/F*7$$\"31]iSmjk>5F2$!3!))*\\(=nsqg*F/7$Fjo$!3;jmm;RFj!*F/ 7$$\"3E+]7GTt%4\"F2$!3y%***\\PuJ0\")F/7$F_p$!3'3LL$e4OZrF/7$Fdp$!3[(** ***\\n\\!*\\F/7$Fip$!3V)*****\\ixCGF/7$F^q$!3Y\"******\\KqP(Fbq7$Fdq$ \"3]om;z*ev:\"F/7$Fiq$\"3GNLL$347T$F/7$F^r$\"3tOLLLjM?`F/7$Fcr$\"3$3++ D\"o7TvF/7$Fhr$\"3Qmm\"HK5S_)F/7$F]s$\"3$>LLL$Q*o]*F/7$$\"3_ek.HW#)))> F2$\"39q\"H2e)[w(*F/7$Fbs$\"3/)\\PfLep(**F/7$$\"3*=/wP!Ry:?F2$\"3:\"eR A'4;U)*F/7$Fgs$\"3Ek;a)ejtq*F/7$F\\t$\"3[Ie9T)oxV*F/7$Fat$\"3q'**\\P4u \"o\"*F/7$Fft$\"3e****\\(=QZl)F/7$F[u$\"3*z**\\7G-89)F/7$F`u$\"3%)GL$3 Fp)pqF/7$Feu$\"3YKL3-$ff3'F/7$Fju$\"38nm;z%zY-&F/7$F_v$\"35kmT5!3B#RF/ 7$Fdv$\"3C***\\iS!piHF/7$Fiv$\"3lim;/rFE>F/7$F^w$\"35JL3F*=4R\"F/7$Fcw $\"3Q&******\\2cb)Fbq7$Fiw$\"3af*\\P%)*)o$fFbq7$F^x$\"37o**\\(=s\"=LFb q7$$\"3G]iSm=\"*zHF2$\"3Ts\\PfL\")3?Fbq7$Fcx$\"3'pw*\\7`a%*pFgw7$$\"3A D\"G)[Ab**HF2$\"37%)yu=<^xW!#@7$$\"3>](oH/*41IF2$\"3.Q+v$f3)>7Fbq7$$\" 3F27$$\"3Y\\7`p2_#)RF2$\"3!*)\\i!R:/l>F27$$ \"3*G$e*[$zV4SF2$\"\"#F)7$$\"3L;/E+^NOSF2F[]m7$F_\\lF[]m7$$\"3ym\"zW7@ ^6%F2F[]m7$Fd\\lF[]m7$F^]lF[]m7$Fh]lF[]m7$Fb^lF[]m7$F\\_lF[]m7$Fa_lF[] m7$Ff_lF[]m7$F[`lF[]m7$F``lF[]m-Fe`l6&Fg`lFh`lF(F(-%&TITLEG6#Q\\oNatur al~cubic~spline~compared~with~piecewise~linear~interpolation6\"-%+AXES LABELSG6$Q\"xFc^mQ!6\"-%*THICKNESSG6#\"\"$-%%VIEWG6$;F(F``l%(DEFAULTG " 1 2 0 1 10 3 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "C urve 2" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 138 "As you can see the spline does not oscillate much betwee n the data points. Let's try the same comparison for the interpolation polynomial." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "ex7:=interp(XX6,YY6,x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$ex7G,.*$)%\"xG\"\"&\"\"\"#!#6\"#g*&#\"#b\"#CF*)F(\" \"%F*F**&\"#5F*)F(\"\"$F*!\"\"*&#\"$D%F1F*)F(\"\"#F*F**&#\"$*eF-F*F(F* F8F*F8" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 133 "plot([ex7,ex6lin ],x=0..5,thickness=3,color=[blue,red],title=\"Interpolation polynomial compared with piecewise linear interpolation\");" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6(-%'CURVESG6$7fo7$$\"\"!F)$!\"\"F) 7$$\"3WmmmT&)G\\a!#>$!31q3az?&R[\"!#<7$$\"3GLLL3x&)*3\"!#=$!3;>Xj@R;s= F27$$\"3))**\\i!R(*Rc\"F6$!3q&za*Ra1R@F27$$\"3umm\"H2P\"Q?F6$!3Q1q&Q)3 %fM#F27$$\"3!***\\PMnNrDF6$!3hp1U'ykN^#F27$$\"3MLL$eRwX5$F6$!3i>;d@4L> EF27$$\"3_LLe*[`HP$F6$!3Gd!Q8AX8l#F27$$\"3rLLL$eI8k$F6$!31]mN[DFqEF27$ $\"3_L$3-8>bx$F6$!3!4?F27$$\"3>MLe*)>VB$)F6$!3uAk(oyTBk\"F27$$\"3Y++DJbw!Q*F 6$!3n1SkQ'*=R7F27$$\"3%ommTIOo/\"F2$!3>LkQqCq2#)F67$$\"3YLL3_>jU6F2$!3 aGk3W^vcYF67$$\"37++]i^Z]7F2$!3Oy\\ma,2(e*F/7$$\"33++](=h(e8F2$\"31`:$ oDt^J#F67$$\"3/++]P[6j9F2$\"3mX!4>+m.(\\F67$$\"3UL$e*[z(yb\"F2$\"3j@7B >Q$R\"pF67$$\"3wmm;a/cq;F2$\"3cMy![M'R@')F67$$\"3%ommmJF2$\"3_G(G()4;B,\"F27$$\"3gmmm\"pW`(>F2$\"3(*y!> yD-m+\"F27$$\"3dLe9TOEH?F2$\"3=i8JI`\"3*)*F67$$\"3K+]i!f#=$3#F2$\"3M^# 30[_qg*F67$$\"3?+](=xpe=#F2$\"3M!*yV>H8/))F67$$\"37nm\"H28IH#F2$\"3%oC K)*\\i@n(F67$$\"3um;zpSS\"R#F2$\"3!GKzay*y]kF67$$\"3GLL3_?`(\\#F2$\"3: :^.,URL]F67$$\"3fL$e*)>pxg#F2$\"3]utE\\Y(*fNF67$$\"33+]Pf4t.FF2$\"3W*p o>RSfO#F67$$\"3uLLe*Gst!GF2$\"3i*RUrKO$f7F67$$\"30+++DRW9HF2$\"3U;As`t uxSF/7$$\"3:++DJE>>IF2$!3!zR5;^$H\\`!#?7$$\"3A+v$4^n)pIF2$!3`I*f2)Q`B7 F/7$$\"3F+]i!RU07$F2$!3YHaM5f:m!)F[y7$$\"39+vo/#3o<$F2$\"3]O!)oL3\"y- \"F/7$$\"3+++v=S2LKF2$\"3%H]134[fP%F/7$$\"3Jmmm\"p)=MLF2$\"3A*)eN\\GEP 9F67$$\"3B++](=]@W$F2$\"3[oq,))[d'3$F67$$\"35L$e*[$z*RNF2$\"3i&or*Q+w+ ^F67$$\"3e++]iC$pk$F2$\"3UQ7uG]TWyF67$$\"3[m;H2qcZPF2$\"3/)4wFfG&*3\"F 27$$\"3O+]7.\"fF&QF2$\"3g')eDd88\\9F27$$\"3Ymm;/OgbRF2$\"3kYQ'z0T+$=F2 7$$\"3w**\\ilAFjSF2$\"3\"yw3Dz,WC#F27$$\"3yLLL$)*pp;%F2$\"3%p#yu_LyREF 27$$\"3)RL$3xe,tUF2$\"3g:_Qt(yo,$F27$$\"3Cn;HdO=yVF2$\"31?#f=#*eeL$F27 $$\"3MMe9\"z-lU%F2$\"3#**)=6CtFbMF27$$\"3a+++D>#[Z%F2$\"3#)pALcM#Gb$F2 7$$\"3ZnT5::^-XF2$\"3Imd$=:Iuf$F27$$\"3TM$3_5,-`%F2$\"3QWfnHR*Gj$F27$$ \"3O,DJ&p!*yb%F2$\"3;qU3Z5XeOF27$$\"3SnmT&G!e&e%F2$\"3W)GP`]6Ln$F27$$ \"3fLe*[=Y.h%F2$\"3i?S$3&Q(on$F27$$\"3m+]P%37^j%F2$\"3%)yWF&pB1n$F27$$ \"3unT&Q)z()fYF2$\"3cjL-:5%Rl$F27$$\"3#RLLL)Qk%o%F2$\"3VQ%*\\:\")=EOF2 7$$\"3-n\"z>6but%F2$\"3!olG@U&*p_$F27$$\"37+]iSjE!z%F2$\"3'>n9ewDyO$F2 7$$\"3L+++DM\"3%[F2$\"3O'HfcaJG:$F27$$\"3a+]P40O\"*[F2$\"31d;^O,)*pGF2 7$$\"3s+voa-oX\\F2$\"3slP&>#4d#[#F27$$\"\"&F)$\"35s)***********>F2-%'C OLOURG6&%$RGBGF(F($\"*++++\"!\")-F$6$7goF'7$F4F*7$F?F*7$FIF*7$FaoF*7$F [pF*7$FepF*7$FjpF*7$F_qF*7$$\"3wmmTg()4_))F6F*7$FdqF*7$$\"3=nT&)3\\m_' *F6F*7$$\"3+N$ekGkX#**F6F*7$$\"31]iSmjk>5F2$!3!))*\\(=nsqg*F67$Fiq$!3; jmm;RFj!*F67$$\"3E+]7GTt%4\"F2$!3y%***\\PuJ0\")F67$F^r$!3'3LL$e4OZrF67 $Fcr$!3[(*****\\n\\!*\\F67$Fhr$!3V)*****\\ixCGF67$F]s$!3Y\"******\\KqP (F/7$Fbs$\"3]om;z*ev:\"F67$Fgs$\"3GNLL$347T$F67$F\\t$\"3tOLLLjM?`F67$F ft$\"3$3++D\"o7TvF67$F[u$\"3Qmm\"HK5S_)F67$F`u$\"3$>LLL$Q*o]*F67$$\"3_ ek.HW#)))>F2$\"39q\"H2e)[w(*F67$$\"3?]iSmTI-?F2$\"3/)\\PfLep(**F67$$\" 3*=/wP!Ry:?F2$\"3:\"eRA'4;U)*F67$Feu$\"3Ek;a)ejtq*F67$$\"3'pT&)e6Bi0#F 2$\"3[Ie9T)oxV*F67$Fju$\"3q'**\\P4u\"o\"*F67$$\"3/++D\"=EX8#F2$\"3e*** *\\(=QZl)F67$F_v$\"3*z**\\7G-89)F67$Fdv$\"3%)GL$3Fp)pqF67$Fiv$\"3YKL3- $ff3'F67$F^w$\"38nm;z%zY-&F67$Fcw$\"35kmT5!3B#RF67$Fhw$\"3C***\\iS!piH F67$F]x$\"3lim;/rFE>F67$$\"3)om\"H2\"34'GF2$\"35JL3F*=4R\"F67$Fbx$\"3Q &******\\2cb)F/7$$\"3S+Dc,6jSHF2$\"3af*\\P%)*)o$fF/7$$\"3K+]7y#=o'HF2$ \"37o**\\(=s\"=LF/7$$\"3G]iSm=\"*zHF2$\"3Ts\\PfL\")3?F/7$$\"3C+voaa+$* HF2$\"3'pw*\\7`a%*pF[y7$$\"3AD\"G)[Ab**HF2$\"37%)yu=<^xW!#@7$$\"3>](oH /*41IF2$\"3.Q+v$f3)>7F/7$$\"3F27$$\"3Y\\7`p2 _#)RF2$\"3!*)\\i!R:/l>F27$$\"3*G$e*[$zV4SF2$\"\"#F)7$$\"3L;/E+^NOSF2F \\`m7$Fd\\lF\\`m7$$\"3ym\"zW7@^6%F2F\\`m7$Fi\\lF\\`m7$F^]lF\\`m7$Fc]lF \\`m7$F]^lF\\`m7$Fa_lF\\`m7$Fe`lF\\`m7$F_alF\\`m7$FialF\\`m7$FcblF\\`m -Fhbl6&FjblF[clF(F(-%&TITLEG6#Q`oInterpolation~polynomial~compared~wit h~piecewise~linear~interpolation6\"-%+AXESLABELSG6$Q\"xFdamQ!6\"-%*THI CKNESSG6#\"\"$-%%VIEWG6$;F(Fcbl%(DEFAULTG" 1 2 0 1 10 3 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "Curve 2" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 66 "As you can see \+ the oscillation is worse than for the cubic spline." }}{PARA 0 "" 0 " " {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}} {MARK "61 2 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }