{VERSION 6 1 "Windows XP" "6.1" } {USTYLETAB {PSTYLE "Warning" -1 7 1 {CSTYLE "" -1 -1 "Courier" 1 12 0 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Dash Item" -1 16 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 1 0 -1 3 3 1 0 1 0 2 2 -1 3 }{PSTYLE "Heading 4" -1 20 1 {CSTYLE "" -1 -1 "MS Serif" 1 12 0 0 0 0 1 0 0 2 2 2 0 0 0 1 }1 1 0 -1 0 0 1 0 1 0 2 2 -1 1 }{PSTYLE "Heading 3" -1 5 1 {CSTYLE "" -1 -1 " MS Serif" 1 14 0 0 0 0 1 1 0 2 2 2 0 0 0 1 }1 1 0 -1 0 0 1 0 1 0 2 2 -1 1 }{PSTYLE "Error" -1 8 1 {CSTYLE "" -1 -1 "Courier" 1 12 255 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "A uthor" -1 19 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }3 1 0 -1 8 8 1 0 1 0 2 2 -1 1 }{PSTYLE "Heading 2" -1 4 1 {CSTYLE "" -1 -1 "MS Serif" 1 16 0 0 0 0 0 1 0 2 2 2 0 0 0 1 }1 1 0 -1 8 2 1 0 1 0 2 2 -1 1 }{PSTYLE "Text Output" -1 2 1 {CSTYLE "" -1 -1 "Courier" 1 12 0 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Heading 1" -1 3 1 {CSTYLE "" -1 -1 "MS Serif" 1 18 0 0 0 0 0 1 0 2 2 2 0 0 0 1 }1 1 0 -1 8 4 1 0 1 0 2 2 -1 1 } {PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Maple Plot" -1 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Line Printed Output" -1 6 1 {CSTYLE "" -1 -1 "Courier" 1 12 0 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Title" -1 18 1 {CSTYLE "" -1 -1 "Times" 1 18 0 0 0 0 0 1 1 2 2 2 0 0 0 1 }3 1 0 -1 12 12 1 0 1 0 2 2 -1 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 3 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "List Item" -1 14 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 1 0 -1 3 3 1 0 1 0 2 2 -1 5 }{PSTYLE "Bullet Item" -1 15 1 {CSTYLE "" -1 -1 "Ti mes" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 1 0 -1 3 3 1 0 1 0 2 2 -1 2 } {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 2 1 2 0 0 0 1 }{CSTYLE "2D Input" -1 19 "Times" 1 12 255 0 0 1 0 0 0 2 1 2 0 0 0 1 }{CSTYLE "Hyperlink" -1 17 "MS Serif" 1 12 0 128 128 1 0 0 1 2 2 2 0 0 0 1 }{CSTYLE "Text" -1 200 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }{CSTYLE "2D Math" -1 2 "Times" 1 12 0 0 0 1 0 0 0 2 2 2 0 0 0 1 } {CSTYLE "Dictionary Hyperlink" -1 45 "MS Serif" 1 12 147 0 15 1 0 0 1 2 2 2 0 0 0 1 }{CSTYLE "Maple Input Placeholder" -1 201 "Courier" 1 12 200 0 200 1 0 1 0 2 1 2 0 0 0 1 }{CSTYLE "2D Output" -1 20 "Times" 0 1 0 0 255 1 0 0 2 2 2 2 0 0 0 1 }{CSTYLE "Page Number" -1 33 "Times" 1 10 0 0 0 0 0 0 2 2 2 2 0 0 0 1 }{PSTYLE "_pstyle1" -1 200 1 {CSTYLE "" -1 -1 "Courier" 0 1 255 0 0 1 0 1 0 2 1 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "_pstyle2" -1 201 1 {CSTYLE "" -1 -1 "Times" 0 1 0 0 255 1 0 0 2 2 2 2 0 0 0 1 }3 3 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "_pstyle3" -1 202 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }3 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{CSTYLE "_cstyle 1" -1 202 "Times" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{PSTYLE "_pstyle4" -1 203 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{CSTYLE "_cstyle2" -1 203 "Times" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }} {SECT 0 {EXCHG {PARA 200 "> " 0 "" {MPLTEXT 1 0 20 "#Andrew Sommese 20 04" }}}{EXCHG {PARA 200 "> " 0 "" {MPLTEXT 1 0 8 "restart;" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 200 "> " 0 "" {MPLTEXT 1 0 32 "g := x -> 4/x- 2/x^2+ln(x)-3/2; " }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 16 "p:= x -> -4/x; " }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 19 "q := x -> -2/x^2; " }{MPLTEXT 1 0 0 "" }}{PARA 200 " > " 0 "" {MPLTEXT 1 0 24 "r := x -> 2*ln(x)/x^2; " }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 56 "expand(diff(g(x),x$2)-p(x)*dif f(g(x),x)-q(x)*g(x)-r(x));" }{MPLTEXT 1 0 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#>I\"gG6\"f*6#I\"xGF%F%6$I)operatorGF%I&arrowGF%F%,**$9$ !\"\"\"\"%*$F.!\"#F2-I#lnGF%6#F.\"\"\"#!\"$\"\"#F6F%F%F%" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#>I\"pG6\"f*6#I\"xGF%F%6$I)operat orGF%I&arrowGF%F%,$*$9$!\"\"!\"%F%F%F%" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#>I\"qG6\"f*6#I\"xGF%F%6$I)operatorGF%I&arrowGF%F%, $*$9$!\"#F/F%F%F%" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#> I\"rG6\"f*6#I\"xGF%F%6$I)operatorGF%I&arrowGF%F%,$*&-I#lnGF%6#9$\"\"\" F1!\"#\"\"#F%F%F%" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6# \"\"!" }{TEXT 20 0 "" }}}{EXCHG {PARA 200 "> " 0 "" {MPLTEXT 1 0 51 "a := 1; alpha := 1/2; b := 2; beta := evalf(ln(2));" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 61 "N := 4; h := (b-a)/(N+1); #we \+ assume N >= 2 in this worksheet" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 39 "with(LinearAlgebra): A:= Matrix(N,N,0):" } {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 48 "xx:= Vector(N ): w := Vector(N): bb := Vector(N):" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 67 "for j from 1 to N do xx[j] := a+ j*h; A[j,j] \+ := 2+h^2*q(xx[j]); od:" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 22 "for j from 1 to N-1 do" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 30 "A[j,j+1] := -1 + p(xx[j])*h/2;" } {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 3 "od:" } {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 20 "for j from 2 \+ to N do" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 29 "A[j ,j-1] := -1- p(xx[j])*h/2;" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 3 "od:" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 46 "bb[1] := -h^2*r(xx[1])+(1+h/2*p(xx[1]))*alpha:" } {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 45 "bb[N] := -h^2 *r(xx[N])+(1-h/2*p(xx[N]))*beta:" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 22 "for j from 2 to N-1 do" }{MPLTEXT 1 0 0 "" }} {PARA 200 "> " 0 "" {MPLTEXT 1 0 23 "bb[j] := -h^2*r(xx[j]);" } {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 3 "od:" } {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 23 "w := LinearSo lve(A,bb):" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 34 " u := vector(N+2): v:= vector(N+2):" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> \+ " 0 "" {MPLTEXT 1 0 23 "for j from 1 to N+2 do " }{MPLTEXT 1 0 0 "" }} {PARA 200 "> " 0 "" {MPLTEXT 1 0 19 "u[j] := a+ (j-1)*h;" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 3 "od:" }{MPLTEXT 1 0 0 "" } }{PARA 200 "> " 0 "" {MPLTEXT 1 0 30 "v[1] := alpha: v[N+2] := beta:" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 42 "for j from 2 to N+1 do v[j] := w[j-1]; od:" }{MPLTEXT 1 0 0 "" }{MPLTEXT 1 0 31 " \nwith(student):readlib(spline);" }{MPLTEXT 1 0 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#>I\"aG6\"\"\"\"" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#>I&alphaG6\"#\"\"\"\"\"#" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#>I\"bG6\"\"\"#" }{TEXT 20 0 "" }}{PARA 201 "" 1 " " {XPPMATH 20 "6#>I%betaG6\"$\"+1=ZJp!#5" }{TEXT 20 0 "" }}{PARA 201 " " 1 "" {XPPMATH 20 "6#>I\"NG6\"\"\"%" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#>I\"hG6\"#\"\"\"\"\"&" }{TEXT 20 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#f*6&I#XXG6\"I#YYGF&I\"zGF&I\"dGF&6$I\"XGF&I\"YGF& 6#IfnCopyright~(c)~2000~Waterloo~Maple~Inc.~All~rights~reserved.GF&F&C *@$/9#\"\"$O-I'splineG6$I*protectedGF8I(_syslibGF&6&9$9%9&F3@$/9'.I'li nearGF&O-F66&F;F8$-I(convertGF86$F;.I'VectorGF7>F\\oF;@%-Fgn6$F8%-F^o6$F FgoF " 0 "" {MPLTEXT 1 0 2 "A;" } {MPLTEXT 1 0 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#-I'RTABLEG6\"6%\" ''R^&-I'MATRIXGF%6#7&7&#\"#N\"#=#!\"%\"\"$\"\"!F37&#!\"&\"\"(#\"#'*\"# \\#!\"*F7F37&F3#!\"$\"\"%#\"#j\"#K#F6F@7&F3F3#!\"(\"\"*#\"$g\"\"#\")I' MatrixG6$I*protectedGFNI(_syslibGF%" }{TEXT 20 0 "" }}}{EXCHG {PARA 200 "> " 0 "" {MPLTEXT 1 0 40 "sltn := makeproc(spline(u,v,x,cubic),x) :" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 28 "plot(\{sltn(x),g(x)\},x=a..b);" }{MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 200 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 202 "" 1 "" {GLPLOT2D 400 400 400 {PLOTDATA 2 "6&-%'CURVESG6$7U7$$\"\"\"\"\"!$\"3++++++++]!#=7$$\"3> LL3x&)*3,\"!#<$\"3u5N)G_rg5&F-7$$\"3hmm;arz@5F1$\"3;55T&pGl?&F-7$$\"35 +D\"y%*z7.\"F1$\"3K*)G,ktg*G&F-7$$\"3OL$e9ui2/\"F1$\"3y11qT'e)o`F-7$$ \"3smm\"z_\"4i5F1$\"3i5r;bf/MbF-7$$\"3qmmT&phN3\"F1$\"3cGiSwAf$o&F-7$$ \"3UL$e*=)H\\5\"F1$\"3g>`)oj]u\"eF-7$$\"3sm;z/3uC6F1$\"3q.ERH@_HfF-7$$ \"3-+]7LRDX6F1$\"3*omiQ!3aMgF-7$$\"3em;zR'ok;\"F1$\"3fan!)*4xC8'F-7$$ \"3-+]i5`h(=\"F1$\"3RK&G_1U.A'F-7$$\"3YLL$3En$47F1$\"3Pvt!e7`:I'F-7$$ \"3cmmT!RE&G7F1$\"3%Q#)=;!*4hO'F-7$$\"3)*****\\K]4]7F1$\"3ON0'3946izn$y TlF-7$$\"3bm;z*ev:J\"F1$\"3-M4_#>8Oe'F-7$$\"3ELL$347TL\"F1$\"3,TNa8zFG mF-7$$\"3=LLLjM?`8F1$\"3)[ID8C$>imF-7$$\"3#***\\7o7Tv8F1$\"3A=')e08b(p 'F-7$$\"3ALLLQ*o]R\"F1$\"3];+e]>^DnF-7$$\"3-+]7=lj;9F1$\"3C$4;40FHv'F- 7$$\"3&***\\PaRY2a\"F1$\"3AAJ\\Eh:foF -7$$\"3imm\"zXu9c\"F1$\"3A9(*[MZNqoF-7$$\"3'******\\y))Ge\"F1$\"3))p([ l7Z/)oF-7$$\"3!****\\i_QQg\"F1$\"3#*['\\P1;!*)oF-7$$\"3#***\\7y%3Ti\"F 1$\"3mdHUq^?'*oF-7$$\"3#****\\P![hY;F1$\"3#=HUI(H0.pF-7$$\"3ELLLQx$om \"F1$\"3+!=wuE(H3pF-7$$\"3')****\\P+V)o\"F1$\"3lET=6Q08pF-7$$\"3im;zpe *zq\"F1$\"3iiT7#>&p;pF-7$$\"3)*****\\#\\'QH\\UOH\"HpF-7$$\"3_mmTvJga=F1$\"3fGiu=X&*HpF-7$$ \"3KL$e9tOc(=F1$\"3vk!zr/Y0$pF-7$$\"3'******\\Qk\\*=F1$\"3\"3VgvuE4$pF -7$$\"3@LL3dg6<>F1$\"3(*)f#3\")3@JpF-7$$\"3_mmmw(Gp$>F1$\"3[!G(ec%f8$p F-7$$\"3-+]7oK0e>F1$\"3%f(y_^&R9$pF-7$$\"3-+](=5s#y>F1$\"3?#)z?DuYJpF- 7$$\"\"#F*$\"3(RX*f0=ZJpF--%&COLORG6&%$RGBG$\"#5!\"\"$F*Fi[lFj[l-F$6$7 S7$F($\"3W**************\\F-7$F5$\"33*p,*e(zC;&F-7$F?$\"3+$))*yC&GEI&F -7$FD$\"3#R(=e\\#HvX&F-7$FI$\"3m\">_:c@\"4cF-7$FN$\"3Y\"\\6aJ)=adF-7$F S$\"3Edae5n@#)eF-7$FX$\"3dGzN+S$p+'F-7$Fgn$\"3o?4e-%Qh7'F-7$F\\o$\"3*y %pmn+kLiF-7$Fao$\"3O&zT`*=&4L'F-7$Ffo$\"3g)*yy]SH0kF-7$F[p$\"3w]6u/MRx kF-7$F`p$\"35Nhs.ZvQlF-7$Fep$\"3]:-*Qe6()e'F-7$Fjp$\"3Ay*pR'HKFmF-7$F_ q$\"3K#\\jA%e=mmF-7$Fdq$\"3'Rn**fq.Up'F-7$Fiq$\"3Tl)z#3'eCs'F-7$F^r$\" 3$z25DU6Zu'F-7$Fcr$\"3[Ll\"*GrJnnF-7$Fhr$\"3]5e#GjIuy'F-7$F]s$\"3)Ho6# Gd%p!oF-7$Fbs$\"3!3!yD-*QN#oF-7$Fgs$\"3$G7)*R,2+%oF-7$F\\t$\"3]:&)*\\@ Ob&oF-7$Fat$\"3)Q#=:&=Ux'oF-7$Fft$\"3!3;EN`[&zoF-7$F[u$\"3Es8!=kR-*oF- 7$F`u$\"3[&4D7&H@**oF-7$Feu$\"3BYTzHo_1pF-7$Fju$\"3\"[1;')>2K\"pF-7$F_ v$\"3DFNW1.0=pF-7$Fdv$\"3!>>#H=@:ApF-7$Fiv$\"3UmH?:k.DpF-7$F^w$\"3<$)) p[B=u#pF-7$Fcw$\"3O`0Id&[!HpF-7$Fhw$\"3>QOfDCDIpF-7$F]x$\"3=-P%3]k5$pF -7$Fbx$\"3AN.ZwSmJpF-7$Fgx$\"3dRA\\,%o?$pF-7$F\\y$\"3A$pF-7$Fez$\"3e3`C;:+KpF-7$Fjz$\"3F(33T#fvJpF-7$F_[l$\"3cAo01=ZJ pF--Fd[l6&Ff[lFj[lFg[lFj[l-%+AXESLABELSG6$Q\"x6\"Q!Fgel-%%VIEWG6$;Fg[l $\"#?Fi[l;$\"2cu\"p93Nh\\F1$\"1y4s]%36(p!#;" 1 2 2 0 10 1 2 6 1 4 2 1.0 45.0 45.0 1 0 "Curve 1" "Curve 2" }}{TEXT 202 0 "" }}}{EXCHG {PARA 200 "> " 0 "" {MPLTEXT 1 0 0 "" }}{PARA 201 "" 1 "" {XPPMATH 20 "6#-I'RTABLEG6\"6%\")SO%)=-I'MATRIXGF%6#7$7$\"\"!F-F,I'MatrixG6$I*prot ectedGF0I(_syslibGF%" }{TEXT 20 0 "" }}}{EXCHG {PARA 200 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{PARA 203 "" 0 "" {TEXT 203 0 "" }}{PARA 203 "" 0 "" {TEXT -1 0 "" }}}{MARK "0 0 0" 0 }{VIEWOPTS 1 1 0 2 1 1805 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }{RTABLE_HANDLES 551396 }{RTABLE M7R0 I3RTABLE_SAVE/551396X,%)anythingG6"6"[gl!"%!!!#1"%"%#"#N"#=#!"&""(""!F-#!"%""$# "#'*"#\#!"$""%F-F-#!"*F,#"#j"#K#!"(""*F-F-#F+F6#"$g""#")6" }