{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 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 "" {MPLTEXT 1 0 968 "restart;\nphi:=proc
(x,L)\n       if x>L then \n          0 \n       else\n          if x>
=0 then 1-x else\n             phi(-x,L)\n          end if;\n        e
nd if;\nend proc;\n\nphiLocal:=proc(i,x,h,L)\n     phi(x/h-(i-1),L)\ne
nd proc;\n\nyapprox:=proc(x,c,h,n,L)\n    local i;\n    add(c[i]*phiLo
cal(i,x,h,L),i=1..n)\nend proc;\n\nL:=1; nPoints:=100; nElements:=nPoi
nts-1;\nq:=4; f:=4;\nnShapeFunctions:=nPoints; \nh:=evalf(L/nElements)
;\n#plot([''phiLocal(i,x,h,L)'' $i=1..nPoints],x=0..L);\nwith(LinearAl
gebra):\n\nmakeMatrix:=proc(q,h,f,nPoints)\n local A,i,j;\n A:=Matrix(
nPoints);\n for i from 1 to nPoints do\n   for j from 1 to nPoints do
\n     if i=j then\n        A[i,j]:=2+2/3 * q*h^2;\n     else\n       \+
 if j=i-1 or j=i+1 then\n           A[i,j]:= -1+1/6 *q *h^2;\n        \+
end if;\n     end if;\n   end do;\n end do;\n A\nend proc:\n\n\nA:=mak
eMatrix(q,h,f,nPoints):\nb:=Vector(nPoints):\nb[1..-1]:=h^2*f:\nc:=Lin
earSolve(A,b):\nwith(plots):\n\nlistplot(['yapprox(i*h,c,h,nPoints,L)'
 $i=0..nPoints-1]);\n\n" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$phiGf*6$
%\"xG%\"LG6\"F)F)@%29%9$\"\"!@%1F.F-,&\"\"\"F2F-!\"\"-F$6$,$F-F3F,F)F)
F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%)phiLocalGf*6&%\"iG%\"xG%\"hG%
\"LG6\"F+F+-%$phiG6$,(*&9%\"\"\"9&!\"\"F29$F4F2F29'F+F+F+" }}{PARA 11 
"" 1 "" {XPPMATH 20 "6#>%(yapproxGf*6'%\"xG%\"cG%\"hG%\"nG%\"LG6#%\"iG
6\"F.-%$addG6$*&&9%6#8$\"\"\"-%)phiLocalG6&F69$9&9(F7/F6;F79'F.F.F." }
}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"LG\"\"\"" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%(nPointsG\"$+\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%
*nElementsG\"#**" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"qG\"\"%" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fG\"\"%" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#>%0nShapeFunctionsG\"$+\"" }}{PARA 11 "" 1 "" {XPPMATH 
20 "6#>%\"hG$\"+55555!#6" }}{PARA 13 "" 1 "" {GLPLOT2D 369 369 369 
{PLOTDATA 2 "6#-%'CURVESG6#7`q7$$\"\"\"\"\"!$\"+<7>N:!#67$$\"\"#F*$\"+
$3%>IIF-7$$\"\"$F*$\"+$y=c[%F-7$$\"\"%F*$\"+Z$f?!fF-7$$\"\"&F*$\"+&*Q4
!G(F-7$$\"\"'F*$\"+q[G?')F-7$$\"\"(F*$\"+q#zJ#**F-7$$\"\")F*$\"+()3$*=
6!#57$$\"\"*F*$\"+W!>>C\"FQ7$$\"#5F*$\"+'fK8O\"FQ7$$\"#6F*$\"+\"G?sZ\"
FQ7$$\"#7F*$\"+*RH'*e\"FQ7$$\"#8F*$\"+Heg)p\"FQ7$$\"#9F*$\"+^S>/=FQ7$$
\"#:F*$\"+frV1>FQ7$$\"#;F*$\"+%)oP0?FQ7$$\"#<F*$\"+4O0,@FQ7$$\"#=F*$\"
+%Q1N>#FQ7$$\"#>F*$\"+VHx#G#FQ7$$\"#?F*$\"+?(*))oBFQ7$$\"#@F*$\"+l=*=X
#FQ7$$\"#AF*$\"+aK\"=`#FQ7$$\"#BF*$\"+2lo3EFQ7$$\"#CF*$\"++Ia#o#FQ7$$
\"#DF*$\"+xGT`FFQ7$$\"#EF*$\"+l]K@GFQ7$$\"#FF*$\"+!G2j)GFQ7$$\"#GF*$\"
+YgQ[HFQ7$$\"#HF*$\"++ne2IFQ7$$\"#IF*$\"+0M$R1$FQ7$$\"#JF*$\"+f\"\\u6$
FQ7$$\"#KF*$\"+.e:oJFQ7$$\"#LF*$\"+NS2;KFQ7$$\"#MF*$\"+6MAhKFQ7$$\"#NF
*$\"+gBi.LFQ7$$\"#OF*$\"+'=)GVLFQ7$$\"#PF*$\"+zqB!Q$FQ7$$\"#QF*$\"+@T[
9MFQ7$$\"#RF*$\"+)GVgW$FQ7$$\"#SF*$\"+hu#\\Z$FQ7$$\"#TF*$\"+I%[6]$FQ7$
$\"#UF*$\"+(*orCNFQ7$$\"#VF*$\"+![Uca$FQ7$$\"#WF*$\"+AP$Rc$FQ7$$\"#XF*
$\"+(3)fzNFQ7$$\"#YF*$\"+p>k#f$FQ7$$\"#ZF*$\"+#pqIg$FQ7$$\"#[F*$\"+7&)
)3h$FQ7$$\"#\\F*$\"+?')4;OFQ7$$\"#]F*$\"+UJq=OFQ7$$\"#^F*F\\[l7$$\"#_F
*Fgz7$$\"#`F*Fbz7$$\"#aF*F]z7$$\"#bF*Fhy7$$\"#cF*Fcy7$$\"#dF*F^y7$$\"#
eF*Fix7$$\"#fF*Fdx7$$\"#gF*F_x7$$\"#hF*Fjw7$$\"#iF*Few7$$\"#jF*F`w7$$
\"#kF*F[w7$$\"#lF*Ffv7$$\"#mF*Fav7$$\"#nF*F\\v7$$\"#oF*Fgu7$$\"#pF*Fbu
7$$\"#qF*F]u7$$\"#rF*Fht7$$\"#sF*Fct7$$\"#tF*F^t7$$\"#uF*Fis7$$\"#vF*F
ds7$$\"#wF*F_s7$$\"#xF*Fjr7$$\"#yF*Fer7$$\"#zF*F`r7$$\"#!)F*F[r7$$\"#
\")F*Ffq7$$\"##)F*Faq7$$\"#$)F*F\\q7$$\"#%)F*Fgp7$$\"#&)F*Fbp7$$\"#')F
*F]p7$$\"#()F*Fho7$$\"#))F*Fco7$$\"#*)F*F^o7$$\"#!*F*Fin7$$\"#\"*F*FZ7
$$\"##*F*FU7$$\"#$*F*FO7$$\"#%*F*FJ7$$\"#&*F*FE7$$\"#'*F*F@7$$\"#(*F*F
;7$$\"#)*F*F67$$\"#**F*F17$$\"$+\"F*F+" 1 2 0 1 10 0 2 9 1 4 2 
1.000000 45.000000 45.000000 0 0 "Curve 1" }}}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 0 "" }}}}{MARK "1 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 
1 1 }{PAGENUMBERS 0 1 2 33 1 1 }