{VERSION 2 3 "IBM INTEL NT" "2.3" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "R3 Fo nt 0" -1 256 1 {CSTYLE "" -1 -1 "Helvetica" 1 12 255 0 0 1 2 1 2 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 12 0 0 0 1 2 2 2 0 0 0 0 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }} {SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 79 "PHYSICS RESOURCE PACKETS P ROJECT File: wvbend.MWS" }}{PARA 0 "" 0 "" {TEXT -1 85 "Rose-Hulman Institute of Technology \+ Author: mjm" }}{PARA 0 "" 0 "" {TEXT -1 97 "http:/ /www.rose-hulman.edu/~moloney Softwa re: Maple V Release 4" }}{PARA 0 "" 0 "" {TEXT -1 0 "" }}{PARA 0 "" 0 "" {TEXT -1 79 "bending of wavefronts due to velocity increasing with \+ height. Use a sequence of" }}{PARA 0 "" 0 "" {TEXT -1 17 "parametric p lots." }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "restart; with(plots):" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 24 "N:=20: radius:=5; h:=10;" }} }{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for k from 1 to N do" }} {PARA 0 "> " 0 "" {MPLTEXT 1 0 66 "circ.k:=[(radius+k/h)*cos(t),k+(rad ius+k/h)*sin(t),t=-Pi/2..Pi/2];" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "o d:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "plot(\{circ.(1..N)\}, color=black,title=`wave speed increases with height`);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "for k from 1 to N do" }}{PARA 0 "> \+ " 0 "" {MPLTEXT 1 0 66 "circ.k:=[(radius-k/h)*cos(t),k+(radius-k/h)*si n(t),t=-Pi/2..Pi/2];" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 3 "od:" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 73 "plot(\{circ.(1..N)\},color=b lack,title=`wave speed decreases with height`);" }}}}{MARK "5 0 0" 73 }{VIEWOPTS 1 1 0 1 1 1803 }