{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 }
{CSTYLE "" -1 256 "" 0 1 0 0 0 0 1 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 
257 "" 0 1 0 0 0 0 1 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 258 "" 0 1 0 0 
0 0 0 1 0 0 0 0 0 0 0 0 }{CSTYLE "" -1 259 "" 0 1 0 0 0 0 0 1 0 0 0 0 
0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 14 0 0 0 
0 2 2 2 0 0 0 0 0 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 0 0 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 0 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 
0 -1 0 }{PSTYLE "R3 Font 0" -1 256 1 {CSTYLE "" -1 -1 "Helvetica" 1 
14 0 0 0 0 2 1 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }
{PSTYLE "R3 Font 2" -1 257 1 {CSTYLE "" -1 -1 "Courier" 1 14 0 0 0 0 
2 2 2 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 
258 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }3 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 259 1 {CSTYLE "" -1 -1 "" 0 1 0 
0 0 0 0 1 0 0 0 0 0 0 0 0 }3 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "
" 0 260 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }0 0 0 
-1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 0 261 1 {CSTYLE "" -1 -1 "" 0 1 
0 0 0 0 0 1 0 0 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE 
"" 0 262 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 }0 0 0 
-1 -1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 262 "" 0 "" {TEXT -1 172 "Este ejemplo resuelve u
n problema de frontera en dimension 1 sin hacer una discretizacion del
 dominio. (La aproximacion se hace en el espacion generado por tres fu
nciones)." }}{PARA 260 "" 0 "" {TEXT -1 21 "Resolver la ecuacion " }}
{PARA 258 "" 0 "" {TEXT 256 19 "-u\"+cu=f en [0, 1 ]" }}{PARA 259 "" 
0 "" {TEXT -1 28 "con condiciones de frontera " }{TEXT 257 11 "u(0)=u(
1)=0" }}{PARA 261 "" 0 "" {TEXT -1 110 "El espacio de elementos finito
s esta dado por la base dada en forma de  un vector de dim N de funcio
nes phi[i]" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "with(linalg):Digits:=
10:" }}}{EXCHG {PARA 0 "" 0 "" {TEXT 258 15 "solucion exacta" }}{PARA 
0 "> " 0 "" {MPLTEXT 1 0 41 "c:=4:u:=4*x**2*(x-1);f:=-diff(u,x$2)+c*u;
" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"uG,$*&)%\"xG\"\"#\"\"\",&F(F*F
*!\"\"F*\"\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fG,(%\"xG!#C\"\")
\"\"\"*(\"#;F))F&\"\"#F),&F&F)F)!\"\"F)F)" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 154 "N:=3:phi:=vector(N,0):F:=vector(N,0):A:=matrix(N,N
,0):\nphi[1]:=x*sin(evalf(Pi)*x):phi[2]:=x*(1-cos(x-1.)):phi[3]:=x*(x-
1.):for i to 3 do print(phi[i]):od:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6
#*&%\"xG\"\"\"-%$sinG6#,$F$$\"+aEfTJ!\"*F%" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#*&%\"xG\"\"\",&F%F%-%$cosG6#,&F$F%$F%\"\"!!\"\"F-F%" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#*&%\"xG\"\"\",&F$F%$F%\"\"!!\"\"F%" }}
}{EXCHG {PARA 0 "" 0 "" {TEXT 259 46 "Construccion de la matriz y del \+
segundo miembr" }{TEXT -1 1 "o" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }
{MPLTEXT 1 0 134 "for i to N do for j to N do A[i,j]:=value(int(c*phi[
i]*phi[j]+diff(phi[i],x)*diff(phi[j],x),x=0..1))od: F[i]:=int(f*phi[i]
,x=0..1):od:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 16 "print(A,det
(A));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6$-%'matrixG6#7%7%$\"+^&z-Y#!\"
*$\"+m*zwC\"!#5$!+w/KY*)F-7%F+$\"(!f@R!\")$!*e6'H6F*7%F.F4$\"+nmmmYF-$
\")m6t>!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 26 "sol:=evalf(li
nsolve(A,F));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$solG-%'vectorG6#7%
$!*[(*ow\"!\"*$\"+B_b;sF+$\"+l-63MF+" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 37 "solcal:=sum(sol[ii]*phi[ii],ii=1..3);" }}{PARA 11 "" 
1 "" {XPPMATH 20 "6#>%'solcalG,(*&%\"xG\"\"\"-%$sinG6#,$F'$\"+aEfTJ!\"
*F($!*[(*ow\"F/*($\"+B_b;sF/F(F'F(,&F(F(-%$cosG6#,&F'F($F(\"\"!!\"\"F<
F(F(*($\"+l-63MF/F(F'F(F9F(F(" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 100 "p1 := plot(solcal, x=0..1, style=point):\np2 := plot(u, x=0..1,
 style=line):\nplots[display](\{p1,p2\});" }}{PARA 13 "" 1 "" 
{GLPLOT2D 399 316 316 {PLOTDATA 2 "6&-%'CURVESG6%7W7$$\"\"!F)F(7$$\"3d
mmm;arz@!#>$!3$3>*Q=!R!f=!#?7$$\"3[LL$e9ui2%F-$!3K.g/-(yaP'F07$$\"3nmm
m\"z_\"4iF-$!3TCU07\"*QY9F-7$$\"3[mmmT&phN)F-$!3pTG$y0L'fDF-7$$\"3BLLe
*=)H\\5!#=$!3LbKPPY)>%RF-7$$\"3gmm\"z/3uC\"FC$!3PQs#=>2xW&F-7$$\"3%)**
*\\7LRDX\"FC$!3]Ce9wSh8sF-7$$\"3^mm\"zR'ok;FC$!3JCTz<WYR#*F-7$$\"3w***
\\i5`h(=FC$!3z0RI]:#Q9\"FC7$$\"3YLLL3En$4#FC$!30MR6tVG'Q\"FC7$$\"3qmm;
/RE&G#FC$!3]6M@7re6;FC7$$\"3\")*****\\K]4]#FC$!3,Iu_'*z=w=FC7$$\"3$***
***\\PAvr#FC$!3w:sqvGA^@FC7$$\"3(******\\nHi#HFC$!3oGz$>\"p&GU#FC7$$\"
3jmm\"z*ev:JFC$!35jP&)><FtEFC7$$\"3?LLL347TLFC$!39%y/uPYL(HFC7$$\"3,LL
LLY.KNFC$!3PD$)*HA%eFKFC7$$\"3v***\\7o7Tv$FC$!3Ps2lX?-@NFC7$$\"3'GLLLQ
*o]RFC$!3ptnt$o#pwPFC7$$\"3A++D\"=lj;%FC$!3oVo1g!\\00%FC7$$\"31++vV&R<
P%FC$!3@q(Q6c;FI%FC7$$\"3XLL$e9Ege%FC$!3KBc[b!*eaXFC7$$\"3HLLeR\"3Gy%F
C$!3ZFxN'\\$ytZFC7$$\"3emm;/T1&*\\FC$!3WtKOO#f]*\\FC7$$\"3&em;zRQb@&FC
$!36xJu3l%e?&FC7$$\"3\\***\\(=>Y2aFC$!3d\\Lf<4br`FC7$$\"39mm;zXu9cFC$!
3:Q]h_'p)HbFC7$$\"3l******\\y))GeFC$!3+)o-D<(pocFC7$$\"3'*)***\\i_QQgF
C$!3n&*=gm6&zx&FC7$$\"3@***\\7y%3TiFC$!3f8/4/7cceFC7$$\"36****\\P![hY'
FC$!3!Gy\\NLl,\"fFC7$$\"3jKLL$Qx$omFC$!3yf]5bZ#f#fFC7$$\"3!)*****\\P+V
)oFC$!3&yEa'Hyc1fFC7$$\"3?mm\"zpe*zqFC$!3uZI?B!yZ&eFC7$$\"3%)*****\\#
\\'QH(FC$!3Y;OhKeqedFC7$$\"3GKLe9S8&\\(FC$!3;)4kZlP'GcFC7$$\"3R***\\i?
=bq(FC$!3bym$QD&R\\aFC7$$\"3\"HLL$3s?6zFC$!3gk(ygkn#H_FC7$$\"3a***\\7`
Wl7)FC$!32t[xsV(*[\\FC7$$\"3#pmmm'*RRL)FC$!3&GuiKN='GYFC7$$\"3Qmm;a<.Y
&)FC$!3->pJ')HgZUFC7$$\"3<LLe9tOc()FC$!3@$Gri8qT\"QFC7$$\"3t******\\Qk
\\*)FC$!3PXZ:,$y^O$FC7$$\"3DLL$3dg6<*FC$!3@C^@\\Nb)y#FC7$$\"3x***\\(oT
Aq#*FC$!3%\\1dJY\"f3DFC7$$\"3HmmmmxGp$*FC$!3P>g#pVYY@#FC7$$\"3sK$eRA5
\\Z*FC$!36\\Up!yub)=FC7$$\"3B++D\"oK0e*FC$!3Qn()Ho#f+a\"FC7$$\"3C+++]o
i\"o*FC$!3;XKhgUp$>\"FC7$$\"3B++v=5s#y*FC$!3/V\"pU6\"e<$)F-7$$\"37+]P4
0O\"*)*FC$!3+-1M-An^UF-7$$\"\"\"F)F(-%'COLOURG6&%$RGBG$\"#5!\"\"F(F(-%
&STYLEG6#%%LINEG-F$6%7WF'7$F+$!3,(oz<U+b[$F07$F2$!3%>^d_+.%3*)F07$F7$!
3[$R!z+Gua<F-7$F<$!3axoT*G'o#)GF-7$FA$!3()fe:E@#zC%F-7$FG$!3h1N6H$)Q;d
F-7$FL$!3VYYGe$z\"GuF-7$FQ$!3yXg:.TA(Q*F-7$FV$!3'**[`zXj8:\"FC7$Fen$!3
?7,v-.E'Q\"FC7$Fjn$!3p!)4<ez60;FC7$F_o$!3!z1H6)y+j=FC7$Fdo$!3c\\4\"H,O
@8#FC7$Fio$!3m'Gu*e00*R#FC7$F^p$!3pq@j/46YEFC7$Fcp$!3HKc?[sUVHFC7$Fhp$
!3I2ZYX!)R'>$FC7$F]q$!3#[?99'4e*[$FC7$Fbq$!3K$3G\\Psgu$FC7$Fgq$!3Gzdt8
L(=-%FC7$F\\r$!3=!*Gqt&)zwUFC7$Far$!39WG,&y\\B`%FC7$Ffr$!3H)H&)oX.bv%F
C7$F[s$!39&GI5S3:)\\FC7$F`s$!3i'[96?>v>&FC7$Fes$!34sp)pjmyO&FC7$Fjs$!3
%RE<,eK6`&FC7$F_t$!3BI#GS*=#[n&FC7$Fdt$!3TfA#)z$>%)y&FC7$Fit$!3KMm8Y+q
qeFC7$F^u$!3Y:ABn+iFfFC7$Fcu$!3.w9!*4&pb%fFC7$Fhu$!3(p![fe*Rw#fFC7$F]v
$!3.+ZmkaGweFC7$Fbv$!3rONFZ)Q(zdFC7$Fgv$!3?@h%)GmL[cFC7$F\\w$!3B9E\"*y
!RoY&FC7$Faw$!3Lu<^[avV_FC7$Ffw$!3@Y[Tlqsf\\FC7$F[x$!3W3z4nDNNYFC7$F`x
$!3M=YxUc0]UFC7$Fex$!3CtUd+:^7QFC7$Fjx$!3'>!*Hmcp,O$FC7$F_y$!3UUv]okj!
y#FC7$Fdy$!3W'30,PJ)*\\#FC7$Fiy$!3)f<o)R3S0AFC7$F^z$!38dk7A\"pi(=FC7$F
cz$!3!Q.1saO7`\"FC7$Fhz$!3oW4'4O7f=\"FC7$F][l$!3U#G.WPOlD)F-7$Fb[l$!3_
@`D%>#e;UF-7$Fg[l$\"3yu>m>Y$zC(!#GFi[l-Fa\\l6#%&POINTG-%+AXESLABELSG6%
Q\"x6\"Q!F[gl-%%FONTG6#%(DEFAULTG-%%VIEWG6$;F(Fg[lF`gl" 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 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "
" {MPLTEXT 1 0 0 "" }}}}{MARK "0 0 0" 172 }{VIEWOPTS 1 1 0 1 1 1803 1 
1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }