%!PS-Adobe-2.0
%%Creator: dvips 5.47 Copyright 1986-91 Radical Eye Software
%%Title: airfoil.dvi
%%Pages: 6 1
%%BoundingBox: 0 0 612 792
%%DocumentFonts: Times-Bold Times-Roman Courier
%%EndComments
%%BeginProcSet: tex.pro
/TeXDict 200 dict def TeXDict begin /N /def load def /B{bind def}N /S /exch
load def /X{S N}B /TR /translate load N /isls false N /vsize 10 N /@rigin{
isls{[0 1 -1 0 0 0]concat}if 72 Resolution div 72 VResolution div neg scale
Resolution VResolution vsize neg mul TR matrix currentmatrix dup dup 4 get
round 4 exch put dup dup 5 get round 5 exch put setmatrix}N /@letter{/vsize 10
N}B /@landscape{/isls true N /vsize -1 N}B /@a4{/vsize 10.6929133858 N}B /@a3{
/vsize 15.5531 N}B /@ledger{/vsize 16 N}B /@legal{/vsize 13 N}B /@manualfeed{
statusdict /manualfeed true put}B /@copies{/#copies X}B /FMat[1 0 0 -1 0 0]N
/FBB[0 0 0 0]N /nn 0 N /IE 0 N /ctr 0 N /df-tail{/nn 8 dict N nn begin
/FontType 3 N /FontMatrix fntrx N /FontBBox FBB N string /base X array
/BitMaps X /BuildChar{CharBuilder}N /Encoding IE N end dup{/foo setfont}2
array copy cvx N load 0 nn put /ctr 0 N[}B /df{/sf 1 N /fntrx FMat N df-tail}
B /dfs{div /sf X /fntrx[sf 0 0 sf neg 0 0]N df-tail}B /E{pop nn dup definefont
setfont}B /ch-width{ch-data dup length 5 sub get}B /ch-height{ch-data dup
length 4 sub get}B /ch-xoff{128 ch-data dup length 3 sub get sub}B /ch-yoff{
ch-data dup length 2 sub get 127 sub}B /ch-dx{ch-data dup length 1 sub get}B
/ch-image{ch-data dup type /stringtype ne{ctr get /ctr ctr 1 add N}if}B /id 0
N /rw 0 N /rc 0 N /gp 0 N /cp 0 N /G 0 N /sf 0 N /CharBuilder{save 3 1 roll S
dup /base get 2 index get S /BitMaps get S get /ch-data X pop /ctr 0 N ch-dx 0
ch-xoff ch-yoff ch-height sub ch-xoff ch-width add ch-yoff setcachedevice
ch-width ch-height true[1 0 0 -1 -.1 ch-xoff sub ch-yoff .1 add]{ch-image}
imagemask restore}B /D{/cc X dup type /stringtype ne{]}if nn /base get cc ctr
put nn /BitMaps get S ctr S sf 1 ne{dup dup length 1 sub dup 2 index S get sf
div put}if put /ctr ctr 1 add N}B /I{cc 1 add D}B /bop{userdict /bop-hook
known{bop-hook}if /SI save N @rigin 0 0 moveto}N /eop{clear SI restore
showpage userdict /eop-hook known{eop-hook}if}N /@start{userdict /start-hook
known{start-hook}if /VResolution X /Resolution X 1000 div /DVImag X /IE 256
array N 0 1 255{IE S 1 string dup 0 3 index put cvn put}for}N /p /show load N
/RMat[1 0 0 -1 0 0]N /BDot 260 string N /rulex 0 N /ruley 0 N /v{/ruley X
/rulex X V}B /V statusdict begin /product where{pop product dup length 7 ge{0
7 getinterval(Display)eq}{pop false}ifelse}{false}ifelse end{{gsave TR -.1 -.1
TR 1 1 scale rulex ruley false RMat{BDot}imagemask grestore}}{{gsave TR -.1
-.1 TR rulex ruley scale 1 1 false RMat{BDot}imagemask grestore}}ifelse B /a{
moveto}B /delta 0 N /tail{dup /delta X 0 rmoveto}B /M{S p delta add tail}B /b{
S p tail}B /c{-4 M}B /d{-3 M}B /e{-2 M}B /f{-1 M}B /g{0 M}B /h{1 M}B /i{2 M}B
/j{3 M}B /k{4 M}B /w{0 rmoveto}B /l{p -4 w}B /m{p -3 w}B /n{p -2 w}B /o{p -1 w
}B /q{p 1 w}B /r{p 2 w}B /s{p 3 w}B /t{p 4 w}B /x{0 S rmoveto}B /y{3 2 roll p
a}B /bos{/SS save N}B /eos{clear SS restore}B end
%%EndProcSet
%%BeginProcSet: texps.pro
TeXDict begin /rf{655360 div mul Resolution mul 7227 div /PixPerEm X findfont
dup length 1 add dict /nn X{1 index /FID ne{nn 3 1 roll put}{pop pop}ifelse}
forall 256 dict begin nn /Encoding get 0 1 255{2 copy get 3 index 2 index get
1000 mul PixPerEm div def pop}for pop pop nn /Metrics currentdict put end
/fontname X /nn dup nn definefont[PixPerEm 0 0 PixPerEm neg 0 0]makefont N
fontname{/foo setfont}2 array copy cvx N fontname load 0 nn put}N
/ObliqueSlant{dup sin S cos div neg}B /SlantFont{/foo X[1 0 foo 1 0 0]
TransFont}N /ExtendFont{/foo X 3 2 roll[S{foo div}forall]3 1 roll[foo 0 0 1 0
0]TransFont}N /TransFont{S findfont S makefont dup length dict /nn X{1 index
/FID ne{nn 3 1 roll put}{pop pop}ifelse}forall dup nn definefont pop}N end
%%EndProcSet
%%BeginProcSet: special.pro
TeXDict begin /SDict 200 dict N SDict begin /@SpecialDefaults{/hs 612 N /vs
792 N /ho 0 N /vo 0 N /hsc 1 N /vsc 1 N /ang 0 N /CLIP false N /BBcalc false N
/p 3 def}B /@scaleunit 100 N /@hscale{@scaleunit div /hsc X}B /@vscale{
@scaleunit div /vsc X}B /@hsize{/hs X /CLIP true N}B /@vsize{/vs X /CLIP true
N}B /@hoffset{/ho X}B /@voffset{/vo X}B /@angle{/ang X}B /@rwi{10 div /rwi X}
B /@llx{/llx X}B /@lly{/lly X}B /@urx{/urx X}B /@ury{/ury X /BBcalc true N}B
/magscale true def end /@MacSetUp{userdict /md known{userdict /md get type
/dicttype eq{md begin /letter{}N /note{}N /legal{}N /od{txpose 1 0 mtx
defaultmatrix dtransform S atan/pa X newpath clippath mark{transform{
itransform moveto}}{transform{itransform lineto}}{6 -2 roll transform 6 -2
roll transform 6 -2 roll transform{itransform 6 2 roll itransform 6 2 roll
itransform 6 2 roll curveto}}{{closepath}}pathforall newpath counttomark array
astore /gc xdf pop ct 39 0 put 10 fz 0 fs 2 F/|______Courier fnt invertflag{
PaintBlack}if}N /txpose{pxs pys scale ppr aload pop por{noflips{pop S neg S TR
pop 1 -1 scale}if xflip yflip and{pop S neg S TR 180 rotate 1 -1 scale ppr 3
get ppr 1 get neg sub neg ppr 2 get ppr 0 get neg sub neg TR}if xflip yflip
not and{pop S neg S TR pop 180 rotate ppr 3 get ppr 1 get neg sub neg 0 TR}if
yflip xflip not and{ppr 1 get neg ppr 0 get neg TR}if}{noflips{TR pop pop 270
rotate 1 -1 scale}if xflip yflip and{TR pop pop 90 rotate 1 -1 scale ppr 3 get
ppr 1 get neg sub neg ppr 2 get ppr 0 get neg sub neg TR}if xflip yflip not
and{TR pop pop 90 rotate ppr 3 get ppr 1 get neg sub neg 0 TR}if yflip xflip
not and{TR pop pop 270 rotate ppr 2 get ppr 0 get neg sub neg 0 S TR}if}
ifelse scaleby96{ppr aload pop 4 -1 roll add 2 div 3 1 roll add 2 div 2 copy
TR .96 dup scale neg S neg S TR}if}N /cp{pop pop showpage pm restore}N end}if}
if}N /normalscale{Resolution 72 div VResolution 72 div neg scale magscale{
DVImag dup scale}if}N /psfts{S 65536 div N}N /startTexFig{/psf$SavedState save
N userdict maxlength dict begin /magscale false def normalscale currentpoint
TR /psf$ury psfts /psf$urx psfts /psf$lly psfts /psf$llx psfts /psf$y psfts
/psf$x psfts currentpoint /psf$cy X /psf$cx X /psf$sx psf$x psf$urx psf$llx
sub div N /psf$sy psf$y psf$ury psf$lly sub div N psf$sx psf$sy scale psf$cx
psf$sx div psf$llx sub psf$cy psf$sy div psf$ury sub TR /showpage{}N
/erasepage{}N /copypage{}N /p 3 def @MacSetUp}N /doclip{psf$llx psf$lly
psf$urx psf$ury currentpoint 6 2 roll newpath 4 copy 4 2 roll moveto 6 -1 roll
S lineto S lineto S lineto closepath clip newpath moveto}N /endTexFig{end
psf$SavedState restore}N /@beginspecial{SDict begin /SpecialSave save N gsave
normalscale currentpoint TR @SpecialDefaults}N /@setspecial{CLIP{newpath 0 0
moveto hs 0 rlineto 0 vs rlineto hs neg 0 rlineto closepath clip}if ho vo TR
hsc vsc scale ang rotate BBcalc{rwi urx llx sub div dup scale llx neg lly neg
TR}if /showpage{}N /erasepage{}N /copypage{}N newpath}N /@endspecial{grestore
clear SpecialSave restore end}N /@defspecial{SDict begin}N /@fedspecial{end}B
/li{lineto}B /rl{rlineto}B /rc{rcurveto}B /np{/SaveX currentpoint /SaveY X N 1
setlinecap newpath}N /st{stroke SaveX SaveY moveto}N /fil{fill SaveX SaveY
moveto}N /ellipse{/endangle X /startangle X /yrad X /xrad X /savematrix matrix
currentmatrix N TR xrad yrad scale 0 0 1 startangle endangle arc savematrix
setmatrix}N end
%%EndProcSet
TeXDict begin 1000 300 300 @start /Fa [ 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 10 14 17 21 21
35 32 14 14 14 21 23 10 14 10 12 21 21 21 21 21 21 21 21 21
21 12 12 23 23 23 18 38 30 28 28 30 25 23 30 30 14 16 30 25
37 30 30 23 30 28 23 25 30 30 39 30 30 25 14 12 14 19 21 14
18 21 18 21 18 14 21 21 12 12 21 12 32 21 21 21 21 14 16 12
21 21 30 21 21 18 20 8 20 22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 14 21 21 7 21 21 21 21
7 18 21 14 14 23 23 0 21 21 21 10 0 19 15 14 18 18 21 42 42
0 18 0 14 14 14 14 14 14 14 14 0 14 14 0 14 14 14 42 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 37 0 11 0 0 0 0 25 30 37 13 0 0 0
0 0 28 0 0 0 12 0 0 12 21 30 21 0 0 0 0 ] /Times-Roman 1000
655360 rf /Fb [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 20 20 20 20 20 20 20 20 20 20 20 20 20
20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20
20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20
20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20
20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20
20 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 20 20 20 20 20 20 20 20 20 20 20 20 20 0 0 0
20 20 20 20 0 20 20 20 20 20 20 20 0 0 20 0 20 20 20 20 20
20 20 20 0 20 20 0 20 20 20 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 20 0 0 0 0 20 20 0 20 0 0 0 0 0 0 0 0 0 20 0 0 20 20
0 20 0 0 0 0 ] /Courier 1000 524288 rf /Fc [ 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 17 20
25 25 41 39 17 17 17 25 28 12 17 12 14 25 25 25 25 25 25 25
25 25 25 14 14 28 28 28 22 46 36 33 33 36 30 28 36 36 17 19
36 30 44 36 36 28 36 33 28 30 36 36 47 36 36 30 17 14 17 23
25 17 22 25 22 25 22 17 25 25 14 14 25 14 39 25 25 25 25 17
19 14 25 25 36 25 25 22 24 10 24 27 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 17 25 25 8 25 25
25 25 9 22 25 17 17 28 28 0 25 25 25 12 0 23 17 17 22 22 25
50 50 0 22 0 17 17 17 17 17 17 17 17 0 17 17 0 17 17 17 50
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 44 0 14 0 0 0 0 30 36 44 15
0 0 0 0 0 33 0 0 0 14 0 0 14 25 36 25 0 0 0 0 ] /Times-Roman
1000 786432 rf /Fd [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 18 23 29 35 35 59 55 23 23 23 35
40 18 23 18 20 35 35 35 35 35 35 35 35 35 35 20 20 40 40 40
31 65 51 47 47 51 43 39 51 51 23 27 51 43 63 51 51 39 51 47
39 43 51 51 67 51 51 43 23 20 23 33 35 23 31 35 31 35 31 23
35 35 20 20 35 20 55 35 35 35 35 23 27 20 35 35 51 35 35 31
34 14 34 38 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 23 35 35 12 35 35 35 35 13 31 35 23 23
39 39 0 35 35 35 18 0 32 25 23 31 31 35 71 71 0 31 0 23 23
23 23 23 23 23 23 0 23 23 0 23 23 23 71 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 63 0 19 0 0 0 0 43 51 63 22 0 0 0 0 0 47 0 0 0
20 0 0 20 35 51 35 0 0 0 0 ] /Times-Roman 1000 1114112 rf /Fe
[ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 18 23 39 35 35 71 59 23 23 23 35 40 18 23 18 20 35 35
35 35 35 35 35 35 35 35 23 23 40 40 40 35 66 51 47 51 51 47
43 55 55 27 35 55 47 67 51 55 43 55 51 39 47 51 51 71 51 51
47 23 20 23 41 35 23 35 39 31 39 31 23 35 39 20 23 39 20 59
39 35 39 39 31 27 23 39 35 51 35 35 31 28 16 28 37 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
23 35 35 12 35 35 35 35 20 35 35 23 23 39 39 0 35 35 35 18
0 38 25 23 35 35 35 71 71 0 35 0 23 23 23 23 23 23 23 23 0
23 23 0 23 23 23 71 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 71 0 21
0 0 0 0 47 55 71 23 0 0 0 0 0 51 0 0 0 20 0 0 20 35 51 39 0
0 0 0 ] /Times-Bold 1000 1114112 rf end
%%EndProlog
%%BeginSetup
%%Feature: *Resolution 300
TeXDict begin 
%%EndSetup
%%Page: 1 1
bop 367 482 a Fe(gnuplot)19 b Fd(airfoil)f(and)h(complex)f(variable)h(demo)
664 602 y Fc(Nicolas)13 b(Brouard)e(and)i(Alex)f(W)l(oo)262
863 y Fb(This)19 b(demo)h(shows)f(how)h(to)g(use)g(bezier)f(splines)g(to)h
(define)g(NACA)f(four)262 903 y(series)g(airfoils)g(and)h(complex)f
(variables)h(to)f(define)h(Joukowski)262 942 y(Airfoils.)39
b(It)19 b(will)h(be)g(expanded)f(after)h(overplotting)f(in)h(implemented)262
982 y(to)f(plot)h(Coefficient)f(of)h(Pressure)f(as)h(well.)1098
1021 y(Alex)g(Woo,)g(Dec.)f(1992)262 1061 y(The)g(definitions)g(below)h
(follows:)f("Bezier)h(presentation)f(of)g(airfoils",)262 1100
y(by)g(Wolfgang)h(Boehm,)f(Computer)g(Aided)h(Geometric)f(Design)h(4)g
(\(1987\))f(pp)h(17-22.)1098 1139 y(Gershon)g(Elber,)f(Nov.)h(1992)281
1179 y(m)40 b(=)20 b(percent)f(camber)281 1218 y(p)40 b(=)20
b(percent)f(chord)h(with)g(maximum)f(camber)281 1258 y(mm)h(=)g(0.6)g(#)g
(NACA6xxx)262 1297 y(thick)f(=)h(0.09)39 b(#)20 b(nine)g(percent)39
b(NACAxx09)281 1337 y(pp)20 b(=)g(0.4#)g(NACAx4xx)281 1376
y(Combined)g(this)f(implies)h(NACA6409)f(airfoil)281 1416 y(Airfoil)h
(thickness)f(function.)262 1534 y(pause)g(0)40 b("NACA)19 b(four)h(series)g
(airfoils)f(by)h(bezier)f(splines")262 1573 y(pause)g(0)40
b("Will)19 b(add)h(pressure)f(distribution)g(later)h(with)g(Overplotting")262
1613 y(mm)f(=)h(0.6)262 1652 y(thick)f(=)h(0.09)262 1692 y(pp)f(=)h(0.4)262
1731 y(set)f(xlabel)h("NACA6409)f(--)h(9\045)g(thick,)f(40\045)h(max)f
(camber,)h(6\045)g(camber")262 1770 y(x0)f(=)h(0.0)262 1810
y(y0)f(=)h(0.0)262 1849 y(x1)f(=)h(0.0)262 1889 y(y1)f(=)h(0.18556)262
1928 y(x2)f(=)h(0.03571)262 1968 y(y2)f(=)h(0.34863)262 2007
y(x3)f(=)h(0.10714)262 2046 y(y3)f(=)h(0.48919)262 2086 y(x4)f(=)h(0.21429)
262 2125 y(y4)f(=)h(0.58214)262 2165 y(x5)f(=)h(0.35714)262
2204 y(y5)f(=)h(0.55724)262 2244 y(x6)f(=)h(0.53571)262 2283
y(y6)f(=)h(0.44992)262 2323 y(x7)f(=)h(0.75000)262 2362 y(y7)f(=)h(0.30281)
262 2401 y(x8)f(=)h(1.00000)262 2441 y(y8)f(=)h(0.01050)967
2574 y Fa(1)p eop
%%Page: 2 2
bop 281 347 a Fb(Directly)20 b(defining)f(the)h(order)f(8)h(Bezier)g(basis)f
(function)h(for)f(a)h(faster)g(evaluation.)262 425 y(bez_d4_i0\(x\))f(=)99
b(\(1)20 b(-)g(x\)**4)262 465 y(bez_d4_i1\(x\))f(=)g(4)h(*)g(\(1)g(-)g
(x\)**3)f(*)h(x)262 504 y(bez_d4_i2\(x\))f(=)g(6)h(*)g(\(1)g(-)g(x\)**2)f(*)h
(x**2)262 544 y(bez_d4_i3\(x\))f(=)g(4)h(*)g(\(1)g(-)g(x\)**1)f(*)h(x**3)262
583 y(bez_d4_i4\(x\))f(=)358 b(x**4)262 662 y(bez_d8_i0\(x\))19
b(=)119 b(\(1)20 b(-)g(x\)**8)262 702 y(bez_d8_i1\(x\))f(=)39
b(8)20 b(*)g(\(1)g(-)g(x\)**7)f(*)h(x)262 741 y(bez_d8_i2\(x\))f(=)g(28)h(*)g
(\(1)g(-)g(x\)**6)f(*)h(x**2)262 780 y(bez_d8_i3\(x\))f(=)g(56)h(*)g(\(1)g(-)
g(x\)**5)f(*)h(x**3)262 820 y(bez_d8_i4\(x\))f(=)g(70)h(*)g(\(1)g(-)g(x\)**4)
f(*)h(x**4)262 859 y(bez_d8_i5\(x\))f(=)g(56)h(*)g(\(1)g(-)g(x\)**3)f(*)h
(x**5)262 899 y(bez_d8_i6\(x\))f(=)g(28)h(*)g(\(1)g(-)g(x\)**2)f(*)h(x**6)262
938 y(bez_d8_i7\(x\))f(=)39 b(8)20 b(*)g(\(1)g(-)g(x\))79 b(*)20
b(x**7)262 978 y(bez_d8_i8\(x\))f(=)378 b(x**8)262 1096 y(m0)19
b(=)h(0.0)262 1135 y(m1)f(=)h(0.1)262 1175 y(m2)f(=)h(0.1)262
1214 y(m3)f(=)h(0.1)262 1254 y(m4)f(=)h(0.0)262 1293 y(mean_y\(t\))f(=)h(m0)f
(*)h(mm)g(*)g(bez_d4_i0\(t\))f(+)h(\\)341 1333 y(m1)g(*)g(mm)g(*)f
(bez_d4_i1\(t\))h(+)f(\\)341 1372 y(m2)h(*)g(mm)g(*)f(bez_d4_i2\(t\))h(+)f
(\\)341 1411 y(m3)h(*)g(mm)g(*)f(bez_d4_i3\(t\))h(+)f(\\)341
1451 y(m4)h(*)g(mm)g(*)f(bez_d4_i4\(t\))262 1530 y(p0)g(=)h(0.0)262
1569 y(p1)f(=)h(pp)g(/)g(2)262 1609 y(p2)f(=)h(pp)262 1648
y(p3)f(=)h(\(pp)g(+)g(1\))f(/)h(2)262 1687 y(p4)f(=)h(1.0)262
1727 y(mean_x\(t\))f(=)h(p0)f(*)h(bez_d4_i0\(t\))f(+)h(\\)341
1766 y(p1)g(*)g(bez_d4_i1\(t\))f(+)h(\\)341 1806 y(p2)g(*)g(bez_d4_i2\(t\))f
(+)h(\\)341 1845 y(p3)g(*)g(bez_d4_i3\(t\))f(+)h(\\)341 1885
y(p4)g(*)g(bez_d4_i4\(t\))262 1963 y(z_x\(x\))f(=)h(x0)g(*)f(bez_d8_i0\(x\))h
(+)f(x1)h(*)g(bez_d8_i1\(x\))f(+)h(x2)g(*)f(bez_d8_i2\(x\))h(+)f(\\)281
2003 y(x3)h(*)g(bez_d8_i3\(x\))f(+)h(x4)g(*)g(bez_d8_i4\(x\))f(+)g(x5)h(*)g
(bez_d8_i5\(x\))f(+)h(\\)281 2042 y(x6)g(*)g(bez_d8_i6\(x\))f(+)h(x7)g(*)g
(bez_d8_i7\(x\))f(+)g(x8)h(*)g(bez_d8_i8\(x\))262 2121 y(z_y\(x,)f(tk\))h(=)g
(\\)321 2161 y(y0)g(*)g(tk)g(*)g(bez_d8_i0\(x\))f(+)g(y1)h(*)g(tk)g(*)g
(bez_d8_i1\(x\))f(+)h(y2)f(*)h(tk)g(*)g(bez_d8_i2\(x\))f(+)h(\\)321
2200 y(y3)g(*)g(tk)g(*)g(bez_d8_i3\(x\))f(+)g(y4)h(*)g(tk)g(*)g
(bez_d8_i4\(x\))f(+)h(y5)f(*)h(tk)g(*)g(bez_d8_i5\(x\))f(+)h(\\)321
2240 y(y6)g(*)g(tk)g(*)g(bez_d8_i6\(x\))f(+)g(y7)h(*)g(tk)g(*)g
(bez_d8_i7\(x\))f(+)h(y8)f(*)h(tk)g(*)g(bez_d8_i8\(x\))281
2358 y(Given)g(t)g(value)f(between)h(zero)f(and)h(one,)g(the)f(airfoild)h
(curve)f(is)h(defined)g(as)262 2437 y(c\(t\))f(=)h(mean\(t1\(t\)\))f(+/-)h
(z\(t2\(t\)\))f(n\(t1\(t\)\),)967 2574 y Fa(2)p eop
%%Page: 3 3
bop 281 347 a Fb(where)20 b(n)g(is)g(the)f(unit)h(normal)f(to)h(the)g(mean)g
(line.)f(See)h(the)f(above)h(paper)g(for)f(more.)281 425 y(Unfortunately,)g
(the)h(parametrization)f(of)h(c\(t\))f(is)h(not)g(the)f(same)h(for)g
(mean\(t1\))281 465 y(and)g(z\(t2\).)g(The)f(mean)h(line)g(\(and)f(its)h
(normal\))f(can)h(assume)f(linear)h(function)f(t1)h(=)g(t,)1318
504 y(-1)281 544 y(but)g(the)g(thickness)f(z_y)h(is,)g(in)f(fact,)h(a)g
(function)f(of)h(z_x)39 b(\(t\).)20 b(Since)g(it)f(is)281 583
y(not)h(obvious)g(how)f(to)h(compute)f(this)h(inverse)f(function)h
(analytically,)f(we)h(instead)281 623 y(replace)g(t)g(in)f(c\(t\))h(equation)
f(above)h(by)g(z_x\(t\))f(to)h(get:)262 702 y(c\(z_x\(t\)\))f(=)h
(mean\(z_x\(t\)\))f(+/-)g(z\(t\))h(n\(z_x\(t\)\),)281 780 y(and)g(compute)g
(and)f(display)h(this)f(instead.)h(Note)f(we)h(also)g(ignore)f(n\(t\))h(and)f
(assumes)281 820 y(n\(t\))h(is)g(constant)f(in)h(the)g(y)g(direction,)262
938 y(airfoil_y1\(t,)e(thick\))i(=)g(mean_y\(z_x\(t\)\))f(+)h(z_y\(t,)f
(thick\))262 978 y(airfoil_y2\(t,)f(thick\))i(=)g(mean_y\(z_x\(t\)\))f(-)h
(z_y\(t,)f(thick\))262 1017 y(airfoil_y\(t\))g(=)g(mean_y\(z_x\(t\)\))262
1056 y(airfoil_x\(t\))g(=)g(mean_x\(z_x\(t\)\))262 1096 y(set)g(nogrid)262
1135 y(set)g(nozero)262 1175 y(set)g(parametric)262 1214 y(set)g(xrange)h
([-0.1:1.1])262 1254 y(set)f(yrange)h([-0.1:.7])262 1293 y(set)f(trange)h([)g
(0.0:1.0])262 1333 y(set)f(title)h("NACA6409)f(Airfoil")262
1411 y(plot)g(airfoil_x\(t\),)g(airfoil_y\(t\))g(title)h("mean)f(line")h(w)g
(l)f(2,)h(\\)361 1451 y(airfoil_x\(t\),)f(airfoil_y1\(t,)g(thick\))h(title)f
("upper)h(surface")f(w)h(l)g(1,)f(\\)361 1490 y(airfoil_x\(t\),)g
(airfoil_y2\(t,)g(thick\))h(title)f("lower)h(surface")f(w)h(l)g(1)262
2277 y @beginspecial 50 @llx 50 @lly 410 @urx 302 @ury 2405
@rwi @setspecial
%%BeginDocument: airfoil1.eps
/gnudict 40 dict def
gnudict begin
/Color false def
/Solid false def
/gnulinewidth 5.000 def
/vshift -60 def
/dl {10 mul} def
/hpt 31.5 def
/vpt 31.5 def
/M {moveto} bind def
/L {lineto} bind def
/R {rmoveto} bind def
/V {rlineto} bind def
/vpt2 vpt 2 mul def
/hpt2 hpt 2 mul def
/Lshow { currentpoint stroke M
  0 vshift R show } def
/Rshow { currentpoint stroke M
  dup stringwidth pop neg vshift R show } def
/Cshow { currentpoint stroke M
  dup stringwidth pop -2 div vshift R show } def
/DL { Color {setrgbcolor Solid {pop []} if 0 setdash }
 {pop pop pop Solid {pop []} if 0 setdash} ifelse } def
/BL { stroke gnulinewidth 2 mul setlinewidth } def
/AL { stroke gnulinewidth 2 div setlinewidth } def
/PL { stroke gnulinewidth setlinewidth } def
/LTb { BL [] 0 0 0 DL } def
/LTa { AL [1 dl 2 dl] 0 setdash 0 0 0 setrgbcolor } def
/LT0 { PL [] 0 1 0 DL } def
/LT1 { PL [4 dl 2 dl] 0 0 1 DL } def
/LT2 { PL [2 dl 3 dl] 1 0 0 DL } def
/LT3 { PL [1 dl 1.5 dl] 1 0 1 DL } def
/LT4 { PL [5 dl 2 dl 1 dl 2 dl] 0 1 1 DL } def
/LT5 { PL [4 dl 3 dl 1 dl 3 dl] 1 1 0 DL } def
/LT6 { PL [2 dl 2 dl 2 dl 4 dl] 0 0 0 DL } def
/LT7 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 1 0.3 0 DL } def
/LT8 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 0.5 0.5 0.5 DL } def
/P { stroke [] 0 setdash
  currentlinewidth 2 div sub M
  0 currentlinewidth V stroke } def
/D { stroke [] 0 setdash 2 copy vpt add M
  hpt neg vpt neg V hpt vpt neg V
  hpt vpt V hpt neg vpt V closepath stroke
  P } def
/A { stroke [] 0 setdash vpt sub M 0 vpt2 V
  currentpoint stroke M
  hpt neg vpt neg R hpt2 0 V stroke
  } def
/B { stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
  0 vpt2 neg V hpt2 0 V 0 vpt2 V
  hpt2 neg 0 V closepath stroke
  P } def
/C { stroke [] 0 setdash exch hpt sub exch vpt add M
  hpt2 vpt2 neg V currentpoint stroke M
  hpt2 neg 0 R hpt2 vpt2 V stroke } def
/T { stroke [] 0 setdash 2 copy vpt 1.12 mul add M
  hpt neg vpt -1.62 mul V
  hpt 2 mul 0 V
  hpt neg vpt 1.62 mul V closepath stroke
  P  } def
/S { 2 copy A C} def
end
gnudict begin
gsave
50 50 translate
0.050 0.050 scale
0 setgray
/Times-Roman findfont 180 scalefont setfont
newpath
LTa
LTb
864 451 M
63 0 V
5994 0 R
-63 0 V
756 451 M
(-0.1) Rshow
864 991 M
63 0 V
5994 0 R
-63 0 V
756 991 M
(0) Rshow
864 1531 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.1) Rshow
864 2070 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.2) Rshow
864 2610 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.3) Rshow
864 3150 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.4) Rshow
864 3690 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.5) Rshow
864 4229 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.6) Rshow
864 4769 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.7) Rshow
1369 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0) Cshow
2378 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.2) Cshow
3388 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.4) Cshow
4397 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.6) Cshow
5407 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.8) Cshow
6416 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(1) Cshow
864 451 M
6057 0 V
0 4318 V
-6057 0 V
864 451 L
3892 91 M
(NACA6409 -- 9% thick, 40% max camber, 6% camber) Cshow
3892 4949 M
(NACA6409 Airfoil) Cshow
LT0
LT1
6318 4526 M
(mean line) Rshow
6426 4526 M
324 0 V
1369 991 M
1 0 V
2 1 V
3 1 V
4 1 V
5 1 V
5 2 V
6 2 V
7 2 V
8 3 V
9 2 V
9 3 V
10 4 V
12 3 V
11 3 V
13 4 V
14 4 V
14 4 V
16 5 V
16 4 V
17 5 V
17 5 V
19 5 V
19 5 V
21 6 V
21 5 V
22 6 V
23 6 V
23 5 V
25 6 V
25 6 V
27 6 V
27 6 V
28 7 V
29 6 V
30 6 V
31 6 V
32 6 V
32 6 V
34 7 V
35 6 V
35 6 V
37 6 V
37 6 V
39 6 V
40 5 V
40 6 V
42 6 V
43 5 V
44 6 V
45 5 V
46 5 V
48 5 V
48 4 V
50 5 V
51 4 V
53 4 V
53 4 V
55 3 V
57 4 V
57 3 V
59 3 V
61 2 V
62 3 V
63 1 V
65 2 V
67 1 V
68 1 V
69 1 V
72 0 V
73 0 V
74 -1 V
77 -1 V
78 -1 V
79 -2 V
82 -3 V
83 -3 V
85 -3 V
87 -4 V
89 -5 V
91 -5 V
92 -6 V
95 -7 V
96 -7 V
98 -9 V
100 -9 V
102 -9 V
104 -11 V
105 -12 V
108 -13 V
109 -14 V
111 -15 V
112 -16 V
115 -17 V
116 -19 V
117 -21 V
119 -21 V
121 -24 V
121 -25 V
LT0
6318 4346 M
(upper surface) Rshow
6426 4346 M
324 0 V
1369 991 M
0 7 V
1 8 V
2 8 V
3 7 V
4 9 V
5 8 V
5 9 V
6 8 V
7 9 V
8 9 V
9 10 V
9 9 V
10 10 V
12 9 V
11 10 V
13 10 V
14 10 V
14 11 V
16 10 V
16 10 V
17 11 V
17 10 V
19 11 V
19 10 V
21 11 V
21 10 V
22 11 V
23 11 V
23 10 V
25 11 V
25 10 V
27 10 V
27 10 V
28 11 V
29 10 V
30 9 V
31 10 V
32 9 V
32 10 V
34 9 V
35 9 V
35 8 V
37 9 V
37 8 V
39 8 V
40 7 V
40 8 V
42 6 V
43 7 V
44 6 V
45 6 V
46 6 V
48 5 V
48 5 V
50 4 V
51 4 V
53 3 V
53 4 V
55 2 V
57 2 V
57 2 V
59 1 V
61 1 V
62 0 V
63 -1 V
65 0 V
67 -2 V
68 -2 V
69 -3 V
72 -3 V
73 -4 V
74 -4 V
77 -6 V
78 -5 V
79 -7 V
82 -7 V
83 -8 V
85 -9 V
87 -10 V
89 -10 V
91 -11 V
92 -13 V
95 -13 V
96 -14 V
98 -15 V
100 -16 V
102 -18 V
104 -18 V
105 -20 V
108 -21 V
109 -22 V
111 -24 V
112 -26 V
115 -27 V
116 -28 V
117 -31 V
119 -32 V
121 -35 V
121 -36 V
6318 4166 M
(lower surface) Rshow
6426 4166 M
324 0 V
1369 991 M
0 -7 V
1 -7 V
2 -7 V
3 -6 V
4 -6 V
5 -5 V
5 -5 V
6 -5 V
7 -5 V
8 -4 V
9 -4 V
9 -3 V
10 -4 V
12 -3 V
11 -2 V
13 -3 V
14 -2 V
14 -2 V
16 -1 V
16 -1 V
17 -1 V
17 -1 V
19 0 V
19 0 V
21 0 V
21 0 V
22 1 V
23 1 V
23 1 V
25 1 V
25 2 V
27 1 V
27 2 V
28 2 V
29 3 V
30 2 V
31 3 V
32 3 V
32 3 V
34 3 V
35 4 V
35 3 V
37 4 V
37 4 V
39 4 V
40 4 V
40 4 V
42 4 V
43 4 V
44 5 V
45 4 V
46 5 V
48 4 V
48 5 V
50 4 V
51 5 V
53 5 V
53 4 V
55 5 V
57 4 V
57 5 V
59 4 V
61 5 V
62 4 V
63 5 V
65 4 V
67 4 V
68 4 V
69 4 V
72 3 V
73 4 V
74 3 V
77 3 V
78 3 V
79 3 V
82 2 V
83 2 V
85 2 V
87 1 V
89 1 V
91 1 V
92 0 V
95 0 V
96 -1 V
98 -1 V
100 -2 V
102 -2 V
104 -3 V
105 -4 V
108 -5 V
109 -5 V
111 -6 V
112 -7 V
115 -8 V
116 -9 V
117 -10 V
119 -12 V
121 -12 V
121 -14 V
stroke
grestore
end
showpage
%%EndDocument
 @endspecial 967 2574 a Fa(3)p eop
