{VERSION 4 0 "SUN SPARC SOLARIS" "4.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 "" 0 1 0 0 0 0 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 "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 " " 0 256 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 "" 0 257 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 "" 0 258 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 }} {SECT 0 {EXCHG {PARA 256 "" 0 "" {TEXT -1 11 "Math 2270-1" }}{PARA 257 "" 0 "" {TEXT -1 20 "Wednesday October 25" }}{PARA 258 "" 0 "" {TEXT -1 22 "Fourier Series example" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 25 "with(linalg):with(plots):" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 45 "Dot:=(f,g)->(1/Pi) *int(f(x)*g(x),x=-Pi..Pi);" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%$DotGR6$%\"fG%\"gG6\"6$%)operatorG%& arrowGF)*&-%$intG6$*&-9$6#%\"xG\"\"\"-9%F4F6/F5;,$%#PiG!\"\"F " 0 "" {MPLTEXT 1 0 242 "Mag:=f->sqrt(Dot(f ,f)); #computes the magnitude of\n #a vector\nDist:=(f,g)->Mag(f-g ); #computes the ``distance'' between\n #two vectors. \nangle:=( f,g)->arccos(Dot(f,g)/(Mag(f)*Mag(g)));\n #computes the ``angle'' \+ between two vectors" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$MagGR6#%\"fG 6\"6$%)operatorG%&arrowGF(-%%sqrtG6#-%$DotG6$9$F2F(F(F(" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%%DistGR6$%\"fG%\"gG6\"6$%)operatorG%&arrowGF)- %$MagG6#,&9$\"\"\"9%!\"\"F)F)F)" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%& angleGR6$%\"fG%\"gG6\"6$%)operatorG%&arrowGF)-%'arccosG6#*&-%$DotG6$9$ 9%\"\"\"*&-%$MagG6#F4F6-F96#F5F6!\"\"F)F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "Dot(x->sin(3*x),x->sin(4*x));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#\"\"!" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 211 "f: =x->abs(x); #we'll do Fourier for the absolute value function\nn:=4; \+ #order 10 expansion\na:=vector(n); #sin coefficients\nb:=vector(n); \+ #cos coefficients\nc:=Dot(f,1/sqrt(2)); #order zero Fourier coeffici ent" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"fG%$absG" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%\"nG\"\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"aG -%&arrayG6$;\"\"\"\"\"%7\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"bG-% &arrayG6$;\"\"\"\"\"%7\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"cG,$*& %#PiG\"\"\"-%%sqrtG6#\"\"#F(#F(F," }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 115 "for i from 1 to n do #compute the projection coeffi cients\na[i]:=Dot(f,x->sin(i*x));\nb[i]:=Dot(f,x->cos(i*x));\nod:\n" } }{PARA 0 "> " 0 "" {MPLTEXT 1 0 88 "evalm(a);#why will an even functio n have Fourier sine\n #coefficients all equal to zero?" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#-%'vectorG6#7&\"\"!F'F'F'" }}}{EXCHG {PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 10 "evalm(b); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6# -%'vectorG6#7&,$*&\"\"\"F)%#PiG!\"\"!\"%\"\"!,$F(#F,\"\"*F-" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 97 "approx:=x->sum(a[k]*sin(k*x) ,k=1..n)\n +sum(b[m]*cos(m*x),m=1..n)\n +c/sqrt(2) ;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'approxGR6#%\"xG6\"6$%)operator G%&arrowGF(,(-%$sumG6$*&&%\"aG6#%\"kG\"\"\"-%$sinG6#*&F4F59$F5F5/F4;F5 %\"nGF5-F.6$*&&%\"bG6#%\"mGF5-%$cosG6#*&FDF5F:F5F5/FDF " 0 "" {MPLTEXT 1 0 49 "approx(x); #order 4 Fourier expansion for abs(x):" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#,(*&-%$cosG6#%\"xG\"\"\"%#PiG!\"\"!\"%*&#\"\"%\"\"*F) *&-F&6#,$F(\"\"$F)F*F+F)F+*&#F)\"\"#F)F*F)F)" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 118 "almost:=plot(approx(x),x=-Pi..Pi,color=`black`) :\nexact:=plot(abs(x),x=-Pi..Pi,color=`black`):\ndisplay(\{almost,exac t\});" }}{PARA 13 "" 1 "" {GLPLOT2D 400 300 300 {PLOTDATA 2 "6&-%'CURV ESG6$7]o7$$!3)****4tk#fTJ!#<$\"3n)=c?$p]&)HF*7$$!3E[2g9v#F*$\"3$y\" o&fG1N!GF*7$$!3)4577.flh#F*$\"3A%4[]cI>n#F*7$$!3u\"QM::*H#[#F*$\"3'e_g (>Q=@DF*7$$!3x\\AhcI#yN#F*$\"3!4Ux1YKIP#F*7$$!3&e=d-FN*GAF*$\"3(=*3g\" \\k%>AF*7$$!3u%*GBP$Rc4#F*$\"3G:HI'\\F*$\"3^92Btt gE>F*7$$!3'or4AN*4E=F*$\"3%oU*G5=F%z\"F*7$$!3QQQ*3**=dq\"F*$\"3,Q]6yBN 'o\"F*7$$!3d!>jg@*>q:F*$\"3_XPP$[*Gq:F*7$$!3_;m,%)H7M9F*$\"3%z\\cM'Rr` 9F*7$$!3v(4F,K))HI\"F*$\"3nxl1jfpN8F*7$$!3R*4!R#[0R=\"F*$\"30\"f,pCG,A \"F*7$$!3oY@QqWIU5F*$\"39='3?iW-2\"F*7$$!3Wp3w$R)\\B#*!#=$\"3')H2b1?)[ L*F^r7$$!3EJc8o39GyF^r$\"3r3$**\\&=8uwF^r7$$!3'=[BP-8If'F^r$\"3K)o7fdz T?'F^r7$$!3CB3<@?)yB&F^r$\"3P<03NaR$o%F^r7$$!3_*)RF^r$\"3pxh+- /N`?F^r7$$!3pTJY)ocYO\"F^r$\"3kn/M<09'z\"F^r7$$!3-w^X*zX7.\"F^r$\"3'4K %y**Qm&p\"F^r7$$!3b.@Z/\"\\$yp!#>$\"3gZ<$R(\\tA;F^r7$$!3%oW#R9-CWOF\\u $\"3XIK$G5dxd\"F^r7$$!38,z7VKJ,J!#?$\"3N6g&onz4c\"F^r7$$\"3zdsCy22`JF \\u$\"3l!enkL5Nd\"F^r7$$\"3p0t!3)GF;mF\\u$\"3&>H#\\A=\\;;F^r7$$\"3ONnL )\\Zz+\"F^r$\"3g'\\=-'pm*o\"F^r7$$\"39SFf3xEa8F^r$\"3wR)R&*)*)f#z\"F^r 7$$\"302**yIK@d>F^r$\"3=s*>$\\6*3/#F^r7$$\"3'R2()Hve,c#F^r$\"3+h1l,R+t BF^r7$$\"3IcwRz^ YF^r7$$\"3q]1XIqOClF^r$\"39&Q'RRs+ChF^r7$$\"3%HoBlmlzz(F^r$\"3'**oP2%) 4zj(F^r7$$\"3^u'f#4)z?@*F^r$\"3D2:^Q(=:K*F^r7$$\"38HvkBm C$o2\"F*7$$\"3Nlj%G%3%R=\"F*$\"3]!4Ts&R;?7F*7$$\"3\\-)REfwoI\"F*$\"3Qt +Z`,KR8F*7$$\"3cM2vQyFT9F*$\"3^BV!Q_X*f9F*7$$\"3z(32s$*Qxc\"F*$\"3>CZ'eYk>F*$\"3]()4xigKG>F*7$$\"3kYuyfix%4#F*$\"3&p$ 3/r`Om?F*7$$\"3S4q))fu.GAF*$\"30]K`I1T=AF*7$$\"3y?w'**=&>gBF*$\"3&y'pC &)*zeP#F*7$$\"3!*>3a=Wj\"[#F*$\"3nO\"=]>1/_#F*7$$\"3A?c*)yu\"3i#F*$\"3 o]u$*QcWwEF*7$$\"3,i-HnWIXFF*$\"3iFlsM\")3)z#F*7$$\"3u=RWhN.yGF*$\"3gP Sd\\^d**GF*7$$\"3tX=#4!HbTHF*$\"3xf46rtRNHF*7$$\"3ss(*RSA20IF*$\"3([Bk r')f>'HF*7$$\"3;HB\"HM-#RIF*$\"3R1+D;n@sHF*7$$\"3.')[UXCLtIF*$\"3>vA`J eezHF*7$$\"3#HWPzaiu5$F*$\"3_ayj+X-%)HF*7$$\"3!)***\\/l#fTJF*F+-%'COLO URG6&%$RGBG\"\"!F\\`lF\\`l-F$6$7U7$F($\"3)****4tk#fTJF*7$F=$\"3w\"*pr) 3PY+$F*7$FG$\"3?3E[*ysa)GF*7$FL$\"3fX&))>2g9v#F*7$FQ$\"3)4577.flh#F*7$ FV$\"3u\"QM::*H#[#F*7$Fen$\"3x\\AhcI#yN#F*7$Fjn$\"3&e=d-FN*GAF*7$F_o$ \"3u%*GBP$Rc4#F*7$Fdo$\"3lCj&f)3xi>F*7$Fio$\"3'or4AN*4E=F*7$F^p$\"3QQQ *3**=dq\"F*7$Fcp$\"3d!>jg@*>q:F*7$Fhp$\"3_;m,%)H7M9F*7$F]q$\"3v(4F,K)) HI\"F*7$Fbq$\"3R*4!R#[0R=\"F*7$Fgq$\"3oY@QqWIU5F*7$F\\r$\"3Wp3w$R)\\B# *F^r7$Fbr$\"3EJc8o39GyF^r7$Fgr$\"3'=[BP-8If'F^r7$F\\s$\"3CB3<@?)yB&F^r 7$Fas$\"3_*)R " 0 "" {MPLTEXT 1 0 0 "" }}}} {MARK "11 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }