{VERSION 2 3 "IBM INTEL NT" "2.3" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 114 32 105 0 0 0 2 0 0 
0 0 0 0 }{CSTYLE "2D Comment" 2 18 "" 0 1 200 136 1 0 0 0 0 0 0 0 0 0 
0 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }
{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 194 136 40 0 0 0 0 0 
0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 
2 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 0 }1 0 
0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Warning" 2 7 1 {CSTYLE "" -1 -1 
"" 0 1 0 0 255 1 0 0 0 0 0 0 1 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }
{PSTYLE "Maple Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 24 0 0 0 0 0 0 
0 0 0 0 0 0 }3 3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple Plot" 0 
13 1 {CSTYLE "" -1 -1 "" 0 1 202 136 1 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg):" }}
{PARA 7 "" 1 "" {TEXT -1 32 "Warning, new definition for norm" }}
{PARA 7 "" 1 "" {TEXT -1 33 "Warning, new definition for trace" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 59 "De
monstration of finite difference method for\n             " }{XPPEDIT 
18 0 "d^2y/dx^2+y=0" "/,&*(%\"dG\"\"#%\"yG\"\"\"*$%#dxGF&!\"\"F(F'F(\"
\"!" }{TEXT -1 7 ",      " }{XPPEDIT 18 0 "x" "I\"xG6\"" }{TEXT -1 9 "
  in  [0," }{XPPEDIT 18 0 "Pi/4" "*&%#PiG\"\"\"\"\"%!\"\"" }{TEXT -1 
1 "]" }}{PARA 0 "" 0 "" {TEXT -1 12 "            " }{XPPEDIT 18 0 "y(0
)=1" "/-%\"yG6#\"\"!\"\"\"" }{TEXT -1 6 ",     " }{XPPEDIT 18 0 "y(Pi/
4)=1" "/-%\"yG6#*&%#PiG\"\"\"\"\"%!\"\"F(" }}{PARA 0 "" 0 "" {TEXT -1 
25 "(The exact solution is   " }{XPPEDIT 18 0 "y=cos(x)+(sqrt(2)-1)*si
n(x)" "/%\"yG,&-%$cosG6#%\"xG\"\"\"*&,&-%%sqrtG6#\"\"#F)F)!\"\"F)-%$si
nG6#F(F)F)" }{TEXT -1 3 "  )" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 
0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 6 "n:=20:" }}
}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 22 "a:=0:  b:=evalf(Pi/4):" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 18 "h:=evalf((b-a)/n):" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "m:=n-1:" }}}{EXCHG {PARA 0 ">
 " 0 "" {MPLTEXT 1 0 15 "A:=matrix(m,m):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 59 "for i to m do \n   for j to m do\n      A[i,j]:=0.0\n
   od\nod:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "for i to m do
 A[i,i]:=-2.0+h^2 od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "fo
r i to m-1 do A[i,i+1]:=1.0 od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 33 "for i to m-1 do A[i+1,i]:=1.0 od:" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 13 "x:=vector(m):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 29 "for i to m do x[i]:=a+i*h od:" }}}{EXCHG {PARA 0 "> \+
" 0 "" {MPLTEXT 1 0 13 "b:=vector(m):" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 25 "for i to m do b[i]:=0 od:" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 11 "b[1]:=-1.0:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 11 "b[m]:=-1.0:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "y:=l
insolve(A,b):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "P:=[seq([x
[i],y[i]],i=1..m)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "plot(
P);" }}{PARA 13 "" 1 "" {INLPLOT "6#-%'CURVESG6$757$$\"1+++=3*p#R!#<$
\"1+++%4$\\:5!#:7$$\"1+++O;)R&yF*$\"1+++q,UH5F-7$$\"1+++Xs4y6!#;$\"1++
+](f</\"F-7$$\"1+++Fjzq:F6$\"1+++0G\\_5F-7$$\"1+++4a\\j>F6$\"1+++%y-;1
\"F-7$$\"1+++\"\\%>cBF6$\"1+++Qc2p5F-7$$\"1+++tN*)[FF6$\"1+++U)**[2\"F
-7$$\"1+++aEfTJF6$\"1+++:k1z5F-7$$\"1+++O<HMNF6$\"1+++K*o:3\"F-7$$F)F6
$\"1+++LNS#3\"F-7$$\"1++++**o>VF6FU7$$\"1+++#)*)Q7ZF6FP7$$\"1+++j!)30^
F6FK7$$\"1+++Xry(\\&F6FF7$$\"1+++Fi[!*eF6FA7$$\"1+++4`=$G'F6F<7$$\"1++
+\"R%)en'F6F77$$\"1+++sMeoqF6F17$$\"1,++aDGhuF6F+-%'COLOURG6&%$RGBG$\"
#5!\"\"\"\"!Fgp" 2 350 350 350 2 0 1 0 2 9 0 4 2 1.000000 45.000000 
45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 8276 8254 0 0 0 
0 0 0 }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "f:=x->cos(x)+(sqrt(
2)-1)*sin(x);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fG:6#%\"xG6\"6$%)
operatorG%&arrowGF(,&-%$cosG6#9$\"\"\"*&,&-%%sqrtG6#\"\"#F1!\"\"F1F1-%
$sinGF/F1F1F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "z:=vecto
r(m):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 38 "for i to m do z[i]
:=evalf(f(x[i])) od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "Q:=
[seq([x[i],y[i]-z[i]],i=1..m)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 8 "plot(Q);" }}{PARA 13 "" 1 "" {INLPLOT "6#-%'CURVESG6$757$$\"1++
+=3*p#R!#<$\"1++++++5@!#@7$$\"1+++O;)R&yF*$\"1++++++:SF-7$$\"1+++Xs4y6
!#;$\"1++++++4dF-7$$\"1+++Fjzq:F6$\"1++++++(=(F-7$$\"1+++4a\\j>F6$\"1+
+++++Y%)F-7$$\"1+++\"\\%>cBF6$\"1++++++#[*F-7$$\"1+++tN*)[FF6$\"1+++++
+H5!#?7$$\"1+++aEfTJF6$\"1+++++!p3\"FM7$$\"1+++O<HMNF6$\"1+++++q@6FM7$
$F)F6$\"1+++++IL6FM7$$\"1++++**o>VF6FV7$$\"1+++#)*)Q7ZF6FQ7$$\"1+++j!)
30^F6FK7$$\"1+++Xry(\\&F6FF7$$\"1+++Fi[!*eF6$\"1,+++++Z%)F-7$$\"1+++4`
=$G'F6$\"1++++++)=(F-7$$\"1+++\"R%)en'F6F77$$\"1+++sMeoqF6F17$$\"1,++a
DGhuF6F+-%'COLOURG6&%$RGBG$\"#5!\"\"\"\"!F\\q" 2 350 350 350 2 0 1 0 
2 9 0 4 2 1.000000 45.000000 45.000000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
0 0 1 1 0 0 0 16365 0 0 0 0 0 0 0 }}}}{MARK "0 0 0" 13 }{VIEWOPTS 1 1 
0 1 1 1803 }