%%Page: 4 4
bop 262 347 a Fb(mm)19 b(=)h(0.0)262 386 y(pp)f(=)h(.5)262
425 y(thick)f(=)h(.12)262 465 y(set)f(title)h("NACA0012)f(Airfoil")262
504 y(set)g(xlabel)h("12\045)f(thick,)h(no)g(camber)f(--)h(classical)f(test)h
(case")262 583 y(plot)f(airfoil_x\(t\),)g(airfoil_y\(t\))g(title)h("mean)f
(line")h(w)g(l)f(2,)h(\\)361 623 y(airfoil_x\(t\),)f(airfoil_y1\(t,)g
(thick\))h(title)f("upper)h(surface")f(w)h(l)g(1,)f(\\)361
662 y(airfoil_x\(t\),)g(airfoil_y2\(t,)g(thick\))h(title)f("lower)h(surface")
f(w)h(l)g(1)262 1449 y @beginspecial 50 @llx 50 @lly 410 @urx
302 @ury 2405 @rwi @setspecial
%%BeginDocument: airfoil2.eps
/gnudict 40 dict def
gnudict begin
/Color false def
/Solid false def
/gnulinewidth 5.000 def
/vshift -60 def
/dl {10 mul} def
/hpt 31.5 def
/vpt 31.5 def
/M {moveto} bind def
/L {lineto} bind def
/R {rmoveto} bind def
/V {rlineto} bind def
/vpt2 vpt 2 mul def
/hpt2 hpt 2 mul def
/Lshow { currentpoint stroke M
  0 vshift R show } def
/Rshow { currentpoint stroke M
  dup stringwidth pop neg vshift R show } def
/Cshow { currentpoint stroke M
  dup stringwidth pop -2 div vshift R show } def
/DL { Color {setrgbcolor Solid {pop []} if 0 setdash }
 {pop pop pop Solid {pop []} if 0 setdash} ifelse } def
/BL { stroke gnulinewidth 2 mul setlinewidth } def
/AL { stroke gnulinewidth 2 div setlinewidth } def
/PL { stroke gnulinewidth setlinewidth } def
/LTb { BL [] 0 0 0 DL } def
/LTa { AL [1 dl 2 dl] 0 setdash 0 0 0 setrgbcolor } def
/LT0 { PL [] 0 1 0 DL } def
/LT1 { PL [4 dl 2 dl] 0 0 1 DL } def
/LT2 { PL [2 dl 3 dl] 1 0 0 DL } def
/LT3 { PL [1 dl 1.5 dl] 1 0 1 DL } def
/LT4 { PL [5 dl 2 dl 1 dl 2 dl] 0 1 1 DL } def
/LT5 { PL [4 dl 3 dl 1 dl 3 dl] 1 1 0 DL } def
/LT6 { PL [2 dl 2 dl 2 dl 4 dl] 0 0 0 DL } def
/LT7 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 1 0.3 0 DL } def
/LT8 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 0.5 0.5 0.5 DL } def
/P { stroke [] 0 setdash
  currentlinewidth 2 div sub M
  0 currentlinewidth V stroke } def
