{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 "Warning" -1 7 1 {CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 2 2 2 2 2 1 1 1 3 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 Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 3 0 0 0 0 1 0 1 0 2 2 0 1 }{PSTYLE "Maple P lot" -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 2 "Q1" }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "restart:w:=1/(1+X): A:=1.: \+ B:=2.:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 171 "for i from 0 to 2 do ps i[i]:=X^i-add(int(psi[i-k]*X^i*w,X=A..B)/alpha[i-k]*psi[i-k],k=1..i): \npsi[i]:=expand(psi[i]):alpha[i]:=int(w*psi[i]^2,X=A..B): print(psi[i ]) end do:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&%\"XG\"\"\"$\"+iMIm9!\"*!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,(*$)%\"XG\"\"#\"\"\"F($\"+@mCE@!\"*F(*&$\"+/L0tHF+F(F& F(!\"\"" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 141 "xn:=[-coeff(psi [2],X,1)/2-sqrt(coeff(psi[2],X,1)^2/4-coeff(psi[2],X,0)),\n-coeff(psi[ 2],X,1)/2+sqrt(coeff(psi[2],X,1)^2/4-coeff(psi[2],X,0))];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#xnG7$$\"+ls`(>\"!\"*$\"+Rg^v " 0 "" {MPLTEXT 1 0 68 "lagra:=proc(xn,k,n);mul(x-xn[i],i=1 ..k-1)*mul(x-xn[i],i=k+1..n) end;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#> %&lagraGf*6%%#xnG%\"kG%\"nG6\"F*F**&-%$mulG6$,&%\"xG\"\"\"&9$6#%\"iG! \"\"/F5;F1,&9%F1F1F6F1-F-6$F//F5;,&F:F1F1F19&F1F*F*F*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 140 "lagrpol:=proc(xn) local n,pol,aux, i;n:=nops(xn);pol:=0;for i to n do aux:=lagra(xn,i,n);\npol:=pol+aux/e val(aux,x=xn[i])*y(xn[i]) end do end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "lagrpol(xn);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&*($ \"+Vt;I\"F'F(F- *($\"+Vt;I\"F'F-F(-F/6#$\"+Rg^v " 0 "" {MPLTEXT 1 0 28 "rule:=int(%/(1.+x),x=1..2.);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%%ruleG,&*&$\"+#y&>p@!#5\"\"\"-%\"yG6 ##\"*`u]R#\"*++++#F*F**&$\"+/`X&)=F)F*-F,6##\"+Rg^v<\"+++++5F*F*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 14 "y:=x->exp(-x):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "evalf(rule);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+B$)[V(*!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 72 "evalf(int(exp(-x)/(1+x),x=1..2)); # 'exact' value, just for co mparison" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+S#>cu*!#6" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q4:" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 94 "F:=[x[1]^2/(x[2]+sin(x[1]+x[2]))-ln(x[1]^2+x[2]^4),\n (x[2]-x[1]^2)*exp(-x[1])-3*x[1]+4*x[2]^2];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"FG7$,&*&&%\"xG6#\"\"\"\"\"#,&&F)6#F,F+-%$sinG6#,&F( F+F.F+F+!\"\"F+-%#lnG6#,&*$)F(F,F+F+*$)F.\"\"%F+F+F4,(*&,&F.F+F9F4F+-% $expG6#,$F(F4F+F+*&\"\"$F+F(F+F4*&F=F+)F.F,F+F+" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}{PARA 7 "" 1 "" {TEXT -1 50 "W arning, the name changecoords has been redefined\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 63 "implicitplot(\{F[1],F[2]\},x[1]=-4..4,x[2]= -4..4,numpoints=2000);" }}{PARA 13 "" 1 "" {GLPLOT2D 494 494 494 {PLOTDATA 2 "6$-%'CURVESG6]u7$7$$!3(>AAAAAA#Q!#<$\"3uh<#4-4:P$F*7$$!3/ gEzM\\6?RF*$\"3ojmmmmmmMF*7$7$F.$\"37kmmmmmmMF*7$$!\"%\"\"!$\"3,IiT$)Q xTNF*7$7$$!3#[WWWWWWk$F*$\"376P$QAwL?$F*7$$!3&oQx))fQKt$F*$\"33')))))) ))))))G$F*7$FB7$$!3TAAAAAAAQF*F+7$7$$!3!ommmmmmY$F*$\"3g#4871Zu.$F*7$$ !3jE#)*zkDQa$F*$\"3]3666666JF*7$7$$!33F#)*zkDQa$F*FTF=7$7$$!3?*))))))) )))))G$F*$\"3kVF?))y\"Q(GF*7$$!3Yp&eG8=i7FF*7$$!3%y!>5%yOn:$F*$\"3K`bbbbbbFF*7$7$ $!3S2>5%yOn:$F*Fio7$$!3k*))))))))))))G$F*Fhn7$7$$!3YMLLLLLLHF*$\"3Q^#H <3\\Sb#F*7$$!3usmH7aneHF*$\"3svxxxxxxDF*7$FhpFao7$7$$!3)obbbbbbv#F*$\" 3KJ8QIsP)R#F*7$$!3]pb!\\A$HdFF*$\"37)************R#F*7$7$$!31pb!\\A$Hd FF*Fgq7$Fdp$\"3%4DH<3\\Sb#F*7$F_q7$$!3v6]&3+$QVFF*$\"3Y`%*HXu#yQ#F*7$7 $$!3GzxxxxxxDF*$\"30aa(orXKC#F*Far7$7$$!3o,++++++CF*$\"3u*ox8$*)o!4#F* 7$$!3!)HnvQyq^DF*$\"3a?AAAAAAAF*7$FbsFgr7$7$$!35CAAAAAAAF*$\"3]yPt.GIT >F*7$$!3-Sf5lnnTBF*$\"3&HWWWWWW/#F*7$F^tF]s7$7$$!3^YWWWWWW?F*$\"3!e+z) QP!ez\"F*7$$!3iOF\"p3+l7#F*$\"3#\\mmmmmm'=F*7$7$F[u$\"39lmmmmmm=F*Fis7 $7$$!3qommmmmm=F*$\"3_R7!Rwral\"F*7$$!3@T=sC340>F*$\"3K())))))))))))o \"F*7$FiuFet7$7$Feu$\"3uR7!Rwral\"F*7$$!3ytQ2,mnHNI>Z0ZFiy7$7$$!3/#4D]-\\cy(Fiy$!3'>YWWWWWW%FiyFg\\l7 $F]]l7$$!39V8rn;4=xFiy$!3!*e`&*)*\\d[HFiy7$7$$!35mr;r*o,m(Fiy$!3%Qommm mmm#FiyFc]l7$Fi]l7$$!3?Nf$o())[kuFiy$!3s))H07+SC9Fiy7$7$$!3>19dp\"y%3s Fiy$!3%*f!*))))))))))))!#>F_^l7$Fe^l7$$!3723.P!f--(Fiy$!3YvQ.3u?&3*!#? 7$7$$!35::27\"o0F'Fiy$\"3a=()))))))))))))Fj^lF\\_l7$Fc_l7$$!3#ets,n^wD 'Fiy$\"3%[m$RoL=V#*Fj^l7$7$Fe[l$\"3m\"[wv_*z)[*Fj^lFi_l7$F^[l7$$!3eiUr ([/^J'Fiy$\"3U%ypK/g1()*Fiy7$7$Fe[l$\"33e\"Qob#y\\)*FiyFc`l7$Fd[l7$$!3 $p6hn$oYThFiy$!3ES6Y&Qb23)Fiy7$7$F`]l$!3UCRG,w95\"*FiyF]al7$F_`l7$$!33 \"*)oVItl=&Fiy$\"3(4I8)[x,J;Fiy7$7$F`]l$\"3sY\\XCR^\"f\"FiyFgal7$Fi`l7 $$!3w,*>\"y`a?XFiy$\"3;,KX6(yQ&)*Fiy7$7$F`]l$\"3q%R@N2xU%)*FiyFabl7$7$ F`]l$!3`DRG,w95\"*Fiy7$$!3))\\d!*=`erFFiy$!3%3tQb7fGn*Fiy7$7$F\\^l$!3$ 4%*\\*[N/:(*FiyF^cl7$7$F`]l$\"3WY\\XCR^\"f\"Fiy7$$!3l*pQN0jwB$Fiy$\"37 ()3wv_))f9Fiy7$7$F\\^l$\"3E,'[/SR4?\"FiyF[dl7$7$F`]l$\"3f$R@N2xU%)*Fiy 7$$!3rAMTc)\\k\"GFiy$\"3N,X_n4cF**Fiy7$7$F\\^l$\"3Fi;[;HG<**FiyFhdl7$7 $Fh^l$!3g**o7gI?i**Fiy7$$!3GufLv:LgAFiyFgy7$FeelFdcl7$7$Fh^l$\"31+&z**FiyFhfl7$Fbel7$$!3!fdqb4Yp)pFj^l$!3=W'4r0sz'**Fiy7$7$$\"3#*>( )))))))))))))Fj^l$!31;gm!R<$e**FiyFbgl7$Fjel7$$\"3(3(o'*Qi*[j\"Fa_l$!3 u\"Hq*Qi*[j\"Fa_l7$7$Ff_l$!30M(Rg\"po9XFj^lF^hl7$7$Fh^l$\"3%pS/J>+&z** Fiy7$$\"3%f&f%p[)RwnFj^l$\"3uQ?(z\"o-*)**Fiy7$7$Ff_l$\"31&o9`=>E)**Fiy F[il7$7$Ff_l$!3'\\,m1R<$e**Fiy7$$\"3av#4];$3g5Fiy$!3s+#)*Q0s*[**Fiy7$7 $$\"3h\"p%>#)f#*3AFiyFgyFhil7$F^jl7$$\"3+6\\QEM<.EFiy$!3Zeg\\PXG9(*Fiy 7$7$$\"3U\\mmmmmmEFiy$!3?@O-R7)))o*FiyFbjl7$7$Fijl$!3CE1\"R!)\\gZ\"Fiy 7$$\"3iIdB0_eA>FiyFh^l7$7$$\"3MIdB0_eA>FiyFh^lFdhl7$7$Figl$\"3='o9`=>E )**Fiy7$$\"3TLS$Rq==`#Fiy$\"3#oP50uDE\"**Fiy7$7$Fijl$\"3&Qdn*=h34**Fiy F\\\\m7$Fhjl7$$\"39%*=A#32n$RFiy$!3#HEbbTS+F*Fiy7$7$$\"3bFWWWWWWWFiy$! 3\\#Qn!Qb09*)FiyFf\\m7$7$F]]m$!3Qm%)*\\9P,'GFiy7$$\"3A[r+KJw+VFiyF\\^l 7$7$$\"3x[r+KJw+VFiyF\\^lF^[m7$7$F]]m$\"3y()32Ys^g(*Fiy7$$\"3')>!z3XFc F%Fiy$\"3-jxxxxxx(*Fiy7$7$$\"3I>!z3XFcF%Fiy$\"3#>wxxxxxx*FiyFb\\m7$F\\ ]m7$$\"3_;s-fpyW]Fiy$!3>2Ge9DM+')Fiy7$7$$\"3c,%4\"Q%**)ycFiyF_[lF\\_m7 $Fb_m7$$\"3q,ND[(*pgeFiy$!3888.EvZQwFiy7$7$$\"307?v\\pD@iFiyFe[lFf_m7$ F\\`m7$$\"3^Ph))GG\"3A'Fiy$!3$><'))GG\"3A'Fiy7$7$$\"3#psEm<\")Hv&FiyF` ]lF``m7$Ff`mFb]m7$F]^m7$$\"3MxM21]inWFiy$\"3m7([h@(fa(*Fiy7$7$$\"3n0AA AAAAiFiy$\"31&owe4#*=_*FiyF[am7$7$Fbam$\"3<'owe4#*=_*Fiy7$$\"3Wd2qIH'y c'Fiy$\"3D6#*Hpq8K%*Fiy7$7$$\"3m#)************zFiy$\"31xnHb3kn#*FiyFja m7$7$$\"3y$)************zFiyFcbm7$$\"3(\\L3k@FOj)Fiy$\"3j4%p8c]T9*Fiy7 $7$Fc^m$\"3U\"fd\\g7&>\"*FiyFibm7$F_cm7$$\"3G)eCsi0m/\"F*$\"33T'4LG*\\ *3*Fiy7$7$$\"37abbbbbb6F*$\"3C6Sq!z8fG*FiyFccm7$7$F_zF\\dm7$$\"3A>\\8[ \\v&>\"F*$\"3Y2T)>&Qyv$*Fiy7$7$$\"3#328\"['=.F\"F*Fi^mF`dm7$Ffdm7$$\"3 k(e \"F*$\"3;P,M$HasD\"F*7$7$$\"36y'RNWUyi\"F*FcxFhfm7$7$$\"3Ly'RNWUyi\"F* Fcx7$$\"3ofc5CEr`;F*$\"3Ofl6)f4&o8F*7$7$F\\v$\"3>_j.px+O9F*Fegm7$7$F\\ v$\"3T_j.px+O9F*7$$\"3#>9#>o1c9UeOnt?x\"F*$\"3%HN>T5/dg\"F*7$7$$\"3]xG#z#o-1 =F*F\\vF\\im7$Fbim7$$\"3*4\"p$[\\Oh#=F*$\"3FT'=21>%H3y_=F*7$7$$\"3U&=%Q!\\!Q ')=F*FauF`jm7$Ffjm7$$\"3s@Pe8'3'Q>F*$\"3Q'QFv\\-D(>F*7$7$$\"3K([`.Kid( >F*FatFjjm7$F`[n7$$\"3u^%z`3?6+#F*$\"3sL%4N!)ox3#F*7$7$Fat$\"3/0;&fF(e `@F*Fd[n7$7$Fat$\"3g/;&fF(e`@F*7$$\"3=[vK())G62#F*$\"3I:\"R$zx`&>#F*7$ 7$$\"3/Y6i[94\"4#F*FesFa\\n7$7$$\"3[Y6i[94\"4#F*Fes7$$\"3ud,c#f-0:#F*$ \"3M$G%)=&=%RH#F*7$7$Fes$\"3#f&zf*)[BkBF*F^]n7$Fd]n7$$\"3YEeYV9xRAF*$ \"3A#Rc(y2X#Q#F*7$7$$\"3A&*>wE\\6hAF*FgqFh]n7$F^^n7$$\"3;R))Qw.:SBF*$ \"3mc6hB'\\)fCF*7$7$Fgq$\"3qu9u09V%\\#F*Fb^n7$Fh^n7$$\"3Ic-#Q7A5X#F*$ \"3+=v&Rlbn_#F*7$7$$\"3Gt8-g&o!eDF*F[qF\\_n7$Fb_n7$$\"3CVyJlBqrDF*$\"3 m3xB!>`Qe#F*7$7$F[q$\"3+/zp')=!ee#F*Ff_n7$F\\`n7$$\"36=Fk[*)\\,FF*$\"3 P61p%QM=j#F*7$7$Fio$\"3stfPmjDYEF*F``n7$Ff`n7$$\"3v'\\pKnZ&QGF*$\"3))4 ;%yVjDn#F*7$7$F]o$\"3iD;(*fZL%p#F*Fj`n7$F`an7$$\"3[S@#>u;7)HF*$\"3>Wn' p9swq#F*7$7$FT$\"3T_Dh!GoYt#F*Fdan7$Fjan7$$\"3%=7B*>:vFJF*$\"3)*RNuY^ \"*QFF*7$7$$\"3fM'\\-QH4@$F*FioF^bn7$Fdbn7$$\"3#*=\\&or\\kF$F*$\"3[?&* eFZ*zw#F*7$7$FE$\"3/!QX2.#*)pFF*Fhbn7$7$$\"3k&))))))))))))G$F*$\"3gz`u I?*)pFF*7$$\"30/+eY_hDMF*$\"3]7Akvpg'z#F*7$7$F0$\"3cc2!>llO!GF*Fgcn7$7 $$\"3CjmmmmmmMF*$\"37c2!>llO!GF*7$$\"3eo>vQ3_tNF*$\"3,E![7;zk#GF*7$7$$ \"3FTWWWWWWOF*$\"3!pPg&ot'4%GF*Ffdn7$F\\en7$$\"3%Q@x0Mg'=PF*$\"3Ke0?Pu 6fGF*7$7$$\"3T=AAAAAAQF*$\"3')=Ye!>]])GF*Fben7$Fhen7$$\"3]0>]IkmfQF*$ \"3#Qk`]7*)e*GF*7$7$$\"3SPkniIuzRF*F]oF^fn7$Fdfn7$$\"3#\\y!)pFU`*RF*$ \"3(Ga_j0\"*z$HF*7$7$$\"3W'*************RF*$\"3\"p_7;Xl&RHF*Fhfn-%'COL OURG6&%$RGBG\"\"\"F9F9-F$6fx7$7$F`q$\"3)))GnJdi,m$F*7$$!3Sz;r)))G$*z#F *$\"3I>AAAAAAQF*7$F^hn7$$!3=e\"HC)>ze4\\'eW)QF*7$7$$!3%>j@U2 o0&GF*F_gnFdhn7$7$Fhr$\"37tD-))\\E8JF*7$$!3#G;l\"yGXIEF*FE7$Fain7$$!3I cR*[.'[cEF*$\"3mi]+YrfnLF*7$7$$!3926'[18zo#F*FbdnFein7$7$F\\jnF07$$!3S cS:e/(ou#F*$\"3CTH/Z$fdj$F*7$7$$!3%op7A$>y\\FF*F]enF`jn7$7$Fgjn$\"3rTW WWWWWOF*7$$!3UcbbbbbbFF*F\\hn7$7$F^s$\"3;pNYnq!zh#F*7$$!3k+^?&3$3XCF*F io7$Fd[o7$$!3#*fU>vR%3Z#F*$\"3a6)\\2`*REGF*7$7$$!3!R?;@3=#3DF*F]oFh[o7 $F^\\o7$$!3SM#o>ppjd#F*$\"3ij:IDIq4JF*7$7$$!3?]ei(=%*od#F*FTFb\\o7$7$$ !3m]ei(=%*od#F*$\"3%*3666666JF*7$Fhr$\"3ctD-))\\E8JF*7$7$Fjs$\"3[D\"33 TJ1<#F*7$$!3Reav_#=-C#F*Fes7$Fh]o7$$!3e3!)Qopz^AF*$\"3/0!)Qopz^AF*7$7$ $!3;&=[G@>wI#F*$\"3o(************R#F*F\\^o7$7$Fc^oFgq7$$!3C-qP\"fc&pBF *$\"3IwZ:pVLZDF*7$7$$!3KQcM.L\">Q#F*F[qFi^o7$7$F`_o$\"3;wxxxxxxDF*7$F^ s$\"3soNYnq!zh#F*7$7$$!3GZ=y\"fm-6#F*F77$$!3YW^#)R^)))4#F*$!3[-$>YIfb% RF*7$7$$!3e69LBb%42#F*FIF]`o7$Fc`o7$$!3!fY]I)zbm?F*$!3-.ih$o3,!QF*7$7$ Fft$!3#GKh[4@xp$F*Fg`o7$7$Fft$\"3#4jpZn'Hm$Q3#F*$\"3Olm mmmmm=F*7$Fdao7$$!3M%o^%)4y`6#F*$\"3?.Rn?.gP>F*7$7$$!3!yXy6>c7;#F*FatF jao7$F`boFe]o7$F]ao7$$!3:)*>^XBxK?F*$!3u#*oPVl6cOF*7$7$$!3s[hb?cDH?F*F >Febo7$F[co7$$!3[[RFiZ,&*>F*$!3/lr$)[j4;NF*7$7$$!3-?&*)4x'z!)>F*$!3Cnm mmmmmMF*F_co7$Feco7$$!3))o#**R6%yd>F*$!3#o1M$>#\\bP$F*7$7$$!33#y@$>)zS $>F*F`pF[do7$Fado7$$!3p+_))\\f7@>F*$!3Ud.n0'HWB$F*7$7$$!3ejLGXP>*)=F*F boFedo7$F[eo7$$!3A]%*oC74&)=F*$!3sI$)3`lo#4$F*7$7$$!3#*ommmmmm=F*$!3#* =hVe\"Gr,$F*F_eo7$7$Feu$\"3sUkp#R_0S\"F*7$$!3+v'3r\"od9>F*$\"3u4666666 :F*7$7$$!3yu'3r\"od9>F*Fgw7$$!3qC#4\\])=i>F*$\"3_lON\\Hj1;F*7$7$$!3Jb` ;$p()=+#F*F\\vFgfo7$F]goFaao7$7$FeuFheo7$$!3'GuH$4hMY=F*$!3Ig-n!*Ql`HF *7$7$$!3<\"R:B4?$R=F*FdpFcgo7$Figo7$$!3!H:/@G#f/=F*$!3Ws!=,%*Hw\"GF*7$ 7$$!3%)o')zjOO%y\"F*F`qF]ho7$Fcho7$$!3)4Zi&\\e^j)[fG4o#F*7$7$ $!3aN'3q-N?t\"F*FhrFgho7$7$$!3xN'3q-N?t\"F*Fhr7$$!3g5X^)[$>B>nj?[jU\"F*FjvF_`p7$7$$!3%H+#zB%>jU\"F*Fjv7$$!354 Hss/q-9F*$!3/<$*\\\\<_>;F*7$7$$!3A\"))[?!>Bh8F*F`wF\\ap7$7$$!3+\"))[?! >Bh8F*F[^p7$$!3_r=Xqt%\\N\"F*$!3!od#*R2(\\*[\"F*7$7$F\\x$!3t+rt/\"*\\C 9F*Fiap7$7$F\\x$\"3[z]e591xZFiy7$$!3!zVB&)fRET\"F*Fbam7$FfbpF`^p7$F_bp 7$$!3\")3L$R2FWI\"F*$!3qhLt#fRAO\"F*7$7$$!3#p!>K&HU$*G\"F*F\\xF[cp7$Fa cp7$$!3CyQ!3Ej;D\"F*$!3Y9]3GcAP7F*7$7$$!3k*R&H(>#*\\@\"F*$!3Adbbbbbb6F *Fecp7$7$F\\dpFay7$$!3BDb\\6Bz+7F*$!3m*e:'*z=.6\"F*7$7$Fay$!3My;AO#[/) )*FiyFbdp7$7$Fay$\"3U\"F*Fijl7$7$$!3ih`vc#R 2>\"F*Fijl7$$!3WcJPJj7d7F*$\"3KRE%[UuBo$Fiy7$7$$!3E5I*p^D_I\"F*F]]mFfe p7$F\\fpFcbp7$Fhdp7$$!3!f$R*)*za>:\"F*$!3-9SRM`y8)*Fiy7$7$$!3i8k/lZq\\ 6F*FgyFafp7$Fgfp7$$!3c**p[%*Q(*)4\"F*$!3'ff&o5m\"ec)Fiy7$7$$!33'[e;a;1 2\"F*F_[lF[gp7$Fagp7$$!3)>>:(3d9]5F*$!3#f*ei!p?jF(Fiy7$7$$!3!zPXuao>, \"F*$!3?TAAAAAAiFiyFegp7$7$$!3ox`WZ&o>,\"F*Fe[l7$$!39Dp$y&\\-35F*$!3\" RyI;U](>fFiy7$7$$!3ea>RsR,y(*FiyF`]lFdhp7$Fjhp7$$!3o_co\"QC!y(*Fiy$!30 0m`Sy>WWFiy7$7$Fgy$!3tW!3M5rWV%FiyF^ip7$7$Fgy$!3Xl*G7#e%3S#Fiy7$$!3Z$* p!oM0_+\"F*Fh^l7$F[jp7$$!3tw)eVZQ\"H5F*$!3QS$z2B>Gv$Fj^l7$7$$!33NFDw?] #3\"F*Ff_lF_jp7$Fejp7$Fay$\"3qsL\"Fiy$\"3j7J(eAO$=AFiy7$7$$\"3SxGV6RP=; FiyFijlF\\_q7$7$$\"3oxGV6RP=;FiyFijl7$$\"3[5*RTAb?;#Fiy$\"3#yQ$>4\"y7< $Fiy7$7$Fijl$\"39()G(HY'yzOFiyFi_q7$7$F]]m$!3\\Lvn8&ep$oFiy7$$\"3yAi$G ug%*R$FiyFe[l7$Ff`qFa^q7$7$Fijl$\"3e')G(HY'yzOFiy7$$\"3oYN$[_3*\\IFiy$ \"3GIvF'e-71%Fiy7$7$$\"3as5Hl`?4MFiyF]]mF^aq7$Fdaq7$$\"3*\\a[J%=f/SFiy $\"3?6.uXqH%)[Fiy7$7$F]]m$\"3h!42'RfoG_FiyFhaq7$7$F]]m$!3gMvn8&ep$oFiy 7$$\"3MZymHT#*\\fFiy$!3xecW2>qFxFiy7$7$Fbam$!3u:4*)Rb1\"*yFiyFebq7$F^b q7$$\"3M1'*pwWw4]Fiy$\"3wDq'**=-pl&Fiy7$7$$\"3poudMpC?dFiyFbamF_cq7$Fe cq7$$\"3\"o=90#4CLgFiy$\"3_C-$R_.7T'Fiy7$7$$\"3c/AAAAAAiFiy$\"3Ha2#3v \"RIlFiyFicq7$7$Fgbm$!33_=C8i3J()Fiy7$$\"3-7a].+#fS'FiyF_[l7$FhdqF[cq7 $7$FbamFbdq7$$\"3&zREN_b*QrFiy$\"3Q!z&p)pmK3(Fiy7$7$$\"3)[)*********** *zFiy$\"33\\P3Z]@AwFiyF^eq7$Fedq7$$\"3u^>G%QT5L*Fiy$!3:')>G%QT5L*Fiy7$ 7$Fi^m$!3dZ'zq8q]`*FiyFjeq7$Fdeq7$$\"3.C:P2ZjO#)Fiy$\"3WU%GEHlLw(Fiy7$ 7$$\"3=0?*Q^N4i)FiyFabmFdfq7$Fjfq7$$\"3y#*RT#y!G&R*Fiy$\"3#Hvj`*p\\#Q) Fiy7$7$Fi^m$\"3'**=*Q38-x&)FiyF^gq7$7$F_z$!31)QVl;)*G-\"F*7$$\"3GG`sn& o@.\"F*Fgy7$F[hq7$$\"3!3wxxxxxx*Fiy$!3YY'zq8q]`*Fiy7$7$F`hqFegq7$$\"3W s(3Q(Rte5F*$\"3]*znu\"e@o*)Fiy7$7$F_z$\"3k-:Jz?#zW*FiyFfhq7$Fhgq7$$\"3 &*omu\"*QxD7F*$!3q%*)oR6'*z/\"F*7$7$Fcx$!3[$=d4Tpo3\"F*F`iq7$F\\iq7$$ \"3g\"Q(pm7(z<\"F*$\"3#e[fjm?Ob*Fiy7$7$$\"3R,&R4'*)\\C7F*Fi^mFjiq7$7$F ajqFc^m7$$\"3PSd_7B:-8F*$\"3Tn`e)ze*35F*7$7$Fcx$\"37kumRSt@5F*Fejq7$Ff iq7$$\"3qn&\\&))o%>]\"F*$!3'e,%*HL\"RY6F*7$7$Fgw$!3'R*yo)*pl\\6F*F_[r7 $7$Fcx$\"3MkumRSt@5F*7$$\"3tnXjB*y&H9F*$\"3'yJa_'*4$f5F*7$7$$\"3I46666 66:F*$\"3Izug*\\(p\"4\"F*F\\\\r7$7$F\\v$!3W(Qe$oqA.7F*7$$\"3m$[.H,V!G: F*F^dp7$F[]rFe[r7$7$FgwFe\\r7$$\"3qYA+7EQd:F*$\"3r;WmaSG46F*7$7$$\"3q5 bLy2Lx;F*F_zFa]r7$Fg]r7$$\"3Gk8\"=\"z,'o\"F*$\"3;xIjKlUe6F*7$7$F\\v$\" 3#fzZ\"*)))Qf6F*F[^r7$Fh\\r7$$\"31Qyus&3hv\"F*$!3;3XTR_xA7F*7$7$F]u$!3 W8w9_t;b7F*Fe^r7$7$$\"3a())))))))))))o\"F*Fb^r7$$\"3X@9M)>!\\?=F*$\"3e (z!)Q-KF*$\"3t?= ]OrjW7F*7$7$Fat$\"3st>gTm\"HF\"F*Ff`r7$7$Fes$!3Y^t3CP?a8F*7$$\"3?qWM<% y%Q@F*F\\x7$FcarFb`r7$F\\ar7$$\"3Oy_]XKk!4#F*$\"33'\\sA`MrG\"F*7$7$Fes $\"3f]D'[0*)zK\"F*Fhar7$7$Fes$!3o^t3CP?a8F*7$$\"3E$y8NvB*\\AF*$!3)z*[i k[.h8F*7$7$Fgq$!3%>&fxm,A)R\"F*Febr7$F^br7$$\"3Jw9M%\\DjA#F*$\"351[qO\" F*7$7$Fgq$\"3_Dr)*)z*4w8F*Ficr7$F[cr7$$\"3w*=n.XFg[#F*$!3KF0q$yg$>9F*7 $7$F[q$!3tx#Q$G#)*>W\"F*Fcdr7$7$$\"3c)************R#F*$\"3wDr)*)z*4w8F *7$$\"3wJ6pCB%p]#F*$\"36w*>kyoTS\"F*7$7$F[q$\"3))yX\"4AwGU\"F*Fber7$Fi dr7$$\"36ujNa&\\\"*4^Ix9F*7$7$Fio$!37)=7+TBc[\"F*F\\fr7$Fh er7$$\"33o6Obl%yk#F*$\"3QpQc:F:F*7$$\"3gwD#pUn&fGF*$!3H8666666:F*7$7$FdgrF`wFbfr7$F\\gr7$ $\"3skV23;(*)y#F*$\"3M)H#fe]px9F*7$7$$\"3V\\Nm62_=HF*Fc\\rF[hr7$7$Fbhr Fgw7$$\"3(f#oa\"F*7$7$FT$\"3 ?v)G!)*Q(\\b\"F*F]jr7$Ffir7$$\"3sxaRc*\\5=$F*$!3_#[&Rc*\\5e\"F*7$7$FE$ !3'HPBE=cXg\"F*Fgjr7$Fcjr7$$\"38x?rwLK?KF*$\"3q=zGBmnz:F*7$7$FE$\"3))e `a!)HB&f\"F*Fa[s7$F][s7$$\"3X_TG?+S3MF*$!3mzj]UAiI;F*7$7$Fbdn$!3#*4%zN ?_Lk\"F*F[\\s7$Fg[s7$$\"3%*Q%oB]?`O$F*$\"3[M$4aFdCh\"F*7$7$Fbdn$\"3DGJ J%*\\NN;F*Fe\\s7$Fa\\s7$$\"3%)G[>!)f%fj$F*$!3Yy#RYU!R!o\"F*7$7$F]en$!3 7#p1o__Ao\"F*F_]s7$7$F0F\\]s7$$\"31MC67DP5NF*$\"3];JWVI=X;F*7$7$F]en$ \"3Q,.PYmRv;F*Fj]s7$7$$\"3')=AAAAAAQF*$!3c$HeYPM!=w]KnYn$F *Fjv7$7$$\"3a#>w]KnYn$F*Fjv7$F[[oFf]s7$F`^s7$$\"3M$yh)z:YbOF*$\"3#[arM vryn\"F*7$7$$\"3/Z@8\">7Wq$F*F\\vFb_s7$7$$\"3fY@8\">7Wq$F*F\\v7$$\"33& 4'eg.Q-QF*$\"3*3,D0vI(3vceQF*$!3+!4\"z&=M_s\"F*7$7$F_gn$!3Ms'fhmzKv\"F*Fi`s7$7$Fien$\"3?_ !zz_cFr\"F*7$$\"3waTlE^@]RF*$\"3\")GZBiPnQ " 0 "" {MPLTEXT 1 0 15 "J:=matrix (2,2):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 67 "for i to 2 do for j to 2 do J[i,j]:=diff(F[i],x[j]) end do end do:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "with(linalg):" }}{PARA 7 "" 1 "" {TEXT -1 80 "Warning, the protected names norm and trace have been redefined an d unprotected\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 11 "x:=[1.,1 .]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "for i to 7 do x:=eva lm(x-linsolve(J,F)) end do;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG- %'vectorG6#7$$\"+huc56!\"*$\"+U=d)H*!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+>\"f-6\"!\"*$\"+>5ai#*!#5" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+R$f-6\"!\"*$\"+&QvCE*!#5 " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+R$f-6\"!\" *$\"+(QvCE*!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$ $\"+S$f-6\"!\"*$\"+!RvCE*!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG -%'vectorG6#7$$\"+R$f-6\"!\"*$\"+)QvCE*!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+R$f-6\"!\"*$\"+'QvCE*!#5" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "x:=[1.2,1.2]:" }}{PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 47 "for i to 7 do x:=evalm(x-linsolve(J,F)) end do ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+]Bf%e\"! \"*$\"+^[%48\"F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6# 7$$\"+Wvjd9!\"*$\"+*4kh2\"F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG -%'vectorG6#7$$\"+%fEeS\"!\"*$\"+^SGa5F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+A)ybR\"!\"*$\"+\"Q/*\\5F+" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+Rp;&R\"!\"*$\"+=vs\\5F+ " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$\"+vi;&R\"!\" *$\"+Lss\\5F+" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"xG-%'vectorG6#7$$ \"+vi;&R\"!\"*$\"+Kss\\5F+" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q3: " }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "resta rt:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "restart:w:=1/sqrt(1- X^2): A:=-1: B:=1:" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 171 "for i from 0 to 3 do psi[i]:=X^i-add(int(psi[i-k]*X^i*w,X=A..B)/alpha[i-k]*psi[i -k],k=1..i):\npsi[i]:=expand(psi[i]):alpha[i]:=int(w*psi[i]^2,X=A..B): print(psi[i]) end do:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#%\"XG" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,&*$)%\"XG\"\"#\"\"\"F(#F(F'!\"\"" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#,&*$)%\"XG\"\"$\"\"\"F(*&#F'\"\"%F(F&F(!\"\"" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 29 "f:=eval(tan(x),x=Pi/8*(1+X));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fG-%$tanG6#,$*(\"\")!\"\"%#PiG\"\"\",&%\"XGF-F -F-F-F-" }}}{EXCHG {PARA 13 "" 1 "" {TEXT -1 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 43 "int(f*psi[0]/sqrt(1-X^2),X=-1..1)/alpha[0];" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#,$-%$intG6$*&-%$tanG6#,$*(\"\")!\"\" %#PiG\"\"\",&%\"XGF0F0F0F0F0F0,&F0F0*$)F2\"\"#F0F.#F.F6/F2;F.F0*$F/F. " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "a:=evalf(%);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"aG$\"+.<_cX!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "b:=evalf(int(f*psi[1]/sqrt(1-X^2),X=-1..1)/alpha [1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bG$\"+Oko%*[!#5" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 54 "c:=evalf(int(f*psi[2]/sqrt(1 -X^2),X=-1..1.)/alpha[2]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cG$ \"+z(p'p&)!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 53 "d:=evalf(i nt(f*psi[3]/sqrt(1-X^2),X=-1..1)/alpha[3]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"dG$\"+Ttt(4%!#6" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 75 "res:=evalf(expand(eval(a*psi[0]+b*psi[1]+c*psi[2]+d*p si[3],X=8./Pi*x-1.)));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$resG,*$\" *snC@\"!#6!\"\"*&$\"+B]vW5!\"*\"\"\"%\"xGF.F.*&$\"+F " 0 "" {MPLTEXT 1 0 27 "plot(tan(x)-res,x=0..Pi/4);" }}{PARA 13 "" 1 "" {GLPLOT2D 494 494 494 {PLOTDATA 2 "6%-%'CURVESG6$7cp7$$\"\"!F)$\"3++++?xY77!#?7$$\"3 )H$**)Q?h)zUF,$\"3YIYu/$=`-\"F,7$$\"3%f')zxSA(f&)F,$\"3Y:Dz*zH&o%)!#@7 $$\"3)*zp;h$eRG\"!#>$\"3s@mM!)G5pnF77$$\"3>tfb\"[W>r\"F;$\"3eZ&\\?<)G` ^F77$$\"3k)>)R^8scCF;$\"30uw\\b&>g`#F77$$\"3WC/C@#)\\,KF;$\"3wZY5#e)p& e\"!#A7$$\"3G0#f^qx!RSF;$!3fj\"Qt>r\"RAF77$$\"3W&)z2*=dm([F;$!3K-*yqOu jN%F77$$\"3)f()3Np)y>dF;$!341FU&3Mk@'F77$$\"3%ewRz>?Hc'F;$!3BK^VkOg;yF 77$$\"3^Fc]XW/-uF;$!3c([AG]&=j\"*F77$$\"3>*[rIpo6C)F;$!3q#f*[+GjF5F,7$ $\"3ah&>)y)>rz*F;$!3+qyqotHw6F,7$$\"3%Hc8V&oEg5!#=$!3KNEAkI\"eA\"F,7$$ \"3sp^u?<#39\"F^p$!3)eRS%4f&zD\"F,7$$\"3AF^p$!3kwG)Qw6![#)F77$$\"31V([!yqLM@F^p$!3%yXk+2'\\ $='F77$$\"3k>'GySb#)H#F^p$!3\"*zc[c!H0+%F77$$\"3c$QD\\&*3rW#F^p$!3f\"R #Q$\\V^\">F77$$\"3Tcn%)>-6CEF^p$\"3Cw-TG%4N='FM7$$\"3m,$H&3N0uFF^p$\"3 !py9oQ]@v#F77$$\"39!=X\"*>t%[HF^p$\"3E)**pDv4K:&F77$$\"3%RW>'zT'G5$F^p $\"3MK6pz[J_rF77$$\"3ulP&[bbAF$F^p$\"33N;[R3c`\"*F77$$\"3%ylZ;?cNV$F^p $\"3c\"=!Q25V#3\"F,7$$\"3ah!G]]c=g$F^p$\"3$QdZ5\\@(3 kv$F^p$\"3E*GW)4))\\J8F,7$$\"3sY-TV9wRQF^p$\"3a-!RhlwYP\"F,7$$\"3RI0il T6BRF^p$\"3%>$3@1wS39F,7$$\"3>379=Up4SF^p$\"38X\"3*o\")4L9F,7$$\"3)f)= mqUF'4%F^p$\"3Q`mV$=WpW\"F,7$$\"3'*=))*3YU;<%F^p$\"3!*=(fJ(G(*\\9F,7$$ \"3%>vN6l5qC%F^p$\"3$[dZA$oYW9F,7$$\"3e#oq8O5%GVF^p$\"3a\"RP$Q#G)G9F,7 $$\"3n7cgr+\")4WF^p$\"34bm7`'fIS\"F,7$$\"3[tcCPR!R\\%F^p$\"3s1yqmMyl8F ,7$$\"3yLd)G!y*zd%F^p$\"3KGh$)[Cs<8F,7$$\"3DY2d&oODu%F^p$\"3/1&o[2iH> \"F,7$$\"3,9-+:lt,\\F^p$\"3UceW5/'[.\"F,7$$\"3;jPv#y+&y]F^p$\"3)HbGL*) H:>)F77$$\"3%)Q(Q#R8LP_F^p$\"3=s-Bz+7IfF77$$\"3WE**)*fo\"pS&F^p$\"3)[ \\CKfu4B$F77$$\"3#*H[:ZlegbF^p$\"3)o8TMEmk)eFM7$$\"3#*G%4Z5)eGdF^p$!3e QGJkf&eV#F77$$\"3is2qxWm')eF^p$!3\"43'fx;&QK&F77$$\"3[KK6N)**=0'F^p$!3 #3.n69WpF)F77$$\"366vD*fZM@'F^p$!3ULtSrQX(4\"F,7$$\"334!\\o8tDQ'F^p$!3 CtX(>%\\9W8F,7$$\"3KplfI6YXlF^p$!3*>iv`n\"*)H:F,7$$\"3\\?rh!Q\\(GmF^p$ !35=UYp1,*f\"F,7$$\"3nrwjIw.7nF^p$!3m`rn/W'pk\"F,7$$\"3w'>n7#pL`nF^p$! 3mB.im>(>m\"F,7$$\"3%Hs'*=@OYz'F^p$!3%elB(fRlq;F,7$$\"30'[6s&eG:oF^p$! 3unuyZ%*\\s;F,7$$\"3:\\i_-b$f$oF^p$!3!yR*H8Rhs;F,7$$\"3C75%y9&ecoF^p$! 33cm'R(e%4n\"F,7$$\"3Cud:$zMs(oF^p$!3=0HxgIWn;F,7$$\"3V4f1_X=:pF^p$!31 F%fXhVgl\"F,7$$\"3iWg(4JMJ&pF^p$!3`-GUz4$zj\"F,7$$\"3!)zh))pS3\"*pF^p$ !3W*y+DDfFh\"F,7$$\"3+:jzGQ.HqF^p$!31!=O4$Ho&*Q@'etF^p$!3EU76[iXw#*F77$$\"3cg/ZJ#*4+uF^p$!3G)zQ ]rWd!zF77$$\"3nAFPtqdTuF^p$!3Mfk-8Fn#R'F77$$\"3!f)\\F:\\0$[(F^p$!3I)Rv Bp$)3t%F77$$\"37\\s$Rg(F^p$\"3O7=JRQ:A5F77$$\"36n(Q'=:jVwF^p$\"3h, ()=<<1GKF77$$\"3!og#z06L$o(F^p$\"3GG:RHK0,cF77$$\"3ab>\"Rt$*fs(F^p$\"3 ?>HM>]#\\M)F77$$\"3;.8.ijloxF^p$\"3eeSKULxH6F,7$$\"3z]1:!**=8\"yF^p$\" 3KrV3D5\"oW\"F,7$$\"3_***p#=;)R&yF^p$\"3G)=:c$Q\\'y\"F,-%'COLOURG6&%$R GBG$\"#5!\"\"F(F(-%+AXESLABELSG6$Q\"x6\"Q!Fdgl-%%VIEWG6$;F($\"+N;)R&y! #5%(DEFAULTG" 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 36 "evalf(tan(Pi/4 .)-eval(res,x=Pi/4.));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\")S\\'y\"! #5" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q2:" }{MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 149 "LUdec:=proc(a,b,c) local L,U,k,n; n:=nops(b ); L:=[b[1]]; U:=[];\nfor k to n-1 do U:=[op(U),c[k]/L[k]];L:=[op(L),b [k+1]-a[k]*U[k]] end do;\n[a,L,U] end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 224 "tridset:=proc(M,r) local n,k,a,L,U,y,x; a:=M[1]; L:= M[2]; U:=M[3]; n:=nops(L);\ny:=[r[1]/L[1]]; for k to n-1 do y:=[op(y), (r[k+1]-a[k]*y[k])/L[k+1]] end do;\nx:=[y[n]]; for k to n-1 do x:=[y[n -k]-U[n-k]*x[1],op(x)] end do end:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "n:=4: A:=-1.: B:=3: h:=(B-A )/n:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "p:=z->sin(z):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "q:=z->-1./(1.+z^3):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "r:=z->exp(-z*1.):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "x:=[seq(A+h*k,k=1..n-1)]:" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "a:=[seq(1-h/2*p(x[k]),k=2. .n-1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "b:=[seq(-2+h^2*q (x[k]),k=1..n-1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "c:=[s eq(1+h/2*p(x[k]),k=1..n-2)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "R:=[seq(h^2*r(x[k]),k=1..n-1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "R[1]:=R[1]-0.5*(1-h/2*p(x[1])):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "R[n-1]:=R[n-1]-4*(1+h/2*p(x[n-1])):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 28 "y4:=tridset(LUdec(a,b,c),R); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#y4G7%$\"+4?L2T!#5$\"+.'*>K " 0 "" {MPLTEXT 1 0 19 "n:=12: h:= (B-A)/n:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "x:=[seq(A+h*k,k =1..n-1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "a:=[seq(1-h/2 *p(x[k]),k=2..n-1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 34 "b:= [seq(-2+h^2*q(x[k]),k=1..n-1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 33 "c:=[seq(1+h/2*p(x[k]),k=1..n-2)]:" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 31 "R:=[seq(h^2*r(x[k]),k=1..n-1)]:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "R[1]:=R[1]-0.5*(1-h/2*p(x[1])):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 35 "R[n-1]:=R[n-1]-4*(1+h/2*p(x[n-1])): " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "y12:=tridset(LUdec(a,b, c),R);" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$y12G7-$\"+&*3AmD!#5$\"+o' 4xV#F($\"+h;*>B%F($\"+fgg2wF($\"+)H-i@\"!\"*$\"+k.wKPF1" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 12 "with(plots):" }}{PARA 7 "" 1 "" {TEXT -1 50 "W arning, the name changecoords has been redefined\n" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 55 "pointplot([[A,0.5],seq([x[k],y12[k]],k=1..n -1),[B,4]]);" }}{PARA 13 "" 1 "" {GLPLOT2D 494 494 494 {PLOTDATA 2 "6# -%'POINTSG6/7$$!\"\"\"\"!$\"\"&F(7$$!+nmmmm!#5$\"+&*3AmDF/7$$!+MLLLLF/ $\"+o'4xV#F/7$$F(F/$\"+h;*>B%F/7$$\"*LLLL$!\"*$\"+fgg2wF/7$$\"*mmmm'F> $\"+)H-i@\"F>7$$\"+++++5F>$\"+k.wK7$$\"+LLLL8F>$\"+7?_UAF>7$$\"+mmm m;F>$\"+*R$)*)p#F>7$$\"+++++?F>$\"+Q.B*3$F>7$$\"+LLLLBF>$\"+X\"QHU$F>7 $$\"+mmmmEF>$\"+PF>7$$\"\"$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 19 "(y12[3]*9-y4[1])/8;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+=TdZU!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "(y12[ 6]*9-y4[2])/8;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+g/$Gt\"!\"*" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 19 "(y12[9]*9-y4[3])/8;" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+&oTH3$!\"*" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 3 "Q5:" }{MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 8 "restart:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 50 "f:=x->exp(-x)/(2.+cos(x)): A:=-1: B:=2.: n:=2;" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%\"nG\"\"#" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 150 "for k from 1 to 3 do In[k]:=evalf( (B-A)*(\nf(A)+f(B)+\n2*sum( f(A+2*(B-A)*i/n),i=1..n/2-1)+\n4*sum(f(A+(2*i-1)*(B-A)/n),i=1..n/2)\n) /3/n); n:=n*3 end do;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%#InG6#\"\" \"$\"+lH5$***!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"\"'" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>&%#InG6#\"\"#$\"+X)>HY*!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"nG\"#=" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>&%#InG6#\"\"$$\"+289^%*!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"n G\"#a" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 57 "for i from 1 to 2 \+ do J[i]:=(81*In[i+1]-In[i])/80 end do; " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"JG6#\"\"\"$\"+cDHc%*!#5" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >&%\"JG6#\"\"#$\"+vS*4X*!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 56 "for i from 1 to 1 do K[i]:=(729*J[i+1]-J[i])/728 end do;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"KG6#\"\"\"$\"+(z')4X*!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 80 "evalf(int(exp(-x)/(2+cos(x)),x=-1.. 2)); # 'exact' value, just for comparison" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+p3)4X*!#5" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "73 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 } {PAGENUMBERS 0 1 2 33 1 1 }