{VERSION 5 0 "IBM INTEL NT" "5.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 1 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "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 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 } } {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q1:" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 8 "phi0:=1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%% phi0G\"\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "phi1:=x-c;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%phi1G,&%\"xG\"\"\"%\"cG!\"\"" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "int(sin(x)*phi0*phi1,x=0..Pi /2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&\"\"\"F$%\"cG!\"\"" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "solve(%,c);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "phi1:=x-1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%phi1G,&%\"xG\"\" \"F'!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "phi2:=x^2-c1*x +c0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%phi2G,(*$)%\"xG\"\"#\"\"\"F **&%#c1GF*F(F*!\"\"%#c0GF*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "int(sin(x)*phi2*phi0,x=0..Pi/2);" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#,*%#c1G!\"\"%#PiG\"\"\"\"\"#F%%#c0GF'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "int(sin(x)*phi2*phi1,x=0..Pi/2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,,%#PiG!\"\"*&#\"\"$\"\"%\"\"\")F$\"\"#F*F**&%#c1GF*F$F *F%F)F%*&F(F*F.F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "solv e(\{%,%%\},\{c0,c1\});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<$/%#c0G,$*& ,(%#PiG!#;*$)F)\"\"#\"\"\"F.\"#SF.F.,&F)F.\"\"$!\"\"F2#F2\"\"%/%#c1G,$ *&,(F)!\"%*&F1F.F,F.F.\"#;F2F.F0F2#F.F4" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "assign(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 5 "phi2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,(*$)%\"xG\"\"#\"\"\"F(*&# F(\"\"%F(*(,(%#PiG!\"%*&\"\"$F()F.F'F(F(\"#;!\"\"F(,&F.F(F1F4F4F&F(F(F 4*&#F(F+F(*&,(F.!#;*$F2F(F(\"#SF(F(F5F4F(F4" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "evalf(phi2);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,( *$)%\"xG\"\"#\"\"\"F(*&$\"+ghcS=!\"*F(F&F(!\"\"$\"+X^t*)p!#5F(" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q2:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 48 "Note: x = (X + 3) / 2 and X = 2 * x - 3" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "a:=int(ln((X+3)/2),X=-1..1)/2;" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%\"aG,&-%#lnG6#\"\"#F)\"\"\"!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "b:=int(ln((X+3)/2)*X,X=-1..1)/(2/3);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bG,&-%#lnG6#\"\"#!\"'#\"\"*F)\"\" \"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 45 "c:=int(ln((X+3)/2)*(X ^2-1/3),X=-1..1)/(8/45);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cG,&-% #lnG6#\"\"#\"#X#\"$D\"\"\"%!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "p:=a+b*X+c*(X^2-1/3);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"pG,*-%#lnG6#\"\"#F)\"\"\"!\"\"*&,&F&!\"'#\"\"*F)F*F*%\"XGF*F **&,&F&\"#X#\"$D\"\"\"%F+F*,&*$)F1F)F*F*#F*\"\"$F+F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "p:=subs(X=2*x-3,p);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"pG,*-%#lnG6#\"\"#F)\"\"\"!\"\"*&,&F&!\"'#\"\"* F)F*F*,&%\"xGF)\"\"$F+F*F**&,&F&\"#X#\"$D\"\"\"%F+F*,&*$)F1F)F*F*#F*F3 F+F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "p:=evalf(expand(p ));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"pG,($!)$*)H9\"!\"(\"\"\"*&$ \")jv#Q\"F(F)%\"xGF)F)*&$\"(v]L#F(F))F-\"\"#F)!\"\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 14 "Just to check:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "plot(\{ln(x),p\},x=1..2);" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6&-%'CURVESG6$7S7$$\"\"\"\"\"!$\"3q c++++]fi!#?7$$\"3hmm;arz@5!#<$\"3M=9*f664h#!#>7$$\"3OL$e9ui2/\"F1$\"3k (\\\"y`H'*>VF47$$\"3smm\"z_\"4i5F1$\"3-H)zA:->A'F47$$\"3qmmT&phN3\"F1$ \"33]$=wh\"*\\6)F47$$\"3UL$e*=)H\\5\"F1$\"3![$p>#*)3x(**F47$$\"3sm;z/3 uC6F1$\"3Qs?utOco6!#=7$$\"3-+]7LRDX6F1$\"3+()fr'>yMM\"FN7$$\"3em;zR'ok ;\"F1$\"3KG`(f$zIA:FN7$$\"3-+]i5`h(=\"F1$\"3_<,*3Ss%)p\"FN7$$\"3YLL$3E n$47F1$\"33VO*4G+v(=FN7$$\"3cmmT!RE&G7F1$\"3#=Vj^&yNL?FN7$$\"3)*****\\ K]4]7F1$\"3?!)\\DR`w1AFN7$$\"3))****\\PAvr7F1$\"3y%ym#o))pyBFN7$$\"3/+ +]nHi#H\"F1$\"3iL'z\\4;Ba#FN7$$\"3bm;z*ev:J\"F1$\"31:,b(zL\"*o#FN7$$\" 3ELL$347TL\"F1$\"3?p)f`!4`hGFN7$$\"3=LLLjM?`8F1$\"3$p!)>t5=d+$FN7$$\"3 #***\\7o7Tv8F1$\"3!Qg.QX+8<$FN7$$\"3ALLLQ*o]R\"F1$\"3gcijsr%fJ$FN7$$\" 3-+]7=lj;9F1$\"3#[,ueOrDZ$FN7$$\"3&***\\PaROFN7$$\"3 GLLe9Ege9F1$\"3]6\"Hb\\.UFN7$$\"3%)** \\(=>Y2a\"F1$\"3cR&*))H-kJVFN7$$\"3imm\"zXu9c\"F1$\"3m&z\"Qxx5oWFN7$$ \"3'******\\y))Ge\"F1$\"3&\\3`l![)pg%FN7$$\"3!****\\i_QQg\"F1$\"3(z/0& okxSZFN7$$\"3#***\\7y%3Ti\"F1$\"39LgAwdFo[FN7$$\"3#****\\P![hY;F1$\"3Y k*3n&Rf2]FN7$$\"3ELLLQx$om\"F1$\"3C%fEUxf28&FN7$$\"3')****\\P+V)o\"F1$ \"3=Tkpks:g_FN7$$\"3im;zpe*zq\"F1$\"3#\\$4`.-`v`FN7$$\"3)*****\\#\\'QH *GU9cFN7$$\"3%***\\i? =bq:>!*[yQ'FN7$$\"3@LL3dg6<>F1$\"3Q3-'Gdpp\\'FN7$$\"3_mmmw(G p$>F1$\"3S_8O`pi#f'FN7$$\"3-+]7oK0e>F1$\"3Efi]^$)f#p'FN7$$\"3-+](=5s#y >F1$\"3+rG89=L'y'FN7$$\"\"#F*$\"3'3++++L\\)oFN-%'COLOURG6&%$RGBG$\"#5! \"\"$F*F*Fb[l-F$6$7S7$F(Fb[l7$F/$\"3o+;Hz#*Hc@F47$F6$\"3Cvs9f\\Q&*RF47 $F;$\"3V!)[]e.,CgF47$F@$\"3Nzg%G7[`-)F47$FE$\"3ifDyp?=y**F47$FJ$\"3UnY qLh_v6FN7$FP$\"3/'RE2)QEc8FN7$FU$\"3Z-#))4G4)R:FN7$FZ$\"3(>7m*eNZ>FN7$F^o$\"3aT,pXR:e?FN7$Fco$\"3p)zVUu&>KAFN7$Fho$\"3 S9'e5kcRS#FN7$F]p$\"3aY)*R@YtmDFN7$Fbp$\"3s7O*za\"H7FFN7$Fgp$\"3i*z`2q fE)GFN7$F\\q$\"32ef'G\"RB%*f80%FN7$Fds$\"3'G**Rsv?t>%FN7$Fis$\"3 &)HU_#RoEK%FN7$F^t$\"3c<90-aIcWFN7$Fct$\"3_A[(eA:Df%FN7$Fht$\"3Fqhc\\$ )*Rs%FN7$F]u$\"3qUn$QO!f\\[FN7$Fbu$\"3%*3MPja@()\\FN7$Fgu$\"36%e&)\\h# G4^FN7$F\\v$\"33Voz_7*zB&FN7$Fav$\"3wCp3sn?``FN7$Ffv$\"3QTP,qrmxaFN7$F [w$\"31c2hHpP$f&FN7$F`w$\"3fNqi5E\"Hr&FN7$Few$\"30')Rd\\_TGeFN7$Fjw$\" 3g#)HL'>Bz%fFN7$F_x$\"3cEz8@*)ohgFN7$Fdx$\"3I:km:vqwhFN7$Fix$\"3)=9e)= >[*G'FN7$F^y$\"3[e_iT/+#R'FN7$Fcy$\"3-i%QA.?#3lFN7$Fhy$\"3N)4(HQh.6mFN 7$F]z$\"3l@k8\"\\2&>nFN7$Fbz$\"3V^^#)yyBAoFN7$Fgz$\"3'GX*f0=ZJpFN-F\\[ l6&F^[lFb[lF_[lFb[l-%+AXESLABELSG6$Q\"x6\"Q!F]el-%%VIEWG6$;F(Fgz%(DEFA ULTG" 1 2 0 1 10 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Curve 1 " "Curve 2" }}}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q3:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 137 "solve(\{c0+c1+c2=int(1/x^(1/3),x=0..2),\n c1 +2*c2=int(x/x^(1/3),x=0..2),\n c1+4*c2=int(x^2/x^(1/3),x=0..2)\} ,\n \{c0,c1,c2\});" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#<%/%#c0G,$ *$)\"\"##F)\"\"$\"\"\"#\"\"*\"#?/%#c1G,$F'#F.\"#5/%#c2G,$F'#F+F/" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 10 "assign(%);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 27 "The approximate formula is:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "c0*y(0)+c1*y(1)+c2*y(2);" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#,(*&)\"\"##F&\"\"$\"\"\"-%\"yG6#\"\"!F)#\"\"*\"#?*(#F /\"#5F)F%F)-F+6#F)F)F)*(#F(F0F)F%F)-F+6#F&F)F)" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "or:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 9 "eval f(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,(-%\"yG6#\"\"!$\"+MZIVr!#5*& $\"+Z4mG9!\"*\"\"\"-F%6#F/F/F/*&$\"+y:5\"Q#F*F/-F%6#\"\"#F/F/" }}} {EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q4:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "y:=x->cos(sqrt(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"yGf*6#%\"x G6\"6$%)operatorG%&arrowGF(-%$cosG6#-%%sqrtG6#9$F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "I2:=evalf(y(0)+4*y(1/2)+y(1))/6;" } }{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#I2G$\"+DyYNw!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "I4:=evalf((y(0)+4*y(1/4)+2*y(1/2)+4*y(3/4 )+y(1))/12);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#I4G$\"+RfYNw!#5" }} }{EXCHG {PARA 0 "" 0 "" {TEXT -1 8 "Romberg:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "(16*I4-I2)/15;" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#$\"+9eYNw!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 9 "To check:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "int(y(x),x=0..1.);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+9eYNw!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 141 "Our 'approximate' answer is amazingly accurate! This is \+ because cos ( sqrt ( x ) ) is, in the 0 to 1 interval, very 'smo oth' - plot it!" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q5:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 11 "A:=0; B:=1;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\" AG\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"BG\"\"\"" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 17 "The nodes are at:" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 76 "x1:=(A+B)/2-(B-A)/2*sqrt(3/5); x2:=(A+B)/2; \+ x3:=(A+B)/2+(B-A)/2*sqrt(3/5);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#x 1G,&#\"\"\"\"\"#F'*&#F'\"#5F'*$-%%sqrtG6#\"#:F'F'!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#x2G#\"\"\"\"\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#x3G,&#\"\"\"\"\"#F'*&#F'\"#5F'-%%sqrtG6#\"#:F'F'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "y:=x->cos(sqrt(x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"yGf*6#%\"xG6\"6$%)operatorG%&arrowGF(-%$cosG6# -%%sqrtG6#9$F(F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 36 "evalf ((5*y(x1)+8*y(x2)+5*y(x3))/18);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\" +8eYNw!#5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 27 "Again, practically ' exact'." }}}}{MARK "49 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 } {PAGENUMBERS 0 1 2 33 1 1 }