/D { stroke [] 0 setdash 2 copy vpt add M
  hpt neg vpt neg V hpt vpt neg V
  hpt vpt V hpt neg vpt V closepath stroke
  P } def
/A { stroke [] 0 setdash vpt sub M 0 vpt2 V
  currentpoint stroke M
  hpt neg vpt neg R hpt2 0 V stroke
  } def
/B { stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
  0 vpt2 neg V hpt2 0 V 0 vpt2 V
  hpt2 neg 0 V closepath stroke
  P } def
/C { stroke [] 0 setdash exch hpt sub exch vpt add M
  hpt2 vpt2 neg V currentpoint stroke M
  hpt2 neg 0 R hpt2 vpt2 V stroke } def
/T { stroke [] 0 setdash 2 copy vpt 1.12 mul add M
  hpt neg vpt -1.62 mul V
  hpt 2 mul 0 V
  hpt neg vpt 1.62 mul V closepath stroke
  P  } def
/S { 2 copy A C} def
end
gnudict begin
gsave
50 50 translate
0.050 0.050 scale
0 setgray
/Times-Roman findfont 180 scalefont setfont
newpath
LTa
LTb
864 451 M
63 0 V
5994 0 R
-63 0 V
756 451 M
(-0.1) Rshow
864 991 M
63 0 V
5994 0 R
-63 0 V
756 991 M
(0) Rshow
864 1531 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.1) Rshow
864 2070 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.2) Rshow
864 2610 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.3) Rshow
864 3150 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.4) Rshow
864 3690 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.5) Rshow
864 4229 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.6) Rshow
864 4769 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.7) Rshow
1369 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0) Cshow
2378 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.2) Cshow
3388 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.4) Cshow
4397 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.6) Cshow
5407 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.8) Cshow
6416 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(1) Cshow
864 451 M
6057 0 V
0 4318 V
-6057 0 V
864 451 L
3892 91 M
(12% thick, no camber -- classical test case) Cshow
3892 4949 M
(NACA0012 Airfoil) Cshow
LT0
LT1
6318 4526 M
(mean line) Rshow
6426 4526 M
324 0 V
1369 991 M
1 0 V
2 0 V
3 0 V
4 0 V
5 0 V
5 0 V
6 0 V
7 0 V
8 0 V
9 0 V
9 0 V
10 0 V
12 0 V
11 0 V
13 0 V
14 0 V
14 0 V
16 0 V
16 0 V
17 0 V
17 0 V
19 0 V
19 0 V
21 0 V
21 0 V
22 0 V
23 0 V
23 0 V
25 0 V
25 0 V
27 0 V
27 0 V
28 0 V
29 0 V
30 0 V
31 0 V
32 0 V
32 0 V
34 0 V
35 0 V
35 0 V
37 0 V
37 0 V
39 0 V
40 0 V
40 0 V
42 0 V
43 0 V
44 0 V
45 0 V
46 0 V
48 0 V
48 0 V
50 0 V
51 0 V
53 0 V
53 0 V
55 0 V
57 0 V
57 0 V
59 0 V
61 0 V
62 0 V
63 0 V
65 0 V
67 0 V
68 0 V
69 0 V
72 0 V
73 0 V
74 0 V
77 0 V
78 0 V
79 0 V
82 0 V
83 0 V
85 0 V
87 0 V
89 0 V
91 0 V
92 0 V
95 0 V
96 0 V
98 0 V
100 0 V
102 0 V
104 0 V
105 0 V
108 0 V
109 0 V
111 0 V
112 0 V
115 0 V
116 0 V
117 0 V
119 0 V
121 0 V
121 0 V
LT0
6318 4346 M
(upper surface) Rshow
6426 4346 M
324 0 V
1369 991 M
0 9 V
1 10 V
2 10 V
3 9 V
4 9 V
5 10 V
5 9 V
6 9 V
7 9 V
8 9 V
9 8 V
9 9 V
10 9 V
12 8 V
11 8 V
13 9 V
14 8 V
14 8 V
16 8 V
16 7 V
17 8 V
17 7 V
19 8 V
19 7 V
21 7 V
21 7 V
22 6 V
23 7 V
23 6 V
25 6 V
25 6 V
27 6 V
27 5 V
28 5 V
29 6 V
30 4 V
31 5 V
32 4 V
32 4 V
34 4 V
35 4 V
35 3 V
37 4 V
37 2 V
39 3 V
40 2 V
40 3 V
42 1 V
43 2 V
44 1 V
45 1 V
46 1 V
48 0 V
48 0 V
50 0 V
51 0 V
53 -1 V
53 -1 V
55 -1 V
57 -2 V
57 -2 V
59 -2 V
61 -3 V
62 -3 V
63 -3 V
65 -3 V
67 -4 V
68 -4 V
69 -4 V
72 -5 V
73 -5 V
74 -5 V
77 -5 V
78 -6 V
79 -6 V
82 -7 V
83 -6 V
85 -7 V
87 -8 V
89 -7 V
91 -8 V
92 -9 V
95 -8 V
96 -9 V
98 -9 V
100 -10 V
102 -10 V
104 -10 V
105 -11 V
108 -11 V
109 -11 V
111 -12 V
112 -13 V
115 -12 V
116 -13 V
117 -14 V
119 -14 V
121 -14 V
121 -15 V
6318 4166 M
(lower surface) Rshow
6426 4166 M
324 0 V
1369 991 M
0 -10 V
1 -10 V
2 -9 V
3 -9 V
4 -10 V
5 -9 V
5 -9 V
6 -9 V
7 -9 V
8 -9 V
9 -9 V
9 -9 V
10 -8 V
12 -9 V
11 -8 V
13 -8 V
14 -8 V
14 -8 V
16 -8 V
16 -8 V
17 -7 V
17 -8 V
19 -7 V
19 -7 V
21 -7 V
21 -7 V
22 -7 V
23 -6 V
23 -7 V
25 -6 V
25 -6 V
27 -5 V
27 -6 V
28 -5 V
29 -5 V
30 -5 V
31 -4 V
32 -5 V
32 -4 V
34 -4 V
35 -4 V
35 -3 V
37 -3 V
37 -3 V
39 -3 V
40 -2 V
40 -2 V
42 -2 V
43 -1 V
44 -2 V
45 -1 V
46 0 V
48 -1 V
48 0 V
50 0 V
51 1 V
53 0 V
53 2 V
55 1 V
57 2 V
57 1 V
59 3 V
61 2 V
62 3 V
63 3 V
65 4 V
67 3 V
68 4 V
69 5 V
72 4 V
73 5 V
74 5 V
77 6 V
78 6 V
79 6 V
82 6 V
83 7 V
85 7 V
87 7 V
89 8 V
91 8 V
92 8 V
95 9 V
96 9 V
98 9 V
100 10 V
102 10 V
104 10 V
105 11 V
108 11 V
109 11 V
111 12 V
112 12 V
115 13 V
116 13 V
117 13 V
119 14 V
121 15 V
121 15 V
stroke
grestore
end
showpage
%%EndDocument
 @endspecial 967 2574 a Fa(4)p eop
