{VERSION 6 0 "IBM INTEL NT" "6.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 Comment" 2 18 "" 0 1 0 0 0 0 0 0 0 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 "Blue Emphasis" -1 256 "Times" 0 0 0 0 255 1 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "Green Emphasis" -1 257 "Times" 1 12 0 128 0 1 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "Red Emphasis" -1 258 "Times" 1 12 255 0 0 1 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "Maroon Emphasis" -1 259 "Times" 1 12 128 0 128 1 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "Dark Red Emphasis" -1 260 "Times" 1 12 128 0 0 1 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 261 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 262 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "Grey Emphasis" -1 263 "Times" 1 12 96 52 84 1 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "Purple Emphasis" -1 264 "Times" 1 12 115 0 230 1 0 1 0 0 0 0 0 0 0 1 }{CSTYLE "" -1 265 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 266 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 267 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 268 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 269 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 270 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "Purple Emphasis" -1 271 "Times" 1 12 102 0 230 1 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 272 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 273 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 274 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 275 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 276 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 277 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 278 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 279 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 280 "" 1 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 281 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 282 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 } {CSTYLE "" -1 283 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 284 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" 259 285 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" 260 286 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" 259 287 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 } {CSTYLE "" 259 288 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" 260 289 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" 260 290 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" 260 291 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 }{CSTYLE "" 259 292 "" 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 } {CSTYLE "" -1 293 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 294 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{CSTYLE "Orange Emphasis" -1 295 "Times" 1 12 225 100 10 1 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "Magenta Em phasis" -1 296 "Times" 1 12 200 0 200 1 0 1 0 0 0 0 0 0 0 0 }{CSTYLE " " -1 297 "" 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 }{PSTYLE "Normal" -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 "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "Times " 1 18 0 0 128 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 4 3 0 3 0 2 2 0 1 } {PSTYLE "Heading 2" -1 4 1 {CSTYLE "" -1 -1 "Times" 1 14 128 0 0 1 2 1 2 2 2 2 1 1 1 1 }1 1 0 0 8 2 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 "Maple 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 "Bulle t Item" -1 15 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 3 3 1 0 1 0 2 2 15 2 }{PSTYLE "Normal" -1 256 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 }} {SECT 0 {PARA 3 "" 0 "" {TEXT -1 44 "Using Taylor series to approximat e integrals" }}{PARA 0 "" 0 "" {TEXT -1 37 "by Peter Stone, Nanaimo, B .C., Canada" }}{PARA 0 "" 0 "" {TEXT -1 20 "Version: 26.3.2007 " }} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 72 "Introductory example of using a Taylor \+ series to approximate an integral" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 22 "Consider the function " }{XPPEDIT 18 0 "f (x) = sqrt(1+2*sin(x));" "6#/-%\"fG6#%\"xG-%%sqrtG6#,&\"\"\"F,*&\"\"#F ,-%$sinG6#F'F,F," }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 88 "f := x -> sqrt(1+2*sin(x)): \n'f(x)'=f(x);\nplot(f(x),x=-0.6..Pi,y=0..1.8,tickmarks=[4,4]);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*$,&\"\"\"F**&\"\"#F*-%$ sinGF&F*F*#F*F," }}{PARA 13 "" 1 "" {GLPLOT2D 428 255 255 {PLOTDATA 2 "6&-%'CURVESG6$7Y7$$!3%*******>=TN_!#=$\"3?H?YOu#y)**!#?7$$!3o******\\ b#*4_F*$\"3&[7v)*>L7s'!#>7$$!3?+++w#RW=&F*$\"3[m'o21fgX*F37$$!3_+++?b. 2]F*$\"3?\"R\">3m&y*>F*7$$!3))*****fwJ'H[F*$\"3Cv(f3#z7oEF*7$$!3G+++6! GAl%F*$\"3ABYc4(fb?$F*7$$!37+++cU#[Z%F*$\"353>7iJ')oOF*7$$!3!)*****H51 e2%F*$\"3a^w&4$z:_XF*7$$!3-+++]zywOF*$\"3O$H_Q5!)=I&F*7$$!3')*****f%[7 vKF*$\"3iMS%y6)yrfF*7$$!3=+++T@2FL\"F*$\"3uo-Z>1#)o&)F*7$$!3!)******* 4&*=l&F3$\"3&4yz(pf==%*F*7$$\"37+++]Ry&G#F3$\"3/')4se/gA5!#<7$$\"3-+++ w1!)>5F*$\"3a?T&o9!4(4\"F_p7$$\"3%******f/qO$=F*$\"33*[e4='>o6F_p7$$\" 37+++Im_]DF*$\"3YDvE:%=mA\"F_p7$$\"3++++dt`dLF*$\"3w9yP>p+)G\"F_p7$$\" 3/+++Y<'y;%F*$\"3[_i^i9BX8F_p7$$\"31+++Y%f([\\F*$\"3S_Y,#\\ojR\"F_p7$$ \"3e*****fO*)yl&F*$\"3!)fL)fG,&R9F_p7$$\"3-+++OM6,lF*$\"3!)>ypv=z'[\"F _p7$$\"33+++E[V:sF*$\"3;?B'*>>^B:F_p7$$\"3=+++;/OY!)F*$\"3M>k=I,Ui:F_p 7$$\"3)******fNq=y)F*$\"39['3&[dZ$f\"F_p7$$\"3?+++O8%))e*F*$\"3^(><;p: Ri\"F_p7$$\"31+++bosN5F_p$\"3AM&zqk$Q\\;F_p7$$\"3.+++sT!f6\"F_p$\"35o^ BL)3Bn\"F_p7$$\"31+++z>`*=\"F_p$\"3O*>o)\\d3!p\"F_p7$$\"3++++9&\\*o7F_ p$\"3G+>mo&[dq\"F_p7$$\"33+++9?W^8F_p$\"3!=G?Stg\"=DB9 F_p$\"3JgI^LlwDz<$f;F_p$\"3m$H^Wt)yHp$>=F_p$\"3&3`P7#Q@ 9F_p$\"3 Oa8#))e!p%o\"F_p7$$\"3!)*****H9K!\\?F_p$\"3%*Q1QI\"=gm\"F_p7$$\"3++++R r1H@F_p$\"3T+lIK]0U;F_p7$$\"34+++MQP/AF_p$\"3+\"*H(4^:hh\"F_p7$$\"3%** *****G54$G#F_p$\"3u`d$\\8uae\"F_p7$$\"3/+++u90gBF_p$\"3U]qsAV,_:F_p7$$ \"3>+++H>iSCF_p$\"3N-;%G0[K^\"F_p7$$\"3z*****>&3A=DF_p$\"3!yNxp!4Es9F_ p7$$\"3/+++gpd(f#F_p$\"3!euN%HweE9F_p7$$\"3=+++efFwEF_p$\"3#o**>*ycYx8 F_p7$$\"3#)*****\\;#f[FF_p$\"3I>nrS1))G8F_p7$$\"32+++)pu9$GF_p$\"3Ip)p !)[#**o7F_p7$$\"3y*****R#eg0HF_p$\"3/AIjXvW67F_p7$$\"37+++m]k%)HF_p$\" 3'z%R['4\"pX6F_p7$$\"3y*****fq&HgIF_p$\"3!GM[aW`\"y5F_p7$$\"3!)***\\/l #fTJF_p$\"3J$zRJ++++\"F_p-%'COLOURG6&%$RGBG$\"#5!\"\"$\"\"!F`]lF_]l-%+ AXESLABELSG6$Q\"x6\"Q\"yFe]l-%*AXESTICKSG6$\"\"%Fj]l-%%VIEWG6$;$!\"'F^ ]l$\"+aEfTJ!\"*;F_]l$\"#=F^]l" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 44.000000 0 0 "Curve 1" }}}}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 40 "Maple's explicit analytical formula for \+ " }{XPPEDIT 18 0 "Int(sqrt(1+2*sin(x)),x);" "6#-%$IntG6$-%%sqrtG6#,&\" \"\"F**&\"\"#F*-%$sinG6#%\"xGF*F*F0" }{TEXT -1 55 " is rather complica ted and involves elliptic functions." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "Int(sqrt(1+2*sin(x)),x);\n value(%);\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$*$-%%sqrtG6#,& \"\"\"F+-%$sinG6#%\"xG\"\"#F+F/" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$* &*.-%%sqrtG6#*$)-%$cosG6#%\"xG\"\"#\"\"\"F0-%%csgnG6#F+F0-F'6#,&-%$sin GF-F/F/F0F0-F'6#,&F7!\"#F/F0F0-F'6#,&F7F " 0 "" {MPLTEXT 1 0 109 "f := x -> sqrt(1+2*sin(x)): \n'f(x)'=f(x);\nconvert(taylor(f(x),x=0,8),polynom):\np := unapply(%,x ):\n'p(x)'=p(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*$,& \"\"\"F**&\"\"#F*-%$sinGF&F*F*#F*F," }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #/-%\"pG6#%\"xG,2\"\"\"F)F'F)*&#F)\"\"#F)*$)F'F,F)F)!\"\"*&#F)\"\"$F)* $)F'F2F)F)F)*&#\"#6\"#CF)*$)F'\"\"%F)F)F/*&#\"#>\"#IF)*$)F'\"\"&F)F)F) *&#\"$h'\"$?(F)*$)F'\"\"'F)F)F/*&#\"$P%\"$:$F)*$)F'\"\"(F)F)F)" }}} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 22 "The Taylo r polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 34 " \+ can be compared graphically with " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#% \"xG" }{TEXT -1 2 ". " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "plot([f(x),p(x)],x=-0.7..1,y=-.1..1.8,col or=[red,blue],legend=[`f(x)`,`p(x)`]);" }}{PARA 13 "" 1 "" {GLPLOT2D 479 340 340 {PLOTDATA 2 "6&-%'CURVESG6%7P7$$!3I+++&*[aF_!#=$\"3p<*=VeH X#Q!#>7$$!3()*****z2$>;_F*$\"3o)[48PVq&eF-7$$!3m*******p*zZ]F*$\"32c7) e6c-\"=F*7$$!3C+++=jSz[F*$\"3eL/!y>&o(\\#F*7$$!3.+++PJoIXF*$\"3nEx*yLE \"HNF*7$$!3%)*****H7L+<%F*$\"3?!f#p\"R&QeVF*7$$!35+++>(R0\"QF*$\"3c;,m p!R;1&F*7$$!3$******\\cc2W$F*$\"30'>S9(p\"Rq&F*7$$!3&)*****HO^]6$F*$\" 3\"Hs$[6h1@iF*7$$!3=+++ZWQ[FF*$\"3cpdtqIyhnF*7$$!33+++i>@!Q#F*$\"3!))G Nw$oQpsF*7$$!37+++_&4a-#F*$\"3.&)pT!p*)4t(F*7$$!38+++t\\@.:2=dU&*Q6F^q7$$\"3/ +++x_Tm=F*$\"39j!pX;[4<\"F^q7$$\"3))*****HE&o#>#F*$\"3-:s+zz#z>\"F^q7$ $\"3C+++&yl]a#F*$\"3(RL)*e^(=E7F^q7$$\"3/+++Y$4\"4HF*$\"33ptou?Xa7F^q7 $$\"3++++]\\DlKF*$\"3)QHpk_87G\"F^q7$$\"3/+++IT%)4OF*$\"3`Pw8?*)G18F^q 7$$\"3%)******f;X#*RF*$\"3Um,K)*)3KL\"F^q7$$\"36+++]:COVF*$\"3)HMSLj&e c8F^q7$$\"3@+++S1J.ZF*$\"3M**Ri&*4r!Q\"F^q7$$\"3'********yHf.&F*$\"3gy /#y[Q=S\"F^q7$$\"3s******p.d*R&F*$\"3jr[wl#[TU\"F^q7$$\"3l******HysTdF *$\"3ft_!=N'RW9F^q7$$\"3c******\\4Q*4'F*$\"3%oB#e(3)zk9F^q7$$\"3)***** **fA0\\kF*$\"3/y&p)Hk*R[\"F^q7$$\"3%)*******pD^\"oF*$\"3\\de\"F^q7$$\"3K++++#*yF*)F*$\"3e@ *)Rt4E*f\"F^q7$$\"3c+++gb!pG*F*$\"3Z&*=B%>qHh\"F^q7$$\"3x******HdiI'*F *$\"3WS4\"Rn)QD;F^q7$$\"\"\"\"\"!$\"3$y(z*[&)ozj\"F^q-%'COLOURG6&%$RGB G$\"*++++\"!\")$FiyFiyFcz-%'LEGENDG6#%%f(x)G-F$6%7S7$$!3a************* *pF*$!3?kmm\"p?3)\\F*7$$!3&emm\"z$[%HmF*$!3g/)zL^LX0$F*7$$!3RK$3_RLqI' F*$!3$>:3j(y9=;F*7$$!3RmmTDSWWfF*$!3zW=dDg*zA#F-7$$!3clm\"z<^%zbF*$\"3 )=C5v-I)z)*F-7$$!3=L$3x2$>;_F*$\"3\"H0([&)RsN?F*7$$!3Sm;a=jSz[F*$\"3)4 KwNB)o#*GF*7$$!3-+](o8$oIXF*$\"3mF?QBb:&o$F*7$$!3+m;aBJ.qTF*$\"3F%*3P7 (R0\"QF*$\"3us?trB@(3&F*7$$!3yKL$ecc2W$F*$\"3:WU!o(o e8dF*7$$!3^mm\"HO^]6$F*$\"3wO!HprX\\A'F*7$$!35****\\ZWQ[FF*$\"3#H**3j4 ^Iw'F*7$$!3c****\\i>@!Q#F*$\"3W=Qt&e[(psF*7$$!3K****\\_&4a-#F*$\"3^LG* [d!3JxF*7$$!3%omTN(\\@.-8 ;(zC35F^q7$$\"3#G+]PCs&>VF-$\"3.cs3Z\")GU5F^q7$$\"3'[nm\"zWWizF-$\"3wu v,\"41m2\"F^q7$$\"3en;HPQxI6F*$\"3$**HrLF,r5\"F^q7$$\"3'HL$3x*3;\\\"F* $\"3ID@R&oa*Q6F^q7$$\"3/L$ekF:k'=F*$\"33P3R?1&4<\"F^q7$$\"3M+](=E&o#># F*$\"3A(Qx!\\l$z>\"F^q7$$\"3aLLe%yl]a#F*$\"3&*Q_JHY@E7F^q7$$\"3S***** \\M4\"4HF*$\"3!Gq@T(z_a7F^q7$$\"3C++DY\\DlKF*$\"3%p72xk(R\"G\"F^q7$$\" 3L**\\7GT%)4OF*$\"3d!*p'*ebo18F^q7$$\"3s++vj;X#*RF*$\"3&fzx$yS1M8F^q7$ $\"3Ommm^:COVF*$\"3i#[ziT)=e8F^q7$$\"3+++]P1J.ZF*$\"3u^Rv)ozOQ\"F^q7$$ \"3+L$ekyHf.&F*$\"3N6h*GSF^q7$$\"3S+]7eb!pG*F*$\"3#>)zOC,g)4#F^q7$$\"3$***\\(=tD1j*F*$\"3[th_ \"\\')4E#F^q7$Fgy$\"3%\\Oz]OzvZ#F^q-F]z6&F_zFczFczF`z-Fez6#%%p(x)G-%+A XESLABELSG6$Q\"x6\"Q\"yFgjl-%%VIEWG6$;$!\"(!\"\"Fgy;$F_[mF_[m$\"#=F_[m " 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 1 "f(x)" "p(x) " }}}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 39 "Now we integrate the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6 #%\"xG" }{TEXT -1 22 " to give a polynomial " }{XPPEDIT 18 0 "q(x)" "6 #-%\"qG6#%\"xG" }{TEXT -1 89 " of degree 8 which will be an approximat ion for the (complicated) indefinite integral of " }{XPPEDIT 18 0 "f(x )" "6#-%\"fG6#%\"xG" }{TEXT -1 19 " given previously. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "Int(p(x ),x);\nq := unapply(value(%),x):\n'q(x)'=q(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$,2\"\"\"F'%\"xGF'*&#F'\"\"#F'*$)F(F+F'F'!\"\"* &#F'\"\"$F'*$)F(F1F'F'F'*&#\"#6\"#CF'*$)F(\"\"%F'F'F.*&#\"#>\"#IF'*$)F (\"\"&F'F'F'*&#\"$h'\"$?(F'*$)F(\"\"'F'F'F.*&#\"$P%\"$:$F'*$)F(\"\"(F' F'F'F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"qG6#%\"xG,2F'\"\"\"*&#F )\"\"#F)*$)F'F,F)F)F)*&#F)\"\"'F)*$)F'\"\"$F)F)!\"\"*&#F)\"#7F)*$)F'\" \"%F)F)F)*&#\"#6\"$?\"F)*$)F'\"\"&F)F)F5*&#\"#>\"$!=F)*$)F'F1F)F)F)*&# \"$h'\"%S]F)*$)F'\"\"(F)F)F5*&#\"$P%\"%?DF)*$)F'\"\")F)F)F)" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 6 "Since " }{XPPEDIT 18 0 "q(0)=0" "6#/-%\"qG6#\"\"!F'" }{TEXT -1 28 ", the value of the function " } {XPPEDIT 18 0 "q(x);" "6#-%\"qG6#%\"xG" }{TEXT -1 13 " at a number " } {TEXT 261 1 "x" }{TEXT -1 20 " is an estimate for " }{XPPEDIT 18 0 "In t(sqrt(1+2*sin(t)),t = 0 .. x);" "6#-%$IntG6$-%%sqrtG6#,&\"\"\"F**&\" \"#F*-%$sinG6#%\"tGF*F*/F0;\"\"!%\"xG" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 47 "We should obtain a reasonably good estimate if " } {TEXT 262 1 "x" }{TEXT -1 15 " is close to 0." }}{PARA 0 "" 0 "" {TEXT -1 28 "For example, let's estimate " }{XPPEDIT 18 0 "Int(sqrt(1+ 2*sin(t)),t = 0 .. 1/8);" "6#-%$IntG6$-%%sqrtG6#,&\"\"\"F**&\"\"#F*-%$ sinG6#%\"tGF*F*/F0;\"\"!*&F*F*\"\")!\"\"" }{TEXT -1 1 "." }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "q(1/8 );\nevalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"*D:\\C\"\"*'4C&R* " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+s([]K\"!#5" }}}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 264 4 "Note" }{TEXT -1 47 ": W e could simply set up the definite integral " }{XPPEDIT 18 0 "Int(p(x) ,x=0..1/8)" "6#-%$IntG6$-%\"pG6#%\"xG/F);\"\"!*&\"\"\"F.\"\")!\"\"" } {TEXT -1 83 ", evaluate this integral analytically, and then obtain a \+ numerical result by using " }{TEXT 0 5 "evalf" }{TEXT -1 3 ". " }} {PARA 0 "" 0 "" {TEXT -1 109 "When this shorter method is used, it is \+ important to obtain the analytical expression first by including the \+ " }{TEXT 258 8 "value(%)" }{TEXT -1 10 " command." }}{PARA 0 "" 0 "" {TEXT -1 87 "If this is not done, Maple will perform numerical integra tion of the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG " }{TEXT -1 77 " which defeats the intention of avoiding such methods \+ in the current context." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "Int(p(x),x=0..1/8);\nvalue(%);\neva lf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$,2\"\"\"F'%\"xGF'*& #F'\"\"#F'*$)F(F+F'F'!\"\"*&#F'\"\"$F'*$)F(F1F'F'F'*&#\"#6\"#CF'*$)F( \"\"%F'F'F.*&#\"#>\"#IF'*$)F(\"\"&F'F'F'*&#\"$h'\"$?(F'*$)F(\"\"'F'F'F .*&#\"$P%\"$:$F'*$)F(\"\"(F'F'F'/F(;\"\"!#F'\"\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"*D:\\C\"\"*'4C&R*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #$\"+s([]K\"!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 99 "We can compare this value with that obtained from the exp licit analytical formula for the integral." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 62 "Int(sqrt(1+2*sin(x )),x=0..1/8);\nvalue(%);\nevalf(evalf[15](%));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$*$,&\"\"\"F(*&\"\"#F(-%$sinG6#%\"xGF(F(#F(F*/F .;\"\"!#F(\"\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,*-%*EllipticFG6$,$ *&#\"\"#\"\"$\"\"\"*&F+#F,F**&,&*&F*F,-%$sinG6##F,\"\")F,F,F,F,F,,&*&F *F,F2F,F,F*F,!\"\"F.F,F,,$*&F*F9F+F.F,F9-F%6$,$*(F+F9F*F.F+F.F,F:F,-%+ EllipticPiG6%F'#F+\"\"%F:F,-FA6%F>FCF:F9" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+d([]K\"!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 96 "Maple can also evaluate integrals numerically, without first constructing an analytical formula." }}{PARA 0 "" 0 "" {TEXT -1 29 "This is achieved by applying " }{TEXT 0 5 "evalf" }{TEXT -1 41 " to the unevaluated from of the integral." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 41 "Int(sqrt(1+2*sin(x )),x=0..1/8);\nevalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$* $,&\"\"\"F(*&\"\"#F(-%$sinG6#%\"xGF(F(#F(F*/F.;\"\"!#F(\"\")" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+d([]K\"!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 232 "Looking back at the value obtaine d from the truncated series, we see that the last two digits in the va lue obtained differ from the last two more recent values. If we includ e extra terms of the Taylor series we get better agreement." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 119 "f := x -> sqrt(1+2*sin(x));\nconvert(taylor(f(x),x=0,11),polynom);\nq : = unapply(int(%,x),x);\nevalf(evalf[15](q(0.125)));" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%\"fGf*6#%\"xG6\"6$%)operatorG%&arrowGF(-%%sqrtG6#,& \"\"\"F0*&\"\"#F0-%$sinG6#9$F0F0F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,8\"\"\"F$%\"xGF$*&#F$\"\"#F$*$)F%F(F$F$!\"\"*&#F$\"\"$F$*$)F%F. F$F$F$*&#\"#6\"#CF$*$)F%\"\"%F$F$F+*&#\"#>\"#IF$*$)F%\"\"&F$F$F$*&#\"$ h'\"$?(F$*$)F%\"\"'F$F$F+*&#\"$P%\"$:$F$*$)F%\"\"(F$F$F$*&#\"&\">()\"& ?.%F$*$)F%\"\")F$F$F+*&#\"&,$y\"&!oAF$*$)F%\"\"*F$F$F$*&#\")@XQ?\"(+)G OF$*$)F%\"#5F$F$F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"qGf*6#%\"xG6 \"6$%)operatorG%&arrowGF(,89$\"\"\"*&#F.\"\"#F.*$)F-F1F.F.F.*&#F.\"\"' F.*$)F-\"\"$F.F.!\"\"*&#F.\"#7F.*$)F-\"\"%F.F.F.*&#\"#6\"$?\"F.*$)F-\" \"&F.F.F:*&#\"#>\"$!=F.*$)F-F6F.F.F.*&#\"$h'\"%S]F.*$)F-\"\"(F.F.F:*&# \"$P%\"%?DF.*$)F-\"\")F.F.F.*&#\"&\">()\"'!)GOF.*$)F-\"\"*F.F.F:*&#\"& ,$y\"'+oAF.*$)F-\"#5F.F.F.*&#\")@XQ?\")+o\"*RF.*$)F-FCF.F.F:F(F(F(" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+d([]K\"!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 164 "An error arising from a \+ mathematical approximation, such as occurred when the previous degree \+ 8 Taylor polynomial was used to approximate the integral, is called a \+ " }{TEXT 264 16 "truncation error" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 62 "More examples of using Taylor series to approximate integrals " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 "; " }}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 10 "Example 1 " }}{PARA 0 "" 0 " " {TEXT 292 8 "Question" }{TEXT -1 2 ": " }}{PARA 0 "" 0 "" {TEXT -1 13 "(a) Estimate " }{XPPEDIT 18 0 "Int((exp(x)-1)/x,x = -1 .. 1);" "6# -%$IntG6$*&,&-%$expG6#%\"xG\"\"\"F,!\"\"F,F+F-/F+;,$F,F-F," }{TEXT -1 50 " using a Taylor polynomial of degree 7 centred at " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 16 " to approximate " }{XPPEDIT 18 0 "f(x) = (exp(x)-1)/x;" "6#/-%\"fG6#%\"xG*&,&-%$expG6#F'\"\"\"F-!\"\" F-F'F." }{TEXT -1 89 ", and compare this value with the value given by Maple's numerical integration procedure " }{TEXT 0 9 "evalf/Int" } {TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 42 "(b) Find the Taylor poly nomial centred at " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 5 " for " }{XPPEDIT 18 0 "f(x) = (exp(x)-1)/x;" "6#/-%\"fG6#%\"xG*&,&-%$e xpG6#F'\"\"\"F-!\"\"F-F'F." }{TEXT -1 49 " of minimal degree which, wh en used to calculate " }{XPPEDIT 18 0 "Int((exp(x)-1)/x,x = -1 .. 1); " "6#-%$IntG6$*&,&-%$expG6#%\"xG\"\"\"F,!\"\"F,F+F-/F+;,$F,F-F," } {TEXT -1 47 ", gives a value which is correct to 10 digits. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 291 8 "Solution" } {TEXT -1 2 ": " }}{PARA 0 "" 0 "" {TEXT -1 55 "First we construct a Ta ylor polynomial of degree 7 for " }{XPPEDIT 18 0 "f(x)=(exp(x)-1)/x" " 6#/-%\"fG6#%\"xG*&,&-%$expG6#F'\"\"\"F-!\"\"F-F'F." }{TEXT -1 7 " abou t " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT 264 4 "Note" }{TEXT -1 70 ": In order for Maple to construc t a Taylor polynomial of degree 7 for " }{XPPEDIT 18 0 "f(x)" "6#-%\"f G6#%\"xG" }{TEXT -1 48 ", a Taylor polynomial of degree 8 is needed fo r " }{XPPEDIT 18 0 "exp(x)" "6#-%$expG6#%\"xG" }{TEXT -1 44 ". For thi s reason it is necessary to supply " }{TEXT 0 6 "taylor" }{TEXT -1 34 " with the third (order) argument \"" }{TEXT 263 1 "9" }{TEXT -1 3 "\" . " }}{PARA 0 "" 0 "" {TEXT 258 7 "Warning" }{TEXT -1 224 ": If you re turn to re-execute these commands after constructing a higher degree T aylor polynomial from a higher order Maple series structure, you will \+ get a degree 8 Taylor polynomial with the third (order) argument set t o \"" }{TEXT 263 1 "9" }{TEXT -1 96 "\". A degree 7 Taylor polynomial \+ can then be constructed with the third (order) argument set to \"" } {TEXT 263 1 "8" }{TEXT -1 2 "\"." }}{PARA 0 "" 0 "" {TEXT -1 63 "This \+ happens because Maple \"remembers\" the higher order series." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 105 "f := x -> (exp(x)-1)/x:\n'f(x)'=f(x);\ntaylor(f(x),x=0,9):\np := unappl y(convert(%,polynom),x):\n'p(x)'=p(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*&,&-%$expGF&\"\"\"F,!\"\"F,F'F-" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"pG6#%\"xG,2\"\"\"F)*&#F)\"\"#F)F'F)F)*&#F)\"\"' F)*$)F'F,F)F)F)*&#F)\"#CF)*$)F'\"\"$F)F)F)*&#F)\"$?\"F)*$)F'\"\"%F)F)F )*&#F)\"$?(F)*$)F'\"\"&F)F)F)*&#F)\"%S]F)*$)F'F/F)F)F)*&#F)\"&?.%F)*$) F'\"\"(F)F)F)" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 70 "Now we integrate this polynomial to obtain a polynomial a pproximation " }{XPPEDIT 18 0 "q(x)" "6#-%\"qG6#%\"xG" }{TEXT -1 30 " \+ for the indefinite integral " }{XPPEDIT 18 0 "Int((exp(x)-1)/x,x);" " 6#-%$IntG6$*&,&-%$expG6#%\"xG\"\"\"F,!\"\"F,F+F-F+" }{TEXT -1 2 ". " } }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "Int('p(x)',x)=Int(p(x),x); \nvalue(rhs(%)):\nq := unapply(%,x):\n'q(x)'=q(x);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#/-%$IntG6$-%\"pG6#%\"xGF*-F%6$,2\"\"\"F.*&#F.\"\"#F.F *F.F.*&#F.\"\"'F.*$)F*F1F.F.F.*&#F.\"#CF.*$)F*\"\"$F.F.F.*&#F.\"$?\"F. *$)F*\"\"%F.F.F.*&#F.\"$?(F.*$)F*\"\"&F.F.F.*&#F.\"%S]F.*$)F*F4F.F.F.* &#F.\"&?.%F.*$)F*\"\"(F.F.F.F*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-% \"qG6#%\"xG,2F'\"\"\"*&#F)\"\"%F)*$)F'\"\"#F)F)F)*&#F)\"#=F)*$)F'\"\"$ F)F)F)*&#F)\"#'*F)*$)F'F,F)F)F)*&#F)\"$+'F)*$)F'\"\"&F)F)F)*&#F)\"%?VF )*$)F'\"\"'F)F)F)*&#F)\"&!GNF)*$)F'\"\"(F)F)F)*&#F)\"'gDKF)*$)F'\"\")F )F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 15 "The polynomial " } {XPPEDIT 18 0 "q(x)" "6#-%\"qG6#%\"xG" }{TEXT -1 46 " is an approximat ion for an antiderivative of " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG " }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 27 "To obtain an estimate for " }{XPPEDIT 18 0 "Int((exp(x)-1)/x,x = -1 .. 1);" "6#-%$IntG6$*& ,&-%$expG6#%\"xG\"\"\"F,!\"\"F,F+F-/F+;,$F,F-F," }{TEXT -1 21 ", we ne ed to compute " }{XPPEDIT 18 0 "q(1)-q(-1)" "6#,&-%\"qG6#\"\"\"F'-F%6# ,$F'!\"\"F+" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 39 "q(1)-q(-1);\nest := evalf(evalf[14] (%));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"'*\\'=\"&+#))" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$estG$\"+M6]9@!\"*" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 61 "We can obtain a more accurate v alue by numerical integration." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 68 "f := x -> (exp(x)-1)/x:\nare a := evalf(evalf[14](Int(f(x),x=-1..1)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%areaG$\"+^<]9@!\"*" }}}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT -1 27 "The absolute error is . . ." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "ab serr := abs(est-area);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'abserrG$ \"$<'!\"*" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 108 "This error can be reduced by increasing the degree of the Tayl or polynomial used to approximate thefunction " }{XPPEDIT 18 0 "f(x)=( exp(x)-1)/x" "6#/-%\"fG6#%\"xG*&,&-%$expG6#F'\"\"\"F-!\"\"F-F'F." } {TEXT -1 2 ". " }}{PARA 0 "" 0 "" {TEXT -1 43 "Using a degree 8 polyno mial to approximate " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 12 " gives . . ." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 231 "f := x -> (exp(x)-1)/x:\n'f(x)'=f(x);\ntay lor(f(x),x=0,10):\np := unapply(convert(%,polynom),x):\n'p(x)'=p(x);\n Int('p(x)',x);\nvalue(%):\nq := unapply(%,x):\n'q(x)'=q(x);\nq(1)-q(-1 );\nest2 := evalf(evalf[14](%));\nabserr2 := abs(est2-area);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*&,&-%$expGF&\"\"\"F,!\"\"F,F' F-" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"pG6#%\"xG,4\"\"\"F)*&#F)\" \"#F)F'F)F)*&#F)\"\"'F)*$)F'F,F)F)F)*&#F)\"#CF)*$)F'\"\"$F)F)F)*&#F)\" $?\"F)*$)F'\"\"%F)F)F)*&#F)\"$?(F)*$)F'\"\"&F)F)F)*&#F)\"%S]F)*$)F'F/F )F)F)*&#F)\"&?.%F)*$)F'\"\"(F)F)F)*&#F)\"'!)GOF)*$)F'\"\")F)F)F)" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$-%\"pG6#%\"xGF)" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#/-%\"qG6#%\"xG,4F'\"\"\"*&#F)\"\"%F)*$)F'\"\"#F) F)F)*&#F)\"#=F)*$)F'\"\"$F)F)F)*&#F)\"#'*F)*$)F'F,F)F)F)*&#F)\"$+'F)*$ )F'\"\"&F)F)F)*&#F)\"%?VF)*$)F'\"\"'F)F)F)*&#F)\"&!GNF)*$)F'\"\"(F)F)F )*&#F)\"'gDKF)*$)F'\"\")F)F)F)*&#F)\"(?fE$F)*$)F'\"\"*F)F)F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"*(Q^37\")+O:d" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%est2G$\"+Y<]9@!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%(abserr2G$\"\"&!\"*" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT -1 149 "There is no change in the estimated value for the integral if the degree of the Taylor polynomial is increased to 9 bec ause the extra terms added to " }{XPPEDIT 18 0 "q(1)" "6#-%\"qG6#\"\" \"" }{TEXT -1 5 " and " }{XPPEDIT 18 0 "q(-1)" "6#-%\"qG6#,$\"\"\"!\" \"" }{TEXT -1 9 ", namely " }{XPPEDIT 18 0 "``(1/36288000)*`.`*1^10 = \+ 1/36288000;" "6#/*(-%!G6#*&\"\"\"F)\")+!)GO!\"\"F)%\".GF)F)\"#5*&F)F)F *F+" }{TEXT -1 5 " and " }{XPPEDIT 18 0 "``(1/36288000)*`.`*(-1)^10 = \+ 1/36288000" "6#/*(-%!G6#*&\"\"\"F)\")+!)GO!\"\"F)%\".GF),$F)F+\"#5*&F) F)F*F+" }{TEXT -1 48 " respectively, are identical and so cancel out ( " }{XPPEDIT 18 0 "q(x)" "6#-%\"qG6#%\"xG" }{TEXT -1 23 " is an even fu nction). " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 231 "f := x -> (exp(x)-1)/x:\n'f(x)'=f(x);\ntaylor(f(x),x =0,11):\np := unapply(convert(%,polynom),x):\n'p(x)'=p(x);\nInt('p(x)' ,x);\nvalue(%):\nq := unapply(%,x):\n'q(x)'=q(x);\nq(1)-q(-1);\nest3 : = evalf(evalf[14](%));\nabserr3 := abs(est3-area);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#/-%\"fG6#%\"xG*&,&-%$expGF&\"\"\"F,!\"\"F,F'F-" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"pG6#%\"xG,6\"\"\"F)*&#F)\"\"#F)F' F)F)*&#F)\"\"'F)*$)F'F,F)F)F)*&#F)\"#CF)*$)F'\"\"$F)F)F)*&#F)\"$?\"F)* $)F'\"\"%F)F)F)*&#F)\"$?(F)*$)F'\"\"&F)F)F)*&#F)\"%S]F)*$)F'F/F)F)F)*& #F)\"&?.%F)*$)F'\"\"(F)F)F)*&#F)\"'!)GOF)*$)F'\"\")F)F)F)*&#F)\"(+)GOF )*$)F'\"\"*F)F)F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$-%\"pG6# %\"xGF)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"qG6#%\"xG,6F'\"\"\"*&# F)\"\"%F)*$)F'\"\"#F)F)F)*&#F)\"#=F)*$)F'\"\"$F)F)F)*&#F)\"#'*F)*$)F'F ,F)F)F)*&#F)\"$+'F)*$)F'\"\"&F)F)F)*&#F)\"%?VF)*$)F'\"\"'F)F)F)*&#F)\" &!GNF)*$)F'\"\"(F)F)F)*&#F)\"'gDKF)*$)F'\"\")F)F)F)*&#F)\"(?fE$F)*$)F' \"\"*F)F)F)*&#F)\")+!)GOF)*$)F'\"#5F)F)F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"*(Q^37\")+O:d" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%% est3G$\"+Y<]9@!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(abserr3G$\"\" &!\"*" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 34 "Taking a degree 10 polynomial for " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6# %\"xG" }{TEXT -1 16 " to approximate " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG 6#%\"xG" }{TEXT -1 10 ", so that " }{XPPEDIT 18 0 "Int(p(t),t=0..x)=q( x)" "6#/-%$IntG6$-%\"pG6#%\"tG/F*;\"\"!%\"xG-%\"qG6#F." }{TEXT -1 84 " has degree 11, leads to an estimate for the integral which is correct to 10 digits." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 231 "f := x -> (exp(x)-1)/x:\n'f(x)'=f(x);\ntaylor(f (x),x=0,12):\np := unapply(convert(%,polynom),x):\n'p(x)'=p(x);\nInt(' p(x)',x);\nvalue(%):\nq := unapply(%,x):\n'q(x)'=q(x);\nq(1)-q(-1);\ne st4 := evalf(evalf[14](%));\nabseer4 := abs(est4-area);" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*&,&-%$expGF&\"\"\"F,!\"\"F,F'F-" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"pG6#%\"xG,8\"\"\"F)*&#F)\"\"#F) F'F)F)*&#F)\"\"'F)*$)F'F,F)F)F)*&#F)\"#CF)*$)F'\"\"$F)F)F)*&#F)\"$?\"F )*$)F'\"\"%F)F)F)*&#F)\"$?(F)*$)F'\"\"&F)F)F)*&#F)\"%S]F)*$)F'F/F)F)F) *&#F)\"&?.%F)*$)F'\"\"(F)F)F)*&#F)\"'!)GOF)*$)F'\"\")F)F)F)*&#F)\"(+)G OF)*$)F'\"\"*F)F)F)*&#F)\")+o\"*RF)*$)F'\"#5F)F)F)" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#-%$IntG6$-%\"pG6#%\"xGF)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"qG6#%\"xG,8F'\"\"\"*&#F)\"\"%F)*$)F'\"\"#F)F)F)*&# F)\"#=F)*$)F'\"\"$F)F)F)*&#F)\"#'*F)*$)F'F,F)F)F)*&#F)\"$+'F)*$)F'\"\" &F)F)F)*&#F)\"%?VF)*$)F'\"\"'F)F)F)*&#F)\"&!GNF)*$)F'\"\"(F)F)F)*&#F) \"'gDKF)*$)F'\"\")F)F)F)*&#F)\"(?fE$F)*$)F'\"\"*F)F)F)*&#F)\")+!)GOF)* $)F'\"#5F)F)F)*&#F)\"*+[3R%F)*$)F'\"#6F)F)F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\",%%est4G$\"+^<]9@!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(abseer4G$ \"\"!F&" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 9 "Example 2" }}{PARA 0 "" 0 "" {TEXT 285 8 "Question" }{TEXT -1 2 ": " }}{PARA 0 " " 0 "" {TEXT -1 34 "(a) Construct a Taylor polynomial " }{XPPEDIT 18 0 "p(x);" "6#-%\"pG6#%\"xG" }{TEXT -1 66 " of degree 10 as a truncated Maclaurin series, which approximates " }{XPPEDIT 18 0 "f(x)=exp(x)/(x ^2+2*x+2)" "6#/-%\"fG6#%\"xG*&-%$expG6#F'\"\"\",(*$F'\"\"#F,*&F/F,F'F, F,F/F,!\"\"" }{TEXT -1 6 " when " }{TEXT 283 1 "x" }{TEXT -1 11 " is n ear 0." }}{PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 264 4 "Note" }{TEXT -1 60 ": You will need an order 11 Taylor series approximation for " } {XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 3 ". )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 33 "(b) Plot a graph of \+ the function " }{XPPEDIT 18 0 "f(x)=exp(x)/(x^2+2*x+2)" "6#/-%\"fG6#% \"xG*&-%$expG6#F'\"\"\",(*$F'\"\"#F,*&F/F,F'F,F,F/F,!\"\"" }{TEXT -1 6 " for " }{TEXT 284 1 "x" }{TEXT -1 9 " between " }{XPPEDIT 18 0 "-2 ;" "6#,$\"\"#!\"\"" }{TEXT -1 5 " and " }{XPPEDIT 18 0 "2;" "6#\"\"#" }{TEXT -1 60 " along with the graph of the Taylor polynomial found in \+ (a)." }}{PARA 0 "" 0 "" {TEXT -1 13 "(c) Estimate " }{XPPEDIT 18 0 "In t(exp(x)/(x^2+2*x+2),x=-1/2..1/2)" "6#-%$IntG6$*&-%$expG6#%\"xG\"\"\", (*$F*\"\"#F+*&F.F+F*F+F+F.F+!\"\"/F*;,$*&F+F+F.F0F0*&F+F+F.F0" }{TEXT -1 30 " using the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG 6#%\"xG" }{TEXT -1 15 " found in (a). " }}{PARA 0 "" 0 "" {TEXT -1 13 "(d) Find the " }{TEXT 263 8 "absolute" }{TEXT -1 5 " and " }{TEXT 263 14 "relative error" }{TEXT -1 107 " in the estimate found in (c) b y comparing with the value given by Maple's numerical integration proc edure " }{TEXT 0 9 "evalf/Int" }{TEXT -1 3 ". " }}{PARA 0 "" 0 "" {TEXT -1 49 "(e) Increase the degree of the Taylor polynomial " } {XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 20 ", used to estimat e " }{XPPEDIT 18 0 "Int(exp(x)/(x^2+2*x+2),x=-1/2..1/2)" "6#-%$IntG6$ *&-%$expG6#%\"xG\"\"\",(*$F*\"\"#F+*&F.F+F*F+F+F.F+!\"\"/F*;,$*&F+F+F. F0F0*&F+F+F.F0" }{TEXT -1 166 " until the value of the integral agree s to 10 digits with that obtained by using Maple's numerical integrati on. What is the minimum degree needed for this to occur?" }}{PARA 0 " " 0 "" {TEXT -1 2 "( " }{TEXT 264 4 "Note" }{TEXT -1 87 ": You may nee d to use a few guard digits to get agreement in the last decimal place . ) " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 286 8 "S olution" }{TEXT -1 2 ": " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT -1 41 "(a) We can construct a Taylor polynomial " } {XPPEDIT 18 0 "p(x);" "6#-%\"pG6#%\"xG" }{TEXT -1 25 " of degree 10 ce ntred at " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 16 " to appr oximate " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 6 " when \+ " }{TEXT 294 1 "x" }{TEXT -1 12 " is near 0. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 113 "f := x -> e xp(x)/(x^2+2*x+2):\n'f(x)'=f(x);\ntaylor(f(x),x,11);\nconvert(%,polyno m):\np := unapply(%,x):\n'p(x)'=p(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*&-%$expGF&\"\"\",(*$)F'\"\"#F+F+*&F/F+F'F+F+F/F+!\" \"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#+7%\"xG#\"\"\"\"\"#\"\"!#F&\"#7 \"\"$#!\"\"\"#;\"\"%#F&\"#S\"\"&#F&\"$W\"\"\"'#!#8\"$s'\"\"(#\"#h\"%SQ \"\")#!$h\"\"&?f#\"\"*#!$\\$\"'+;?\"#5-%\"OG6#F&\"#6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"pG6#%\"xG,4#\"\"\"\"\"#F**&#F*\"#7F**$)F'\"\"$F *F*F**&#F*\"#;F**$)F'\"\"%F*F*!\"\"*&#F*\"#SF**$)F'\"\"&F*F*F**&#F*\"$ W\"F**$)F'\"\"'F*F*F**&#\"#8\"$s'F**$)F'\"\"(F*F*F8*&#\"#h\"%SQF**$)F' \"\")F*F*F**&#\"$h\"\"&?f#F**$)F'\"\"*F*F*F8*&#\"$\\$\"'+;?F**$)F'\"#5 F*F*F8" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 67 "(b) The following picture shows the graph of the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 11 " (drawn in " }{TEXT 256 4 "blue" }{TEXT -1 26 ") along with the graph of " } {XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 11 " (drawn in " } {TEXT 258 3 "red" }{TEXT -1 3 "). " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 74 "plot([f(x),p(x)],x=-2..2,y=- 1..2,color=[red,blue],legend=[`f(x)`,`p(x)`]);" }}{PARA 13 "" 1 "" {GLPLOT2D 449 385 385 {PLOTDATA 2 "6&-%'CURVESG6%7S7$$!\"#\"\"!$\"37N1 $=;knw'!#>7$$!3MLLL$Q6G\">!#<$\"3G(GzLm8\\0)F-7$$!3bmm;M!\\p$=F1$\"3e# )*49Q#3o$*F-7$$!3MLLL))Qj^7$$!3wmm;C2G!e\"F1$\"3m5Ara-YS:F>7$$!3OLL$3yO5]\"F1$\"3cz7W: Ar\"y\"F>7$$!3&*****\\nU)*=9F1$\"39x@@9:Fe?F>7$$!3SLL$3WDTL\"F1$\"3s;; -(Ro$pBF>7$$!35++]d(Q&\\7F1$\"3/%G_%>\\M)p#F>7$$!3gmmmc4`i6F1$\"3SoS\" zOhk/$F>7$$!3KLLLQW*e3\"F1$\"3[*fE?pr7N$F>7$$!3w++++()>'***F>$\"3Ex](e c#>!o$F>7$$!3E++++0\"*H\"*F>$\"3=:\"GQ\\mI)RF>7$$!35++++83&H)F>$\"3%)3 *\\BD4%RUF>7$$!3\\LLL3k(p`(F>$\"3:[*z***y0PWF>7$$!3Anmmmj^NmF>$\"3^*\\ Ku\"H[EYF>7$$!3)zmmmYh=(eF>$\"3o%>x7$$!3+,++v#\\N)\\F>$\"3eaM Q/d$Q&[F>7$$!3commmCC(>%F>$\"3BQgT=>s;\\F>7$$!39*****\\FRXL$F>$\"37J\\ +X`bg\\F>7$$!3t*****\\#=/8DF>$\"3SnT)pZjS)\\F>7$$!3=mmm;a*el\"F>$\"3/` ><\\ur&*\\F>7$$!3komm;Wn(o)F-$\"3#4#4DtnT**\\F>7$$!3IqLLL$eV(>!#?$\"3K spd$*******\\F>7$$\"3)Qjmm\"f`@')F-$\"3\\RJ*)>2]+]F>7$$\"3%z****\\nZ)H ;F>$\"3P7^G^k>.]F>7$$\"3ckmm;$y*eCF>$\"3<1,7$$\"3f)******R^bJ$ F>$\"3jMY_#)*[Q-&F>7$$\"3'e*****\\5a`TF>$\"3o5hi=REW]F>7$$\"3'o****\\7 RV'\\F>$\"3Q<`q\\Zgr]F>7$$\"3Y'*****\\@fkeF>$\"3S'[)H)yA96&F>7$$\"3_IL LL&4Nn'F>$\"3)Q.!=e>=c^F>7$$\"3A*******\\,s`(F>$\"3)=<#=')3x8_F>7$$\"3 %[mm;zM)>$)F>$\"3O\\f;,.*\\F&F>7$$\"3M*******pfa<*F>$\"3a#oO$Q1/_`F>7$ $\"39HLLeg`!)**F>$\"3(eqL1?]WV&F>7$$\"3w****\\#G2A3\"F1$\"3VL%*R0s8JbF >7$$\"3;LLL$)G[k6F1$\"3t_$*))4YOOcF>7$$\"3#)****\\7yh]7F1$\"3!f$GfzP=e dF>7$$\"3xmmm')fdL8F1$\"35/AVW,7()eF>7$$\"3bmmm,FT=9F1$\"3joFXW)R6.'F> 7$$\"3FLL$e#pa-:F1$\"3G_0$*f@`'='F>7$$\"3!*******Rv&)z:F1$\"3DMe?@![2M 'F>7$$\"3ILLLGUYo;F1$\"3'R03e+v9`'F>7$$\"3_mmm1^rZ7$$\"34++]sI@K=F1$\"3ACExAeRDpF>7$$\"34++]2%)38>F1$\"3`Zc5\"yi49(F>7$ $\"\"#F*$\"3'3lI*)4c!*Q(F>-%'COLOURG6&%$RGBG$\"*++++\"!\")$F*F*Fb[l-%' LEGENDG6#%%f(x)G-F$6%7en7$F($\"3fan^@M=GkF17$$!3SLL$e%G?y>F1$\"3)G;7%o w8oeF17$$!3ymmm\"p0k&>F1$\"31xBKmzw]`F17$$!3&*****\\P&3Y$>F1$\"370;u& \\PM([F17$F/$\"3R8K+_njLWF17$$!31++v3-)[(=F1$\"3C)*Q,J8v]PF17$F5$\"3%p He'R4iiJF17$$!3#)***\\7Y\"H%z\"F1$\"3tXP!4,E6g#F17$F:$\"3f)R#oL#z@8#F1 7$$!3+?9F17$FE$\"3lF,.!\\Dq]* F>7$FJ$\"33FMTx#))3r'F>7$FO$\"3P#e'H>baT\\F>7$FT$\"33+<-/R:[RF>7$FY$\" 3A]J\"Qu0\"=NF>7$Fhn$\"3Whr#[e%oPMF>7$F]o$\"3l?zL9pzUNF>7$Fbo$\"3#*\\. ugNHfPF>7$Fgo$\"321$z=1xD,%F>7$F\\p$\"3]p@r39l\\UF>7$Fap$\"3B(z>[jv0W% F>7$Ffp$\"3GHtn?QKFYF>7$F[q$\"3qR9z\\ir\\ZF>7$F`q$\"3![JrL&)oQ&[F>7$Fe q$\"3za>X&oEn\"\\F>7$Fjq$\"3'[^h>qb0'\\F>7$F_r$\"3A(QI>\\jS)\\F>7$Fdr$ \"33\"*RJ\\ur&*\\F>7$Fir$\"3SJ5DtnT**\\F>7$F^s$\"3yrpd$*******\\F>7$Fd s$\"3c\\I*)>2]+]F>7$Fis$\"3oyQ>^k>.]F>7$F^t$\"3ribIhnL5]F>7$Fct$\"34qN $ey[Q-&F>7$Fht$\"3Q$)R/A7$F]u$\"375zgu+fr]F>7$Fbu$\"3rZSz0qL6^F> 7$Fgu$\"3AKA$=+Ze:&F>7$F\\v$\"3;nvzH&pD@&F>7$Fav$\"3sh!H(*><;F&F>7$Ffv $\"36lWb4gpU`F>7$F[w$\"3MUi$o%e47aF>7$F`w$\"3/hra\\+`zaF>7$Few$\"3?**R %oN\"eEbF>7$Fjw$\"3[j&Qy[:(HbF>7$F_x$\"3qGx3g<3YaF>7$Fdx$\"3/ifc'3(4._ F>7$Fix$\"3'3Lmy4%3(p%F>7$F^y$\"3_o,YPgBgQF>7$Fcy$\"31#e-AfIn@#F>7$Fhy $!3_%)Qy_;)3(=F-7$F]z$!3IUq[#*)3u>%F>7$$\"33+++S2ls=F1$!3=mq#Rl4H$oF>7 $Fbz$!3S#R]$GSj05F17$$\"3/++v.Uac>F1$!3Aa$eM:--V\"F17$Fgz$!3We8pC!e8&> F1-F\\[l6&F^[lFb[lFb[lF_[l-Fd[l6#%%p(x)G-%+AXESLABELSG6$Q\"x6\"Q\"yF^h l-%%VIEWG6$;F(Fgz;$!\"\"F*Fgz" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 1 "f(x)" "p(x)" }}}}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT -1 4 "(c) " }{XPPEDIT 18 0 "Int(p(x),x = -1/ 2 .. 1/2);" "6#-%$IntG6$-%\"pG6#%\"xG/F);,$*&\"\"\"F.\"\"#!\"\"F0*&F.F .F/F0" }{TEXT -1 26 " provides an estimate for " }{XPPEDIT 18 0 "Int(f (x),x = -1/2 .. 1/2);" "6#-%$IntG6$-%\"fG6#%\"xG/F);,$*&\"\"\"F.\"\"#! \"\"F0*&F.F.F/F0" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "Int('p(x)',x=-1/2..1/2);\nva lue(%);\narea_est := evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$I ntG6$-%\"pG6#%\"xG/F);#!\"\"\"\"##\"\"\"F." }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"+$)G1,M\"++sY7o" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %)area_estG$\"+?*4C*\\!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 41 "(d) We can obtain a numerical value for " } {XPPEDIT 18 0 "Int(f(x),x = -1/2 .. 1/2);" "6#-%$IntG6$-%\"fG6#%\"xG/F );,$*&\"\"\"F.\"\"#!\"\"F0*&F.F.F/F0" }{TEXT -1 113 " by Maple's numer ical integration. Using a few guard digits will ensure that the result is correct to 10 digits. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 86 "f := x -> exp(x)/(x^2+2*x+2):\naccu rate_val := evalf(evalf[14](Int(f(x),x=-1/2..1/2)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%-accurate_valG$\"+,#4C*\\!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 114 "We can now calculate the absolute and relative errors in the value obtained by integrating the Taylor polynomial. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 82 "abserr := abs(area_est-accurate_val);\nre lerr := evalf(abserr/abs(accurate_val),4);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'abserrG$\"$>(!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%'relerrG$\"%S9!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 " " {TEXT -1 31 "(e) Taking a Taylor polynomial " }{XPPEDIT 18 0 "q(x)" "6#-%\"qG6#%\"xG" }{TEXT -1 25 " of degree 16 centred at " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 16 " to approximate " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 63 " gives a value for the int egral which is correct to 10 digits. " }}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 158 "f := x -> exp(x)/(x^2+2 *x+2):\n'f(x)'=f(x);\ntaylor(f(x),x,17):\nconvert(%,polynom):\nq := un apply(%,x):\n'q(x)'=q(x);\nInt('q(x)',x=-1/2..1/2);\nvalue(%);\nevalf( %);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*&-%$expGF&\"\"\" ,(*$)F'\"\"#F+F+*&F/F+F'F+F+F/F+!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#/-%\"qG6#%\"xG,@#\"\"\"\"\"#F**&#F*\"#7F**$)F'\"\"$F*F*F**&#F*\"#;F **$)F'\"\"%F*F*!\"\"*&#F*\"#SF**$)F'\"\"&F*F*F**&#F*\"$W\"F**$)F'\"\"' F*F*F**&#\"#8\"$s'F**$)F'\"\"(F*F*F8*&#\"#h\"%SQF**$)F'\"\")F*F*F**&# \"$h\"\"&?f#F**$)F'\"\"*F*F*F8*&#\"$\\$\"'+;?F**$)F'\"#5F*F*F8*&#\"%\" e)\"(!3u " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 9 "Example 3" }}{PARA 0 "" 0 "" {TEXT 287 8 "Question" }{TEXT -1 2 ": " }}{PARA 0 "" 0 "" {TEXT -1 34 "(a) Const ruct a Taylor polynomial " }{XPPEDIT 18 0 "p(x);" "6#-%\"pG6#%\"xG" } {TEXT -1 66 " of degree 10 as a truncated Maclaurin series, which appr oximates " }{XPPEDIT 18 0 "f(x)=(1+x*tanh(x))/(x^2-x+1)" "6#/-%\"fG6#% \"xG*&,&\"\"\"F**&F'F*-%%tanhG6#F'F*F*F*,(*$F'\"\"#F*F'!\"\"F*F*F2" } {TEXT -1 6 " when " }{TEXT 274 1 "x" }{TEXT -1 11 " is near 0." }} {PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 271 4 "Note" }{TEXT -1 60 ": Yo u will need an order 11 Taylor series approximation for " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 3 ". )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 33 "(b) Plot a graph of the f unction " }{XPPEDIT 18 0 "f(x)=(1+x*tanh(x))/(x^2-x+1)" "6#/-%\"fG6#% \"xG*&,&\"\"\"F**&F'F*-%%tanhG6#F'F*F*F*,(*$F'\"\"#F*F'!\"\"F*F*F2" } {TEXT -1 5 " for " }{TEXT 275 1 "x" }{TEXT -1 9 " between " }{XPPEDIT 18 0 "-2;" "6#,$\"\"#!\"\"" }{TEXT -1 5 " and " }{XPPEDIT 18 0 "2;" "6 #\"\"#" }{TEXT -1 60 " along with the graph of the Taylor polynomial f ound in (a)." }}{PARA 0 "" 0 "" {TEXT -1 13 "(c) Estimate " }{XPPEDIT 18 0 "Int((1+x*tanh(x))/(x^2-x+1),x=-1/3..1/3)" "6#-%$IntG6$*&,&\"\"\" F(*&%\"xGF(-%%tanhG6#F*F(F(F(,(*$F*\"\"#F(F*!\"\"F(F(F1/F*;,$*&F(F(\" \"$F1F1*&F(F(F6F1" }{TEXT -1 29 " using the Taylor polynomial " } {XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 15 " found in (a). " }}{PARA 0 "" 0 "" {TEXT -1 13 "(d) Find the " }{TEXT 263 14 "relative \+ error" }{TEXT -1 107 " in the estimate found in (c) by comparing with \+ the value given by Maple's numerical integration procedure " }{TEXT 0 9 "evalf/int" }{TEXT -1 3 ". " }}{PARA 0 "" 0 "" {TEXT -1 49 "(e) Inc rease the degree of the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#- %\"pG6#%\"xG" }{TEXT -1 20 ", used to estimate " }{XPPEDIT 18 0 "Int( (1+x*tanh(x))/(x^2-x+1),x=-1/3..1/3)" "6#-%$IntG6$*&,&\"\"\"F(*&%\"xGF (-%%tanhG6#F*F(F(F(,(*$F*\"\"#F(F*!\"\"F(F(F1/F*;,$*&F(F(\"\"$F1F1*&F( F(F6F1" }{TEXT -1 124 " until the value of the integral agrees to10 di gits with that obtained by using Maple's numerical integration. What \+ is the " }{TEXT 263 14 "minimum degree" }{TEXT -1 26 " needed for this to occur?" }}{PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 264 4 "Note" } {TEXT -1 86 ": You may need to use a few guard digits to get agreement in the last decimal place. )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT 289 8 "Solution" }{TEXT -1 1 ":" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 41 "(a) We can construct a Taylor polynomial " }{XPPEDIT 18 0 "p(x);" "6#-%\"pG6#%\"xG" }{TEXT -1 25 " of degree 10 centred at " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"! " }{TEXT -1 16 " to approximate " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\" xG" }{TEXT -1 6 " when " }{TEXT 293 1 "x" }{TEXT -1 12 " is near 0. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 118 "f := x -> (1+x*tanh(x))/(x^2-x+1):\n'f(x)'=f(x);\ntaylor(f(x),x ,11);\nconvert(%,polynom):\np := unapply(%,x):\n'p(x)'=p(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\"xG*&,&\"\"\"F**&F'F*-%%tanhGF&F* F*F*,(*$)F'\"\"#F*F*F'!\"\"F*F*F2" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#+ 9%\"xG\"\"\"\"\"!F%F%F%\"\"##!\"%\"\"$\"\"%F(\"\"&#F'\"#:\"\"'#\"#AF. \"\"(#\"$.%\"$:$\"\")#!#fF5\"\"*#!%'4%\"%NG\"#5-%\"OG6#F%\"#6" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"pG6#%\"xG,6\"\"\"F)F'F)*$)F'\"\"# F)F)*&#\"\"%\"\"$F)*$)F'F/F)F)!\"\"*&#F/F0F)*$)F'\"\"&F)F)F3*&#F,\"#:F )*$)F'\"\"'F)F)F)*&#\"#AF;F)*$)F'\"\"(F)F)F)*&#\"$.%\"$:$F)*$)F'\"\")F )F)F)*&#\"#fFHF)*$)F'\"\"*F)F)F3*&#\"%'4%\"%NGF)*$)F'\"#5F)F)F3" }}} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 67 "(b) The f ollowing picture shows the graph of the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 11 " (drawn in " }{TEXT 271 6 "purple" }{TEXT -1 26 ") along with the graph of " }{XPPEDIT 18 0 "f(x )" "6#-%\"fG6#%\"xG" }{TEXT -1 11 " (drawn in " }{TEXT 295 6 "orange" }{TEXT -1 3 "). " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 90 "plot([f(x),p(x)],x=-2..2,y=-1..2,color=[coral, COLOR(RGB,.4,0,.9)],legend=[`f(x)`,`p(x)`]);" }}{PARA 13 "" 1 "" {GLPLOT2D 554 334 334 {PLOTDATA 2 "6&-%'CURVESG6%7Y7$$!\"#\"\"!$\"3#Gw tI%f$H=%!#=7$$!3MLLL$Q6G\">!#<$\"3rc#o\\%y:3VF-7$$!3bmm;M!\\p$=F1$\"3y Prs$H*)4U%F-7$$!3MLLL))Qj^'***F-$\"3I/q::b rseF-7$$!3E++++0\"*H\"*F-$\"3`nc$3,]G/'F-7$$!35++++83&H)F-$\"3\"[k;,(Q C8iF-7$$!3\\LLL3k(p`(F-$\"3C5$z%>*pgP'F-7$$!3Anmmmj^NmF-$\"3!z*y*>RtZe 'F-7$$!3)zmmmYh=(eF-$\"3#=/2\\O4/y'F-7$$!3+,++v#\\N)\\F-$\"3DzO>)z8)Rq F-7$$!3commmCC(>%F-$\"3WM\"oGe7%4tF-7$$!39*****\\FRXL$F-$\"3WODef3[kwF -7$$!3t*****\\#=/8DF-$\"3S,Fm'[)Gy!)F-7$$!3=mmm;a*el\"F-$\"3%)yBq'[B*4 ')F-7$$!3komm;Wn(o)!#>$\"31*yW]v9g?*F-7$$!3IqLLL$eV(>!#?$\"3=[FP(R&H!) **F-7$$\"3)Qjmm\"f`@')Fjr$\"3Uj'4)Q&oN4\"F17$$\"3%z****\\nZ)H;F-$\"3%Q /(*=kVT+8F17$$\"3f)******R^bJ$F-$ \"3F7[.a;+@9F17$$\"3'e*****\\5a`TF-$\"3!R)*eSL2j`\"F17$$\"3'o****\\7RV '\\F-$\"3t]35$*pKP;F17$$\"3Y'*****\\@fkeF-$\"3[i*>\"G/YG$)F-$\"3IR6g*Hw;#=F17$$\"33KL$eCZwu)F-$\"3.%oVn@?W\"=F17$ $\"3M*******pfa<*F-$\"3>/aBr9D,=F17$$\"3Ckm;zy*zd*F-$\"3Q/(RnGOSy\"F17 $$\"39HLLeg`!)**F-$\"3GTi&o+>Fw\"F17$$\"3w****\\#G2A3\"F1$\"3gfIy([\"R 22V;F17$$\"3#)****\\7yh]7F1$\"3CF8Us (\\#p:F17$$\"3xmmm')fdL8F1$\"3k)yF:UB_\\\"F17$$\"3bmmm,FT=9F1$\"3'R$e+ &zB\">9F17$$\"3FLL$e#pa-:F1$\"3EShkv81X8F17$$\"3!*******Rv&)z:F1$\"31e z*QM$Gz7F17$$\"3ILLLGUYo;F1$\"3QKBMa;N27F17$$\"3_mmm1^rZF1$\"3Z?2V2H#3 .\"F17$$\"\"#F*$\"318@F1$!3PkhnGG\"[9\"F]^l7$$!3SLL$e%G?y>F1$!3Kfsa(=m>3\"F]^l7$$!3)* ***\\(oUIn>F1$!3ohwwh#RA-\"F]^l7$$!3ymmm\"p0k&>F1$!3-R`q8r$\\l*!#:7$$! 3&*****\\P&3Y$>F1$!3t$3)[(f,Ug)Fb_l7$F/$!3/`qEkWVdwFb_l7$$!3em;/'zXQ*= F1$!38!*\\!oyJ5\"pFb_l7$$!31++v3-)[(=F1$!3utd_7MpIiFb_l7$$!3IL$e9i9f&= F1$!3pu#3HMs6h&Fb_l7$F5$!3[+:7(e*eZ]Fb_l7$$!3#)***\\7Y\"H%z\"F1$!35I@V DqVhRFb_l7$F:$!3wI*)oA<3!4$Fb_l7$$!3\"esL1\"Fb _l7$FI$!3?o$3NoU)=i!#;7$FN$!3%Rgxeh.hW$Fhbl7$FS$!3%p2H[Rd&)y\"Fhbl7$FX $!3]ua``DId()F17$Fgn$!3a/3HC%[X\"QF17$F\\o$!3Gb*Q(H%yrc\"F17$Fao$!3qBi .y,)z2$F-7$Ffo$\"3N?L^gRN/EF-7$F[p$\"3q,h(G\"G-o\\F-7$F`p$\"30E5MmWaEf F-7$Fep$\"3'o<2%ybOpkF-7$Fjp$\"3m4#)e23H\\nF-7$F_q$\"3wl0\\9M\\MqF-7$F dq$\"3D%*Q;=Be3tF-7$Fiq$\"3i%39nt7Wm(F-7$F^r$\"3c@YDy`Gy!)F-7$Fcr$\"37 '*H`aM#*4')F-7$Fhr$\"3y7z,bZ,1#*F-7$F^s$\"3I\\FP(R&H!)**F-7$Fds$\"3#32 7)Q&oN4\"F17$Fis$\"3aw>Xu@Y)=\"F17$F^t$\"32Tb%[e9/I\"F17$Fct$\"3=_MP'* p+@9F17$Fht$\"3(o12Ntij`\"F17$F]u$\"3!4nz')[awj\"F17$Fbu$\"3q<'=sx0+t \"F17$Fgu$\"3GQM7i85!z\"F17$Fav$\"39V3[!Hct#=F17$Few$\"3lYD7*Q$3J=F17$ F_x$\"3%e4MWWw^x\"F17$Fix$\"3['e%\\5,(ze\"F17$F^y$\"3)\\xtI\\)R^5F17$F cy$!34-t[z?&=J#F-7$Fhy$!3S\"*R0a$ojE$F17$F]z$!3'o[W1D!4\"Q*F17$Fbz$!3# )R&>=D#4h@Fhbl7$Fgz$!3vW))fS!GoT%Fhbl7$F\\[l$!3([A./5=P+)Fhbl7$Fa[l$!3 Sd<&>TM#*\\\"Fb_l7$$\"3\"*****\\n'*33Fb_l7$Ff[l$!3[D< +HV)z@KuE%Fb_ l7$$\"35++D1>V_=F1$!3hV%[k)4![\"[Fb_l7$$\"33+++S2ls=F1$!3p7wn9\"HMU&Fb _l7$$\"33++vt&pG*=F1$!3UI'GUi[#*4'Fb_l7$F`\\l$!3Ct.iF5o[oFb_l7$$\"31+] i0j\"[$>F1$!3Z-\"*=CP0WxFb_l7$$\"3/++v.Uac>F1$!3!Q]+AG1=u)Fb_l7$$\"3/+ D\"G:3u'>F1$!3EV.aW3@#G*Fb_l7$$\"3-+](=5s#y>F1$!35'H0FDy?&)*Fb_l7$$\"3 -+v$40O\"*)>F1$!31A)Ro-\"GX5F]^l7$Fe\\l$!3rH*>hf'e36F]^l-%&COLORG6&F\\ ]l$\"\"%!\"\"Fb]l$\"\"*F^^m-Fd]l6#%%p(x)G-%+AXESLABELSG6$Q\"x6\"Q\"yFh ^m-%%VIEWG6$;F(Fe\\l;$F^^mF*Fe\\l" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 1 "f(x)" "p(x)" }}}}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT -1 4 "(c) " }{XPPEDIT 18 0 "Int(p(x),x=-1/3. .1/3)" "6#-%$IntG6$-%\"pG6#%\"xG/F);,$*&\"\"\"F.\"\"$!\"\"F0*&F.F.F/F0 " }{TEXT -1 26 " provides an estimate for " }{XPPEDIT 18 0 "Int(f(x),x = -1/3 .. 1/3);" "6#-%$IntG6$-%\"fG6#%\"xG/F);,$*&\"\"\"F.\"\"$!\"\"F 0*&F.F.F/F0" }{TEXT -1 2 ". " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "Int(p(x),x=-1/3..1/3);\nvalu e(%);\narea_est := evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$Int G6$,6\"\"\"F'%\"xGF'*$)F(\"\"#F'F'*&#\"\"%\"\"$F'*$)F(F.F'F'!\"\"*&#F. F/F'*$)F(\"\"&F'F'F2*&#F+\"#:F'*$)F(\"\"'F'F'F'*&#\"#AF:F'*$)F(\"\"(F' F'F'*&#\"$.%\"$:$F'*$)F(\"\")F'F'F'*&#\"#fFGF'*$)F(\"\"*F'F'F2*&#\"%'4 %\"%NGF'*$)F(\"#5F'F'F2/F(;#F2F/#F'F/" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\"+ACL2Q\"+&>HV_&" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)area_estG$ \"+vh$>*o!#5" }}}{PARA 0 "" 0 "" {TEXT -1 40 "(d) We can obtain a nume rical value for " }{XPPEDIT 18 0 "Int(f(x),x = -1/3 .. 1/3);" "6#-%$In tG6$-%\"fG6#%\"xG/F);,$*&\"\"\"F.\"\"$!\"\"F0*&F.F.F/F0" }{TEXT -1 113 " by Maple's numerical integration. Using a few guard digits will \+ ensure that the result is correct to 10 digits. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "accurate_val := evalf(Int(f(x),x=-1/3..1/3));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %-accurate_valG$\"+^k$>*o!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {PARA 0 "" 0 "" {TEXT -1 114 "We can now calculate the absolute and re lative errors in the value obtained by integrating the Taylor polynomi al. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 82 "abserr := abs(area_est-accurate_val);\nrelerr := eval f(abserr/abs(accurate_val),5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'a bserrG$\"$w#!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'relerrG$\"&Z+%!# 7" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 31 "(e) Taking a Taylor polynomial " }{XPPEDIT 18 0 "q(x)" "6#-%\"qG6#%\"xG" }{TEXT -1 25 " of degree 16 centred at " }{XPPEDIT 18 0 "x=0" "6#/%\"x G\"\"!" }{TEXT -1 16 " to approximate " }{XPPEDIT 18 0 "f(x)" "6#-%\"f G6#%\"xG" }{TEXT -1 63 " gives a value for the integral which is corre ct to 10 digits. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 148 "f := x -> (1+x*tanh(x))/(x^2-x+1);\ntaylor(f (x),x,17):\nconvert(%,polynom):\np := unapply(%,x);\nInt('p(x)',x=-1/3 ..1/3);\nevalf[14](value(%));\nevalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fGf*6#%\"xG6\"6$%)operatorG%&arrowGF(*&,&\"\"\"F.*&9$F.-%%ta nhG6#F0F.F.F.,(*$)F0\"\"#F.F.F0!\"\"F.F.F8F(F(F(" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%\"pGf*6#%\"xG6\"6$%)operatorG%&arrowGF(,B\"\"\"F-9$F- *$)F.\"\"#F-F-*&#\"\"%\"\"$F-*$)F.F4F-F-!\"\"*&#F4F5F-*$)F.\"\"&F-F-F8 *&#F1\"#:F-*$)F.\"\"'F-F-F-*&#\"#AF@F-*$)F.\"\"(F-F-F-*&#\"$.%\"$:$F-* $)F.\"\")F-F-F-*&#\"#fFMF-*$)F.\"\"*F-F-F8*&#\"%'4%\"%NGF-*$)F.\"#5F-F -F8*&#\"$8(\"$n&F-*$)F.\"#6F-F-F8*&#\"&By#\"'Df:F-*$)F.\"#7F-F-F-*&#\" ')*QAFboF-*$)F.\"#8F-F-F-*&#\"(p(ow\"(v53'F-*$)F.\"#9F-F-F-*&#\"(`K1\" F_pF-*$)F.F@F-F-F8*&#\"*z=z<*\"*vG^Q'F-*$)F.\"#;F-F-F8F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$-%\"pG6#%\"xG/F);#!\"\"\"\"$#\"\"\" F." }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"/*z4XO>*o!#9" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+^k$>*o!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 " " {TEXT -1 9 "Example 4" }}{PARA 0 "" 0 "" {TEXT 288 8 "Question" } {TEXT -1 2 ": " }}{PARA 0 "" 0 "" {TEXT -1 34 "(a) Construct a Taylor \+ polynomial " }{XPPEDIT 18 0 "p(x);" "6#-%\"pG6#%\"xG" }{TEXT -1 66 " o f degree 10 as a truncated Maclaurin series, which approximates " } {XPPEDIT 18 0 "f(x) = exp(x)/(1+sin(x^2));" "6#/-%\"fG6#%\"xG*&-%$expG 6#F'\"\"\",&F,F,-%$sinG6#*$F'\"\"#F,!\"\"" }{TEXT -1 6 " when " } {TEXT 276 1 "x" }{TEXT -1 11 " is near 0." }}{PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 271 4 "Note" }{TEXT -1 60 ": You will need an order 11 T aylor series approximation for " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"x G" }{TEXT -1 3 ". )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 " " {TEXT -1 33 "(b) Plot a graph of the function " }{XPPEDIT 18 0 "f(x) = exp(x)/(1+sin(x^2));" "6#/-%\"fG6#%\"xG*&-%$expG6#F'\"\"\",&F,F,-%$ sinG6#*$F'\"\"#F,!\"\"" }{TEXT -1 7 " , for " }{TEXT 277 1 "x" }{TEXT -1 9 " between " }{XPPEDIT 18 0 "-3/2;" "6#,$*&\"\"$\"\"\"\"\"#!\"\"F( " }{TEXT -1 5 " and " }{XPPEDIT 18 0 "3/2;" "6#*&\"\"$\"\"\"\"\"#!\"\" " }{TEXT -1 61 ", along with the graph of the Taylor polynomial found \+ in (a)." }}{PARA 0 "" 0 "" {TEXT -1 13 "(c) Estimate " }{XPPEDIT 18 0 "Int(exp(x)/(1+sin(x^2)),x = -1/3 .. 1/3);" "6#-%$IntG6$*&-%$expG6#%\" xG\"\"\",&F+F+-%$sinG6#*$F*\"\"#F+!\"\"/F*;,$*&F+F+\"\"$F2F2*&F+F+F7F2 " }{TEXT -1 30 " using the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 15 " found in (a). " }}{PARA 0 "" 0 "" {TEXT -1 13 "(d) Find the " }{TEXT 263 14 "relative error" }{TEXT -1 107 " in the estimate found in (c) by comparing with the value given b y Maple's numerical integration procedure " }{TEXT 0 9 "evalf/int" } {TEXT -1 3 ". " }}{PARA 0 "" 0 "" {TEXT -1 49 "(e) Increase the degre e of the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" } {TEXT -1 20 ", used to estimate " }{XPPEDIT 18 0 "Int(exp(x)/(1+sin(x ^2)),x = -1/3 .. 1/3);" "6#-%$IntG6$*&-%$expG6#%\"xG\"\"\",&F+F+-%$sin G6#*$F*\"\"#F+!\"\"/F*;,$*&F+F+\"\"$F2F2*&F+F+F7F2" }{TEXT -1 125 " u ntil the value of the integral agrees to10 digits with that obtained b y using Maple's numerical integration. What is the " }{TEXT 263 14 "m inimum degree" }{TEXT -1 26 " needed for this to occur?" }}{PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 264 4 "Note" }{TEXT -1 86 ": You may need \+ to use a few guard digits to get agreement in the last decimal place. \+ )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT 290 8 "Solu tion" }{TEXT -1 1 ":" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 " " {TEXT -1 41 "(a) We can construct a Taylor polynomial " }{XPPEDIT 18 0 "p(x);" "6#-%\"pG6#%\"xG" }{TEXT -1 25 " of degree 10 centred at \+ " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 16 " to approximate \+ " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 6 " when " }{TEXT 297 1 "x" }{TEXT -1 12 " is near 0. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 114 "f := x -> exp(x)/(1+sin(x ^2)):\n'f(x)'=f(x);\ntaylor(f(x),x,11);\nconvert(%,polynom):\np := una pply(%,x):\n'p(x)'=p(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#% \"xG*&-%$expGF&\"\"\",&F+F+-%$sinG6#*$)F'\"\"#F+F+!\"\"" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#+;%\"xG\"\"\"\"\"!F%F%#!\"\"\"\"#F)#!\"&\"\"'\" \"$#\"#8\"#C\"\"%#\"$,\"\"$?\"\"\"&#!$p#\"$?(F,#!%,M\"%S]\"\"(#\"%TB\" %k!)\"\")#\"'tW>\"'!)GO\"\"*#!'*3c(\"(+)GO\"#5-%\"OG6#F%\"#6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"pG6#%\"xG,8\"\"\"F)F'F)*&#F)\"\"#F)*$)F 'F,F)F)!\"\"*&#\"\"&\"\"'F)*$)F'\"\"$F)F)F/*&#\"#8\"#CF)*$)F'\"\"%F)F) F)*&#\"$,\"\"$?\"F)*$)F'F2F)F)F)*&#\"$p#\"$?(F)*$)F'F3F)F)F/*&#\"%,M\" %S]F)*$)F'\"\"(F)F)F/*&#\"%TB\"%k!)F)*$)F'\"\")F)F)F)*&#\"'tW>\"'!)GOF )*$)F'\"\"*F)F)F)*&#\"'*3c(\"(+)GOF)*$)F'\"#5F)F)F/" }}}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 67 "(b) The following pictu re shows the graph of the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6 #-%\"pG6#%\"xG" }{TEXT -1 11 " (drawn in " }{TEXT 260 5 "brown" } {TEXT -1 26 ") along with the graph of " }{XPPEDIT 18 0 "f(x)" "6#-%\" fG6#%\"xG" }{TEXT -1 11 " (drawn in " }{TEXT 296 7 "magenta" }{TEXT -1 2 ")." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 83 "plot([f(x),p(x)],x=-1.5..1.5,y=-5..8,color=[magenta,b rown],legend=[`f(x)`,`p(x)`]);" }}{PARA 13 "" 1 "" {GLPLOT2D 320 373 373 {PLOTDATA 2 "6&-%'CURVESG6%7U7$$!3++++++++:!#<$\"3!HX`-i)*[D\"!#=7 $$!3&*****\\P&3YV\"F*$\"3764D*yYYE\"F-7$$!3!***\\iv_8F-7$$!3#****\\P\"\\J\\7F*$\"3j)Q^FiTN V\"F-7$$!3'***\\7V0@&=\"F*$\"3&=M&G#3'**Q:F-7$$!33+]i&exd7\"F*$\"3YJ.5 ![5=F-7$$!3****\\i!3%f+5F*$\" 3!G'[:pe'e*>F-7$$!3;++D\"oS:P*F-$\"3z!QGCP#p8AF-7$$!3h*****\\<#)*=()F- $\"3[T]aYFmvCF-7$$!3#*****\\(G3U9)F-$\"3N$z/5x!>TFF-7$$!3Y*****\\-\\r \\(F-$\"3!3)))3S'=B3$F-7$$!3?+++vGVZoF-$\"3Y;J\"=\\qGZ$F-7$$!3_*****\\ (4J@iF-$\"3+$[l$R[.(*QF-7$$!37++D1Bt_cF-$\"3#3E_A#*HQK%F-7$$!3')***** \\FPm(\\F-$\"3/umPPF-$\"3Hl#3./&GSgF-7$$!3'3+++&=$z9$F-$\"3Rq)ztnhAk'F-7$$!3N*** \\iX/4]#F-$\"3&Ru2!Gg>HtF-7$$!3C***\\(o8y%)=F-$\"3iNu^i07)*zF-7$$!33** **\\i:#>C\"F-$\"3MhH3$*=&zp)F-7$$!3O!***\\7ev:l!#>$\"3_\"[F;g*eH$*F-7$ $!3uF++](o2[\"!#?$\"3%)zD?w8=&)**F-7$$\"3i(***\\P>:mkFir$\"3U\"o\\g0cB 1\"F*7$$\"3d***\\iv&QA7F-$\"3H\"=$>XxQ86F*7$$\"3j++]PPBW=F-$\"3))H>>(H xH;\"F*7$$\"3%*)*****\\Nm'[#F-$\"3v[_i!e%o27F*7$$\"36****\\(yb^6$F-$\" 3v0o4`s([C\"F*7$$\"3')***\\PMaKs$F-$\"3)3nK))pJ\\F\"F*7$$\"3a****\\7TW )R%F-$\"3A!R=ik@@I\"F*7$$\"3*y*****\\@80]F-$\"37^d9JL(=K\"F*7$$\"3_+++ D6!Hl&F-$\"3+ktc1]CR8F*7$$\"3j)**\\P4w)RiF-$\"3I!40mqEGN\"F*7$$\"3s,++ vZf\")oF-$\"3+i*o:UNnO\"F*7$$\"3'z**\\P/-a[(F-$\"3ET[)fzK.Q\"F*7$$\"3R ++v=Yb;\")F-$\"3Q(=Di#*emR\"F*7$$\"3s)****\\i@Ot)F-$\"3t3Ge8&ejT\"F*7$ $\"3g)**\\PfL'z$*F-$\"3hPq%4(*))GW\"F*7$$\"3>+++!*>=+5F*$\"3CuS))>!eiZ \"F*7$$\"3-++DE&4Q1\"F*$\"3$*4iTCG$3_\"F*7$$\"3=+]P%>5p7\"F*$\"3=VcF@r ]y:F*7$$\"39+++bJ*[=\"F*$\"3s33hJ&)fY;F*7$$\"33++Dr\"[8D\"F*$\"3qA;O[k `Zsv_ E?F*7$$\"32+++b!)[/9F*$\"36&fSc[$=@@F*7$$\"31+]i0j\"[V\"F*$\"3]9jVcz]H AF*7$$\"3/+D\"G:3uY\"F*$\"3qj0t?v*QO#F*7$$\"3++++++++:F*$\"3^)z&>d7`?D F*-%'COLOURG6&%$RGBG$\"*++++\"!\")$\"\"!F\\\\lFh[l-%'LEGENDG6#%%f(x)G- F$6%7Y7$F($!3wjIG(>)\\O?!#;7$$!3++]PM@l$[\"F*$!3coj#Gm/:#=Fg\\l7$$!3)* ***\\(oUIn9F*$!3pa$QC\\%*oi\"Fg\\l7$$!3'***\\7.k&4X\"F*$!3cgz[[A%4X\"F g\\l7$F/$!3QllkV-0#H\"Fg\\l7$$!3$**\\il:ghS\"F*$!3G)\\wv#>*=0\"Fg\\l7$ F4$!30\\_7\"p%y>&)F*7$$!3)**\\Pff=dM\"F*$!3cHq0g?&pn'F*7$F9$!3Me()yf=, #>&F*7$$!3))*****\\;?:G\"F*$!3g.zNPT)[*RF*7$F>$!3cCQtQ*pA/$F*7$FC$!3GN =sUY!4q\"F*7$FH$!3E7m^UBj`#*F-7$FM$!3U&yy(4\"f$)G%F-7$FR$!3(o&p#)y<<;7 F-7$FW$\"3Q=nbALYMfFir7$Ffn$\"3sy]/'o'G;QId^3(R&F-7$F^q$\"3LSQth??SgF- 7$Fcq$\"3.(3Pj-\\Ak'F-7$Fhq$\"3MM=j<]>HtF-7$F]r$\"3Au#*f<07)*zF-7$Fbr$ \"3![(*HE*=&zp)F-7$Fgr$\"3#=6F;g*eH$*F-F\\s7$Fcs$\"3G9(\\g0cB1\"F*7$Fh s$\"3H2zAXxQ86F*7$F]t$\"3e2EF+t(H;\"F*7$Fbt$\"3O@zTfYo27F*7$Fgt$\"3JMC 7Y\"y[C\"F*7$F\\u$\"3%z)3?6x$\\F\"F*7$Fau$\"3Nom.No:-8F*7$Ffu$\"3\"*)G 9%H+,A8F*7$F[v$\"3*4x>)3'H(R8F*7$F`v$\"3)p?%=4h;a8F*7$Fev$\"3U7(G2ij.P \"F*7$Fjv$\"3%\\ADv,\"z)Q\"F*7$F_w$\"3-X>)ejzbT\"F*7$Fdw$\"3!*>]!fBB_X \"F*7$Fiw$\"3<@yXdC`?:F*7$F^x$\"3Iyq3al!)>;F*7$Fcx$\"3%=@\"=COtxM#F*7$Fby$\"3Y(zg&zE?&)GF*7$Fgy$\" 3W/ye&=k8e$F*7$F\\z$\"3\"y]8[/eeh%F*7$Ffz$\"3c>'>(*)>\">(fF*7$F`[l$\"3 A;RpKePQzF*-Fe[l6&Fg[l$\")#)eqkFj[l$\"))eqk\"Fj[lFggl-F^\\l6#%%p(x)G-% +AXESLABELSG6$Q\"x6\"Q\"yF`hl-%%VIEWG6$;$!#:!\"\"$\"#:Fhhl;$!\"&F\\\\l $\"\")F\\\\l" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 1 "f(x)" "p(x)" }}}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 4 "(c) " }{XPPEDIT 18 0 "Int(p(x),x=-1/3..1/3)" "6#-%$IntG6$- %\"pG6#%\"xG/F);,$*&\"\"\"F.\"\"$!\"\"F0*&F.F.F/F0" }{TEXT -1 26 " pro vides an estimate for " }{XPPEDIT 18 0 "Int(f(x),x = -1/3 .. 1/3);" "6 #-%$IntG6$-%\"fG6#%\"xG/F);,$*&\"\"\"F.\"\"$!\"\"F0*&F.F.F/F0" }{TEXT -1 2 ". " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "Int(p(x),x=-1/3..1/3);\nvalue(%);\narea_est := evalf( %);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#-%$IntG6$,8\"\"\"F'%\"xGF'*&#F' \"\"#F'*$)F(F+F'F'!\"\"*&#\"\"&\"\"'F'*$)F(\"\"$F'F'F.*&#\"#8\"#CF'*$) F(\"\"%F'F'F'*&#\"$,\"\"$?\"F'*$)F(F1F'F'F'*&#\"$p#\"$?(F'*$)F(F2F'F'F .*&#\"%,M\"%S]F'*$)F(\"\"(F'F'F.*&#\"%TB\"%k!)F'*$)F(\"\")F'F'F'*&#\"' tW>\"'!)GOF'*$)F(\"\"*F'F'F'*&#\"'*3c(\"(+)GOF'*$)F(\"#5F'F'F./F(;#F.F 5#F'F5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6##\".hYz)Q;B\".+[oqb`$" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%)area_estG$\"+W(o;b'!#5" }}}{PARA 0 "" 0 "" {TEXT -1 40 "(d) We can obtain a numerical value for " } {XPPEDIT 18 0 "Int(f(x),x = -1/3 .. 1/3);" "6#-%$IntG6$-%\"fG6#%\"xG/F );,$*&\"\"\"F.\"\"$!\"\"F0*&F.F.F/F0" }{TEXT -1 113 " by Maple's numer ical integration. Using a few guard digits will ensure that the result is correct to 10 digits. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "accurate_val := evalf(Int(f(x),x=-1 /3..1/3));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%-accurate_valG$\"+z)o; b'!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 114 "We can now calculate the absolute and relative errors in the valu e obtained by integrating the Taylor polynomial. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 82 "abserr := ab s(area_est-accurate_val);\nrelerr := evalf(abserr/abs(accurate_val),4) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'abserrG$\"$N\"!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'relerrG$\"%g?!#6" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 31 "(e) Taking a Taylor polynomial \+ " }{XPPEDIT 18 0 "q(x)" "6#-%\"qG6#%\"xG" }{TEXT -1 25 " of degree 14 \+ centred at " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 16 " to ap proximate " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 63 " giv es a value for the integral which is correct to 10 digits. " }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 156 "f := x -> exp(x)/(1+sin(x^2)):\n'f(x)'=f(x);\ntaylor(f(x),x,15):\nconve rt(%,polynom):\nq := unapply(%,x):\n'q(x)'=q(x);\nInt('q(x)',x=-1/3..1 /3);\nevalf(value(%));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#/-%\"fG6#%\" xG*&-%$expGF&\"\"\",&F+F+-%$sinG6#*$)F'\"\"#F+F+!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#/-%\"qG6#%\"xG,@\"\"\"F)F'F)*&#F)\"\"#F)*$)F'F,F)F )!\"\"*&#\"\"&\"\"'F)*$)F'\"\"$F)F)F/*&#\"#8\"#CF)*$)F'\"\"%F)F)F)*&# \"$,\"\"$?\"F)*$)F'F2F)F)F)*&#\"$p#\"$?(F)*$)F'F3F)F)F/*&#\"%,M\"%S]F) *$)F'\"\"(F)F)F/*&#\"%TB\"%k!)F)*$)F'\"\")F)F)F)*&#\"'tW>\"'!)GOF)*$)F '\"\"*F)F)F)*&#\"'*3c(\"(+)GOF)*$)F'\"#5F)F)F/*&#\")H_7;\")+o\"*RF)*$) F'\"#6F)F)F/*&#\")pG(>(\"*+;+z%F)*$)F'\"#7F)F)F)*&#\"*,!*or$\"+gTSX7F) *$)F'F9F)F)F)*&#\"+(>ksD*\",+7Hyr)F)*$)F'\"#9F)F)F/" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#-%$IntG6$-%\"qG6#%\"xG/F);#!\"\"\"\"$#\"\"\"F." }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+z)o;b'!#5" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 6 "Tasks " }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 2 " Q1" }}{PARA 0 "" 0 "" {TEXT -1 37 "(a) Construct a Taylor polynomial o f " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 77 " degree 6, a s a truncated Maclaurin series, which approximates the function " } {XPPEDIT 18 0 "f(x) = sin(sqrt(x^2+1));" "6#/-%\"fG6#%\"xG-%$sinG6#-%% sqrtG6#,&*$F'\"\"#\"\"\"F1F1" }{TEXT -1 7 " when " }{TEXT 268 1 "x" } {TEXT -1 13 " is near 0. " }}{PARA 0 "" 0 "" {TEXT -1 33 "(b) Plot a \+ graph of the function " }{XPPEDIT 18 0 "f(x) = sin(sqrt(x^2+1));" "6#/ -%\"fG6#%\"xG-%$sinG6#-%%sqrtG6#,&*$F'\"\"#\"\"\"F1F1" }{TEXT -1 6 " \+ for " }{TEXT 266 1 "x" }{TEXT -1 9 " between " }{XPPEDIT 18 0 "-3;" "6 #,$\"\"$!\"\"" }{TEXT -1 66 " and 3 along with the graph of the Taylor polynomial found in (a)." }}{PARA 0 "" 0 "" {TEXT -1 66 "(c) Construc t a Taylor polynomial of degree 7, which approximates " }{XPPEDIT 18 0 "Int(sin(sqrt(t^2+1)),t = 0 .. x);" "6#-%$IntG6$-%$sinG6#-%%sqrtG6#, &*$%\"tG\"\"#\"\"\"F0F0/F.;\"\"!%\"xG" }{TEXT -1 7 " when " }{TEXT 267 1 "x" }{TEXT -1 11 " is near 0." }}{PARA 0 "" 0 "" {TEXT -1 13 "(d ) Estimate " }{XPPEDIT 18 0 "Int(sin(sqrt(x^2+1)),x = 0 .. 1/4);" "6#- %$IntG6$-%$sinG6#-%%sqrtG6#,&*$%\"xG\"\"#\"\"\"F0F0/F.;\"\"!*&F0F0\"\" %!\"\"" }{TEXT -1 43 " using the Taylor polynomial found in (b). " }} {PARA 0 "" 0 "" {TEXT -1 49 "(e) Increase the degree of the Taylor pol ynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 19 ", used to estimate " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 71 " \+ until you obtain an estimate for the numerical value of the integral \+ " }{XPPEDIT 18 0 "Int(sin(sqrt(x^2+1)),x = 0 .. 1/4);" "6#-%$IntG6$-%$ sinG6#-%%sqrtG6#,&*$%\"xG\"\"#\"\"\"F0F0/F.;\"\"!*&F0F0\"\"%!\"\"" } {TEXT -1 69 " which is correct to 10 digits. What is the mnimum degree needed for " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 20 " f or this to occur? " }}{PARA 0 "" 0 "" {TEXT -1 34 "___________________ _______________" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" } }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" } }}{PARA 0 "" 0 "" {TEXT -1 34 "__________________________________" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 " " {TEXT -1 2 "Q2" }}{PARA 0 "" 0 "" {TEXT -1 34 "(a) Construct a Taylo r polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 66 " \+ of degree 8, as a truncated Maclaurin series, which approximates " } {XPPEDIT 18 0 "f(x) = sin(x)/x;" "6#/-%\"fG6#%\"xG*&-%$sinG6#F'\"\"\"F '!\"\"" }{TEXT -1 7 " when " }{TEXT 265 1 "x" }{TEXT -1 11 " is near \+ 0." }}{PARA 0 "" 0 "" {TEXT -1 33 "(b) Plot a graph of the function " }{XPPEDIT 18 0 "f(x) = sin(x)/x;" "6#/-%\"fG6#%\"xG*&-%$sinG6#F'\"\"\" F'!\"\"" }{TEXT -1 6 " for " }{TEXT 269 1 "x" }{TEXT -1 9 " between \+ " }{XPPEDIT 18 0 "-6;" "6#,$\"\"'!\"\"" }{TEXT -1 66 " and 6 along wit h the graph of the Taylor polynomial found in (a)." }}{PARA 0 "" 0 "" {TEXT -1 66 "(c) Construct a Taylor polynomial of degree 9, which appr oximates " }{XPPEDIT 18 0 "Int(sin(t)/t,t = 0 .. x);" "6#-%$IntG6$*&-% $sinG6#%\"tG\"\"\"F*!\"\"/F*;\"\"!%\"xG" }{TEXT -1 7 " when " }{TEXT 270 1 "x" }{TEXT -1 11 " is near 0." }}{PARA 0 "" 0 "" {TEXT -1 13 "(d ) Estimate " }{XPPEDIT 18 0 "Int(sin(x)/x,x = 0 .. 1);" "6#-%$IntG6$*& -%$sinG6#%\"xG\"\"\"F*!\"\"/F*;\"\"!F+" }{TEXT -1 130 " using the Tayl or polynomial found in (b), and compare this value with the value give n by Maple's numerical integration procedure " }{TEXT 0 9 "evalf/Int" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 49 "(e) Increase the degree of the Taylor polynomial " } {XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 19 ", used to estimat e " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 71 " until you o btain an estimate for the numerical value of the integral " } {XPPEDIT 18 0 "Int(sin(x)/x,x = 0 .. 1);" "6#-%$IntG6$*&-%$sinG6#%\"xG \"\"\"F*!\"\"/F*;\"\"!F+" }{TEXT -1 69 " which is correct to 10 digits . What is the mnimum degree needed for " }{XPPEDIT 18 0 "p(x)" "6#-%\" pG6#%\"xG" }{TEXT -1 20 " for this to occur? " }}{PARA 0 "" 0 "" {TEXT -1 34 "__________________________________" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 34 "__ ________________________________" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 2 "Q4" }}{PARA 0 "" 0 " " {TEXT -1 34 "(a) Construct a Taylor polynomial " }{XPPEDIT 18 0 "p(x );" "6#-%\"pG6#%\"xG" }{TEXT -1 67 " of degree 10, as a truncated Macl aurin series, which approximates " }{XPPEDIT 18 0 "f(x) = exp(x)/(1+si n(x^2));" "6#/-%\"fG6#%\"xG*&-%$expG6#F'\"\"\",&F,F,-%$sinG6#*$F'\"\"# F,!\"\"" }{TEXT -1 6 " when " }{TEXT 272 1 "x" }{TEXT -1 12 " is near \+ 0. " }}{PARA 0 "" 0 "" {TEXT -1 33 "(b) Plot a graph of the function \+ " }{XPPEDIT 18 0 "f(x) = exp(x)/(1+sin(x^2));" "6#/-%\"fG6#%\"xG*&-%$e xpG6#F'\"\"\",&F,F,-%$sinG6#*$F'\"\"#F,!\"\"" }{TEXT -1 6 " for " } {TEXT 273 1 "x" }{TEXT -1 9 " between " }{XPPEDIT 18 0 "-2;" "6#,$\"\" #!\"\"" }{TEXT -1 5 " and " }{XPPEDIT 18 0 "2;" "6#\"\"#" }{TEXT -1 60 " along with the graph of the Taylor polynomial found in (a)." }} {PARA 0 "" 0 "" {TEXT -1 13 "(c) Estimate " }{XPPEDIT 18 0 "Int(exp(x) /(1+sin(x^2)),x = -1/3 .. 1/3);" "6#-%$IntG6$*&-%$expG6#%\"xG\"\"\",&F +F+-%$sinG6#*$F*\"\"#F+!\"\"/F*;,$*&F+F+\"\"$F2F2*&F+F+F7F2" }{TEXT -1 30 " using the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG 6#%\"xG" }{TEXT -1 15 " found in (a). " }}{PARA 0 "" 0 "" {TEXT -1 13 "(d) Find the " }{TEXT 263 14 "relative error" }{TEXT -1 107 " in the \+ estimate found in (c) by comparing with the value given by Maple's num erical integration procedure " }{TEXT 0 9 "evalf/Int" }{TEXT -1 3 ". \+ " }}{PARA 0 "" 0 "" {TEXT -1 49 "(e) Increase the degree of the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 20 ", \+ used to estimate " }{XPPEDIT 18 0 "Int(exp(x)/(1+sin(x^2)),x = -1/3 . . 1/3);" "6#-%$IntG6$*&-%$expG6#%\"xG\"\"\",&F+F+-%$sinG6#*$F*\"\"#F+! \"\"/F*;,$*&F+F+\"\"$F2F2*&F+F+F7F2" }{TEXT -1 125 " until the value \+ of the integral agrees to 10 digits with that obtained by using Maple' s numerical integration. What is the " }{TEXT 263 14 "minimum degree" }{TEXT -1 26 " needed for this to occur?" }}{PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 264 4 "Note" }{TEXT -1 86 ": You may need to use a few g uard digits to get agreement in the last decimal place. )" }}{PARA 0 " " 0 "" {TEXT -1 34 "__________________________________" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 87 "f := x \+ -> exp(x)/(1+sin(x^2));\ntaylor(f(x),x=0,11):\np := unapply(convert(%, polynom),x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fGf*6#%\"xG6\"6$%) operatorG%&arrowGF(*&-%$expG6#9$\"\"\",&F1F1-%$sinG6#*$)F0\"\"#F1F1!\" \"F(F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"pGf*6#%\"xG6\"6$%)oper atorG%&arrowGF(,8\"\"\"F-9$F-*&#F-\"\"#F-*$)F.F1F-F-!\"\"*&#\"\"&\"\"' F-*$)F.\"\"$F-F-F4*&#\"#8\"#CF-*$)F.\"\"%F-F-F-*&#\"$,\"\"$?\"F-*$)F.F 7F-F-F-*&#\"$p#\"$?(F-*$)F.F8F-F-F4*&#\"%,M\"%S]F-*$)F.\"\"(F-F-F4*&# \"%TB\"%k!)F-*$)F.\"\")F-F-F-*&#\"'tW>\"'!)GOF-*$)F.\"\"*F-F-F-*&#\"'* 3c(\"(+)GOF-*$)F.\"#5F-F-F4F(F(F(" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "plot([f(x),p(x)],x=-2..2,y=- 1..4);" }}{PARA 13 "" 1 "" {GLPLOT2D 365 365 365 {PLOTDATA 2 "6&-%'CUR VESG6$7jn7$$!\"#\"\"!$\"3[*)fHc-$[c&!#=7$$!3MLLL$Q6G\">!#<$\"3(pp'f0m8 @HF-7$$!3bmm;M!\\p$=F1$\"3aa:>M\"G22#F-7$$!3MLLL))Qj^WJ8F-7$$!35++]d(Q&\\7F1$\"3UvNfp o@L9F-7$$!3gmmmc4`i6F1$\"3ce[ac9U#e\"F-7$$!3KLLLQW*e3\"F1$\"33'oWxE8Wv \"F-7$$!3w++++()>'***F-$\"3Q!3i&>L&*)*>F-7$$!3E++++0\"*H\"*F-$\"3hcBPp (=gI#F-7$$!35++++83&H)F-$\"3aJ%*49a=oEF-7$$!3\\LLL3k(p`(F-$\"3S\\jvME( *fIF-7$$!3Anmmmj^NmF-$\"3TuzpEw46OF-7$$!3)zmmmYh=(eF-$\"3I^ilX%QY:%F-7 $$!3+,++v#\\N)\\F-$\"3m&eI$p8dw[F-7$$!3commmCC(>%F-$\"3Q;Mz)=*>#f&F-7$ $!39*****\\FRXL$F-$\"3S>8$\"3w,$pS^C#*4*F- 7$$!3IqLLL$eV(>!#?$\"3s7)H/*pB!)**F-7$$\"3)Qjmm\"f`@')Fjr$\"3[W>\"p`)* >3\"F17$$\"3%z****\\nZ)H;F-$\"3aN&f;&\\cY6F17$$\"3ckmm;$y*eCF-$\"3ya66 Jy*e?\"F17$$\"3f)******R^bJ$F-$\"3U8!y#4bSb7F17$$\"3'e*****\\5a`TF-$\" 3+*)*pe#H&HH\"F17$$\"3'o****\\7RV'\\F-$\"3)=vt$)F-$\"3^%Hc*obl-9F17$$\"3M *******pfa<*F-$\"3a<92F*3D5h`\"F17$$\"3;LLL$)G[k6F1$\"3Zc\"3]z^2i \"F17$$\"3#)****\\7yh]7F1$\"35,m(***)piu\"F17$$\"3xmmm')fdL8F1$\"3)))* z*)3+*y\">F17$$\"3bmmm,FT=9F1$\"3)f&=rn$)3p@F17$$\"3FLL$e#pa-:F1$\"32& R$R_q$Q`#F17$$\"3!*******Rv&)z:F1$\"3\"35VgNs1.$F17$$\"3ILLLGUYo;F1$\" 3q$oP?!GAGRF17$$\"3_mmm1^rZYzF17$$\"35++D1>V_=F1$\"3h! o\\kVQt#*)F17$$\"33+++S2ls=F1$\"3`i@DE-+75!#;7$$\"33++vt&pG*=F1$\"3g)R .!*)*z)e6F`[l7$$\"34++]2%)38>F1$\"3[&R/!*p*[U8F`[l7$$\"33+DccB&R#>F1$ \"3%[hSrj,2Y\"F`[l7$$\"31+]i0j\"[$>F1$\"33=3-fY%ff\"F`[l7$$\"31+voa-oX >F1$\"3S*yJby@;v\"F`[l7$$\"3/++v.Uac>F1$\"3$=(G]Cl-K>F`[l7$$\"3/+D\"G: 3u'>F1$\"3$*Q!G+1PE9#F`[l7$$\"3-+](=5s#y>F1$\"3O*)*en?I0R#F`[l7$$\"39] iSwSq$)>F1$\"3f$>,>CZ7`#F`[l7$$\"3-+v$40O\"*)>F1$\"3kh!=&4U(\\o#F`[l7$ $\"3!*\\(oa-oX*>F1$\"3%*HkNZNN`GF`[l7$$\"\"#F*$\"3U*GACP%HQIF`[l-%'COL OURG6&%$RGBG$\"#5!\"\"$F*F*Fd_l-F$6$7en7$F($!3!zuG63icl$!#:7$$!3SLL$e% G?y>F1$!3)eQl+J7$yKF[`l7$$!3ymmm\"p0k&>F1$!3*flU`9sh$HF[`l7$$!3&***** \\P&3Y$>F1$!35sJkZ8JEEF[`l7$F/$!3gg#p)QM,YBF[`l7$$!3em;/'zXQ*=F1$!3onh 'RG;U7#F[`l7$$!31++v3-)[(=F1$!3\\_>Dd6N@>F[`l7$$!3IL$e9i9f&=F1$!3oXXd^ $yft\"F[`l7$F5$!3W=(y=+anc\"F[`l7$$!3#)***\\7Y\"H%z\"F1$!3g6#\\dT!zQ7F [`l7$F:$!34,`n]e!et*F`[l7$$!3O\"F17$F\\o$!3= 2wyjTogdF-7$Fao$!3_Bl!)4(Q1=\"F-7$Ffo$\"3e'RuB!p(Q2\"F-7$F[p$\"3(Q&Gzt jA>AF-7$F`p$\"3#=!y4HVO(*GF-7$Fep$\"3-I$z9yE$pNF-7$Fjp$\"3[Vf6ihSVTF-7 $F_q$\"3=BZ\\P\\mu[F-7$Fdq$\"3O$HYcu/>f&F-7$Fiq$\"3M?\\\"flU)[kF-7$F^r $\"33\"F17$Fis$\"3wT=Y_\\cY6F17$F^t$\"3\"\\*4#4!z*e?\"F17$Fc t$\"3I&zG;D2aD\"F17$Fht$\"3g%oAB8sHH\"F17$F]u$\"3KT#)4q#*y?8F17$Fbu$\" 3$Qh\"Hs_-X8F17$Fgu$\"3y61=U$)*[O\"F17$F\\v$\"3GB8v9:j!R\"F17$Fav$\"3r $Q&)p*[!oU\"F17$Ffv$\"3iL$*e#H0l\\\"F17$F[w$\"3yJB%)Rdj:;F17$F`w$\"3y7 s$*ybnP=F17$Few$\"3bO-j\"4:U@#F17$Fjw$\"3/$eZ$QI,yGF17$F_x$\"3dFIbtLQ9 RF17$Fdx$\"3#=9*eD()fkbF17$Fix$\"3\"*oH_q^mF!)F17$F^y$\"3>r2!y*[1E6F`[ l7$Fcy$\"3R$=()\\5ZOk\"F`[l7$Fhy$\"3sp*Q*3?\\pAF`[l7$Fbz$\"3Jd8)>$GJMJ F`[l7$Fg[l$\"3i8R?78roTF`[l7$Fi^l$\"3#4RM:RM&*\\&F`[l-F^_l6&F`_lFd_lFa _lFd_l-%+AXESLABELSG6$Q\"x6\"Q\"yFg[m-%%VIEWG6$;F(Fi^l;$Fc_lF*$\"\"%F* " 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1" "C urve 2" }}}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" } }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 " " {TEXT -1 34 "__________________________________" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 2 "Q5 " }}{PARA 0 "" 0 "" {TEXT -1 13 "(a) Estimate " }{XPPEDIT 18 0 "Int(1/ (x+exp(x)),x = 0 .. 1/2);" "6#-%$IntG6$*&\"\"\"F',&%\"xGF'-%$expG6#F)F '!\"\"/F);\"\"!*&F'F'\"\"#F-" }{TEXT -1 54 " by using a Taylor polynom ial of degree 12 centred at " }{XPPEDIT 18 0 "x = 1/4;" "6#/%\"xG*&\" \"\"F&\"\"%!\"\"" }{TEXT -1 16 " to approximate " }{XPPEDIT 18 0 "f(x) = 1/(x+exp(x));" "6#/-%\"fG6#%\"xG*&\"\"\"F),&F'F)-%$expG6#F'F)!\"\" " }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 87 "(b) Compare this valu e with the value given by Maple's numerical integration procedure " } {TEXT 0 9 "evalf/Int" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 13 "( c) Estimate " }{XPPEDIT 18 0 "Int(1/(x+exp(x)),x = 0 .. 1/2);" "6#-%$I ntG6$*&\"\"\"F',&%\"xGF'-%$expG6#F)F'!\"\"/F);\"\"!*&F'F'\"\"#F-" } {TEXT -1 54 " by using a Taylor polynomial of degree 12 centred at " } {XPPEDIT 18 0 "x = 0;" "6#/%\"xG\"\"!" }{TEXT -1 16 " to approximate \+ " }{XPPEDIT 18 0 "f(x) = 1/(x+exp(x));" "6#/-%\"fG6#%\"xG*&\"\"\"F),&F 'F)-%$expG6#F'F)!\"\"" }{TEXT -1 1 "." }}{PARA 0 "" 0 "" {TEXT -1 0 " " }}{PARA 0 "" 0 "" {TEXT 264 4 "Note" }{TEXT -1 2 ": " }}{PARA 15 "" 0 "" {TEXT -1 114 "In part (a) force Maple to use floating point arith metic by entering the centre for the Taylor series in the form " } {XPPEDIT 18 0 "x = 0" "6#/%\"xG\"\"!" }{TEXT -1 4 ".25." }}{PARA 15 " " 0 "" {TEXT -1 66 "It is not be sensible to use a Taylor polynomial a pproximaton for " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 12 " centred at " }{XPPEDIT 18 0 "x=0" "6#/%\"xG\"\"!" }{TEXT -1 63 " \+ because the radius of convergence of the Maclaurin series for " } {XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 25 " is not much larg er than " }{XPPEDIT 18 0 "1/2" "6#*&\"\"\"F$\"\"#!\"\"" }{TEXT -1 2 ". " }}{PARA 0 "" 0 "" {TEXT -1 41 "____________________________________ _____" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" } }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" } }}{PARA 0 "" 0 "" {TEXT -1 41 "_______________________________________ __" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 2 "Q6" }{TEXT 280 1 " " }{TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 34 "(a) Construct a Taylor polynomial " }{XPPEDIT 18 0 " p(x);" "6#-%\"pG6#%\"xG" }{TEXT -1 66 " of degree 10 as a truncated Ma claurin series, which approximates " }{XPPEDIT 18 0 "f(x) = cosh(x)/(x ^2-x+1);" "6#/-%\"fG6#%\"xG*&-%%coshG6#F'\"\"\",(*$F'\"\"#F,F'!\"\"F,F ,F0" }{TEXT -1 6 " when " }{TEXT 278 1 "x" }{TEXT -1 11 " is near 0." }}{PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 271 4 "Note" }{TEXT -1 60 ": \+ You will need an order 11 Taylor series approximation for " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" }{TEXT -1 3 ". )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 33 "(b) Plot a graph of the f unction " }{XPPEDIT 18 0 "f(x) = cosh(x)/(x^2-x+1);" "6#/-%\"fG6#%\"xG *&-%%coshG6#F'\"\"\",(*$F'\"\"#F,F'!\"\"F,F,F0" }{TEXT -1 6 " for " } {TEXT 279 1 "x" }{TEXT -1 9 " between " }{XPPEDIT 18 0 "-2;" "6#,$\"\" #!\"\"" }{TEXT -1 5 " and " }{XPPEDIT 18 0 "2;" "6#\"\"#" }{TEXT -1 60 " along with the graph of the Taylor polynomial found in (a)." }} {PARA 0 "" 0 "" {TEXT -1 13 "(c) Estimate " }{XPPEDIT 18 0 "Int(cosh(x )/(x^2-x+1),x = -1/3 .. 1/3);" "6#-%$IntG6$*&-%%coshG6#%\"xG\"\"\",(*$ F*\"\"#F+F*!\"\"F+F+F//F*;,$*&F+F+\"\"$F/F/*&F+F+F4F/" }{TEXT -1 30 " \+ using the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG " }{TEXT -1 15 " found in (a). " }}{PARA 0 "" 0 "" {TEXT -1 13 "(d) Fi nd the " }{TEXT 263 14 "relative error" }{TEXT -1 107 " in the estimat e found in (c) by comparing with the value given by Maple's numerical \+ integration procedure " }{TEXT 0 9 "evalf/int" }{TEXT -1 3 ". " }} {PARA 0 "" 0 "" {TEXT -1 49 "(e) Increase the degree of the Taylor pol ynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 20 ", used to estimate " }{XPPEDIT 18 0 "Int(cosh(x)/(x^2-x+1),x = -1/3 .. 1/3) ;" "6#-%$IntG6$*&-%%coshG6#%\"xG\"\"\",(*$F*\"\"#F+F*!\"\"F+F+F//F*;,$ *&F+F+\"\"$F/F/*&F+F+F4F/" }{TEXT -1 125 " until the value of the int egral agrees to10 digits with that obtained by using Maple's numerical integration. What is the " }{TEXT 263 14 "minimum degree" }{TEXT -1 26 " needed for this to occur?" }}{PARA 0 "" 0 "" {TEXT -1 2 "( " } {TEXT 264 4 "Note" }{TEXT -1 86 ": You may need to use a few guard dig its to get agreement in the last decimal place. )" }}{PARA 0 "" 0 "" {TEXT -1 29 "_____________________________" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" } }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" } }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 29 "_____________________________" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{SECT 1 {PARA 4 "" 0 "" {TEXT -1 2 "Q7" }} {PARA 0 "" 0 "" {TEXT -1 34 "(a) Construct a Taylor polynomial " } {XPPEDIT 18 0 "p(x);" "6#-%\"pG6#%\"xG" }{TEXT -1 65 " of degree 8 as \+ a truncated Maclaurin series, which approximates " }{XPPEDIT 18 0 "f(x ) = (1+x*sinh(x))/(x^2-x+1);" "6#/-%\"fG6#%\"xG*&,&\"\"\"F**&F'F*-%%si nhG6#F'F*F*F*,(*$F'\"\"#F*F'!\"\"F*F*F2" }{TEXT -1 6 " when " }{TEXT 281 1 "x" }{TEXT -1 11 " is near 0." }}{PARA 0 "" 0 "" {TEXT -1 2 "( \+ " }{TEXT 271 4 "Note" }{TEXT -1 59 ": You will need an order 9 Taylor \+ series approximation for " }{XPPEDIT 18 0 "f(x)" "6#-%\"fG6#%\"xG" } {TEXT -1 3 ". )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 33 "(b) Plot a graph of the function " }{XPPEDIT 18 0 "f(x) = (1+x*sinh(x))/(x^2-x+1);" "6#/-%\"fG6#%\"xG*&,&\"\"\"F**&F'F*-%%sinhG 6#F'F*F*F*,(*$F'\"\"#F*F'!\"\"F*F*F2" }{TEXT -1 5 " for " }{TEXT 282 1 "x" }{TEXT -1 9 " between " }{XPPEDIT 18 0 "-2;" "6#,$\"\"#!\"\"" } {TEXT -1 5 " and " }{XPPEDIT 18 0 "2;" "6#\"\"#" }{TEXT -1 60 " along \+ with the graph of the Taylor polynomial found in (a)." }}{PARA 0 "" 0 "" {TEXT -1 13 "(c) Estimate " }{XPPEDIT 18 0 "Int((1+x*sinh(x))/(x^2- x+1),x = -1/4 .. 1/4);" "6#-%$IntG6$*&,&\"\"\"F(*&%\"xGF(-%%sinhG6#F*F (F(F(,(*$F*\"\"#F(F*!\"\"F(F(F1/F*;,$*&F(F(\"\"%F1F1*&F(F(F6F1" } {TEXT -1 29 " using the Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#- %\"pG6#%\"xG" }{TEXT -1 15 " found in (a). " }}{PARA 0 "" 0 "" {TEXT -1 13 "(d) Find the " }{TEXT 263 14 "relative error" }{TEXT -1 107 " i n the estimate found in (c) by comparing with the value given by Maple 's numerical integration procedure " }{TEXT 0 9 "evalf/int" }{TEXT -1 3 ". " }}{PARA 0 "" 0 "" {TEXT -1 49 "(e) Increase the degree of the \+ Taylor polynomial " }{XPPEDIT 18 0 "p(x)" "6#-%\"pG6#%\"xG" }{TEXT -1 20 ", used to estimate " }{XPPEDIT 18 0 "Int((1+x*sinh(x))/(x^2-x+1), x = -1/4 .. 1/4);" "6#-%$IntG6$*&,&\"\"\"F(*&%\"xGF(-%%sinhG6#F*F(F(F( ,(*$F*\"\"#F(F*!\"\"F(F(F1/F*;,$*&F(F(\"\"%F1F1*&F(F(F6F1" }{TEXT -1 124 " until the value of the integral agrees to10 digits with that obt ained by using Maple's numerical integration. What is the " }{TEXT 263 14 "minimum degree" }{TEXT -1 26 " needed for this to occur?" }} {PARA 0 "" 0 "" {TEXT -1 2 "( " }{TEXT 264 4 "Note" }{TEXT -1 86 ": Yo u may need to use a few guard digits to get agreement in the last deci mal place. )" }}{PARA 0 "" 0 "" {TEXT -1 29 "_________________________ ____" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" } }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 " " {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" } }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 0 "" 0 "" {TEXT -1 29 "_____________________________" }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 1 ";" }}}}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 " " 0 "" {TEXT -1 0 "" }}}{MARK "4 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }