{VERSION 6 0 "IBM INTEL LINUX" "6.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 2 1 2 0 0 0 1 }{CSTYLE "2D Output" -1 20 "Times" 0 1 0 0 255 1 0 0 2 2 2 2 0 0 0 1 }{CSTYLE "_cstyle1" -1 205 "Times" 1 18 0 0 0 0 0 1 0 2 2 2 0 0 0 1 }{CSTYLE "_cstyle2" -1 206 "Courier" 0 1 255 0 0 1 0 1 0 2 1 2 0 0 0 1 }{CSTYLE "_cstyle5" -1 209 "Times" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Error" -1 8 1 {CSTYLE "" -1 -1 "Courier" 1 10 255 0 255 1 0 0 0 2 2 1 0 0 0 1 }1 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Maple Output" -1 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }3 3 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "Left Justified Maple Output" -1 12 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 0 0 0 1 }1 3 0 0 0 0 2 0 2 0 2 2 -1 1 }{PSTYLE "Maple Plot" -1 13 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }3 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "_ps tyle1" -1 201 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 } 0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }{PSTYLE "_pstyle5" -1 205 1 {CSTYLE " " -1 -1 "" 0 1 0 0 0 0 0 0 0 2 2 2 0 0 0 1 }0 0 0 -1 -1 -1 1 0 1 0 2 2 -1 1 }} {SECT 0 {PARA 201 "" 0 "" {TEXT 205 54 "Transverse modes of a free bea m 4/14/2004 for MTH 444" }}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 8 "restart:" }{MPLTEXT 1 206 69 "\nu:=c1*exp(alpha*x)+c2*exp(-alpha*x) +c3*cos(alpha*x)+c4*sin(alpha*x);" }{MPLTEXT 1 206 15 "\nu1:=diff(u,x) ;" }{MPLTEXT 1 206 16 "\nu2:=diff(u1,x);" }{MPLTEXT 1 206 16 "\nu3:=di ff(u2,x);" }{MPLTEXT 1 206 43 "\nde:=simplify(diff(u,x,x,x,x)-alpha^4* u)=0;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"uG,**&%#c1G\"\"\"-%$expG6 #*&%&alphaGF(%\"xGF(F(F(*&%#c2GF(-F*6#,$F,!\"\"F(F(*&%#c3GF(-%$cosGF+F (F(*&%#c4GF(-%$sinGF+F(F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u1G,** (%#c1G\"\"\"%&alphaGF(-%$expG6#*&F)F(%\"xGF(F(F(*(%#c2GF(F)F(-F+6#,$F- !\"\"F(F4*(%#c3GF(-%$sinGF,F(F)F(F4*(%#c4GF(-%$cosGF,F(F)F(F(" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u2G,**(%#c1G\"\"\")%&alphaG\"\"#F(- %$expG6#*&F*F(%\"xGF(F(F(*(%#c2GF(F)F(-F-6#,$F/!\"\"F(F(*(%#c3GF(-%$co sGF.F(F)F(F6*(%#c4GF(-%$sinGF.F(F)F(F6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#u3G,**(%#c1G\"\"\")%&alphaG\"\"$F(-%$expG6#*&F*F(%\"xGF(F(F(* (%#c2GF(F)F(-F-6#,$F/!\"\"F(F6*(%#c3GF(-%$sinGF.F(F)F(F(*(%#c4GF(-%$co sGF.F(F)F(F6" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#deG/\"\"!F&" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 40 "bc1:=simplify(subs(x=-HL ,u2)/alpha^2)=0;" }{MPLTEXT 1 206 41 "\nbc2:=simplify(subs(x= HL,u2)/a lpha^2)=0;" }{MPLTEXT 1 206 41 "\nbc3:=simplify(subs(x=-HL,u3)/alpha^3 )=0;" }{MPLTEXT 1 206 41 "\nbc4:=simplify(subs(x= HL,u3)/alpha^3)=0;" }{MPLTEXT 1 206 26 "\nevalm([[coeff(bc1,c1)]]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$bc1G/,**&%#c1G\"\"\"-%$expG6#,$*&%&alphaGF)%#HLGF)! \"\"F)F)*&%#c2GF)-F+6#F.F)F)*&%#c3GF)-%$cosGF,F)F1*&%#c4GF)-%$sinGF,F) F1\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$bc2G/,**&%#c1G\"\"\"-%$e xpG6#*&%&alphaGF)%#HLGF)F)F)*&%#c2GF)-F+6#,$F-!\"\"F)F)*&%#c3GF)-%$cos GF,F)F5*&%#c4GF)-%$sinGF,F)F5\"\"!" }}{PARA 11 "" 1 "" {XPPMATH 20 "6# >%$bc3G/,**&%#c1G\"\"\"-%$expG6#,$*&%&alphaGF)%#HLGF)!\"\"F)F)*&%#c2GF )-F+6#F.F)F1*&%#c3GF)-%$sinGF,F)F)*&%#c4GF)-%$cosGF,F)F1\"\"!" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%$bc4G/,**&%#c1G\"\"\"-%$expG6#*&%&al phaGF)%#HLGF)F)F)*&%#c2GF)-F+6#,$F-!\"\"F)F5*&%#c3GF)-%$sinGF,F)F)*&%# c4GF)-%$cosGF,F)F5\"\"!" }}{PARA 8 "" 1 "" {TEXT -1 31 "Error, unable \+ to compute coeff\n" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 38 "c c1:=expand(exp(-alpha*HL)*(bc3-bc1));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$cc1G/,,*&\"\"#\"\"\"%#c2GF)!\"\"*(-%$expG6#*&%&alphaGF)%#HLGF)F+ %#c3GF)-%$sinGF/F)F+*(F-F+%#c4GF)-%$cosGF/F)F+*(F-F+F3F)F8F)F)*(F-F+F7 F)F4F)F+\"\"!" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 38 "cc2:=e xpand(exp( alpha*HL)*(bc4-bc2));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% $cc2G/,,*&\"\"#\"\"\"%#c2GF)!\"\"*(-%$expG6#*&%&alphaGF)%#HLGF)F)%#c3G F)-%$sinGF/F)F)*(F-F)%#c4GF)-%$cosGF/F)F+*(F-F)F3F)F8F)F)*(F-F)F7F)F4F )F)\"\"!" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 40 "dc1:=simpli fy((cc1-cc2)*exp( alpha*HL));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$dc 1G/,2*&%#c3G\"\"\"-%$sinG6#*&%&alphaGF)%#HLGF)F)!\"\"*&%#c4GF)-%$cosGF ,F)F0*&F(F)F3F)F)*&F2F)F*F)F0*(-%$expG6#,$*(\"\"#F)F.F)F/F)F)F)F(F)F*F )F0*(F8F)F2F)F3F)F)*(F8F)F(F)F3F)F0*(F8F)F2F)F*F)F0\"\"!" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 38 "cc3:=expand(exp( alpha*HL)*(bc3 +bc1));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$cc3G/,,*&\"\"#\"\"\"%#c1 GF)F)*(-%$expG6#*&%&alphaGF)%#HLGF)F)%#c3GF)-%$sinGF.F)!\"\"*(F,F)%#c4 GF)-%$cosGF.F)F5*(F,F)F2F)F8F)F5*(F,F)F7F)F3F)F)\"\"!" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 38 "cc4:=expand(exp(-alpha*HL)*(bc4 +bc2));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$cc4G/,,*&\"\"#\"\"\"%#c1 GF)F)*(-%$expG6#*&%&alphaGF)%#HLGF)!\"\"%#c3GF)-%$sinGF.F)F)*(F,F2%#c4 GF)-%$cosGF.F)F2*(F,F2F3F)F8F)F2*(F,F2F7F)F4F)F2\"\"!" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 40 "dc2:=simplify((cc3-cc4)*exp( al pha*HL));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$dc2G/,2*(-%$expG6#,$*( \"\"#\"\"\"%&alphaGF.%#HLGF.F.F.%#c3GF.-%$sinG6#*&F/F.F0F.F.!\"\"*(F(F .%#c4GF.-%$cosGF4F.F6*(F(F.F1F.F9F.F6*(F(F.F8F.F2F.F.*&F1F.F2F.F6*&F8F .F9F.F.*&F1F.F9F.F.*&F8F.F2F.F.\"\"!" }}}{EXCHG {PARA 201 "" 0 "" {TEXT -1 0 "" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 20 "a:=arra y(1..2,1..2);" }{MPLTEXT 1 206 34 "\na[1,1]:=lhs(subs(c3=1,c4=0,dc1)); " }{MPLTEXT 1 206 34 "\na[1,2]:=lhs(subs(c3=0,c4=1,dc1));" }{MPLTEXT 1 206 34 "\na[2,1]:=lhs(subs(c3=1,c4=0,dc2));" }{MPLTEXT 1 206 34 "\na [2,2]:=lhs(subs(c3=0,c4=1,dc2));" }{MPLTEXT 1 206 10 "\nevalm(a);" } {MPLTEXT 1 206 41 "\ndeterm:=simplify(linalg[det](evalm(a)));" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"aG-%&arrayG6%;\"\"\"\"\"#F(7\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"aG6$\"\"\"F',*-%$sinG6#*&%&alphaG F'%#HLGF'!\"\"-%$cosGF+F'*&-%$expG6#,$*(\"\"#F'F-F'F.F'F'F'F)F'F/*&F3F 'F0F'F/" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"aG6$\"\"\"\"\"#,*-%$co sG6#*&%&alphaGF'%#HLGF'!\"\"-%$sinGF,F0*&-%$expG6#,$*(F(F'F.F'F/F'F'F' F*F'F'*&F4F'F1F'F0" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"aG6$\"\"#\" \"\",*-%$sinG6#*&%&alphaGF(%#HLGF(!\"\"-%$cosGF,F(*&-%$expG6#,$*(F'F(F .F(F/F(F(F(F*F(F0*&F4F(F1F(F0" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\" aG6$\"\"#F',**&-%$expG6#,$*(F'\"\"\"%&alphaGF/%#HLGF/F/F/-%$cosG6#*&F0 F/F1F/F/!\"\"*&F*F/-%$sinGF4F/F/F2F/F8F/" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#K%'matrixG6#7$7$,*-%$sinG6#*&%&alphaG\"\"\"%#HLGF.!\"\"-%$cosGF+ F.*&-%$expG6#,$*(\"\"#F.F-F.F/F.F.F.F)F.F0*&F4F.F1F.F0,*F1F0F)F0F:F.F3 F07$F(,*F:F0F3F.F1F.F)F.Q(pprint06\"" }}{PARA 11 "" 1 "" {XPPMATH 20 " 6#>%'determG,,*&\"\"%\"\"\"-%$expG6#,$*(\"\"#F(%&alphaGF(%#HLGF(F(F(! \"\"F.F1*&F'F()-%$cosG6#*&F/F(F0F(F.F(F(*&F.F(-F*6#,$*(F'F(F/F(F0F(F(F (F1*(F'F(F9F(F3F(F(" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 35 " determ2:=subs(alpha*HL=ahl,determ);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6 #>%(determ2G,,*&\"\"%\"\"\"-%$expG6#,$*&\"\"#F(%$ahlGF(F(F(!\"\"F.F0*& F'F()-%$cosG6#F/F.F(F(*&F.F(-F*6#,$*&F'F(F/F(F(F(F0*(F'F(F7F(F2F(F(" } }}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 38 "determ3:=simplify(dete rm2/exp(4*ahl));" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(determ3G,,*&\" \"%\"\"\"-%$expG6#,$*&\"\"#F(%$ahlGF(!\"\"F(F0*&F.F(-F*6#,$*&F'F(F/F(F 0F(F0*(F'F(F2F()-%$cosG6#F/F.F(F(F.F0*&F'F(F7F(F(" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 24 "plot(determ3,ahl=0..18);" }}{PARA 13 "" 1 "" {GLPLOT2D 484 152 152 {PLOTDATA 2 "6%-%'CURVESG6$7ex7$$\"\"!F) F(7$$\"3z****\\(ofV!\\!#>$!3'*4YctEZ%f&!#A7$$\"3d*****\\P>(3)*F-$!3>,M i/;#[6)!#@7$$\"32++D1zIr9!#=$!3:T/*3/_Us$!#?7$$\"3#******\\(Quh>F:$!32 L#ffnMq1\"F-7$$\"3:++]7ehUHF:$!3.Tl#)=cgQWF-7$$\"3$)******\\x[BRF:$!3w ![LqQ@A:\"F:7$$\"3w***\\ig!RIcF:$!3Ky_7vj7jMF:7$$\"3\"3++DY$HPtF:$!3Mq t>8NsZqF:7$$\"3)3+DJ&))3(H)F:$!36&GlXYnTW*F:7$$\"3'4+]PC%)oD*F:$!3#HiD %H_c&>\"!#<7$$\"35+vVjzm@5F[o$!3h,%Q#e#pRW\"F[o7$$\"36++]-vk<6F[o$!3NQ 'R*y8Gv;F[o7$$\"32+++5!z3J\"F[o$!3oC0nN?F[o7$$\"3-++]<06/:F[o$!3\"> !ey$\\*e%=#F[o7$$\"3!***\\PHOU'p\"F[o$!3M`sspD%Q2#F[o7$$\"3z***\\7uO() )=F[o$!3SzLv#)=Y,F[o$!3\"Gt&>+*3+X\"F[o7$$\"3* )***\\PhNq1#F[o$!3nfm\"4n#ed6F[o7$$\"3%*******\\]=c@F[o$!3DB3m)[$=E$)F :7$$\"3*****\\i[M`C#F[o$!3%)G/(ehNs%[F:7$$\"3$****\\PcVwL#F[o$!3\\ne)Q AX[6\"F:7$$\"3))***\\7k_*HCF[o$\"30W@)\\?l(HEF:7$$\"3$)***\\(=RL1>F[o7$ $\"3))**\\iDb,WIF[o$\"3-)G)R]h'H&>F[o7$$\"3s*****\\`&f\"4$F[o$\"3*e_&R JEw\")>F[o7$$\"3-+]PWbF[o7$$\"3&)***\\Pbbn=$F[o$ \"3j8#R3'z-&)>F[o7$$\"3******\\sb\">G$F[o$\"3&zg?1X:h\">F[o7$$\"3n*** \\7fvqP$F[o$\"3`>0'z=*fxpsqOF[o$\"3f_0?7fr\"y*F: 7$$\"3Y+++&p5'oPF[o$\"3lWFz*>z%3iF:7$$\"32+]7yn#[&QF[o$\"3:K`LF?heGF:7 $$\"39++DhG/TRF[o$!378bi1m0rdF-7$$\"3?+]PW*es-%F[o$!3I.Z/%H/m*RF:7$$\" 3$)****\\F]Z8TF[o$!3_m.qO?K)H(F:7$$\"3y**\\(o\"R`5UF[o$!3xhU/ISBv5F[o7 $$\"3u***\\i!Gf2VF[o$!3-\"))G^j7.Q\"F[o7$$\"3f+]i&p^YS%F[o$!31J-V\"['f L;F[o7$$\"3a+++&e5<]%F[o$!3(e$Hs*ohb#=F[o7$$\"3K++]2!o\"*f%F[o$!30xUlB wM\\>F[o7$$\"35+++Iai'p%F[o$!37.,>v#R$**>F[o7$$\"30+]i&y*)4s%F[o$!3.ZO Q<9-+?F[o7$$\"3*****\\79a`u%F[o$!3#*GKp!)p&f*>F[o7$$\"3%***\\(o\\=(pZF [o$!3li9\\Ub:()>F[o7$$\"3))****\\_G3%z%F[o$!3ArOGHzjt>F[o7$$\"3x***\\P c6G%[F[o$!3**owK0MfK>F[o7$$\"3m*****\\FS:*[F[o$!3K>JbbC@t=F[o7$$\"3m** ****f&ea)\\F[o$!3`x328$f#4o:qPk Q%F:7$$\"3O*\\P4([x4dF[o$\"3l@lAXTIF#)F:7$$\"3!***\\iD\"*=6eF[o$\"3F#G sLo#3t6F[o7$$\"3V+DJ!Q.E\"fF[o$\"3>0@c,tLv9F[o7$$\"33+++Nw,9gF[o$\"3li 'yl!G5<i>F[o7$$\"3!)*\\i:yF[o7$$\"3(** *\\ivbzGiF[o$\"3$o3@(RJ;))>F[o7$$\"39+vopLF]iF[o$\"3/ufB$)Ql&*>F[o7$$ \"3I++vj6vriF[o$\"3!GhGAzi%**>F[o7$$\"3j+](=v1ZJ'F[o$\"3OE=;8Q,'*>F[o7 $$\"33+++SBmdjF[o$\"3vt_>g;%y(>F[o7$$\"3.+DchufdkF[o$\"3&RKZ%p0az=F[o7 $$\"3)***\\7$eKvb'F[o$\"3XI3Q9US1hI&F:7$$\"3O*\\iSP\"yAqF [o$\"3q'[Oi)3\\H=F:7$$\"3[*******)3C6rF[o$!39_+$=CaUq\"F:7$$\"3#)*\\i! **\\H3sF[o$!3QZ43'\\Wg^&F:7$$\"3;+]73\"\\`I(F[o$!3f@++?Yl?\"*F:7$$\"3 \\+v=F[o 7$$\"3g+vo\\,clzF[o$!3]g*>5u2/&>F[o7$$\"3G+]i?\"*)>1)F[o$!37q%G?k8%H=F [o7$$\"3'3]i:4=%e\")F[o$!3O:\"[Dm'eS;F[o7$$\"3X,+]iq%[D)F[o$!3Yr))Q@$G 4R\"F[o7$$\"3c,vof*)RV$)F[o$!3l4P*4?!)f6\"F[o7$$\"3n,](o&3&>V)F[o$!3l' [`_J*>h!)F:7$$\"3z,D1aF]?&)F[o$!3/NKaX8T5ZF:7$$\"38++D^Y04')F[o$!3'[0V zkjA@\"F:7$$\"3<,DJN)pXq)F[o$\"3\">%4AS*=-g#F:7$$\"3V+]P>]3+))F[o$\"3= !Hh]!45=jF:7$$\"3o*\\PM?+c*))F[o$\"3HOu()o\">h!)*F:7$$\"3'*)***\\(Q:6* *)F[o$\"3oD(o(\\tt$H\"F[o7$$\"3i*\\(o>)G.4*F[o$\"3?=5#>j5!p:F[o7$$\"3G +](=DU&*=*F[o$\"3htwER$4Fy\"F[o7$$\"3'4]iSob()G*F[o$\"3iGk?PqWE>F[o7$$ \"3&)***\\i6pzQ*F[o$\"3Ip$=]&Ge%*>F[o7$$\"3H*\\P%3>:J%*F[o$\"39mbknu$) **>F[o7$$\"3]+]i+ZLu%*F[o$\"3+<\"[O`%=!*>F[o7$$\"3#**\\7G\\F[o7$$\"3N*****\\G+2c*F[o$\"3n,I_UWbE>F[o7$$\"3****\\Ppe1Z'*F [o$\"32J:KPje0=F[o7$$\"3j++v`9VL(*F[o$\"3(=%Gnk/)3j\"F[o7$$\"3B,v$4l3n #)*F[o$\"3H92LjI'yQ\"F[o7$$\"31+]7[e)*>**F[o$\"3g[ot>Wo'4\"F[o7$$\"3*) \\7`/jK,5!#;$\"3!))e>0N*[uwF:7$$\"3&****\\U-a1,\"F]_m$\"3^h*=U[;e6%F:7 $$\"3&)*\\P9Z!H?5F]_m$\"3tczeW3Y3HF-7$$\"3#***\\i=p#*H5F]_m$!36:45PU*[ a$F:7$$\"3++D\"eOj&R5F]_m$!3CzT8.&o$\\sF:7$$\"3*)*****H\")*>\\5F]_m$!3 8trvY,ao5F[o7$$\"3')*\\il>F'e5F]_m$!3Mzly0:[m8F[o7$$\"3$)**\\7!ea!o5F] _m$!36?W'e)y)fh\"F[o7$$\"3!)*\\(oj>[x5F]_m$!37h&>c_:#3=F[o7$$\"3&**** \\sM4p3\"F]_m$!3f0US54NO>F[o7$$\"3![(=U'3q94\"F]_m$!3>$e'*o%\\*Q(>F[o7 $$\"3#)\\PfD3.'4\"F]_m$!3C+AG\"HF]*>F[o7$$\"3&[ilZc\"f+6F]_m$!3K?vVTAd **>F[o7$$\"3))*\\PRI_^5\"F]_m$!3,7a'*3?\\()>F[o7$$\"3%*\\7G#ytU6\"F]_m $!3(>cCq9&*R\">F[o7$$\"3****\\ig_RB6F]_m$!3\\u'**\\%f(px\"F[o7$$\"3-]i :2J_L6F]_m$!3IN<$*=$3gb\"F[o7$$\"30+vo`4lV6F]_m$!3]v?GHnTr7F[o7$$\"33] (=-!)yP:\"F]_m$!3\"e**3k.z$[$*F:7$$\"3$****\\nk1R;\"F]_m$!3OHbMqGM+cF: 7$$\"3y*\\7t'p+t6F]_m$!3)\\u>(pdfK?F:7$$\"3\")**\\(yG2@=\"F]_m$\"3yi_e 3\")H-;F:7$$\"3%)*\\P%3w?\">\"F]_m$\"3=#zxr3gU=&F:7$$\"3()******GzI+7F ]_m$\"3_n!fL4g\\f)F:7$$\"3-+]ijW-57F]_m$\"355Jo4X0#>\"F[o7$$\"3+++D)*4 u>7F]_m$\"3)[j-yDP(z9F[o7$$\"3(***\\(G`d%H7F]_m$\"3%ytZ!HUr6F[o7$$ \"32]7.#pyzC\"F]_m$\"3nbn^:z3q>F[o7$$\"3piS;[)z,D\"F]_m$\"3\"ROxg#\\M$ )>F[o7$$\"38voH/5Q_7F]_m$\"3k4j%[weF*>F[o7$$\"3e(oH/;#ea7F]_m$\"3#*>yT !>6$)*>F[o7$$\"3-+Dc;Lyc7F]_m$\"3%fxKPW\"****>F[o7$$\"3Y7`psW)*e7F]_m$ \"3w#yE*oiz(*>F[o7$$\"3\"\\7G)Gc=h7F]_m$\"3QN\"R'>*H<*>F[o7$$\"3NP4'\\ y'Qj7F]_m$\"3PP/x]T!=)>F[o7$$\"3(*\\P4Tzel7F]_m$\"3Qv[q&>Q!o>F[o7$$\"3 .v$fLD!**p7F]_m$\"3bh]!\\\")*4H>F[o7$$\"3#***\\ilDRu7F]_m$\"3s$>=\"Qk@ v=F[o7$$\"3#*\\P%eM=SG\"F]_m$\"3'**GK2GEvq\"F[o7$$\"3#**\\ig7WOH\"F]_m $\"3'[$>+fluw9F[o7$$\"35]7G1*pKI\"F]_m$\"3c:kJ1TS\">\"F[o7$$\"35++]'o& *GJ\"F]_m$\"3g5%*>htT?')F:7$$\"30]7`&z_>K\"F]_m$\"3oqW&pWz\"G_F:7$$\"3 ***\\iX!*45L\"F]_m$\"3'[\\]p^g[m\"F:7$$\"3$*\\Pf8q1S8F]_m$!3f8fv?u$H&> F:7$$\"3))**\\iAT7\\8F]_m$!3Bkl$)['Ho]&F:7$$\"3%****\\7T\"fe8F]_m$!3@z qT#zKs-*F:7$$\"3#)**\\()*pe!o8F]_m$!3h$Rst*e\\A7F[o7$$\"3))****\\))f_x 8F]_m$!3I9-D*fq&)\\\"F[o7$$\"3$***\\7xK*pQ\"F]_m$!308:tR\"z5s\"F[o7$$ \"3'*\\(=AG\\iR\"F]_m$!3()GQ>-?>z=F[o7$$\"3(**\\7tG0bS\"F]_m$!3%f)4*Qo *3t>F[o7$$\"3[Pfe)G>yS\"F]_m$!3'o.:QO/h)>F[o7$$\"3)\\Pf)*GL,T\"F]_m$!3 ?(=yO)e'[*>F[o7$$\"3\\7G8\"HZCT\"F]_m$!3'**H\\5[b$**>F[o7$$\"3**\\iS#H hZT\"F]_m$!3;[>kTNc**>F[o7$$\"3](ozOHvqT\"F]_m$!3y1I!*>'*[&*>F[o7$$\"3 +DJ&\\H*Q>9F]_m$!3I!zm*RC9()>F[o7$$\"3^ilA'H.F[o 7$$\"3,++](HF[o7$$\"3%*\\7.xuqL9F]_m$!3=QBxGbEU= F[o7$$\"3))*\\ill(RV9F]_m$!3>gtXkK%yl\"F[o7$$\"3)*\\P4Oy3`9F]_m$!3G%o! =TtM69F[o7$$\"3#***\\i:!yFY\"F]_m$!3__-6ir+76F[o7$$\"3.](=-\"36s9F]_m$ !3cN&**[A8?%yF:7$$\"3(**\\7[gV9[\"F]_m$!3qSUu8\"G:H%F:7$$\"33]iS*Rw2\\ \"F]_m$!3ZQ=1h&z&>fF-7$$\"3?+++%>4,]\"F]_m$\"3X,'[nXw\"GJF:7$$\"36+vVC Ll4:F]_m$\"38')H3@&p\">oF:7$$\"3-+]([X(>>:F]_m$\"3E*f$>n^CE5F[o7$$\"36 +DJ&eT(G:F]_m$\"3GcGAPVHL8F[o7$$\"3-++v:dGQ:F]_m$\"3k\">T[p5>f\"F[o7$$ \"3!*\\(of\"3vZ:F]_m$\"3)yv$)\\2'H\"z\"F[o7$$\"3(**\\(=;f@d:F]_m$\"3ur K)*y8oE>F[o7$$\"3+voHm%[>c\"F]_m$\"3Dv1stxwo>F[o7$$\"3/]iS;5om:F]_m$\" 3wcNd6&HK*>F[o7$$\"32Dc^mNTr:F]_m$\"3_)3')zgZ)**>F[o7$$\"35+]i;h9w:F]_ m$\"37#GEp\"Gc))>F[o7$$\"3++D\"[+TNf\"F]_m$\"33_!3')f:mz\"F[o7$$\"3\"* *****H*e$4h\"F]_m$\"3!>>A)[KS*Q\"F[o7$$\"3s*\\Pa9/4i\"F]_m$\"3$\\U*z^I (p2\"F[o7$$\"3*)**\\(yRs3j\"F]_m$\"3kUz!\\A'y=sF:7$$\"31+DJ]1%3k\"F]_m $\"3J-*e_3q=Q$F:7$$\"3()***\\F!*33l\"F]_m$!3o3[-Ym=!*eF-7$$\"3&**\\ili C(f;F]_m$!3h\">)[,5cDTF:7$$\"3.+]P].ko;F]_m$!3?0i'4*RFJvF:7$$\"35+v=ug bx;F]_m$!3**Rw%Gg:)p5F[o7$$\"3=+++)zrko\"F]_m$!3+#y&4Hzd_8F[o7$$\"3#* \\i:>y(fp\"F]_m$!31e$p2R@mg\"F[o7$$\"3-+DJSQ[0i%R$>F[o7$$\"3()**\\i#)e\\CF[o7$$\"3z\\(=-7X!HF[o7$$\"3r*\\7yN%fLU$p)>F[o7$$\"3j\\iS&fV\"QTPV#yz&>F[o7$$\"3\"******H$Gp UF[o7$$\"36+v=38z^ \+ " 0 "" {MPLTEXT 1 206 30 "s[0]:=fsolve(determ3=0,ahl=1);" }{MPLTEXT 1 206 31 "\ns[1]:=fsolve(determ3=0,ahl=2);" }{MPLTEXT 1 206 31 "\ns[2]:= fsolve(determ3=0,ahl=4);" }{MPLTEXT 1 206 31 "\ns[3]:=fsolve(determ3=0 ,ahl=6);" }{MPLTEXT 1 206 31 "\ns[4]:=fsolve(determ3=0,ahl=7);" } {MPLTEXT 1 206 10 "\neps:=0.1:" }{MPLTEXT 1 206 10 "\nnmax:=60;" } {MPLTEXT 1 206 25 "\nfor i from 5 to nmax do " }{MPLTEXT 1 206 52 "\n \+ s[i]:=fsolve(determ3=0,ahl=s[i-1]+eps..s[i-1]+2);" }{MPLTEXT 1 206 19 "\nod: unassign('i'):" }{MPLTEXT 1 206 19 "\nseq(s[i],i=1..40);" } {MPLTEXT 1 206 79 "\n#plot([seq([s[i],0],i=1..nmax)],style=point,symbo lsize=15,view=[0..20,-1..1]);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\" sG6#\"\"!$F'F'" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"sG6#\"\"\"$\"+s .-lB!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"sG6#\"\"#$\"+7BgER!\" *" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"sG6#\"\"$$\"+>R!y\\&!\"*" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>&%\"sG6#\"\"%$\"+YFeoq!\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%%nmaxG\"#g" }}{PARA 12 "" 1 "" {XPPMATH 20 "6J$\"+s.-lB!\"*$\"+7BgERF%$\"+>R!y\\&F%$\"+YFeoqF%$\"+H)z$R')F%$\" +7w,@5!\")$\"+Xs4y6F0$\"+yoF0$\"+T]d?@F0$\"+uYlxAF0$\"+2VtMCF0$\"+RR\"=f#F0$\"+sN*)[FF0$ \"+0K(f!HF0$\"+PG0jIF0$\"+qC8?KF0$\"+.@@xLF0$\"+NVF0$\"+J&pnZ%F0$\"+k\"\\Qj%F 0$\"+(zG4z%F0$\"+H%3![\\F0$\"+i!)30^F0$\"+&pn@E&F0$\"+FtC>aF0$\"+gpKwb F0$\"+$f1Mt&F0$\"+Di[!*eF0$\"+eecZgF0$\"+\"\\XY?'F0$\"+C^shjF0" }}} {EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 47 "frequencies:=seq(evalf(2 *s[i]^2/Pi),i=1..nmax);" }{MPLTEXT 1 206 94 "\nplot([seq([evalf(2*s[i] ^2/Pi),0],i=1..nmax)],style=point,symbolsize=15,view=[0..2600,-1..1]); " }}{PARA 201 "" 0 "" {TEXT -1 0 "" }}{PARA 12 "" 1 "" {XPPMATH 20 "6# >%,frequenciesG6hn$\"+q*=3c$!\"*$\"+8Y`:)*F($\"+UsBC>!\")$\"+4>'3=$F-$ \"+B*e;v%F-$\"+tWhOmF-$\"+P$Hd$))F-$\"+Y.!\\8\"!\"($\"+%oVwT\"F8$\"+]H !=t\"F8$\"+U\"yt2#F8$\"+f#pVX#F8$\"+0jxiGF8$\"+x#*f-LF8$\"+v\"QQx$F8$ \"+)*H\\wUF8$\"+]Pc5[F8$\"+I/0w`F8$\"+HI&H(fF8$\"+f:F,mF8$\"+@g+hsF8$ \"+-k:_zF8$\"+:Fsu')F8$\"+`\\qG%*F8$\"+6.T@5!\"'$\"+?<4.6Fen$\"+AZ\"z= \"Fen$\"+;$zeF\"Fen$\"+.b)pO\"Fen$\"+$GL7Y\"Fen$\"+aEie:Fen$\"+?O:f;Fe n$\"+xh#Gw\"Fen$\"+F.kp=Fen$\"+qgfz>Fen$\"+1Mp#4#Fen$\"+MB$*3AFen$\"+b GJGBFen$\"+o\\$3X#Fen$\"+v')\\wDFen$\"+tRI0FFen$\"+l3DPGFen$\"+\\$RB(H Fen$\"+D%p06$Fen$\"+'4T>D$Fen$\"+dVX'R$Fen$\"+5#4Ta$Fen$\"+gc!\\p$Fen$ \"+)pV)[QFen$\"+JL#f+%Fen$\"+eX9mTFen$\"+wt]HVFen$\"+'y6g\\%Fen$\"+!zd cm%Fen$\"+&QX%Q[Fen$\"+rXP9]Fen$\"+c`W$>&Fen$\"+Gxlv`Fen$\"+%p65c&Fen$ \"+_s]\\dFen" }}{PARA 13 "" 1 "" {GLPLOT2D 400 400 400 {PLOTDATA 2 "6' -%'CURVESG6$7hn7$$\"39+++q*=3c$!#<$\"\"!F,7$$\"3_+++8Y`:)*F*F+7$$\"3)) *****>CPU#>!#;F+7$$\"30+++4>'3=$F3F+7$$\"35+++B*e;v%F3F+7$$\"3?+++tWhO mF3F+7$$\"3O*****pLHd$))F3F+7$$\"3,+++Y.!\\8\"!#:F+7$$\"3&******RoVwT \"FCF+7$$\"3,+++]H!=t\"FCF+7$$\"3.+++U\"yt2#FCF+7$$\"31+++f#pVX#FCF+7$ $\"3()*****\\IwF'GFCF+7$$\"3#******pF*f-LFCF+7$$\"3@+++v\"QQx$FCF+7$$ \"3%******z*H\\wUFCF+7$$\"3.+++]Pc5[FCF+7$$\"3B+++I/0w`FCF+7$$\"3C+++H I&H(fFCF+7$$\"3[+++f:F,mFCF+7$$\"3Y*****4-15E(FCF+7$$\"3`+++-k:_zFCF+7 $$\"3L+++:Fsu')FCF+7$$\"3;+++`\\qG%*FCF+7$$\"3)******4J59-\"!#9F+7$$\" 3!*******><4.6FapF+7$$\"3)******>s9z=\"FapF+7$$\"3-+++;$zeF\"FapF+7$$ \"3))*****H]&)pO\"FapF+7$$\"3.+++$GL7Y\"FapF+7$$\"3,+++aEie:FapF+7$$\" 30+++?O:f;FapF+7$$\"3#******pFapF+7$$\"3))*****fS$p#4#FapF+7$$\"33+++MB$*3AFapF+7$$\"38+++bGJG BFapF+7$$\"37+++o\\$3X#FapF+7$$\"3#)*****\\n)\\wDFapF+7$$\"3/+++tRI0FF apF+7$$\"3)******\\'3DPGFapF+7$$\"3')******[$RB(HFapF+7$$\"39+++D%p06$ FapF+7$$\"3-+++'4T>D$FapF+7$$\"33+++dVX'R$FapF+7$$\"33+++5#4Ta$FapF+7$ $\"3!*******fc!\\p$FapF+7$$\"38+++)pV)[QFapF+7$$\"3'******4LBf+%FapF+7 $$\"3'******zbWh;%FapF+7$$\"3d*****fP2&HVFapF+7$$\"3/+++'y6g\\%FapF+7$ $\"3?+++!zdcm%FapF+7$$\"3*******\\QX%Q[FapF+7$$\"3H+++rXP9]FapF+7$$\"3 ;+++c`W$>&FapF+7$$\"38+++Gxlv`FapF+7$$\"3\")*****Rp65c&FapF+7$$\"3M+++ _s]\\dFapF+-%'COLOURG6&%$RGBG$\"#5!\"\"F+F+-%+AXESLABELSG6$Q!6\"Few-%' SYMBOLG6$%(DEFAULTG\"#:-%&STYLEG6#%&POINTG-%%VIEWG6$;F+$\"%+EF,;$FawF, $\"\"\"F," 1 5 0 1 15 0 2 9 1 4 2 1.000000 45.000000 45.000000 0 0 "Cu rve 1" }}}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 42 "multiply_the_ above_by > sqrt(E*I/rho)/L^2;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#2*&*( %\"EG\"\"\"%$rhoG!\"\"^#F'F'#F'\"\"#%\"LG!\"#%6multiply_the_above_byG " }}}{EXCHG {PARA 201 "" 0 "" {TEXT -1 0 "" }}{PARA 201 "" 0 "" {TEXT -1 0 "" }}}{PARA 201 "" 0 "" {TEXT 209 46 "Our brass bar: (see self_lo aded_brass_bar.mws)" }}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 120 " w:=0.025;\nh:=0.0043;\nII:=w*h^3/12;\narea:=w*h;\nsigma := 8600;\nEE:= 9.1e10;\nrho:=sigma*area; # density per unit length" }{MPLTEXT 1 206 20 "\nEI:=EE*II;\nL:=1.33;" }{MPLTEXT 1 206 23 "\nfac:=sqrt(EI/rho )/L^2;" }{MPLTEXT 1 206 52 "\nhertz:=fac*frequencies;\n#omegas:=evalf( 2*Pi)*hertz;" }}{PARA 201 "" 0 "" {TEXT -1 0 "" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\"wG$\"#D!\"$" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%\" hG$\"#V!\"%" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#IIG$\"+LeRc;!#>" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#>%%areaG$\"%v5!\"(" }}{PARA 11 "" 1 " " {XPPMATH 20 "6#>%&sigmaG\"%+')" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>% #EEG$\"#\"*\"\"*" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$rhoG$\"(+]C*!\" (" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%#EIG$\"+3-K2:!\")" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%\"LG$\"$L\"!\"#" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$facG$\"+!>&o#G#!\"*" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%&her tzG6hn$\"+F(G#G\")!\"*$\"+\\vdSA!\")$\"+ayU#R%F+$\"+rj!4E(F+$\"+XTl%3 \"!\"($\"+d,$\\^\"F2$\"+])=p,#F2$\"+7-i!f#F2$\"+UU.OKF2$\"+Y4;`RF2$\"+ >.+UZF2$\"+eBb-cF2$\"+uq\"[`'F2$\"+eWzQvF2$\"+7X[9')F2$\"+Ks)=w*F2$\"+ j-5)4\"!\"'$\"+qI=F7FK$\"+UrVj8FK$\"+#[io]\"FK$\"+\"4ful\"FK$\"+lpA:=F K$\"+2h;!)>FK$\"+;lF_@FK$\"+!>e:L#FK$\"+K6,=DFK$\"+V`j6FFK$\"+>3V7HFK$ \"+jvR?JFK$\"+ub`NLFK$\"+\\[%yb$FK$\"+'RDty$FK$\"+2s(R-%FK$\"+&G+yE%FK $\"+JYz=XFK$\"+X-'px%FK$\"+DrHU]FK$\"+s_![J&FK$\"+&o%[%f&FK$\"+n`L\")e FK$\"+9tNvhFK$\"+I0bwkFK$\"+7]\"\\y'FK$\"+f2X+rFK$\"+zx:BuFK$\"+fg.`xF K$\"+0c3!4)FK$\"+HkIM%)FK$\"+3&)p&y)FK$\"+g=EW\"*FK$\"+#['**4&*FK$\"+n B!H))*FK$\"+_zHE5!\"&$\"+%zA]1\"F^r$\"+iZY/6F^r$\"+dQiW6F^r$\"+!3+b=\" F^r$\"+GM4F7F^r$\"+/RSp7F^r$\"+0:V78F^r" }}}{EXCHG {PARA 201 "> " 0 " " {MPLTEXT 1 206 69 "semitones_from_A440:=seq(12*log[2](hertz[i]/440), i=1..nops([hertz])):" }{MPLTEXT 1 206 12 "\nevalf(%,3);" }}{PARA 12 " " 1 "" {XPPMATH 20 "6hn$!$\"p!\"\"$!$:&F%$!$*RF%$!$7$F%$!$U#F%$!$&=F%$ !$N\"F%$!$<*!\"#$!$K&F4$F/F4$\"$I\"F4$\"$=%F4$\"$&oF4$\"$K*F4$\"$;\"F% $\"$Q\"F%$\"$e\"F%$\"$y\"F%$\"$'>F%$\"$8#F%$\"$I#F%$\"$X#F%$\"$g#F%$\" $v#F%$\"$*GF%$\"$-$F%$\"$:$F%$\"$F$F%$\"$R$F%$\"$^$F%$\"$i$F%$\"$t$F%$ \"$$QF%$\"$$RF%$\"$.%F%$\"$8%F%$\"$A%F%$\"$J%F%$\"$S%F%$\"$\\%F%$\"$d% F%$\"$m%F%$\"$u%F%$\"$\"[F%$\"$*[F%$\"$(\\F%$\"$/&F%$\"$6&F%$\"$=&F%$ \"$D&F%$\"$K&F%$\"$R&F%$\"$X&F%$\"$_&F%$\"$e&F%$\"$k&F%$\"$q&F%$\"$w&F %$\"$#eF%$\"$)eF%" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 36 "no te_name:=proc(semitones_from_A440)" }{MPLTEXT 1 206 53 "\n local not e_names,offset,octave,soctave,nname,ind;" }{MPLTEXT 1 206 14 "\n off set:=5;" }{MPLTEXT 1 206 49 "\n octave:=floor(offset+semitones_from_ A440/12);" }{MPLTEXT 1 206 36 "\n soctave:=convert(octave,string);" }{MPLTEXT 1 206 71 "\n note_names:=[\"A\",\"Bb\",\"B\",\"C\",\"Db\", \"D\",\"Eb\",\"E\",\"F\",\"Gb\",\"G\",\"Ab\"];" }{MPLTEXT 1 206 51 "\n ind:=1+(floor(0.5+semitones_from_A440) mod 12);" }{MPLTEXT 1 206 27 "\n nname:=note_names[ind];" }{MPLTEXT 1 206 34 "\n cat(note_n ames[ind],soctave);" }{MPLTEXT 1 206 10 "\nend proc:" }{MPLTEXT 1 206 74 "\nseq( note_name(semitones_from_A440[i]),i=1..nops([semitones_from _A440]));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6hnQ$C-16\"Q#F0F$Q#F1F$Q#D2 F$Q#A2F$Q$Eb3F$Q#G3F$Q#C4F$Q#E4F$Q#G4F$Q$Bb5F$Q$Db5F$Q#E5F$Q$Gb5F$Q#A5 F$Q#B6F$Q$Db6F$Q$Eb6F$Q#F6F$Q$Gb6F$Q$Ab6F$Q$Bb7F$Q#B7F$Q#C7F$Q#D7F$Q$E b7F$Q#E7F$Q$Gb7F$Q#G7F$Q$Ab7F$Q#A8F$Q$Bb8F$Q#B8F$Q#C8F$Q$Db8F$Q#D8F$Q$ Eb8F$Q#E8F$Q#F8F$Q$Gb8F$Q#G8F$Q$Ab8F$FMQ#A9F$Q$Bb9F$Q#B9F$FPQ#C9F$Q$Db 9F$Q#D9F$FSQ$Eb9F$Q#E9F$FUQ#F9F$FVQ$Gb9F$Q#G9F$FXQ$Ab9F$" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 48 "libname:=\"/home/ringland/src/M apleSound/siggen\"," }{MPLTEXT 1 206 46 "\n \"/home/ringland/s rc/MapleSound/wav\"," }{MPLTEXT 1 206 18 "\n libname;" } {MPLTEXT 1 206 14 "\nwith(SigGen);" }{MPLTEXT 1 206 11 "\nwith(WAV);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(libnameG6%QE/home/ringland/src/Ma pleSound/siggen6\"QB/home/ringland/src/MapleSound/wavF'Q6/usr/local/ma ple9/libF'" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#73%$ASDG%'AttackG%&Decay G%*FM_SignalG%(FM_SineG%.Integer2BytesG%$LFOG%(LimiterG%*NormalizeG%)Q uantizeG%)SawtoothG%'SignalG%%SineG%-ViewEnvelopeG%+WriteSoundG%*word2 byteG%.writeintegersG" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7$%(ReadWAVG% )WriteWAVG" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "evalf(2*Pi); " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"+3`=$G'!\"*" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 332 "mySine:=proc(Amp,freq,phase,tau,t)\n A mp*evalhf(exp(-t/tau)*sin(2*Pi*freq*t+phase)):\nend proc:\nmySignal:=p roc(amplitudes,frequencs,decayconst,t)\n# print(decayconst,t);\n# \+ print(mySine(amplitudes[1],frequencs[1],0,decayconst,t));\n add( my Sine(amplitudes[i],frequencs[i],0,decayconst,t) \n ,i=1..nops(amplit udes));\nend proc:" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 14 "l owestmode:=5;" }{MPLTEXT 1 206 27 "\nhighestmode:=15;\ntau:=4.0;" } {MPLTEXT 1 206 24 "\nfreqs:=[seq(hertz[i],i=" }{MPLTEXT 1 206 10 "lowe stmode" }{MPLTEXT 1 206 36 "..highestmode)]; # Enter Frequencies" } {MPLTEXT 1 206 83 "\namps:=[seq(evalf(exp(-(2*i/nmax))),i=lowestmode.. highestmode)]; # Enter Amplitudes" }{MPLTEXT 1 206 31 "\nDuration:=2.0 : # 2 sec. sample" }{MPLTEXT 1 206 76 "\nSampleFreq:=22050: # For test ing, and speed, try at half the CD sample rate" }{MPLTEXT 1 206 29 "\n deltaT:=evalf(1/SampleFreq);" }{MPLTEXT 1 206 37 "\nNpoints:=round(Dur ation*SampleFreq);" }{MPLTEXT 1 206 167 "\nSigData:=Array([seq(mySigna l(amps,freqs,tau,i*deltaT),i=1..Npoints)],datatype=float):\n#SigData:= Array([seq(Signal(amps,freqs,i*deltaT),i=1..Npoints)],datatype=float): " }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%+lowestmodeG\"\"&" }}{PARA 11 " " 1 "" {XPPMATH 20 "6#>%,highestmodeG\"#:" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%$tauG$\"#S!\"\"" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>% &freqsG7-$\"+XTl%3\"!\"($\"+d,$\\^\"F($\"+])=p,#F($\"+7-i!f#F($\"+UU.O KF($\"+Y4;`RF($\"+>.+UZF($\"+eBb-cF($\"+uq\"[`'F($\"+eWzQvF($\"+7X[9') F(" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%%ampsG7-$\"+\\s\"[Y)!#5$\"+JvI (=)F($\"+kc*)=zF($\"+$Q$GfwF($\"+2A=3uF($\"+1JJlrF($\"+,iSIpF($\"+g/?. nF($\"+5MW$['F($\"+`3*3F'F($\"+(f1`1'F(" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%'deltaTG$\"+#RZ^`%!#9" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%(Np ointsG\"&+T%" }}}{EXCHG {PARA 201 "> " 0 "" {MPLTEXT 1 206 83 "WriteWA V(\"/home/ringland/src/brassbar_modes_5_to_15.wav\", SigData, SampleFr eq, 16);" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#%41~channel(s)~found.G" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#%@44100~points~per~channel~found.G" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#Q9Quantizing~to~16~bits...6\"" }} {PARA 11 "" 1 "" {XPPMATH 20 "6#QXSaving~to~/home/ringland/src/brassba r_modes_5_to_15.wav6\"" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#Q%Done6\"" } }}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 7 "?evalhf" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "eva" }}}{EXCHG {PARA 201 "" 0 "" {TEXT -1 0 "" }}}{PARA 205 "" 0 "" {TEXT -1 0 "" }}}{MARK "27 0 0" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 }{PAGENUMBERS 0 1 2 33 1 1 }