%%Page: 5 5
bop 262 347 a Fb(set)19 b(title)h("")262 386 y(set)f(xlab)h("")262
425 y(set)f(key)262 465 y(set)g(parametric)262 504 y(set)g(samples)h(100)262
544 y(set)f(isosamples)g(10)262 583 y(set)g(data)h(style)f(lines)262
623 y(set)g(function)h(style)f(lines)262 662 y(pause)g(0)40
b("Joukowski)19 b(Airfoil)g(using)h(Complex)f(Variables")262
702 y(set)g(title)h("Joukowski)f(Airfoil)g(using)h(Complex)f(Variables")g
(0,0)262 741 y(set)g(time)262 780 y(set)g(yrange)h([-.2)f(:)h(1.8])262
820 y(set)f(trange)h([0:)f(2*pi])262 859 y(set)g(xrange)h([-.6:.6])262
899 y(zeta\(t\))f(=)h(-eps)f(+)h(\(a+eps\)*exp\(t*{0,1}\))262
938 y(eta\(t\))f(=)h(zeta\(t\))f(+)h(a*a/zeta\(t\))262 978
y(eps)f(=)h(0.06)262 1017 y(a)f(=.250)262 1056 y(set)g(xlabel)h("eps)f(=)h
(0.06)g(real")262 1135 y(plot)f(real\(eta\(t\)\),imag\(eta\(t\)\))262
1922 y @beginspecial 50 @llx 50 @lly 410 @urx 302 @ury 2405
@rwi @setspecial
%%BeginDocument: airfoil3.eps
/gnudict 40 dict def
gnudict begin
/Color false def
/Solid false def
/gnulinewidth 5.000 def
/vshift -60 def
/dl {10 mul} def
/hpt 31.5 def
/vpt 31.5 def
/M {moveto} bind def
/L {lineto} bind def
/R {rmoveto} bind def
/V {rlineto} bind def
/vpt2 vpt 2 mul def
/hpt2 hpt 2 mul def
/Lshow { currentpoint stroke M
  0 vshift R show } def
/Rshow { currentpoint stroke M
  dup stringwidth pop neg vshift R show } def
/Cshow { currentpoint stroke M
  dup stringwidth pop -2 div vshift R show } def
/DL { Color {setrgbcolor Solid {pop []} if 0 setdash }
 {pop pop pop Solid {pop []} if 0 setdash} ifelse } def
/BL { stroke gnulinewidth 2 mul setlinewidth } def
/AL { stroke gnulinewidth 2 div setlinewidth } def
/PL { stroke gnulinewidth setlinewidth } def
/LTb { BL [] 0 0 0 DL } def
/LTa { AL [1 dl 2 dl] 0 setdash 0 0 0 setrgbcolor } def
/LT0 { PL [] 0 1 0 DL } def
/LT1 { PL [4 dl 2 dl] 0 0 1 DL } def
/LT2 { PL [2 dl 3 dl] 1 0 0 DL } def
/LT3 { PL [1 dl 1.5 dl] 1 0 1 DL } def
/LT4 { PL [5 dl 2 dl 1 dl 2 dl] 0 1 1 DL } def
/LT5 { PL [4 dl 3 dl 1 dl 3 dl] 1 1 0 DL } def
/LT6 { PL [2 dl 2 dl 2 dl 4 dl] 0 0 0 DL } def
/LT7 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 1 0.3 0 DL } def
/LT8 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 0.5 0.5 0.5 DL } def
/P { stroke [] 0 setdash
  currentlinewidth 2 div sub M
  0 currentlinewidth V stroke } def
