# 3.05.HagenPoiseuille.draw # # Praat script that creates figure 3.5 of # Boersma (1998): Functional Phonology. # http://fonsg3.hum.uva.nl/paul/ # Erase all Times 12 Line width... 2 Viewport... 1 5 0 2.3 Create simple Matrix... parabol 2 500 0 # Upper half of parabolic profile. Formula... if row=1 then 80-30*((col-1)/(ncol-1))^1.4142 else 50+50*((col-1)/(ncol-1))^0.7071 fi To Polygon Draw... 0 100 0 100 Remove # Lower half. select Matrix parabol Formula... if row=2 then 100 - self else self fi To Polygon Draw... 0 100 0 100 Remove Dotted line # Upper half of boundary-layer profile. select Matrix parabol Formula... if row=1 then 73-23*((col-1)/(ncol-1))^4 else 50+50*((col-1)/(ncol-1))^0.25 fi To Polygon Draw... 0 100 0 100 Remove # Lower half. select Matrix parabol Formula... if row=2 then 100 - self else self fi To Polygon Draw... 0 100 0 100 Remove Plain line select Matrix parabol Remove Line width... 4 Draw line... 0 0 100 0 Draw line... 0 100 100 100 Line width... 2 Dotted line # horizontal line at y=0 Draw line... 0 50 80 50 #Draw line... 50 0 50 100 Plain line Draw arrow... 80 50 100 50 Text... 100 Left 50 Half %%v__max Text... 0 Right 50 Half %y = 0 Text... 50 Centre 100 Bottom %v = 0 Draw arrow... 140 50 160 50 Text... 160 Left 50 Half %x Draw arrow... 140 50 140 90 Text... 140 Centre 90 Bottom %y Draw arrow... 140 50 130 35 Text... 130 Right 35 Top %z Viewport... 0.5 6 0.2 2.0 Line width... 1