/D { stroke [] 0 setdash 2 copy vpt add M
  hpt neg vpt neg V hpt vpt neg V
  hpt vpt V hpt neg vpt V closepath stroke
  P } def
/A { stroke [] 0 setdash vpt sub M 0 vpt2 V
  currentpoint stroke M
  hpt neg vpt neg R hpt2 0 V stroke
  } def
/B { stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
  0 vpt2 neg V hpt2 0 V 0 vpt2 V
  hpt2 neg 0 V closepath stroke
  P } def
/C { stroke [] 0 setdash exch hpt sub exch vpt add M
  hpt2 vpt2 neg V currentpoint stroke M
  hpt2 neg 0 R hpt2 vpt2 V stroke } def
/T { stroke [] 0 setdash 2 copy vpt 1.12 mul add M
  hpt neg vpt -1.62 mul V
  hpt 2 mul 0 V
  hpt neg vpt 1.62 mul V closepath stroke
  P  } def
/S { 2 copy A C} def
end
gnudict begin
gsave
50 50 translate
0.050 0.050 scale
0 setgray
/Times-Roman findfont 180 scalefont setfont
newpath
LTa
LTb
864 451 M
63 0 V
5994 0 R
-63 0 V
756 451 M
(-0.2) Rshow
864 883 M
63 0 V
5994 0 R
-63 0 V
756 883 M
(0) Rshow
864 1315 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.2) Rshow
864 1746 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.4) Rshow
864 2178 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.6) Rshow
864 2610 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.8) Rshow
864 3042 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1) Rshow
864 3474 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.2) Rshow
864 3905 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.4) Rshow
864 4337 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.6) Rshow
864 4769 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.8) Rshow
864 451 M
0 63 V
0 4255 R
0 -63 V
864 271 M
(-0.6) Cshow
1873 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(-0.4) Cshow
2883 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(-0.2) Cshow
3893 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0) Cshow
4902 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.2) Cshow
5912 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.4) Cshow
6921 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.6) Cshow
864 451 M
6057 0 V
0 4318 V
-6057 0 V
864 451 L
3892 91 M
(eps = 0.06 real) Cshow
3892 4949 M
(Joukowski Airfoil using Complex Variables) Cshow
0 4949 M
(Thu Sep 23 19:41:11 1993) Lshow
LT0
6318 4526 M
(real\(eta\(t\)\), imag\(eta\(t\)\)) Rshow
6426 4526 M
324 0 V
6416 883 M
-8 0 V
-23 0 V
-38 1 V
-54 2 V
-67 3 V
-81 4 V
-93 6 V
-105 7 V
-115 8 V
-125 10 V
-132 11 V
-140 13 V
-147 14 V
-151 14 V
-156 16 V
-159 16 V
-161 16 V
-164 16 V
-164 16 V
-165 15 V
-164 15 V
-164 14 V
-162 13 V
-160 12 V
-157 10 V
-155 9 V
-152 7 V
-148 5 V
-144 4 V
-139 2 V
-135 0 V
-130 -2 V
-124 -4 V
-119 -5 V
-112 -8 V
-107 -9 V
-100 -11 V
-93 -13 V
-86 -14 V
-79 -16 V
-72 -17 V
-64 -18 V
-57 -19 V
-49 -21 V
-41 -21 V
-32 -22 V
-25 -22 V
-16 -23 V
-9 -23 V
0 -23 V
9 -23 V
16 -23 V
25 -22 V
32 -22 V
41 -21 V
49 -20 V
57 -20 V
64 -18 V
72 -17 V
79 -15 V
86 -15 V
93 -12 V
100 -11 V
107 -9 V
112 -8 V
119 -6 V
124 -3 V
130 -2 V
135 0 V
139 2 V
144 3 V
148 6 V
152 7 V
155 9 V
157 10 V
160 12 V
162 12 V
164 14 V
164 15 V
165 16 V
164 16 V
164 16 V
161 16 V
159 16 V
156 15 V
151 15 V
147 14 V
140 12 V
132 12 V
125 10 V
115 8 V
105 7 V
93 6 V
81 4 V
67 3 V
54 1 V
38 1 V
23 1 V
8 0 V
stroke
grestore
end
showpage
%%EndDocument
 @endspecial 967 2574 a Fa(5)p eop
%%Page: 6 6
bop 262 347 a Fb(eps)19 b(=)h(0.06*{1,-1})262 386 y(set)f(xlabel)h("eps)f(=)h
(0.06)g(+)g(i0.06")262 465 y(plot)f(real\(eta\(t\)\),imag\(eta\(t\)\))262
1252 y @beginspecial 50 @llx 50 @lly 410 @urx 302 @ury 2405
@rwi @setspecial
%%BeginDocument: airfoil4.eps
/gnudict 40 dict def
gnudict begin
/Color false def
/Solid false def
/gnulinewidth 5.000 def
/vshift -60 def
/dl {10 mul} def
/hpt 31.5 def
/vpt 31.5 def
/M {moveto} bind def
/L {lineto} bind def
/R {rmoveto} bind def
/V {rlineto} bind def
/vpt2 vpt 2 mul def
/hpt2 hpt 2 mul def
/Lshow { currentpoint stroke M
  0 vshift R show } def
/Rshow { currentpoint stroke M
  dup stringwidth pop neg vshift R show } def
/Cshow { currentpoint stroke M
  dup stringwidth pop -2 div vshift R show } def
/DL { Color {setrgbcolor Solid {pop []} if 0 setdash }
 {pop pop pop Solid {pop []} if 0 setdash} ifelse } def
/BL { stroke gnulinewidth 2 mul setlinewidth } def
/AL { stroke gnulinewidth 2 div setlinewidth } def
/PL { stroke gnulinewidth setlinewidth } def
/LTb { BL [] 0 0 0 DL } def
/LTa { AL [1 dl 2 dl] 0 setdash 0 0 0 setrgbcolor } def
/LT0 { PL [] 0 1 0 DL } def
/LT1 { PL [4 dl 2 dl] 0 0 1 DL } def
/LT2 { PL [2 dl 3 dl] 1 0 0 DL } def
/LT3 { PL [1 dl 1.5 dl] 1 0 1 DL } def
/LT4 { PL [5 dl 2 dl 1 dl 2 dl] 0 1 1 DL } def
/LT5 { PL [4 dl 3 dl 1 dl 3 dl] 1 1 0 DL } def
/LT6 { PL [2 dl 2 dl 2 dl 4 dl] 0 0 0 DL } def
/LT7 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 1 0.3 0 DL } def
/LT8 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 0.5 0.5 0.5 DL } def
/P { stroke [] 0 setdash
  currentlinewidth 2 div sub M
  0 currentlinewidth V stroke } def
/D { stroke [] 0 setdash 2 copy vpt add M
  hpt neg vpt neg V hpt vpt neg V
  hpt vpt V hpt neg vpt V closepath stroke
  P } def
/A { stroke [] 0 setdash vpt sub M 0 vpt2 V
  currentpoint stroke M
  hpt neg vpt neg R hpt2 0 V stroke
  } def
/B { stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
  0 vpt2 neg V hpt2 0 V 0 vpt2 V
  hpt2 neg 0 V closepath stroke
  P } def
/C { stroke [] 0 setdash exch hpt sub exch vpt add M
  hpt2 vpt2 neg V currentpoint stroke M
  hpt2 neg 0 R hpt2 vpt2 V stroke } def
/T { stroke [] 0 setdash 2 copy vpt 1.12 mul add M
  hpt neg vpt -1.62 mul V
  hpt 2 mul 0 V
  hpt neg vpt 1.62 mul V closepath stroke
  P  } def
/S { 2 copy A C} def
end
gnudict begin
gsave
50 50 translate
0.050 0.050 scale
0 setgray
/Times-Roman findfont 180 scalefont setfont
newpath
LTa
LTb
864 451 M
63 0 V
5994 0 R
-63 0 V
756 451 M
(-0.2) Rshow
864 883 M
63 0 V
5994 0 R
-63 0 V
756 883 M
(0) Rshow
864 1315 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.2) Rshow
864 1746 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.4) Rshow
864 2178 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.6) Rshow
864 2610 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(0.8) Rshow
864 3042 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1) Rshow
864 3474 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.2) Rshow
864 3905 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.4) Rshow
864 4337 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.6) Rshow
864 4769 M
63 0 V
5994 0 R
-63 0 V
-6102 0 R
(1.8) Rshow
864 451 M
0 63 V
0 4255 R
0 -63 V
864 271 M
(-0.6) Cshow
1873 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(-0.4) Cshow
2883 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(-0.2) Cshow
3893 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0) Cshow
4902 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.2) Cshow
5912 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.4) Cshow
6921 451 M
0 63 V
0 4255 R
0 -63 V
0 -4435 R
(0.6) Cshow
864 451 M
6057 0 V
0 4318 V
-6057 0 V
864 451 L
3892 91 M
(eps = 0.06 + i0.06) Cshow
3892 4949 M
(Joukowski Airfoil using Complex Variables) Cshow
0 4949 M
(Thu Sep 23 19:41:11 1993) Lshow
LT0
6318 4526 M
(real\(eta\(t\)\), imag\(eta\(t\)\)) Rshow
6426 4526 M
324 0 V
6416 883 M
-7 1 V
-22 4 V
-34 7 V
-46 9 V
-58 12 V
-67 15 V
-77 17 V
-85 18 V
-92 21 V
-100 21 V
-106 23 V
-111 24 V
-117 24 V
-121 24 V
-125 24 V
-129 24 V
-132 24 V
-136 23 V
-137 21 V
-140 21 V
-142 19 V
-143 18 V
-144 16 V
-144 14 V
-145 13 V
-144 10 V
-144 9 V
-143 6 V
-142 4 V
-139 3 V
-138 -1 V
-136 -2 V
-132 -4 V
-129 -7 V
-126 -8 V
-122 -11 V
-118 -12 V
-113 -15 V
-108 -16 V
-102 -17 V
-98 -19 V
-91 -21 V
-86 -22 V
-79 -23 V
-73 -24 V
-65 -25 V
-59 -25 V
-52 -26 V
-44 -26 V
-37 -26 V
-29 -26 V
-22 -26 V
-13 -26 V
-6 -25 V
2 -24 V
10 -23 V
19 -22 V
26 -21 V
35 -19 V
42 -17 V
51 -16 V
59 -13 V
68 -12 V
75 -9 V
84 -7 V
91 -5 V
100 -2 V
108 0 V
116 3 V
125 5 V
132 7 V
141 9 V
148 12 V
156 13 V
164 15 V
171 16 V
179 17 V
184 18 V
191 17 V
195 17 V
199 16 V
202 15 V
204 13 V
202 11 V
200 8 V
196 6 V
190 2 V
180 0 V
170 -3 V
157 -5 V
143 -6 V
127 -8 V
111 -9 V
93 -9 V
75 -8 V
58 -8 V
40 -5 V
24 -4 V
7 -1 V
stroke
grestore
end
showpage
%%EndDocument
 @endspecial 967 2574 a Fa(6)p eop
%%Trailer
end
userdict /end-hook known{end-hook}if
%%EOF