(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     75054,       1963]
NotebookOptionsPosition[     70592,       1829]
NotebookOutlinePosition[     71006,       1847]
CellTagsIndexPosition[     70963,       1844]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Numerical Solution of Plane Irrotational Flow", "Title"],

Cell["\<\
APPH 4200 Physics of Fluids
Columbia University\
\>", "Subsubtitle"],

Cell[TextData[{
 "Solve Laplace's equation for the streamfunction, \[Psi][",
 StyleBox["x,y",
  FontSlant->"Italic"],
 "], for two relatively simple examples. Since ",
 StyleBox["Mathematica",
  FontSlant->"Italic"],
 " can solve thousands of simultaneous equations quickly, the only challenge \
to this problem is defining the computational grid and the boundary \
conditions on the flow. Two dimensional solutions to elliptical PDE's can be \
easily and accurately solved using today's computers."
}], "Text"],

Cell[CellGroupData[{

Cell["Define Laplace's Equation on a Grid", "Section"],

Cell[GraphicsData["PDF", "\<\
9E14ARda;S@:9LO/Sj8:=B0`86mRJPXl?2m<IFiWM6PP=R0`858_AVU/M6Eb82m6
K65dIDAUHfmTIChn2W=dLVEQK@YhW7fmBJh^?JhT>?mGLEM`d]E;:jQaE@eZ0@6l
c0CnB21APmYndH`DBIehlG21NcjSMnYLSAU5ommo_Yobil<ono^?Ookegoj_mNNo
okmooNnoc/nIgnUob/m^IoLoooa[[`jSFoh>bjjO@;m4hOohjooilkonn_klmknn
WcKf7?g?ooa[W;ioC_]CIRloHl]]aC9nVU`eIo_IlllhkO^IAO1^?nL8o^[?M?R?
_lKNL/:94oIHNX>eOfXAG<m?`b?61baG[7=nAQE;oGjjW;7>n_WT4F?;gi[`6Sok
j1EV6NeW=[UWZClK^1AmYRBf1?j7IVcF>>?/bC_<EV4?k?VhURE5lFdiHlRITXne
n6c9n6JZ3B?WHoj<>67gCh3TMXl;bT7Ah5B57DDY@5:o17oCBXbi^QRWco7ClaV[
JM[:n/V`cYnjN<6e=9J6i?WnUSBFYH3W5MKgTPM/:HH]1M]`N2;7goYID`4:Y4WA
U7]/h8kB`=[m;NgPj7V3]Kl11LbOcMZGEX7V`d^dOB6Uo`>IjcmMCQVeX/7]>[A^
:m9c`I3K[HAWoAVC>KRFEGlZ2QZYG6Q<LPK:E^XM/?NO[fUIJ]5<?1HU8^eEV[ZT
L7R33/RVDePmOT;G5]o6cmX1Zca8CeL/]eW8<>ZlgnbT[?kS[on@Ejn>lTU1o<no
IY7G//`cd4;n:E0JJVm^nA/F:D:Tk7`ogoTcRcC0<SEMdT05Mk`NPRGWb:YHZZFl
o9VO]>aaOjnSQKC@L>LWICH=bX>nRBi7W/?42]a/?YmD6Fj3Mj0djCC^8mQF9O^e
oM5L[=k@j0g<mK=i:]7ID_oB4LT=92nk/47;<`X^_iQYHRkS380UQE761I;?ShU@
>0hJ`?bT`0mN@nUgN24K`XDh6cEEdPW=GWAYDjdV;>m6]kALbfA9Ojgo54Wm@RHo
jKG@/`3HgC/jGBVHPL?BOdiRNiXlMGC7^69;;KEd1_[/9Nm?6aM84Ik5L`Vk9F>S
jmJbGbQD_0;7^QTmGA[AH5l[_A2[mS_Jn3DE3oiJ`QPjV;XRGF]OK_VK[MH73>Tn
9a/MBWl>==aF5nY14/enTePcdGiZcgSX6iT/k:<5Bfo7>ab0a_9Kk3P5B[=HmiQT
M5Fm45GOb]H2T4L_WUjUgd3IS]YA9582]AGFo/13olaJF>QUb4/^Xe5PNE??i2f^
ADYGG[O0FlXFEkC5]jFaUl05KN>0F:B[aB?ja`^VM0`]@aWeN?heC;kYdW4M[CAi
BFGhT`LF_V:6DBLb`RC;;>QFd0E=jGT2=Hbk>=d<:2flXVQ=3mjX^X]iohZDaaV;
KET7J;cGZn==h@>[GR6M2o/NjCX6[Y3Go=<NaO8PGDYec6:][?@hP`?6_f9YF?XF
QTG^QH8/PflJaQ;5gloNfQ=KobNm3:j@@@63GV0I7GM=F4Y4ViMKY4O/[8^:e<Rd
@EhaJK`E@ccaa:315]IHN^=[K;b1:m]GeljVKcjBeKeH_EdZKdm]3iPOMLWFBU2J
ajbgOEa;cANP6l0CIFS36f:HE`cf[W76HWnE<FZ[bmn^KKIcO1G;eSHk2^^[Bko0
ed8V?h3BAEIm:c2^;lg5@L78:hLVdfD28@DXN21g;2Qmbj^LFL8Rkl^7UcB`I4=6
@QWZ>ehO_>]BmHe3^c@9bGRKDnnifNO<9W<:c_[TVIbRR4E6NCA]VBO9jcjK_:<L
iaXkd<1BBMY1^4FJa44Z5S?:=j_X8koQV8UJjG3GWV]>3_895gA<>?eJd<6RRa[f
8Xb^NAPLB@GGN01c?JDZKZka[T^_/T[0F_0ZhW`eb7bEXcDNb?J7DA^GbcB5Dn4a
k`>TYF::QShAVJRKDk?0lSJ/iAScU8;gOOoioohZd^MoOoh?VH5Sd98VaC[QQ:Eo
7nHmI_Sk[ookGEQ82F6R:GNm2`^TFnX:W@TG5T<kiR;C3Tjcd?4ka<A:9R1HF=`C
]R`XN8=mM64QAH^5AIW55QJbd^7KFZA_a/9Rbe^<dNe3ke@CePJ5:ji5^T`^;:@O
@1?398;?g9fUJOPOVS4<ofJA4EBG9ZE[On3hi/<]N:Nh/92VQaURgE/G5UFG@hHi
BenH`]`C=XX9_@0G5PC]^`/;@QU0^K30VTeZdXH/^@mcMC4GJi>S^ima^ZJ]kYl<
iOfaQLFeS0nUPCk>OT/J/K00h7VEmEeUbXB5aI62aL92NP7TL7=QDILFXK`;F5SX
/Bfca2[3B;^oBk/;2l0RRdD^;>AiVkEOM66aL8Vf[k^`f8>cKIWNLV7AMAH7660N
GEQL;9EQ2h][TIjK2`^TLZ4a=Bh/9RMIl/1a5aIKRdKV@UaHX;eRHE7fm0HMT4e7
YnMn`]@F;`oV`/9P>kJ`<2cU`HF5=83IKgIBE_72CZ@B_B^^i_/VU/5Q_j2[an0W
oBejK1@P6Zk<meV0^mbI7k[MO89dNR]1FIeSbUk^nYOC4P`c:00<6[9T`=/XLhAj
0JJMMW?2PSXDd3QU;iaLK5gm5^gNbRWj3^=HoGC]QP/acG2<KYYc7KM8IC2A5IdX
n`S=`nL0=iLQ]KAdF:I/V>PD[5JAUUHeSck^bJ^0FBc[[27Y7lMe_`<F=Sdoh]=n
=dj@J]3l=9gPGba[RVJ[/F_QQ>2C6TBGTn3fXcPMjldCanFnJeVOQoi7Z[U__Kl]
2RA5TmEDDNZnS;PEnlV[>gY^2InlLoT<[4d?FkZlab=Q[:Rd:E`;fS9VG]mP?bXM
2a=YSl2:YMTRJ>;Un30lLX7J^HQ0FF36l/Tl7_fDHL/5iZ1nA]LIIn1?7m4o;BFf
kPm[NEf4HEY7jdPHTdfk_ETVF`R_akR=56>@I35gaj`f4R1aQPjOLPOfB84URkZb
L8]DCFGE5Ke3GIKYmA>@QJC_PYn05O[4;GEYl/WbV7TNW9mKcNfQnNI:@FKMj1K^
7MSZIWh6g_ZF4X7gkfRof1;LG6[JndW3`7a0GSSUF`8_M5X6^E;X70SR19d_Oe_I
1[baj:J@em=^RU1:653AZj19YHY7_eIkKQXbI3b=AmXU5fAOin@ml>2TW/gE;=:?
l;fGhA1C8SC3@B93ibKOnE1L_jVKe]W5O4KM=1T]Ng?;gk3LEfj@^VTP472=S7=L
8l^kR`TindMK4flM6`KVL1nFR_Ogd65Xf3cl3T>CZe]iSOQR7QIVaJA3>b/^ScUn
KZ=^l0RM3CAB=l`51Wj`6PY6ENY6TMADddDg6TSk]3>X=PhH][D11PhoXg7blY6j
DJ2cHBjF0C_OMZiK3[7>YK0Z2LPkTc6=4fKEfHQ<N]4j7D/:K=5h;LI6H6J1<J[:
@23gWibk0MSM=J>;]EhGYaJU7W/J6;C^V5=YYAO^6JgbK]Zc:i0RE>Y6HC=BA1HX
c3T6L`aTAiV;n]eN6U>^`F6`:gECfk299O_>1f?mi;SMRFI[fVVJiCMe<c4FHXee
f6ZGm7@H><`RCecBn]MfbecXO/ioP^LUkMbbVYiakg3aOLHbU]dBXO1_PgeaIQ9@
i_Mok3MFc1>]5[m`oISGbUmf[_kVKNi9lA2/<<mn;5eG9cNM@oZJKolWf7?Z5R/;
_l?5maWSCP@d4EjhQ:2WI2hFT?=e0igc:/f]PE]POXcPGZF0Md@?]M1K?F;>n_K:
H/k2F9G4W8G9@QIc5UHeBLbI1c>147?V`FRAaI`9]Sb9>O>dTlFL`5O<2H_:7Q?;
PBCV[>lK[iS3S2DaIlTKVlFL`9j?Jc4aIf7TC6;>`Z3lR3VBdig4W7T`ic8aal0E
LbiTebBPIc57kS=N<FMQoIO4W?FETLBLP5O<2D]SJK@[iS2=8NK8NBf;>@]CaB_V
;4aP@/bAVQl^iZbRAH<IT_fnHXi24g<Fe/e9c15LOXTiT[bIa9aiiW8a9h2:>H6_
V1<F5G>TNHdTiTQn_eO<HM]8HPkKJa9cf:2cV;>`VT]Rc/:d8l@LQbKV>5HaId6V
^6;>TmEO<`:Y^i[4W8EGm15c5/J59>H8o[:H/c0d?F;>`R1iaIb55Ok`gbF9>@]c
S11cEYFTQiRc@9:jV2>GcRCVl159c=5LV9QS`<@L@aAc5TJN9>H/[<4N<FMQZI74
W5D`fC8aal0EL`bZV;?0iB@aIh7:cV;>hTZ]Y1?FcV9>h2_V90]E5<7CaIb5QEj8
>Kak4W<FIk<QiRb@4HnH8bGJ/iPSIK=Mc35`aAb3:^H/<<e9c96:jZnH8`GC7S47
9IG5W0Mo;N4[iRc>L?j]V2=9fJnH8iJCaAcR9>HTK2a8/U3<4KaLc14`TYRc@2MO
<DM0bF:>h>lELaHHaTO<FARO@laI=[YN<BO`5G?2XV9>H1Eci8KW4G=FPcXLHXkP
V/BL15G<2H>:>G;3T/FLeDcKMc5WHChmlaWXEZjH4lS473NX=;>`=3PIUnlELlCB
/iQ3W<@L`O<ELaKd]bCV[=Jo;>I87^X[iXRUIc5W=KTbRCV1[iRC;1Ac5US_9>K`
6HnH/d1T9c4W/8XiPJnH4aHELnBI?H/i<V[^El`AblQRC/:EkN]l[iPSeGfbV<?Z
3c4Wh1Ec7T_=5jPd8gRnHXiHEQIcn<Sb2cmR3Y>Ia1c1;H/iP_/Si[0PTYPS;e[;
HXkPmHXiZl^HTlBLeDO=H/kZF75V<FL=6G6BV2=hIc5WSN_2LlFLaG^5V;=6KEW<
2Gc5W;2XV2=gj4W<hB<O<DNV?2L>mgfbV9>`RCUQdNF[G5jcV2=h_68>7Y3572HQ
a1b7Elbi1Q=cn<0TiZ`>QSF;>J_;TSV9>E8a=H/iDWGWGl@L;?h>g6I6dVoFT1i]
kGlWh:aBOPThHWT474`E/h2c@1UU0DLVfHn0/n0ZT@@LVE3m4W0VB?dTh<Sll15`
0Ul19b`Zh<Qdo15`5Udg/X33S2D1Ih63B@9>H<o7]IR0/l099`5W`K_V4G0Fb;BA
CRPQh1Rh0Xi15G0hXdl2cP9Sn0Ph2jhY?IoA/h0Cl0XhHATOB/<57:Ha19`5FST9
>0]N@5O0TAaV0FO1lnP:>3:KM077OUl1Aj49>59ZSh0Sn;N0/d3PQH0Sdol@L0:X
P1?h2SQQD@57VULFL=IGOPThK1]9`65kC@8>6g@FL1KHe2CP;4Pc8N0h=077/@Xh
^XZeg3eIiN8Mi0gO1BXVohCQT776G9YSTf=>3A:Vl_j?M0G668bYT[m5[H9mQXcc
QnB_HA_iZA[i6Io99eboL:k0>U_/cSRkV:ZP;:iF1COlUGC[i;jEXUN0JmGiRkYo
@WZAHHgc6DcZD4Zh7PCHESE3acnl=:R5Kf;9PaF=38P5m>jTU;i0djg;o6?Vg=D_
cS0Wi@NTB9aAMF:0ECnG=QM;Ye5d9G@]B^4^l7I/>kX`UVLd5M>PZVPb<@=LUZnY
2kY3BYO4NN3ifHSX5RTjOHK>9LYQe[V^Vh7AP?K:UXFNUhT0aajXfeCA3AlWOY3W
cX?EklL`kjmn?gk6E[lOM[mLoZSO3ajXDdGhbi2h1TN?L_UDQJ7^aSbX^h=QUZcj
oL@Ij]Gc;aQ;;e]BVFGJVUFiNW@DQKRYC0VW5gf4=36@[M8GL5bYWbTcVlY1h7U/
6^LFbMR^_[bV1]9e<ChML_E[h]8mHD7F`/^[BI:lD::@EU;e0O[257CPdUd_L<CX
@<4Ql86H69LoAbYhUG=KQf@;34SGeS4a@CiLJ25CQG>ARg65c7JTO_d<<?gk3m`/
2Q_;@AZ?=6Pk_nPTheZhK2lIKU@Mn:afk045/iJSRn05oj_fiic:h`jUnGmfnZ3T
E^QO:oLK>/jFZRoQaklGFf?23>BN@KUaJ25Q/FWhX4?@iZYhTFVDbaMJkE7UCP]X
GLPBWARUo?PQ[g5Tl[4hBeIR3GjWePMW@ZOCmn1:?38A:jo48iH_BcaLBbN9Ag1o
9Aj^OZo48o>hMRFNeM]94PnWkb7ab0c_Ba:?`1dBSmbV9hV7SdPBSnK291h39_4H
X/@S/oVF9Ak>mQn9Qd^<9?78d^Zha6?PBS`6EN;Q0S]9?5coI8U733];??:TVBFN
`5OR2H]a55SNG8V7joN@N7Sg9?78m7QTRFOe/Uj9AdZdIXU7bVJja6?PBS`6Ve4U
oL/BSmAZOBDNi`J^a8>BbQ;?Po52>;hBSkBjnUm8?1/cJbbbd8BXl^`RBaAFYEXf
gl;e[aQL2anA;9aQnQd2Ona]^g=Pj;<`icK;gk1lV3VIHAe<Ngj1keh<J:Vie`GD
c@^6<DOdakOJOfMA;I75gkREGeW4W2]WdO4]A/G<h^CKKAK;8ZI>IUW`eSWo:BjN
Dk=X^^8>SPLcGjfbo`<kNYSkPVkbW`KA7h`olW/H;L31KjnV[giO@eO7NmUDCOXd
_8YkjW28=K/0dcX1i5cB<_dN`k8Jacib4g]nBP[feGAA/K5XGlHhb9g7dMVC9Ne2
BcS?=`/jAIjad8E;hK3QkeUdH@fa^fZ2:UmH;TLTVi^;O[bFd[T`lOHKYiY:`4<D
EF`IPc/_n4GQBHFCKa@;2XhE<i06:f>YV`0ggM;Ln2P^MAIZ04a0>`Xg><;6MN10
NDVJMMZ9FGCGC0b]LCQje3n1IJ;fh4=>>5gAIC2C7R0`5^/5V7ei<BHEEcCEVhId
<`^e?>UR8gSolM@?^o]F4I?784OYQAC`=V@1e4VOaTBa0NW31lXdL=THnQa7l]db
bET5UQi_7Gd6YTF]LYF0:k3WZ6P7cC;MbiJMG7?WlWLnP`BH79<9bnah<MPWB=h;
N4_97oP`6E>`BhbL^n=^Lj6`b=A6YXc`HY4G2FonWj:<77k:RBRaIOK:R@>^`M/P
Ti6j3anjb4G0JFCGo@L4iEa@NB4jHP_N`9YJf_l61k[Qj0Se]3RNF35?9<`]b[Li
K5Q4H3;eLG8giL2Vo0<B1TOVAZ_VHQVGb?ZhL5hFF2MA5e>na>`8LQ=GEFAHh32:
G1Wlaem`>Tb6R^H?5cj8CPn^^PQa`jB/8NTn_@B4W1I7oo5GHj[S^4c3dCQ1bJ4n
6^H0L7IH@aLPE<Tn;ZJAkSnH0RbMGBXi:]e8:AN2H><^b6_017:3bO/5YBEe^R]M
bca`oi;[<7f01loi58>OfOAdh??olALVVKa2^S9IZ_ma;>^NWF2I;>nhH4Qk`T<g
54fI8X77GVR]QFcB`?:G_A1h=kCc0@Vad__AL66[P^lMdSQJeND==XhQLH>M1NNh
Z1a`To<T^7TkKA=VX<aghE2?><kn3N5Tl4;73aOjebB4Ek>ObED<dR8U>W7jd_H9
M@;NfiSSdldA;D0B;nOEZJToK=]8oNVL2O;:dBO74<_Dl8D@IT_b[U`3n15DdLEE
<T1o_<HGAIXXU3OD:AZ1g:UBiJM_?Q=Oh@>>hgdZ<`dbkF]dl^2;SKT]lX@15Gl_
QXZ^GH5K:TVefSBCO8TT3GSR1oW9<=;daM6RWM=2=b/54?RcaJ9K_/WkH<YIb:Pf
c@6HfP4B[lCm=LoHB6]iA_L`^@0]PK6fhSCK;I0I;1EXMA<3;2K/l4Gh^<7=WP6W
G5@<]j9RP/R=G0k1<Bb7_mIKV0S0_GIYDoiW6?0Fo25cQDi>lm=HPlFam2>5?[An
`[OIldohllU;i7RbNl85KVV/6_1FlYK:`YFKbV@j>bTAPB4KJ1`hHVV2^]Sno4Xb
IP]eda5aM6gSoo`;:b@J?Xh<M7P7iB1Yemn=WC9oEkf?M;KLZl1LJ7_BIX]MLaaN
L:A]ODNj]Tc79;Qh]U_XNHS93P[2HMMk:LCY>UchlL:6MK/Ln2_R>O1jfYZHc^8K
mK1eP0H/j9>TOF9QSKHf]:fREeT4UK^h9FmBOAQC0Rm;]E^J3Z6hN4UM37RfCMcJ
dTF7T<=jK6PPc;QTEog/i0IS/Jl3nhQn6hX`YjjXIkB0/KRYRU>e3NOka^/`<AU`
aCL>12f]jWiaVK7QN5=?G>WE<3iMB9m[W=i/PikL3HfZJTLQLmYj>2FJW=iFK:PK
E@n05T1abHR=/WJlUo9KHM3YBifFAlN;Th:;LO^??HNO/KAK7S;F=;JmXR<8m0V>
<IGnGG@6Ien8ETQgMe9aZ0H<GPIAjUDg2Q8gj`:@KOCHPCmf4CSO;59@Z7WDFXOg
e8L50f]MdkneInSHLXj]jY]KQk[fi=Ob]eRf=UNcX=VbB?h5;i/W^DEV7kPJkkCn
]W]c>/[W<YP2OT4?ahQlTLjLn5/:5?W1=1a_n/Eb1oS/gZ?h?OE5dMmj?ce;Wm>F
_JN6M5X5Q3TdGmWlNfdnCd[D;M8Xef1cdmmf_gEC<cPEaJmnh8YM71EF?Wo_HOV1
7lAac1;A8gj6hH<^;D6m=ONRH;AMk2?oZF0G3UbClA8VMV9:?h:aA7M8@Zk_m7bT
4eSAaR^@k@>gaFo^Y58dUViK1bG71lQ[afTaO[O5ZY[Z[RiX=>g2oeIDZHcYcfk]
]:0eGdQ?ONdec?9gF<Y`Rlg2bV0?<I_>2:URhjIM=d8IY0^Ye^dmHM?o`Zo7XQ<m
CKU5kiI81OSjW@dSiF9Hi[JFd6Qg0_TOOl7QIQ@>;?mDl0g^l5]BiY89j<j=b_mQ
ahNEOSNJd0bH[4:@I/623TE;3_cI1=X]Te]5L4=DIYoJ/Bo/<lD33F>6g]WEnaU[
FX=k<HY]ZOAb;DgkISPIXdO_gAjafI`^Y/Z^Pm<mPn>c3/FHk3R^aiY=F0ZK3R>[
X6h7KhlIde3QB5i_50B?8<^/@GG/8GWW4=E`<<W5MQNOiXB5;Yb8`84Nm<iS/:D<
RhZHj>bSDka[fA37<0g3CURYKLLnd`V;_VJKBaVMgb3l1dY6ic[Wdcfa0SaMl;=/
OghWW=C8d<63_S:Hk22<2FHfj2ZToNgnLLc^h2fPKmO2gn3O:_]R7/7^0VV8J57h
2hm1U7QWokF[cQHkV0gY6oX^mYMT1dhHk8YhW7<]jHWT9^?na98JFSc=C=aRGlWO
X=AefDN4EU:hX2NfoQiYAhmj/IcM^JZ=<g@^PIL6@lUV9hegC2La:5S]n37;H=n5
N]EYafB^4H>4K@F[m><@_DjSD^/FL7YbCeKRhDI]Bl;@BSJ<A2hMdNlI[FRZSThA
7J=0F:EQXM_DK]aVkVSlg1T5J`Hc/QI7b@P1HK[6Z]MH:mH5HG?Uh^hgM;hK7Xdn
^l0^EOKoMgJ122dX[W]l=l]kiic>/M:C>?lJV]8/e1Qd]3Q4W:YbBEm/iN6L>FQ3
k_PJF/;o5<^R;7`]AG_D1ibkiR6LdhkIQHk]e_>>0OAFF=<]Oh]UJZB1N`eLSTok
Eab?M4^mInPM7=/c39<b6nXJ>]F_A_88kin_Gl?O<7bj9nMCob7Y6g5Cn>:C=6k<
f^3n<RIST<c1AUamk?RFKMBmV9;AaERWLb1fRc`D[`afgSB]>2BSdV<@C2ENS9]`
;:eWP/gH[8^GKU<`Ckl=TQ=4n4NFD9l=Af8ml:TofA[a6k/kU6bi5ZcL=cW]2@H`
hFDb];eXL^_f@`7eW]67nP`/cL;5G@TRG64Fn9^0[h@024FC_XJEFm;F=6Y3aB@m
hSnhfO:Nd]@]H^]V0>P1^=NIUYZUQ`fSc3ie:;UW37<UPZ^gn_@]Tn[hF^::jfbT
dRV[]^PHQ>i9AY0?ZQCJJ56BZ8?f8LEE5b]JBaGceeKo15`LCQc3kLdBfKm/PM1E
D89PG8?gAKeUhQLHC?`UOQdklN/F8gja@J@?9GPG77NLnEe]>_Fk/5fP>_Fk6Wfe
4oF;hAO4[U>o1K`9m^@KmE_@a/maj_ORX7kMHZjF5ekZ5`Y8bLc_E5m;Igha_NH>
56=n7M^kH_Pb_i;2@GMOIGH;59OIP_T]6:g2D0HZfHUOQlKkGZ`l;YefChJMnb`_
kE/`<V1`l^>KoH3A_PFU1;kWd[j5oLiaf[MP;>1KnEWJYkX66Ne;=f8hD=kS[N/[
nD:WOJoUd[id2>oCJEln02bjdkh5:nFcWOIe[;B_@jMm[nGB__@9aG;Pd[kLnO7`
_WA<7L7k>SKN5i<c=37WOC756behGmAh9WhM>o>K;;Z5iF:ULk5ROLQOn<JfWHjg
c?ibVmW9k2ni/QG/;jJUV>MLmQL=X6Kb5jF`Ei2o:?5bW?a5Uk=W9WnaaNIQOjFV
eTTH[S^IoRd;bnGRm6mQJ9?Ym:oP2CWKjMnbh?c@WOh]R4H85dBSObl>n]L]B_mb
?dQeOYL?Y0QafEm9DT^79k^XBol6__B_Fhcn;ES5UN[d;oIGd5g1jEllP7TfnYM9
J<knG^STkcDHml/7LYndLKmLaYG<oF;KBRIoiiT?nh/1oMoC_`Fa=1kn=baLC0P/
V@4^V3Bf]2hZ6hjCn@A/9nW>03_fME5HU04^R9jR5731S;kNIE5Q>9IYAbaIbP7o
CSOFA@DaI1hFF2`edl22Ro?01Y@8EV1<L=VjRInm?j:A[4/5BhX>Qa/N`HAn>AG/
f:WP/>PVM2`8adj@3YU>1IM=Yj0hg=1H_A?2M?[[`@A;F^R]MIWP/]DC58B_Y9Z_
VS;11@5kS0ZF_=6eiE;1PINWfRa61O=Z;3N=2RiHKWhSDL5lF^:22eJXKCPGG10G
DjKScPFca_^UPRDSEN>h:1EL45<cLl5/LgDj5eb`j4YTL61SPmePM73IfQ^2mIG4
UN23bjkL7ji72_O3GSkHlNF3`j2cVX8PAVeTS4e=@@LG;;PaEOLC:RMBU`h^F6IS
H74jF=XQ_F`^7LbV67``jf:NF;5YlfcY^>jR=D8hl2F4gF:4<>^^Cb>4]NhOA[RL
/H?2DV@LUX2ABJabh4YmP/@B`e83BBb1m7bh95KPBf:iADT/gU3:ke9D1J5hY8T7
RBDFn^c56I>;Sd]RGA`TUU^<a4:><X/UMdCOVeR/PWgCIJ@ch7gJW<Ebk2aFF9C5
HSj<a^;mejFai>T<Kj971/?0FV4cHDbjdePEKPPedeQQ8HeE?j`]^WO65Jk_IJG^
^S88iTUW`<nW>8gUf;_[/2R=EK7[P[@WcTHHS67mMOdd5QF?g8@YSoDkiAC]?_BL
RLNZf0Za;heE?gA/RfAEo@[M_iG6ZXRf:NnM7X4/gYc6Z]RX369mJFh;=`R1a/8m
b6?I3b>b52VCECme^09WEK6;`KV/PSd>2^CWifbF0ZNc32ZOEN2R^if/HQkfBW@F
LgQj>V>Lc6OE3i>fU_P/:Fo^=70nRaG<_W=JmPoGiLIWGAQlUU^<cm;J3;JZdSFo
93j[H^M;[NV<bYf4Ulmbk7aFF<QA2Fb9d:YdcDj<EXEcMlnW[laYJC/HRM>Zf5h@
Y5I5A:9P]B@7c@;Lo5P3WR`d?fiaPhcELWaI[C0X>B:h95Z;4>?^Odi[e@JIJU`9
[GH=lWL]DVGX3ENLDm5l:7CoaWM=kYHikQUj1lOfS>UeaT@Ho5/QijZW><A60VAO
`LJ?`Qnh74ja]?77?A?0kV8WYFO@>oMi:Sb9bJQI<Z_6Q?PG71UeBkUWj1dLfc<<
gdA4f@9R6X9Y`XEMJF@5CMTlI_H2:jmdK:N[0>b>k8<GEh69315;cF@8LB930XM<
O2e:QU@89DJ6E2PoCXID>4dJ6E8QSP@ID/W1I3:THZkoT24EblY4Q]BYc]jG3;Th
b12g:1WRd<R@2Ro7C8ID;Wf33:V/gB130R/IL_4U@b[952Lkj]@9P9<QMF[Ljg]l
B0NH_>0^_4i`QYGMZ>?DA8IDlRZ93:V3Zn=dO6DbY>8=cfB8Y6EW<ZA2oP`bY68Q
VlV@bPfd<hkgU/P@Qdj6G</U@`AoV@cQ0ahbY5ZOM/T@adZ6>7@bi5X^6E:adTaT
B9dJnBc8T3Yemm9MQ`@^eZ@d=82C8MDRJE`bQ3FNb930U`c95Tc90Pl]Z7TN<ZAR
_AaT21mHG_R@8DaR8T<4cdb6E1AX9T=H2XT<ZF0LPPbYZ9273:W[jiT<ZN@]Tb_L
X]M0mXEK^fLbY:iC<QUB4K@STb4E7fm8I4SMR6XCI<S5bAO^F/`G3]a8T25lh4>6
B9;bhEDc6A;HON6^iO[2PMi9I8SPoI8QUNA6T25<@Y0Q5hHWW1V^8a`NV<R@2XKY
8D<ZDYo84:VFTLT@ZKSjGi0Q5C5G7c8T;3XUgNLQ@bYIR3bk9V>ACeP?6N8hI]M^
DC:Tj^:GDfS8LDj6E;87dhjLC8KlC[Mja?GffbF^?fB8h21339QG786ka@dW@b[F
dDj6B8ZL3:U@FQ<Ih]S9T;3`CJU@`8<<4ORB8AEjNd^77c:T8^Y]9T<Zb9A4QUCZ
/TZ6E10SCXID/ToV5kO7@hH4M[nhJkUnLKPjTB6E7kW8I0ROU/P@`@lID_4UT4b6
/<JM3:VPNA8ID_MnbA2f^DB6E4BTCFA8H2=3g61TB0FKHVA89C?ScW6k>AUBMgg8
4<O^7>L6LgGKmB53:ShoTlV@2S8^TB74R@bYX>lN<ZA2cTiT29]RLXi3GF@bA9]W
BlLO<RC`9D?LHV@8jlk94:gkQ`aYAD=SQh]LhiKljIIFB_T=[R^L@WRE8DY[mMmj
Bo[1e@>Zk_ZXe@?eWTDVKaB4^O7nGZXLh[Y[`K_H=L`aOm_mUVFP8P3UmA5CI4iR
[F[/[deWc15^HPf117ZoKV:=fjP9iCML8YNiRCWRd[AQokciRCE/7E7_:=eFh7iR
Q>4XYU0maNbgTRh=VjXOGk66cbVDUYbbg>:nHT`3fYCiN^6FW9^H:eS3aa4NIk46
QnfICQPe>h^eNSVJhANH9I;1OO7k/MBD4`GZ;jHP>HbeeTK@KHZ<KV]`O4idFh?_
MJKK6R8>1MgFh1^Lj;K0UficRm9]_66Rf`B_Ufi[L=]=M5][_FJjkN:PfmaRM5^c
F46GC6]`J7kX]/I81b>MDKi<]cUf^Rd/B[La7dJglOi>]cGhDa_MeQRU`NTf9^bU
faZR`3adFeQ8]cDDOJ;K6/[UXM]J7bOCK@dkOQ?MiSSlhmfRM9_PNNVf1ZK8jKKF
=MBo7[64:Mgf>nFH4DQ5m8M^JmS[iWAKjmnlM5][k?F=KV]=Xj[Y4OCF@KLeKR^n
M5]3R2BSfaZm5gNi?ha^DjAdFd<<9Z?HFU<GOZGKF^</B0nP_AWMY/3Y=X=:]lTe
<m5]c<=3]c67RFi[GKLAG[Z]MLbE<]gF/2D_dff/h4Bg=GG3_GCKQD6g^LGX=Zg=
X=_TVNFUfjAFFZKK91l[dff>WFh;2oVca^m3c0C;@kN9XBNjSC3AKMX><]dV:JZ9
KY=f=1;MeSYmcY:;N]M1:hjOC;Li_WAK69A^JkZgmM9]Q;oX=PAFTZjYhHgYN=Lk
GdihbdS2[^EGH:D2=V2Bif91dQM[H@YCZdI^:0R=MJP[bb1a0[O;NKU5bY@JhVSl
K0TMHk5QNTc^A;dH<WWAf4GgS:WmHM5=;0jah1bZg:/5@BJ6NHO82eKAmgM6Ua]P
OORi6g<=@6_3BYjBbD26S^daAWenjf[;Q4Fo<TPo41W0:n=n<J82OS?1S6N_AlYV
ge0>GA[YSO5QNdGGCK3dUL6Lh><:_I2`n8?@<]@1hNPRbjf;cGD6[h6O`BnPHASZ
>i0nSEh<=5AeY2g65B6N=Yb3IZnfHDI?QS>jk[QAiA^O6=TI<bk0aB`@O/G6caQ;
Wjin=YRlfeh[W0^8J4a;PaQ:NQ42V?[V[iI5<G[XBlTH4?Q/9:>[HBUC]6lXF<JR
i[Vb:H7WOJoMPPSh;7ech^bY/F4kL?O6I^^Q>4?7=307SAVdd=EXl4NcIIJncL>/
<oXD9X2WJI2hdVl^n>k1:@?NJ]Ho`1U<JCJ41E5im^<WRQ1ZRgi^e5PMDQKMJ4Ia
PQCLBQ0Q;1YObn]SX=n1@la?NZo/ZTUTXbChYTIF<_@a5Q;4e4?93COm^>]JiAal
[G1D>h;08EBBe<GYHTcaC=lfRgU;bWd?MU?Snf>F`ND0YgoLonn7ej6cIND72QVU
m31ohGbc=TT0aWO]jSX3n/i_07Xh?h7K^D^L05j3fE7oK/NkNE/bbkXNBHNRgX=g
7>LE3;_/Ija?7ECQKXIJGWBUiS=d]6BZl1YDTnDGe=lAm@YWo]5c@eQS_VM<jZAE
G_Pe0WiC_gIa3@PK/IK6W]d<j6a9e0MP?kQi3S2V;]YohAMIEn<TPnl7XmdFO^kc
H/]3ZnV<CoWgRl5M8Kb1>^]Pdm/O3<`Oo5_1JGl8i2=i:]<a=^kfJRFTUZ:>j[`1
oIG=MhIUg1fS5_KAcKOU^_=XZ=W2B=0G7_gV[5i02d:aCDf`QYBNWVG4@;V@kYIl
5N:4DOAK5nX=2JV24IRJ>D6PcQ3>ik_W@edkjGY`I3<oPLa0BP9N?XAbH>3PRi[^
j[MgTj4;e2]fM`fc79RAXPgbo<;81G62AVJegDelFe7ikM?_G30mg81_WnVE[:eL
iG94G/KL:=InV`dHKP@hT_lCaXi]gFKV5RjXd3BFOSZ^XRMPC6@=_X?0hN?l2a7i
OHdkPbLWj]9EZhQ^U[mQ/EN]<c[@>J[Fh:>BS<>`nMdimXXFM`5O0D7?R>h@nhG_
ChKDZIlj8a@600GjFed[2[]n5R<n6dHAl8OG[Cl<5:@^@kSkK=a1RHj=bEo;i>6_
LO_XoUBkof`7;lLnJA=HV@ikohofo1LcB@`jiVO/DkTY^O@;S7jV^aD83di82_I]
W<6IT=hIKP4>NF<6dH`C2VNQmFR@GX?H?gVUOK=h[8_2Sc0e3457`mPfWKGPeQdC
TJVESJm`PBcYme7jaIb;cBMU]gC6>[[GVYgjd_fgAcdk0K6mFA<a6Al91@mgGSST
H78h_4l6=Hc@AdNg2?g=D]9iAm59K/;ZXfAh=doMdHlRVnEGk:>f/3OcfiL@faY@
i5YTdH@k[SQ7<>G/Ol75ed3GHYkmOPO7mPcO2j29<?Rg@PBYEL7Kh=H=i`ABJ4dg
>7=E^[3:X`doo4b2IY=o`_@<j4cWNBYFlM;5N3:W1/[j5a`I_IIK5?L>S^dIgoJ<
8Q5A]X2HhZ_PKK3ZUVb0NXc8XR_5LIX;QI2>kG@E`G;1^f<UlPSN_NXf[b]h4bO1
>o3UO/>RPWM7h0@C_3/o^gX5kfj13OD8>U@G_3^j/dO`kPP0nPSNWD5V@o3^H5ZC
h7ea2=i^DL7KX@WN/_HMSn3M<GU>PWM7G8`TN0MF`O_R:gQgO?8X1>n>3dMU`K/S
o>I<afGl3L7KX@WN5j^2;FPT`K_c0fA9l>j<i7WR>5bd@o3^WgfHfH>n`=DU2MhM
kUhQN7M41lj2MhNGE`SN0VLB_1fjh7d]Eo3^QM_CGO3V0ak1^o=[EB5h>eK1fj4;
g]Mb1Nn>B2U9l>j82OL8gV8YFO0>G:a9]OT8gQdOaTZ2=f/l2Mj1[n2M;EP111iJ
D:dnPWO7G2H4Kcj`_?0A_9W497Qgc3>Bh2eInQk1VjF@1>o>koFih=dAFoDA_3^F
ZDW`kXQ:TPC_3[kR4K`kY`/QN7Ml8R49gYeOT4j2MnLGYD;`k]D249_POG48gViA
`K/SN6PXfWcP8gQ;TTHjg4HF_0=O`M/]=^?XF7hU`K_3NnTA_3^o3ab2=i<@P_N5
;WQOP`WNO60B_3/R/3j2MnLX7X9gibNVIl9So1N2MnLWL;?P7AJB`[gY1nh^OMga
dM//N7Ml/V3T4e[;P[MSikO3XX9gicB99;KL[8GPgC7W<l7KTjF2mnmdPg3YV0Pn
PWMWD;/@_0Eg5k`=Z>2]`0C_gWBAa]hOe9<;gYgOBVQfA3oiNPE_abihQhDBM/N7
64;`k_PPAQJlaC13l2I:P[OPl@SNDUHW2mhMLFM<l>hPb5c`kXR=HX9gaj@b2Mj1
UjOJ;2IhlnXTN7MlQohA_?Vd97QgQ[@9`K_c6ga9l6J=^n0]6AUIl1Hl7l6KKBh9
gX9g5[`3Vn3]1Q>lNk>?DUNDJndQN?NZgm_T4@C2Bh:ghb]hQj5K5gM:5[`k0YMV
`E/<8`_Na4W`kYQo?H9gAgRF97Rc:HKPcK[8P[LfcaK7fiL5kl1Gl7J;2Mj/>aNl
]NkON32BY90S5ITLfK4Q=lVAdPZoAhk/n2A2b94M/mHTA`JnLZAKE8kT3I<LfJ53
?79Tad`jbI6B_YGUb8]3SWB;bI78DIHS>cBaAhk/f56Ki<R>B;I9SWC/LVAHE8iT
?Tb>i?eMS^b@DTb>U=lWbI5<f2]7M_CRSa`I5/ZAWNY3b962fb]7mS56UR<knY<T
AcZ>k]X]bZKdXMnTHJllaPTi/V?/=cWB4jIbi>nD/klNmG_Tb0i=gNG83[_9TAge
jG9TQiIXLZCl;UV>k7g]T2>ka/jQ78UkD8jd7bI7:U8iDWkG:d5f2ZQGS^a@F^l1
=52C8aFh76U@iLR>IQmb9??`b97<HI8Sne19j<ZA7N?F8dMfO;l]bi6/h2A7M_Ba
8DMN67:TFdb>e=X<>K:SECmbY:CbI3UBVWc;LZASUb?3d[DQfhI4Qn^A8c_f]8LL
BISTB6d7FHk/n7A;b95me9;Tb0i5;/^AWImKBLOkb7:Thb]7QT79oXhG?nA8`Umb
Y;BnVUb2>_Qm^:ihe:`>BS4<OFWLl7o5ee?8;I0cLK65[O9KJm0/?UNSI^5W2Y^U
d1XF`@f=IIAc2Ybe_RlOoR9dUP:?W@EhPgAEgBe/28]^CB0FeV>o_cefEUPHlH[O
Km3OM[_U2K;PFORIXVLY]?1I13Lj5WSj_G?l[6EKOo`<`aHobf74ci93GcQ7:C;W
Z4i6RPiA7FBD>dMek<8<ijR>cI7V7=GQP^W>DHkdcLE6>g>>jQ]k:^cWoY9c564h
AbUDibSkKF>Bk@8<ijP>e_aaSW:;>dLa3LTi2[O<cU6MoYgI>J[3JC>LXoZY9c]7
bH_OGnLX]d@bl3fTacV:iN`iDJ3>D@ZBLeAR7@Nn^?a4SaT<m1OAHhH/27j3^mm3
XDEjl@/3mcMjc00?lDB?6OSlO8XN<l2<Y>PaS^>AK[4DgC/4g[nSaha1cBBRa`c@
5SUjc<2P6m5SQXd@=g[</14jX/L<O=4ZAHlI@knV7[R>=g[<j<XBAOBH0MkRRAhc
/6a9dF=^`RejS4>;7^=HXlN<XMmhAi2H`NPe=g[<P?^DAXlI2>]_dF?l]dN?LH]A
[G9QbM5SA_N?V=S09YJMXlL<3=`YNXaSSakS5X/N<o2U2H/N`nMim9P1WaF;7X?O
mlN<j35b`DSAH`JnXIbRa`a`>bE7SaUh:e;dV37/^`dF?LJaAhmaRdF?HMFVj343
<oTWN/b0AeB:7W=;eJ;71=CX<HhmNPbKhjoX<F?YMlj2?ai8Ln:?RA=o7?P^YL:R
o?70Ic2=?aj<TG3ih`7]g?SSPLm[1Wll447^hHl7B=Z7?ajLXPIo?:Ink_CbaaL7
On`FiHlM6Wll64LSlLL3nfLBOca`:?77PNeM<Gcih`6f=_SSPLe_VClNH8]V>XiU
/O?73Xdo_UP9HD4[lLL3nidbOcch^H0Cad_=o?70Y2[ca`<^>HTo7Pbmh_cah<L?
4WllY^j;Rn<klLL>WCnnU//OSfVKfHdoiP<NoWP`^<m>@L<=gj3Q1U?@L;ELoWQP
0ecRS`LL9aknF2`]llN1RcDYoO2:llL3O63RSeWSRCl>O?WSK66@G<M3dkKk`al?
=;KPSoW0l/:7?fHB4glld3PCOccPCI;iHiI2hXo7DQLPhhl7IfNI?ah8Fi;hh`4G
Ul@O3fb?N_SS/O@;SiLo7_2bBOcaH?S0a1l?MQo17d/=U/`OGacl/E^D?ak/Si`P
iP<OoURB]>;`fR_cah4_On`Fhho7>UoVSlLjkNF?lH3<7c<9`AmOj?ca=AQoc0LV
oWQ`QeKVS`M6d<@O3kaZRClNo4[4_nF?YELl;glL5]d4LkjInN>1JGCVSlMF_j@h
@EYZhXlM>b4A5^F?1FoSSnEV8oSS0AL0hhlmFLXOodhgn8R1X2D?Oc`@:R?aah:W
llL6U3mFH?caf?XM3oKnR1OSo;4l]UonN9`b<Wo/f?WS/>SF`U=ZhXl7W0`bOb`g
gl4O4bGnN230JNJ?1b9Q9_ih<;2Ul/L32`7WS`NfDQUo?;1g:?77PBmok1KSSgUe
hXl7PW@no36OU_SS`C0h`Al?K>g:o35[g?WSPCE9hXl5khLoIY];o?70A[C47`Lf
o]P=aQl?[8><?aiK?h^So?7PQZAQAm[9o;7Sban7XMo0ebgca`>>CITo7XPTV_QS
h/@O3cQ:?Oca@32Ga1nc:@Ioc;[8o;4fcaK7MlglLN3;7k_5n6?FWO?7F_L?Obb5
gX8oEVCllL07l19o?>S;T_SSLOB;>/HO2bbI?`ilnF>g:7o<6bKnN1`=T1CllL17
K1=o;<^T;o?75`Mok9KVhNUKiXmUTEAOoURNdC=o?13F9O77SYdo3X_balb7lLNl
_o?7P_KUSbDU=O77C=S;7dnh<SclLES87lm2cdg_S2No4I;ih`TG]l@OCgX91Wo/
f;_[/5R/LNbG<oih5ZGF]KnNYJc;7g_2U3onWG;dei=A>Q9o?14]aOUS6N?jiHmU
4W62?ijOkW:d8bgcaa<O9GCnN7hV9bSKN<POf`oSSaDYOccaFE_SS2N2USQo;:Sl
g0=XPLHO:g3nf:3baa=QDX8oIQhNoYPiC?baC=IkiXlWX_4lo?74=f0bOl`:C_ca
;6Af;WmlHO37KS7nF6/cn>>9n3@?Oc`A=3Oaam9cS<`O>gKn>2cM^^FMn6>1in6?
9cePd^VSI_iHfd7VSbNlGX8oW_cH/?NATX?ml<O`4Dgl/L2EnF?7Ucl>Pi8Slf?T
X</O4okkj36cYG0WdZJkKUR?22LCLo8D?FKRlf<YNTcP^bIgRhEMlC/4[Vod64f4
ag4QS>Pa2RdVS04/Z?1ULRePi5=]o77?Y9?`gDU3V9h1Mm4WN/c4YeiC>RJn_Y6R
acR>S;[5R^;N8K0n`j?7J2:RK1]iU1/mAZ75PE5P4F:HfHPN`d:8HjAckUD0:G[<
G>/G6B:FQ``QCVA8h3/_28^B8A?CNR=39^9=>QTRb<V@BCm]9d<VG>TO<VARGOf@
8A?MO2939[dCP`bi><P@]bPIh]38T<V5Ib939QJJR@bIF;TW<RB`TR4GGc94kYo9
T3WWBhK<=kK^a5/KI8Q38d<^EWISL]VDhD^6C706@HH@9S94l4^6C:E;h`]Z/fLb
I<kiTR4CC4>@8G?>C8HhM3;TFRhI<^Mlb10nh253iXdgHfB8HbE373XILRfG39TX
^TB6C2a]7c9THW6Kb930aI[DNLV@b@dF@HJ`aQ<I4_RB8MW2kgThe[07lkaTb4CP
fR13n<3b`XL<HA8C6C;9K@@I<T4PI3:4YI3849VeI39TH]?9@hI8:CiTb0B7VLP@
JNBob11iPahbI>;[3HT<VOPV@bI39Z8e939Tk_J@8AL76N8F9D<V@UL66L87?VB8
92W8T<Th_4661;iTR5^<397;7c9THUO6@hKP0IT<HA:23;W@bI1[<3:43daTb5ck
5aTRKnY3Q/cE7S9ThQ>Woih<VMPPl90QHM4Y:KIG939TTXG8/f/b5^T4lQU1QSR>
fKEKU0bI^_SU51Xk7Y`<VF@?YQdYV@ciWFi>[Q6di253iYh?6B8hb103BXHX<39T
H[e/I<S4>][947V/Tb7ce8L<LNaTB5Q8K`S<I8S0U`bAVbLbQ2RA8A=K=c8I<_VI
UB139QP08d<VR14W@bKIYfE7mT>61;iTR5^<3>7ER@bIH38N<XA?BfC8g>LQ@bKR
i6@bQ3G^I<P4cI?8T0TJ:9<QK7>9396L?fA8H2=3g61Tb0BKHVC893=cbI29F396
QTa8bhT<LGc9T33X[6I2=dmTb8BRWLV@2C8^TB74R@bIX>lN<V@RV4/R@mPDP`aQ
GF@bA9]WRn?W8D<2Gc;4;DJ6/>jL3=6jOlP@ZL0JI8PR8d=TS3VI39UHifLbA0`U
TB4Cgk1>I4SPBhJhALT@gS2A8H;kBhJ8IFDbI29lJB93;Phba2e6QR17V@bAC9FG
395W]4b6c6<ACXd<LNaTB5RD364nS0cQoId<VMSdIVC8??b4ebe/9^`U@aJLiQlb
92`T@aIf92HbI74SB2I35PR;A8H/dQU1QSSfkSX/BXH/oBP[NnD5l/G9445>QWS2
U0ciWG;de`]AOC<I/[jEb925X:i6QR`/X9d<4NATR?anb131R@aIS=J[I0S^@C;4
OQPIXTS94?W]I<P2Mn5TR:1;Q/Q/8lP@1Dj665@bI85J2C:4NGS84>H`TB7[>`lI
8Q>?GfC8`^M]<aW22TiTb2XeTb4G1QWR5R=3]3J335WHf_>@8H^QHf/jHcaTR6<W
@l;2cTAP9T<F0PIW<VAQ1dm?Yln73=5fT<T@BE4V@j@MIC9T`BLUTb4;O4n@8@]d
H29375lb90a:QRa/n0TbQ?0G6K;P3oj4f5[`idhQ]UHKhcNhSW<:9@N;o_7UoQhA
HT_VWce2K2elTMa2K<U/]Ml@FokK@fb51EF]je3kgG>8[MGaoLC[AJK8_<P43?<R
Fg0oLBnbQLo<QaOIjP`E@<naABoZjdGVB9//o9o=Rfc16J?HCha2kTE665iT2]F;
c7kO[o?1RBQkTBgj?VL_<[Nh5aWCT;c8L<_/AKKPj?9hTBflY3>ML4Kf8U_`07Zl
b=bBTW7:jdG6LXjLc>UNI0ZB5aVRagB46ZliN/bRWnSnMm5S1]kI9g[<`8MgD_BH
0M8WAHmak=5Sg6:kX@Ol]E?dV86niHTN<a119DF?6E>o[FGAHaajm9Q[/NPa0f_n
53eVP7=hX/L<;1TmNP`l>6kd6?j>j34:=G[<<7h7<F;X8>7AHn@1oDJ?XF]:RQhc
j0bGX/N<6e5>XlN<bD28=cH<7FFNj35dE4SAHnQIhM5S7=gX<MMPdF>^VhE5Sd59
YNPa0i6BDo@H9Si5Sf5=?=5Sj6RCX/O0ULJSaa14m1R5VQmjbJCX<KmK5/EX<6Hi
NX`HCXhN@lnF53g6/DN?LH]5SaV8[9bRak2a?M5S1VRM53eV38eMMJ?7G1cAHmaR
dF<60^=4m1S=aA<mQ_9_hgOB1od5Z41faPaAeOARKN0>nfgPem9dTT;=4H4>6dU^
NEgYEF`@idm6hXgSWhI^/G4B`B`j0hiH>h259Jo[?:a;66bVjiM<XOSRPVmIW9B1
UDa[o2`k>N:]Ni8A7dOUjlkhF^PFm8e3n0Uld@J_:5c>]WkH[W`j24XI5[@CGNS@
/HBQVRjVSa^SC<DIi<Cod<TC?9;QH]l6iaEZ/@PB0jcl@H0PgAbkm:]QS>FR7TI`
=4N04OZQ7oEEHK1lkQ0>N<8mAbeXc?X4<[>eZAbiU9TeR?>EVOGSY5j;;2QWQ/gL
e0b[G3fF<[<1^O@cb7^CVOGSb/`bLTgWX=del/Pd=kFUc2a2?m9]D8UI=35=>YTY
Pba<4[=aG<>=oX;UcVC3@Y4LEC7h5@Y]DG1OI0b?KWh2R6:303_3W:gPPHZj@dQY
_6@GBjjgGD03?RO4T9Ul0/N<aYY0XCZTUa3SCoP9T1ZT_EFT1nnL59NlDG3a:]/l
ZCIN6WC?4ci`^o13UacHnL33JBRR>=TEo78GaQVk0[87P/YhRaf8F><HEjRJ6FLd
bMhO2=0OVgA36^^9lnWWjQIj<YD<^CLGl^kiYSj0RdE/g]HVQTlToL6FGIa`8Ci_
[Zk7JQSlfSM[0EoBoC>;?j1R9lK5eXZXPMlcX:J3^X6>dacSLmGJB8/6Yi:RHgZJ
k^Gk@HbIITFd7;=@jD3/Ib1PkYHKKRef>YC899?OmKg:H]2CGoeBC9V9KjPl<FG4
lZFH<Q<OETTaIBJnQo;4U15;_C5Uj6h`k>Lj4E=V779O5U=6eUIORRU3oh/KDdI^
dE=<6MhmH/YXlRfV3474U55X<FD6@^>TV39d2WURbTRB_QaC1]iP7U=6PLND8KBH
<P=4JHXY<dQ2YYPbPh5NBSZQLH5X<FD^S9PbK[V>d4OSW26>S=c/RiPb_7F::C=>
7cVVc<2gPijH<_@7BC5UY5bVaiAAh35U26m<6O[UA4`INY`l<FGXN_K4U:4oGHXY
Tc4mUbkNhKa@jWlADfISW?Vhb`FSkco3P<QOVen_i_3BlEUR19k3eS6>aa_9j_dN
ZIR?@>4nAB4I[J?3b/Hb3]/E9kFOGM1aLPlJi=J:^/A_W5Rik=4SWhIM@?AYJC7k
fga_>/8RLl6e4L83hB^ThlGg^<VaJ?Xf_XF3gH5<jh8g4EPo[9gYiJ><7f[5LKnB
RE^PU@fN<Gnh0O[3Kof1deCfE[=jhl6_7mBA`fecI;5`dnEDkFccjdWl@;LUUCnA
9@B8f7IPLZS7AHfk=S@Tk:RfTfXS30RfMGE>/SIS/VSb]5P]YK//1/kP0EK`a8i4
E/F5fP3hfO?K13k^MDY]XW2^j4f2dGFREBRdQ^7AMj9]H7O4CXe3K]nnLU^73<B=
GM]6J0G4[S;oPm]0V[]=g4HBIfQ3HEaJO_UL6dXkMe[V3@G/2kd2[:D`FTiZ;HIC
Rc7;KA>=//4>?>iR`RghB=8ZgVZJO@XIMk7OdGKlX;HGa3e0[a?hKYZcm]<da?a]
?^d<JVcFQ2j<IQ@WJ5?RmLNKDP=W_DmZBVfX_lA]BiI0JdjNoXgHi5QAfV5[AFgY
l]MK5IY@dK1[/WiHnY9]S8/:^MegHlRW<Za_blHLI2XKHZUKAMFH@Gic@dgUYP;l
G[J:aRRh^M>Qfi7:ShQR:cHl=JInLh`Z4U[jaTJ[YG[c`1OKVAE;dTFFGYj]U];d
jo93mmLRW3:jh/d]6iEEY_;dITCIPE3Vb=g6ihgUOF[l3NZG?g7NHJa_7USZVhM[
N5O>MM0@3_EM50D:2cG9ODaJS>SSkfm?/Kc9O0aFiQ<USQdckA3]ZB6I4C:2kF=F
[WCJ?4JQKg:Wf^0QH5f<kgDoV=m=bUL`4W<[L@HTPT8<E`<:n:@U=RPAe3<:1_FU
KZP3GKVU7;mahVJ?`R<@/ZLjnj0J44Nb5bFajVdZO2Z9@HNCgcYc6:VnUT:fhT;d
7FQk_3mKmZM;jde?oL890QJ?83/JH^P[bIJ:g;G[oe?noOl@A[`]IFiTLgAbIF5]
2VE^I6mRJPXf830PKf9Z2S4g<CDe2VE^I6mRJPXd830PKf9Z2S`l;eAiL6D_D65W
IBm=IFAYHD9_N21K<20`83Pg=R0a>3=M2Rm@HG9UKW@P<b0`858:;e9ULfmeLV=U
Lc`l;e1bKf=CIGAK;e14AUd:;dEhM4MCM65dIB0h830PDPXn?PX_@fm^M6E^M7<P
=B0`858:?Sh:IFiTKf9Z2S<P<21_HVX:?3`P;eAiL6DP;e1QIfEc82m;JFAc85/:
=20`858:GB0_@fmeKW@P<@Xn?PYUKVA_HVX:<B0`86mRJPXl?2mDNG1U82m3HGAQ
K6mW82m@HFMULb0c830PDPX_CFEdHFAQM64P>B0`858:?Sh:IFiTKf9Z2SLP<21_
HVX:?3`_E7U`IBm5N7A7DgAQM6D:;dm@CB0a?SiUKVA_HVX:>20`86mRJPXl?2mB
=`Xg830PDShn2VE^I6mRJPXi830PKf9Z2S`l;eAiL6D_CFEdHFAQM64:;e=eHWAi
L6D_F4e<;daUKVMdJ20a<c0`?SicM79UHFd:?3mhL65SJfEd869UIfU^?BO_^klW
86UT?BMG=Dd`CG13IFQYB7YbIE=jCUASNV]S>F@W?ch:?3mQI6mRIBehHG0]IVU/
M6EbLb1ULf<m8T=BC4HR?ch:?7PjN6e`KFEdHB1hKFa^LcYh?BMQI6mRICY^LcY]
IGAQ;bLPN3YhKG1dJcdWF4e@87A_Kfa[JG@P<Rhi;S4]<C</86IbHFeUMfmbJb0a
;SHW?PXlLVAV>U94AR1hKFa^LcYbI6Hm9fQdM70j;bmgMgL^Mc<^Kg9W;c4i>CT_
<38_<S8]LVAV;G=iKWAQN2e^Lb<W87Q]K6ic>VUH?BMXM7A`>Rl_KW<^HFA_HVD^
Hfm];fUH;c4^<2lW?PXlLVAV>TAULf=bJG1dJFm^879TISYQHVmeM3dW=f<`IC0c
<c@]IS0aHbda<FDi;C0`<30]HSTe=68`IV=RIV4i9b1hKFa^LcY`I6Hm9fQdM70j
;bm^LbiQI6mRIBiSKfd_L6AV;c4^<blW871TISY@LVmTMF=ULSdWAe1<84MXKg=d
Lf=bJG1d83P^=S<W;ch:?79TISY4IG=SLVU`M6U_KR1bI6HjHF9_MG@m9cMS<6D`
<c<d;FH`<F<]<C5U>Bd`<30`;F8i=CAR<6ISHVIQ>BLPN6e/KW<jN65`?BMXM7A`
>Rl_KW<^HFA_HVD^Hfm];gQQL2la;S0_9b1hHG0jCFmTJFIiA65dICdW<S0`>Bda
<2da<e@`>CXe=3Xc=2d`=3X`<2LPN65`>T=bIF5dIDAQM6Dm9c8`<3T]<C0]<C=D
<3Tj=C@j<c@]<3@j<30W?SahHG0j@g9UHGA_LUA_Kf`nAe1<84MXKg=dLf=bJG1d
83Pf<b1L:71cMg9YM6EL:C`_N65`>T=bIF5dKg9DKfm/?S`_LVAV>TAULf=bJG1d
JFm^?PXlLVAV>TAULf=bJG1dJFm^879TISYQHVmeM3dW=f<`IC0c<c@]IS0aHbda
<FDi;C0`<30]HSTe=68`IV=RIV4i9b1hKFa^LcYhHG1=CCdWJ7AdL3X_;fic;V5T
Kf9U;V=_KBmhHG0_<Bh`;fe];bLPN65`CDdjA6mSMFeUKWA9A3dW=f<`IC0c<c@]
IS0aHbda<FDi;C0`<30]HSTe=68`IV=RIV4i9bln2SabI6HjA6EcHg9YL7AYKfhP
LVAV>V5RKgEd?BLgHc1U<3<c=2eV<35S;C4aICT]<30`<2eR>CDdHS1VHf9VHCTW
87Q]K6ic>VAS?BMXM7A`>Rl_L7EbK2i_LVL_I6<_IFaUKFE^M7<_<Bha;bLPI6<j
IVmbKF5d?BMQL71/JF=QM6U_KRm`I6HW?SaTHcYdJGA/IChlLVAV>T5/M3hlLVAV
>VaY87Q]K3Y/HFiW?BMh;FAUIV5eK7@W?UE^M6UdK6ET?2mbI6HjK6Tn?2mbI6Hj
@Fad?S`_I6<jM6UdK6Dn?2mbI6HjA6EcHg9YL7AYKfhn2S`_LVAV>U94ASh:?2mh
>WQ]L6eUM64n2R0P820P820P820P820P820P820P820P820P820P820P820P820P
820P820P820P820P820P820P820P820P820P820P820P80XP820P820P820P820P
820P820P820P820P820P820P820P820P820P820P820P820P820P820P820P820P
820P820P820P820:?3mhL65SJfEd86E^I3dWMbLo?PYUKVAcM79UHFd:IFiTKf9Z
2S8P<21_HVX:?3`_D79_I7ESIG8XAe1<84MXKg=dLf=bJG1d83P^=S<Y2Rm3LVEQ
M6U_KTAQM6DXA3Xb<30i<C0a<c0i=C@c=2d`=2L`<2LY2Rm=KfA4HGAU:4@j<S0`
>C4`<C<`>CDd<c@]<3@W<30W:@X_@g9UHGA_LRQ7D4`PAfQ_LgAcHg9YL7@P>3Hc
85`XL7=gLVUdIE`Y:ChnIFiTKf9Z2VmRJPXl?0X_AFiSKfAYKVLP;deQHe9_KF5^
AFiSKfAYKVL:;e1bIF5]HVaU82Q5Df5^KVm`AE=cK65cJ6A_HgE]IFidHfaQLg=K
<C1`M5e5DfaUIWARLV5SJf5bM6USK6E5Dg9YIfQdHW9QHf/:AE=cK65cJ7EcIG1Q
Hf]QIfEKMG=UKV5]IG=MAE=/IFIdHW9QHf]SKfa_LTECLVUWJ7ARLV5SJb0UMG=U
I21VKg8PIVm^M21SKfa_LPY5Dg=/HG=XMG=UL65SJf5WIDECK6EVM69bHF=[HFec
LgU]HTECLVUWJ7ARLV5SJb0UKF5dJ7<:AE=cK65cJ7EcIG1QHf]QIfE5DfaUIWAR
LV5SJf5]LfeQM6Q5Dg9YIfQdHW9QHf/P9FeQM6Qc2TECLfaQLfQeLfE`HF=[HFMU
FgEdISQMAE=/IFIdHW9QHf]YKW1eM6E^HdECLVUWJ7ARLV5SJb0UMG=UIWE/87A_
87AiL6DPI6UbIF=dK7TPI6UQHg9YM6US86=XHG9QHgAULW<:AE=cK65cJ6iUMf=_
KFeQKVA5DfaUIWARLV5SJdECLfaQLfQWLV5TAE=bJFMXM69bHF=[AE=/IFIdHW9Q
Hf]5Dg=/HG=XIFicMG9UKF5dJ4ECK6EVM69bHF=[AE=cK65cJ6iQHVaQAE=bJFMX
M69bHF=[AE=bJFMXM69bHF=[2TECLfaQLfQ^IGMSKfe]HFiTAE=/IFIdHW9QHf]5
Dg=/HG=XHg9_Lg=5Dg9YIfQdHW9QHf]5DfaUIWARLV5SJdECLfaQLfQUKW=eLVE]
HGAXAE=/IFIdHW9QHf]5Dg=/HG=XM6U]IG=5Dg9YIfQdHW9QHf]5Dg9YIfQdHW9Q
Hf/:AE=cK65cJ6iUMf=_KFeQKVA5DfaUIWARLV5SJdECLfaQLfQ]HVI5Dg9YIfQd
HW9QHf]K<Ee5DfaUIWARLV5SJdECLfaQLfQUKW=eLVE]HGAXAE=/IFIdHW9QHf]5
Dg=/HG=XKF5dJ69VAE=/IFIdHW9QHf/S<DECLVUWJ7ARLV5SJdECLVUWJ7ARLV5S
JdECLVUWJ7ARLV5SJdECHFi^Kg1UKV@Y2Rm5Lf=QL6ETD79UHFeRK6DP:4ECHFi^
KfE`9CE3I6mSMFeUKWASK65cLbDe@S4`L7@U=D@U=d9QLWAYHfaU9CM49C119CE3
MG=UL65SJf5WIBDe@WEcIFiQKFEc9CE49CM2Hfm/Kg8U=d@U<S0U<SEeLfET9C8`
IVmb9C8`IVm^M2Db<6=_K6mb9C119CE3MG=UL65SJf5WIBDg@V5]Lg=iKF8U=d@U
<S0U<SE]HGAXLbD`@BDe@gEcIG1QHf]QIfDU=d9QKG=]HGAX9CM49C8`9C8eKF5d
J7<U<44U=D=eLfE`HF=[HFMU9CE2MGAV>2DeA2Dg@VU^L7EdIFiS9CM49C8`9C8e
MG=UIWE/9C8`M6lU<S1dNG1U9C8`I6UbIF=dK7TU<S1TJF5SLVUdJF<U<S1SJ65b
HF=dIG9c9C119CE3KVEgHfm]KF5^I2Dg@RDe@fMbHF@U=d@U=d8U=D=UKW=eLVE]
HGAX9CM29CE3KV5RK64U=d@U=d@U<44U=D=^IGMSKfe]HFiT9CM29CE3Hg9_Lg<U
=d@U=d8U=D=UKW=eLVE]HGAX9CM29CE3M6U]IG<U=d@U=d@U<44U=D=^IGMSKfe]
HFiT9CM29CE3KF9V9CM49CE2<BDeA2Dg@RDe@fE^LgEbIFeQM6PU=d8U=D=]HGAX
HVHU=d8U<S<a9CM49CM49CM4AE=QKVi_IG1UKV@Y2RmCMF9ZIF=d82Q5Df5^KVmd
AE=cK65cJ6IbHF=5DfaUIWARLV5SJdECLfaQLfQ`HG9dJF5/GS8PAE=cK65cJ71c
JDECLVUWJ7ARLV5SJdECK6EVM69bHF=[AE=cK65cJ71QLWAYHF`PN5hbAE=bJFMX
M69bHF=[82/PAE=cK65cJ6IbHF=5DfaUIWARLV5SJdECLfaQLfQ`HG9dJF5/GS8P
AE=cK65cJ71cJDECLVUWJ7ARLV5SJdECK6EVM69bHF=[AE=cK65cJ71QLWAYHF`P
NEhbAE=bJFMXM69bHF=[82HP?B0V8315Dg=/HG=XAE=cK65cJ0Y5Dg=/HG=XIW9Q
HdECK6EVM69bHF=[AE=cK65cJ71cJEm5DfaUIWARLV5SJfT[<BaZAE=bJFMXM69b
HF=[82dP<TECLfaQLfQ`LfUOAE=/IFIdHW9QHf]Y;6Y5Dg9YIfQdHW9QHf/P:b15
Dg=/HG=XL7=YGdECK6EVM69bHF=[JBda;6Y5Dg9YIfQdHW9QHf]5Dg9YIfQdHW9Q
Hf]5DfaUIWARLV5SJdECLfaQLfQTIFadHB1hGS95Dg9YIfQdHW9QHf/P:b15Dg=/
HG=XIW9QHdECK6EVM69bHF=[AE=cK65cJ71cJEm5DfaUIWARLV5SJfT/JR/aAE=b
JFMXM69bHF=[82dP<TECLfaQLfQ`LfUOAE=/IFIdHW9QHf]Y;6Y5Dg9YIfQdHW9Q
Hf/P:b15Dg=/HG=XL7=YGdECK6EVM69bHF=[JBaZ;C55Dg9YIfQdHW9QHf]5Dg9Y
IfQdHW9QHf]5DfaUIWARLV5SJdECLfaQLfQTIFadHB1iGS95Dg9YIfQdHW9QHf/P
9R15Dg=/HG=XHG1`LVmh82HP<4ECHFi^KgAUKV@Y2Rm5Lf=QL6ETDgERJVESM20X
AE=QKVi_IG<U=D=VLV5S9CM29CE3L65bM6UQK2DeAC8U<S0U=D=`LfTU=d@U=d8U
=D=`HG9dJF5/9C8`N2DeAC8U=d@U<S0[9C8`9CE3IW9QHbDg@RDe@g1QLWAYHF`U
=DDb9C8`9CE3L7=Y9CM49CM29CE3L65bM6UQK2Db<7TU=DDb9CM49C8`9RDb<3dU
<S0V9C8`<2De@bDe@bD`@BDe@fIbHF<U=d8U=D=`LfUO9CM2JB/a;6XU=d@U<S0]
9C8`<RDe@g1cJElU=d9Y;6XU=d@U<S0[9C8`9CE3L7=YGbDg@VT]<BaZ9CM49CM4
9CM29CE3I6E/M64U<S1h9CE5<RDgA2Db<2/U<S0U=D=VLV5S9CM29CE3L7=YGbDg
@VT/JR/a9CM49C8`;BDb<38U=D=`LfUO9CM2JBaZ9CM49C8`:bDb<2De@g1cJElU
=d9Y;6X]<BDgA2DgA2Dg@RDe@fAUK7AQ9C8`NBDeAC8U=d@U<S0V9C8`9CE3HG1`
LVmh9C8`9RDb<315Df5^KVmULfE^I2T:;eAiL6DP:4E5M7U`IC=5AGAiL6EUKV@Y
2Rm3Kfa_LR0XADESKf``;S0`<30`<20`;S0`<30`<20`;S0`<30`<20a;S0`<30`
<4E5Hfm/IFiT:@X_@T]3Kfa_LR0XADERJf<a;S0`<30`<20a;S0`<30`<20a;S0`
<30`<20a;S0`<30`<4E5HV]SIFiT:@X_E6UdK6DP:4E5M6UdK6E5AGAYM6aUIFiT
:@X_CF5WKVUVJF=QM6U_KR0XADE]HFLc=Rh`<30`<315AFeQIfE^I2T:;d9QLfE/
JFiU82Q5AF9QLc0^<30`<30`ADERHG=UKV@Y2Shn2VE^I6mRJPYhLVEV2S0P<C0:
<30`<30`<30`<20f=CDc=B1V80X`<30`<34g=3Da830`<30`86hP2S0`<30`<CPi
=S4P<30`<30PKR0:<30`<30a=c<i<R0`<30`<21^80X`<30`<34g<SHa830`<30`
86hP2S0`<30`<30`<CDP<30`<30PKR0:<30`<30a=c8d<20`<30`<21^80X`<30`
<34g=C4e830`<30`86hP2S0`<30`<CLe=CHP<30`<30PKR0:<30`<30a=cDh=B0`
<30`<21^80YdLV5YK6Eb2S`l82mCJGYU834`82mBKfmd834P<21B82m9KVI_838P
<21B2Rm9A21K?3U4AC=2=CE4ADI6@S52=3Hd=3Q3ADA1=D4`A4A2@cHa?S`iA4Dc
@SDeA4E6AT8a@S@f=3@h@dE4@CE1<4A4@T<f<CiM2Shn2W=dHG9dN79UIPXb<CPc
<0XU9DE?AP\>"], "Graphics",
 ImageSize->{389, 76}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"General", "::", "spell1"}], "]"}], ";"}]], "Input",
 CellLabel->"In[1]:="],

Cell["Define Laplace's Equation\[Ellipsis]", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"eqn", "[", 
   RowBox[{"i_", ",", "j_"}], "]"}], " ", ":=", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"(", " ", 
      RowBox[{
       RowBox[{"\[Psi]", "[", 
        RowBox[{
         RowBox[{"i", "+", "1"}], ",", "j"}], "]"}], " ", "-", " ", 
       RowBox[{"2", " ", 
        RowBox[{"\[Psi]", "[", 
         RowBox[{"i", ",", "j"}], "]"}]}], " ", "+", " ", 
       RowBox[{"\[Psi]", "[", 
        RowBox[{
         RowBox[{"i", "-", "1"}], ",", "j"}], "]"}]}], ")"}], "/", 
     RowBox[{"\[Delta]x", "^", "2"}]}], " ", "+", " ", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Psi]", "[", 
        RowBox[{"i", ",", 
         RowBox[{"j", "+", "1"}]}], "]"}], " ", "-", " ", 
       RowBox[{"2", " ", 
        RowBox[{"\[Psi]", "[", 
         RowBox[{"i", ",", "j"}], "]"}]}], " ", "+", " ", 
       RowBox[{"\[Psi]", "[", 
        RowBox[{"i", ",", 
         RowBox[{"j", "-", "1"}]}], "]"}]}], ")"}], "/", 
     RowBox[{"\[Delta]y", "^", "2"}]}]}], " ", "\[Equal]", " ", 
   "0"}]}]], "Input",
 CellLabel->"In[2]:="],

Cell[BoxData[
 RowBox[{"eqn", "[", 
  RowBox[{"i", ",", "j"}], "]"}]], "Input",
 CellLabel->"In[3]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Define a Grid", "Section"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"height", " ", "=", " ", "1.0"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox["meter",
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"length", " ", "=", " ", "5.0"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox["meter",
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"wallHeight", " ", "=", " ", "0.5"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox["meter",
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"wallLength", " ", "=", " ", "0.5"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox["meter",
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], ";"}]}], "Input",
 CellLabel->"In[4]:="],

Cell[TextData[{
 "Grid points go from 1 to ",
 StyleBox["numX",
  FontWeight->"Bold"],
 " and from 1 to ",
 StyleBox["numY",
  FontWeight->"Bold"],
 "."
}], "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"numX", " ", "=", " ", "80"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"numY", " ", "=", " ", "60"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{"\"\<Number of Grid Points = \>\"", ",", " ", 
    RowBox[{"numX", " ", "numY"}]}], "]"}], ";"}]}], "Input",
 CellLabel->"In[8]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"dx", " ", "=", " ", 
   RowBox[{"length", "/", 
    RowBox[{"(", 
     RowBox[{"numX", "-", "1"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dy", " ", "=", " ", 
   RowBox[{"height", "/", 
    RowBox[{"(", 
     RowBox[{"numY", " ", "-", "1"}], ")"}]}]}], ";"}]}], "Input",
 CellLabel->"In[11]:="],

Cell[BoxData[
 RowBox[{"wallLength", "/", "dx"}]], "Input",
 CellLabel->"In[13]:="],

Cell[BoxData[
 RowBox[{"wallHeight", "/", "dy"}]], "Input",
 CellLabel->"In[14]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"flow", " ", "=", " ", "5.0"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{"m", "/", "s"}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], ";"}]], "Input",
 CellLabel->"In[15]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"grid", " ", "=", " ", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         StyleBox[
          RowBox[{
           RowBox[{
            RowBox[{"length", "/", "2"}], "-", 
            RowBox[{"wallLength", "/", "2"}]}], "\[LessEqual]", " ", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"ix", "-", "1"}], ")"}], " ", "dx"}], " ", 
           "\[LessEqual]", " ", 
           RowBox[{
            RowBox[{"length", "/", "2"}], " ", "+", 
            RowBox[{"wallLength", "/", "2"}]}]}],
          FontColor->RGBColor[1, 0, 1]], 
         StyleBox[" ",
          FontColor->RGBColor[1, 0, 1]], 
         StyleBox["&&",
          FontColor->GrayLevel[0]], 
         StyleBox["  ",
          FontColor->GrayLevel[0]], 
         StyleBox[
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"iy", "-", "1"}], ")"}], " ", "dy"}], "\[LessEqual]", 
           " ", "wallHeight"}],
          FontColor->GrayLevel[0]]}], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox["||",
         FontColor->GrayLevel[0]], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox[
         RowBox[{"ix", " ", "\[Equal]", " ", "1"}],
         FontColor->GrayLevel[0]], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox["||",
         FontColor->GrayLevel[0]], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox[
         RowBox[{"ix", " ", "==", " ", "numX"}],
         FontColor->GrayLevel[0]], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox["||",
         FontColor->GrayLevel[0]], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox[
         RowBox[{"iy", " ", "\[Equal]", " ", "1"}],
         FontColor->GrayLevel[0]], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox["||",
         FontColor->GrayLevel[0]], 
        StyleBox[" ",
         FontColor->GrayLevel[0]], 
        StyleBox[
         RowBox[{"iy", " ", "==", " ", "numY"}],
         FontColor->GrayLevel[0]]}], ",", " ", "1", ",", " ", "0"}], "]"}], 
     ",", " ", 
     RowBox[{"{", 
      RowBox[{"ix", ",", " ", "1", ",", " ", "numX"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"iy", ",", " ", "1", ",", " ", "numY"}], "}"}]}], "]"}]}], 
  ";"}]], "Input",
 CellLabel->"In[16]:="],

Cell[BoxData[
 RowBox[{"Dimensions", "[", "grid", "]"}]], "Input",
 CellLabel->"In[17]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"Count", "[", 
   RowBox[{
    RowBox[{"Flatten", "[", "grid", "]"}], ",", 
    RowBox[{"x_", "/;", 
     RowBox[{"x", "\[Equal]", "1"}]}]}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    StyleBox["number",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["of",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["known",
     FontColor->RGBColor[1, 0, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["grid",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["values",
     FontColor->RGBColor[0, 1, 0]]}], " ", "*)"}]}]], "Input",
 CellLabel->"In[18]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"Count", "[", 
   RowBox[{
    RowBox[{"Flatten", "[", "grid", "]"}], ",", 
    RowBox[{"x_", "/;", 
     RowBox[{"x", "\[Equal]", "0"}]}]}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    StyleBox["number",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["of",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["unknown",
     FontColor->RGBColor[1, 0, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["grid",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox[" ",
     FontColor->RGBColor[0, 1, 0]], 
    StyleBox["values",
     FontColor->RGBColor[0, 1, 0]]}], " ", "*)"}]}]], "Input",
 CellLabel->"In[19]:="],

Cell[BoxData[
 RowBox[{"Length", "[", 
  RowBox[{"Flatten", "[", "grid", "]"}], "]"}]], "Input",
 CellLabel->"In[20]:="],

Cell[BoxData[
 RowBox[{"ArrayPlot", "[", 
  RowBox[{
   RowBox[{"Transpose", "[", "grid", "]"}], ",", 
   RowBox[{"Mesh", "\[Rule]", "True"}]}], "]"}]], "Input",
 CellLabel->"In[21]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"wallBound", "[", "ix_", "]"}], " ", ":=", " ", 
  RowBox[{"If", "[", 
   RowBox[{
    StyleBox[
     RowBox[{
      RowBox[{
       RowBox[{"length", "/", "2"}], "-", 
       RowBox[{"wallLength", "/", "2"}]}], "\[LessEqual]", " ", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"ix", "-", "1"}], ")"}], " ", "dx"}], " ", "\[LessEqual]", 
      " ", 
      RowBox[{
       RowBox[{"length", "/", "2"}], " ", "+", 
       RowBox[{"wallLength", "/", "2"}]}]}],
     FontColor->RGBColor[1, 0, 1]], ",", " ", 
    RowBox[{
     RowBox[{"Floor", "[", 
      RowBox[{"wallHeight", "/", "dy"}], "]"}], "+", "1"}], ",", " ", "1"}], 
   "]"}]}]], "Input",
 CellLabel->"In[22]:="],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"wallBound", "[", "ix", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"ix", ",", 
       RowBox[{"Floor", "[", 
        FractionBox["numX", "4"], "]"}], ",", 
       RowBox[{"Floor", "[", 
        FractionBox[
         RowBox[{"3", " ", "numX"}], "4"], "]"}]}], "}"}]}], "]"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"PointSize", "[", "0.02`", "]"}]}], ",", 
   RowBox[{"Axes", "\[Rule]", "False"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<x\>\"", ",", "\"\<h\>\"", ",", "\"\<Wall Boundary Index\>\"", ",", 
      "\"\<\>\""}], "}"}]}]}], "]"}]], "Input",
 CellLabel->"In[23]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Define Boundary Conditions", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Clear", "[", "\[Psi]", "]"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox[
    RowBox[{
    "\[Psi]", " ", "is", " ", "the", " ", "variable", " ", "used", " ", "to", 
     " ", "store", " ", "the", " ", 
     RowBox[{"streamfunction", "."}]}],
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}]}]], "Input",
 CellLabel->"In[24]:="],

Cell["Top and bottom\[Ellipsis]", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"ix_", ",", "numY"}], "]"}], " ", "=", " ", 
   RowBox[{"flow", "  ", "height"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"ix_", ",", "1"}], "]"}], " ", "=", " ", "0"}], ";"}]}], "Input",
 CellLabel->"In[25]:="],

Cell["Left and right sides\[Ellipsis]", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"1", ",", "iy_"}], "]"}], " ", "=", " ", 
   RowBox[{"flow", " ", 
    RowBox[{"(", 
     RowBox[{"iy", "-", "1"}], ")"}], " ", "dy"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"numX", ",", "iy_"}], "]"}], " ", "=", " ", 
   RowBox[{"flow", " ", 
    RowBox[{"(", 
     RowBox[{"iy", "-", "1"}], ")"}], " ", "dy"}]}], ";"}]}], "Input",
 CellLabel->"In[27]:="],

Cell["The wall\[Ellipsis]", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"\[Psi]", "[", 
        RowBox[{"ix", ",", "iy"}], "]"}], " ", "=", " ", "0"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"iy", ",", " ", "1", ",", " ", 
        RowBox[{"wallBound", "[", "ix", "]"}]}], "}"}]}], "]"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"ix", ",", " ", "1", ",", " ", "numX"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellLabel->"In[29]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Define the Equations", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{"eqns", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Flatten", "[", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{"eqn", "[", 
           RowBox[{"ix", ",", "iy"}], "]"}], ",", 
          RowBox[{"{", 
           RowBox[{"iy", ",", " ", 
            RowBox[{
             RowBox[{"wallBound", "[", "ix", "]"}], "+", "1"}], ",", " ", 
            RowBox[{"numY", "-", "1"}]}], "}"}]}], "]"}], " ", ",", 
        RowBox[{"{", 
         RowBox[{"ix", ",", " ", "2", ",", " ", 
          RowBox[{"numX", "-", "1"}]}], "}"}]}], "]"}], "]"}], " ", "/.", " ",
      
     RowBox[{"\[Delta]x", " ", "\[Rule]", " ", "dx"}]}], " ", "/.", " ", 
    RowBox[{"\[Delta]y", " ", "\[Rule]", " ", "dy"}]}]}], ";"}]], "Input",
 CellLabel->"In[30]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"eqns", " ", "=", " ", 
   RowBox[{"Simplify", "[", "eqns", "]"}]}], ";"}]], "Input",
 CellLabel->"In[31]:="],

Cell[BoxData[
 RowBox[{"Dimensions", "[", "eqns", "]"}]], "Input",
 CellLabel->"In[32]:="],

Cell[BoxData[
 RowBox[{"eqns", "\[LeftDoubleBracket]", 
  RowBox[{"{", 
   RowBox[{"1", ",", "2", ",", "3"}], "}"}], 
  "\[RightDoubleBracket]"}]], "Input",
 CellLabel->"In[33]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"unknowns", " ", "=", " ", 
   RowBox[{"Flatten", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"\[Psi]", "[", 
         RowBox[{"ix", ",", "iy"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"iy", ",", " ", 
          RowBox[{
           RowBox[{"wallBound", "[", "ix", "]"}], "+", "1"}], ",", " ", 
          RowBox[{"numY", "-", "1"}]}], "}"}]}], "]"}], " ", ",", 
      RowBox[{"{", 
       RowBox[{"ix", ",", " ", "2", ",", " ", 
        RowBox[{"numX", "-", "1"}]}], "}"}]}], "]"}], "]"}]}], ";"}]], "Input",\

 CellLabel->"In[34]:="],

Cell[BoxData[
 RowBox[{"unknowns", "\[LeftDoubleBracket]", 
  RowBox[{"{", 
   RowBox[{"1", ",", "2", ",", "3"}], "}"}], 
  "\[RightDoubleBracket]"}]], "Input",
 CellLabel->"In[35]:="],

Cell[BoxData[
 RowBox[{"Dimensions", "[", "unknowns", "]"}]], "Input",
 CellLabel->"In[36]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Finding the Solution", "Section"],

Cell[BoxData[
 RowBox[{"Timing", "[", 
  RowBox[{
   RowBox[{"sol", " ", "=", " ", 
    RowBox[{"Flatten", "[", 
     RowBox[{"Solve", "[", 
      RowBox[{"eqns", ",", "unknowns"}], "]"}], "]"}]}], ";"}], 
  "]"}]], "Input",
 CellLabel->"In[37]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Psi]Sol", " ", "=", " ", 
   RowBox[{"Interpolation", "[", 
    RowBox[{"Flatten", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"ix", "-", "1"}], ")"}], " ", "dx"}], ",", " ", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"iy", "-", "1"}], ")"}], " ", "dy"}], ",", " ", 
           RowBox[{"\[Psi]", "[", 
            RowBox[{"ix", ",", "iy"}], "]"}]}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"ix", ",", " ", "1", ",", " ", "numX"}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"iy", ",", " ", "1", ",", " ", "numY"}], "}"}]}], "]"}], 
       " ", "/.", " ", "sol"}], ",", "1"}], "]"}], "]"}]}], ";"}]], "Input",
 CellLabel->"In[38]:="],

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{"\[Psi]Sol", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "length"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
   RowBox[{"Contours", "\[Rule]", "18"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
   RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
   " ", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Streamfunction Around Wall\>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.432286002259787*^9, 3.4322860161594954`*^9}},
 CellLabel->"In[39]:="],

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{"\[Psi]Sol", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", 
     RowBox[{
      FractionBox["length", "2"], "+", 
      FractionBox["height", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
   RowBox[{"Contours", "\[Rule]", "18"}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
   RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Streamfunction Around Wall\>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{3.4322860358686333`*^9},
 CellLabel->"In[40]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Computing the Velocity", "Section"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"uX", "[", 
    RowBox[{"xx_", ",", "yy_"}], "]"}], " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"\[Psi]Sol", "[", 
        RowBox[{"x", ",", "y"}], "]"}], ",", "y"}], "]"}], " ", "/.", " ", 
     RowBox[{"x", " ", "\[Rule]", " ", "xx"}]}], " ", "/.", " ", 
    RowBox[{"y", " ", "\[Rule]", " ", "yy"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"uY", "[", 
    RowBox[{"xx_", ",", "yy_"}], "]"}], " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"\[Psi]Sol", "[", 
         RowBox[{"x", ",", "y"}], "]"}], ",", "x"}], "]"}]}], "/.", " ", 
     RowBox[{"x", " ", "\[Rule]", " ", "xx"}]}], " ", "/.", " ", 
    RowBox[{"y", " ", "\[Rule]", " ", "yy"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"u2", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"uX", "[", 
      RowBox[{"x", ",", "y"}], "]"}], "^", "2"}], " ", "+", " ", 
    RowBox[{
     RowBox[{"uY", "[", 
      RowBox[{"x", ",", "y"}], "]"}], "^", "2"}]}]}], ";"}]}], "Input",
 CellLabel->"In[41]:="],

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{
    FractionBox["1", "2"], " ", 
    RowBox[{"u2", "[", 
     RowBox[{"x", ",", "y"}], "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", 
     RowBox[{
      FractionBox["length", "2"], "+", 
      FractionBox["height", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
   RowBox[{"Contours", "\[Rule]", "18"}], ",", 
   RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<U^2/2 Around Wall\>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{3.432286041621207*^9},
 CellLabel->"In[44]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], " ", "=", " ", 
   RowBox[{"1", " ", "-", " ", 
    RowBox[{
     RowBox[{"u2", "[", 
      RowBox[{"x", ",", "y"}], "]"}], "/", 
     RowBox[{"flow", "^", "2"}]}]}]}], ";"}]], "Input",
 CellLabel->"In[45]:="],

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", 
     RowBox[{
      FractionBox["length", "2"], "+", 
      FractionBox["height", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
   RowBox[{"Contours", "\[Rule]", "18"}], ",", 
   RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Normalized Pressure\>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{3.432286045613276*^9},
 CellLabel->"In[46]:="],

Cell[BoxData[
 RowBox[{"Plot3D", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", 
     RowBox[{
      FractionBox["length", "2"], "+", 
      FractionBox["height", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Normalized Pressure\>\""}]}], 
  "]"}]], "Input",
 CellLabel->"In[47]:="],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"%", ",", 
   RowBox[{"ViewPoint", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"1.2`", ",", "1.2`", ",", "1.2`"}], "}"}]}]}], "]"}]], "Input",
 CellLabel->"In[48]:="],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<y\>\"", ",", "\"\<Cp\>\""}], "}"}]}], ",", 
   RowBox[{
   "PlotLabel", "\[Rule]", "\"\<Pressure Along Leading Wall Edge\>\""}]}], 
  "]"}]], "Input",
 CellLabel->"In[49]:="],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{"x", ",", "wallHeight"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "length"}], "}"}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<Cp\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Pressure Along Wall Height\>\""}]}], 
  "]"}]], "Input",
 CellLabel->"In[50]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Flow Past a Cylinder", "Section"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{"\[Psi]", ",", "uX", ",", "uY", ",", "u2"}], "]"}], ";"}]], "Input",\

 CellLabel->"In[51]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"radius", " ", "=", " ", "0.2"}], " ", 
  RowBox[{"(*", " ", 
   StyleBox["meter",
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"centerY", " ", "=", " ", 
   RowBox[{"height", "/", "2"}]}], "  ", 
  RowBox[{"(*", " ", 
   StyleBox["meter",
    FontColor->RGBColor[0, 1, 0]], " ", "*)"}], ";"}]}], "Input",
 CellLabel->"In[52]:="],

Cell[CellGroupData[{

Cell["Define Boundary Conditions", "Subsection"],

Cell["Top and bottom\[Ellipsis]", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"ix_", ",", "numY"}], "]"}], " ", "=", " ", 
   RowBox[{"flow", "  ", "height"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"ix_", ",", "1"}], "]"}], " ", "=", " ", "0"}], ";"}]}], "Input",
 CellLabel->"In[54]:="],

Cell["Left and right sides\[Ellipsis]", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"1", ",", "iy_"}], "]"}], " ", "=", " ", 
   RowBox[{"flow", " ", 
    RowBox[{"(", 
     RowBox[{"iy", "-", "1"}], ")"}], " ", "dy"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"numX", ",", "iy_"}], "]"}], " ", "=", " ", 
   RowBox[{"flow", " ", 
    RowBox[{"(", 
     RowBox[{"iy", "-", "1"}], ")"}], " ", "dy"}]}], ";"}]}], "Input",
 CellLabel->"In[56]:="],

Cell["The cylinder", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      StyleBox[
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"ix", "-", "1"}], ")"}], " ", "dx"}], "  ", "-", " ", 
            RowBox[{"length", "/", "2"}]}], ")"}], "^", "2"}], " ", "+", " ", 
         
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"iy", " ", "-", "1"}], ")"}], "dy"}], " ", "-", " ", 
            "centerY"}], ")"}], "^", "2"}]}], " ", "\[LessEqual]", " ", 
        RowBox[{"radius", "^", "2"}]}],
       FontColor->RGBColor[1, 0, 1]], ",", " ", 
      RowBox[{
       RowBox[{"\[Psi]", "[", 
        RowBox[{"ix", ",", "iy"}], "]"}], " ", "=", " ", 
       RowBox[{"flow", " ", "centerY"}]}]}], "]"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"iy", ",", " ", "1", ",", "numY"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"ix", ",", " ", "1", ",", " ", "numX"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellLabel->"In[58]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Define the Equations", "Subsection"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"eqns", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Flatten", "[", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"If", "[", 
         RowBox[{
          StyleBox[
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"ix", "-", "1"}], ")"}], " ", "dx"}], "  ", "-", 
                " ", 
                RowBox[{"length", "/", "2"}]}], ")"}], "^", "2"}], " ", "+", 
             " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"iy", " ", "-", "1"}], ")"}], "dy"}], " ", "-", " ",
                 "centerY"}], ")"}], "^", "2"}]}], " ", "\[LessEqual]", " ", 
            RowBox[{"radius", "^", "2"}]}],
           FontColor->RGBColor[1, 0, 1]], ",", 
          RowBox[{"-", "1"}], ",", 
          RowBox[{"eqn", "[", 
           RowBox[{"ix", ",", "iy"}], "]"}]}], "]"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"iy", ",", " ", "2", ",", 
          RowBox[{"numY", "-", "1"}]}], "}"}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"ix", ",", " ", "2", ",", " ", 
          RowBox[{"numX", "-", "1"}]}], "}"}]}], "]"}], "]"}], " ", "/.", " ",
      
     RowBox[{"\[Delta]x", " ", "\[Rule]", " ", "dx"}]}], " ", "/.", " ", 
    RowBox[{"\[Delta]y", " ", "\[Rule]", " ", "dy"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"eqns", " ", "=", " ", 
   RowBox[{"DeleteCases", "[", 
    RowBox[{"eqns", ",", " ", 
     RowBox[{"_", "?", "Negative"}]}], "]"}]}], ";"}]}], "Input",
 CellLabel->"In[59]:="],

Cell[BoxData[
 RowBox[{"Dimensions", "[", "eqns", "]"}]], "Input",
 CellLabel->"In[61]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"eqns", " ", "=", " ", 
   RowBox[{"Simplify", "[", "eqns", "]"}]}], ";"}]], "Input",
 CellLabel->"In[62]:="],

Cell[BoxData[
 RowBox[{"eqns", "\[LeftDoubleBracket]", 
  RowBox[{"{", 
   RowBox[{"1", ",", "2", ",", "3"}], "}"}], 
  "\[RightDoubleBracket]"}]], "Input",
 CellLabel->"In[63]:="],

Cell[BoxData[{
 RowBox[{
  RowBox[{"unknowns", " ", "=", " ", 
   RowBox[{"Flatten", "[", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        StyleBox[
         RowBox[{
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"(", 
                RowBox[{"ix", "-", "1"}], ")"}], " ", "dx"}], "  ", "-", " ", 
              
              RowBox[{"length", "/", "2"}]}], ")"}], "^", "2"}], " ", "+", 
           " ", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{
               RowBox[{"(", 
                RowBox[{"iy", " ", "-", "1"}], ")"}], "dy"}], " ", "-", " ", 
              "centerY"}], ")"}], "^", "2"}]}], " ", "\[LessEqual]", " ", 
          RowBox[{"radius", "^", "2"}]}],
         FontColor->RGBColor[1, 0, 1]], ",", 
        RowBox[{"-", "1"}], ",", 
        RowBox[{"\[Psi]", "[", 
         RowBox[{"ix", ",", "iy"}], "]"}]}], "]"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"iy", ",", " ", "2", ",", 
        RowBox[{"numY", "-", "1"}]}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"ix", ",", " ", "2", ",", " ", 
        RowBox[{"numX", "-", "1"}]}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"unknowns", " ", "=", " ", 
   RowBox[{"DeleteCases", "[", 
    RowBox[{"unknowns", ",", " ", 
     RowBox[{"_", "?", "Negative"}]}], "]"}]}], ";"}]}], "Input",
 CellLabel->"In[64]:="],

Cell[BoxData[
 RowBox[{"unknowns", "\[LeftDoubleBracket]", 
  RowBox[{"{", 
   RowBox[{"1", ",", "2", ",", "3"}], "}"}], 
  "\[RightDoubleBracket]"}]], "Input",
 CellLabel->"In[66]:="],

Cell[BoxData[
 RowBox[{"Dimensions", "[", "unknowns", "]"}]], "Input",
 CellLabel->"In[67]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Solution", "Subsection"],

Cell[TextData[{
 "The solution is found quickly using  ",
 StyleBox["Mathematica",
  FontSlant->"Italic"],
 "'s built-in ",
 StyleBox["Solve",
  FontWeight->"Bold"],
 " function."
}], "Text"],

Cell[BoxData[
 RowBox[{"Timing", "[", 
  RowBox[{
   RowBox[{"sol", " ", "=", " ", 
    RowBox[{"Flatten", "[", 
     RowBox[{"Solve", "[", 
      RowBox[{"eqns", ",", "unknowns"}], "]"}], "]"}]}], ";"}], 
  "]"}]], "Input",
 CellLabel->"In[68]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Psi]Sol", " ", "=", " ", 
   RowBox[{"Interpolation", "[", 
    RowBox[{"Flatten", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"ix", "-", "1"}], ")"}], " ", "dx"}], ",", " ", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"iy", "-", "1"}], ")"}], " ", "dy"}], ",", " ", 
           RowBox[{"\[Psi]", "[", 
            RowBox[{"ix", ",", "iy"}], "]"}]}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"ix", ",", " ", "1", ",", " ", "numX"}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"iy", ",", " ", "1", ",", " ", "numY"}], "}"}]}], "]"}], 
       " ", "/.", " ", "sol"}], ",", "1"}], "]"}], "]"}]}], ";"}]], "Input",
 CellLabel->"In[69]:="],

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{"\[Psi]Sol", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "length"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
   RowBox[{"Contours", "\[Rule]", "24"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
   RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
   RowBox[{
   "PlotLabel", "\[Rule]", "\"\<Streamfunction Around Cylinder\>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{3.4322860537699738`*^9},
 CellLabel->"In[70]:="],

Cell[BoxData[
 RowBox[{"plt", "=", 
  RowBox[{"ContourPlot", "[", 
   RowBox[{
    RowBox[{"\[Psi]Sol", "[", 
     RowBox[{"x", ",", "y"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", 
      RowBox[{
       FractionBox["length", "2"], "-", 
       FractionBox["height", "2"]}], ",", 
      RowBox[{
       FractionBox["length", "2"], "+", 
       FractionBox["height", "2"]}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
    RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
    RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
    RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
    RowBox[{"Contours", "\[Rule]", "24"}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
    RowBox[{
    "PlotLabel", "\[Rule]", "\"\<Streamfunction Around Cylinder\>\""}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{3.432286057117145*^9},
 CellLabel->"In[71]:="],

Cell[BoxData[
 RowBox[{"cylinder", "=", 
  RowBox[{"Graphics", "[", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
     RowBox[{"Circle", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         FractionBox["length", "2"], ",", "centerY"}], "}"}], ",", "radius"}],
       "]"}]}], "}"}], "]"}]}]], "Input",
 CellLabel->"In[72]:="],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plt", ",", "cylinder"}], "]"}]], "Input",
 CellLabel->"In[73]:="]
}, Closed]],

Cell[CellGroupData[{

Cell["Computing the Velocity", "Subsection"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"uX", "[", 
    RowBox[{"xx_", ",", "yy_"}], "]"}], " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"\[Psi]Sol", "[", 
        RowBox[{"x", ",", "y"}], "]"}], ",", "y"}], "]"}], " ", "/.", " ", 
     RowBox[{"x", " ", "\[Rule]", " ", "xx"}]}], " ", "/.", " ", 
    RowBox[{"y", " ", "\[Rule]", " ", "yy"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"uY", "[", 
    RowBox[{"xx_", ",", "yy_"}], "]"}], " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"\[Psi]Sol", "[", 
         RowBox[{"x", ",", "y"}], "]"}], ",", "x"}], "]"}]}], "/.", " ", 
     RowBox[{"x", " ", "\[Rule]", " ", "xx"}]}], " ", "/.", " ", 
    RowBox[{"y", " ", "\[Rule]", " ", "yy"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"u2", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"uX", "[", 
      RowBox[{"x", ",", "y"}], "]"}], "^", "2"}], " ", "+", " ", 
    RowBox[{
     RowBox[{"uY", "[", 
      RowBox[{"x", ",", "y"}], "]"}], "^", "2"}]}]}], ";"}]}], "Input",
 CellLabel->"In[74]:="],

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{
    FractionBox["1", "2"], " ", 
    RowBox[{"u2", "[", 
     RowBox[{"x", ",", "y"}], "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", 
     RowBox[{
      FractionBox["length", "2"], "+", 
      FractionBox["height", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
   RowBox[{"Contours", "\[Rule]", "18"}], ",", 
   RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<U^2/2 Around Cylinder\>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{3.432286063149158*^9},
 CellLabel->"In[77]:="],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{"x_", ",", "y_"}], "]"}], " ", "=", " ", 
   RowBox[{"1", " ", "-", " ", 
    RowBox[{
     RowBox[{"u2", "[", 
      RowBox[{"x", ",", "y"}], "]"}], "/", 
     RowBox[{"flow", "^", "2"}]}]}]}], ";"}]], "Input",
 CellLabel->"In[78]:="],

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", 
     RowBox[{
      FractionBox["length", "2"], "+", 
      FractionBox["height", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "False"}], ",", 
   RowBox[{"Contours", "\[Rule]", "18"}], ",", 
   RowBox[{"PerformanceGoal", " ", "\[Rule]", " ", "\"\<Speed\>\""}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Normalized Pressure\>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{3.432286067540759*^9},
 CellLabel->"In[79]:="],

Cell[BoxData[
 RowBox[{"Plot3D", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{"x", ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{
      FractionBox["length", "2"], "-", 
      FractionBox["height", "2"]}], ",", 
     RowBox[{
      FractionBox["length", "2"], "+", 
      FractionBox["height", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", "80"}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "\"\<Normalized Pressure\>\""}]}], 
  "]"}]], "Input",
 CellLabel->"In[80]:="],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"%", ",", 
   RowBox[{"ViewPoint", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"1.2`", ",", "1.2`", ",", "1.2`"}], "}"}]}]}], "]"}]], "Input",
 CellLabel->"In[81]:="],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"%%", ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "Input",
 CellLabel->"In[82]:="],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{
     FractionBox["length", "2"], ",", "y"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"y", ",", "0", ",", "height"}], "}"}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<y\>\"", ",", "\"\<Cp\>\""}], "}"}]}], ",", 
   RowBox[{
   "PlotLabel", "\[Rule]", "\"\<Pressure Profile through Cylinder\>\""}]}], 
  "]"}]], "Input",
 CellLabel->"In[83]:="],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Cp", "[", 
    RowBox[{
     RowBox[{
      FractionBox["length", "2"], "+", 
      RowBox[{"1.1`", " ", "radius", " ", 
       RowBox[{"Cos", "[", "\[Theta]", "]"}]}]}], ",", 
     RowBox[{"centerY", "+", 
      RowBox[{"1.1`", " ", "radius", " ", 
       RowBox[{"Sin", "[", "\[Theta]", "]"}]}]}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "0", ",", 
     RowBox[{"2", " ", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<\[Theta]\>\"", ",", "\"\<Cp\>\""}], "}"}]}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Range", "[", 
       RowBox[{"0", ",", 
        RowBox[{"2", " ", "\[Pi]"}], ",", 
        FractionBox["\[Pi]", "2"]}], "]"}], ",", "Automatic"}], "}"}]}], ",", 
   
   RowBox[{
   "PlotLabel", "\[Rule]", "\"\<Pressure Profile around Cylinder\>\""}]}], 
  "]"}]], "Input",
 CellLabel->"In[84]:="]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{

Cell["Summary", "Section"],

Cell[TextData[{
 "Two dimensional potential flow requires the solution of an elliptic partial \
differential equation. These problems can be solved ",
 StyleBox["very easily ",
  FontSlant->"Italic"],
 "with today's computers."
}], "Text"]
}, Closed]]
}, Open  ]]
},
WindowSize->{653, 828},
WindowMargins->{{19, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{1, Automatic},
ShowSelection->True,
FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (May 21, 2008)",
StyleDefinitions->"TutorialBook.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 62, 0, 103, "Title"],
Cell[655, 25, 78, 3, 43, "Subsubtitle"],
Cell[736, 30, 511, 11, 74, "Text"],
Cell[CellGroupData[{
Cell[1272, 45, 54, 0, 87, "Section"],
Cell[1329, 47, 30110, 464, 84, "Graphics"],
Cell[31442, 513, 133, 4, 27, "Input"],
Cell[31578, 519, 52, 0, 23, "Text"],
Cell[31633, 521, 1083, 33, 68, "Input"],
Cell[32719, 556, 102, 3, 27, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[32858, 564, 32, 0, 55, "Section"],
Cell[32893, 566, 769, 24, 88, "Input"],
Cell[33665, 592, 164, 8, 22, "Text"],
Cell[33832, 602, 356, 9, 68, "Input"],
Cell[34191, 613, 354, 11, 48, "Input"],
Cell[34548, 626, 83, 2, 27, "Input"],
Cell[34634, 630, 83, 2, 27, "Input"],
Cell[34720, 634, 224, 7, 27, "Input"],
Cell[34947, 643, 2470, 77, 168, "Input"],
Cell[37420, 722, 90, 2, 27, "Input"],
Cell[37513, 726, 785, 27, 48, "Input"],
Cell[38301, 755, 787, 27, 48, "Input"],
Cell[39091, 784, 120, 3, 27, "Input"],
Cell[39214, 789, 185, 5, 27, "Input"],
Cell[39402, 796, 711, 22, 68, "Input"],
Cell[40116, 820, 790, 22, 142, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[40943, 847, 45, 0, 55, "Section"],
Cell[40991, 849, 377, 11, 48, "Input"],
Cell[41371, 862, 41, 0, 23, "Text"],
Cell[41415, 864, 339, 10, 48, "Input"],
Cell[41757, 876, 47, 0, 23, "Text"],
Cell[41807, 878, 489, 16, 48, "Input"],
Cell[42299, 896, 35, 0, 22, "Text"],
Cell[42337, 898, 497, 15, 48, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[42871, 918, 39, 0, 55, "Section"],
Cell[42913, 920, 843, 23, 108, "Input"],
Cell[43759, 945, 143, 4, 27, "Input"],
Cell[43905, 951, 90, 2, 27, "Input"],
Cell[43998, 955, 180, 5, 27, "Input"],
Cell[44181, 962, 644, 19, 108, "Input"],
Cell[44828, 983, 184, 5, 27, "Input"],
Cell[45015, 990, 94, 2, 27, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[45146, 997, 39, 0, 55, "Section"],
Cell[45188, 999, 248, 8, 27, "Input"],
Cell[45439, 1009, 876, 24, 88, "Input"],
Cell[46318, 1035, 819, 20, 88, "Input"],
Cell[47140, 1057, 1011, 26, 183, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[48188, 1088, 41, 0, 55, "Section"],
Cell[48232, 1090, 1230, 39, 68, "Input"],
Cell[49465, 1131, 1042, 28, 171, "Input"],
Cell[50510, 1161, 305, 10, 27, "Input"],
Cell[50818, 1173, 995, 26, 183, "Input"],
Cell[51816, 1201, 593, 18, 107, "Input"],
Cell[52412, 1221, 212, 6, 27, "Input"],
Cell[52627, 1229, 519, 16, 107, "Input"],
Cell[53149, 1247, 432, 12, 68, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[53618, 1264, 39, 0, 55, "Section"],
Cell[53660, 1266, 154, 5, 27, "Input"],
Cell[53817, 1273, 421, 13, 48, "Input"],
Cell[CellGroupData[{
Cell[54263, 1290, 48, 0, 35, "Subsection"],
Cell[54314, 1292, 41, 0, 23, "Text"],
Cell[54358, 1294, 339, 10, 48, "Input"],
Cell[54700, 1306, 47, 0, 23, "Text"],
Cell[54750, 1308, 489, 16, 48, "Input"],
Cell[55242, 1326, 28, 0, 22, "Text"],
Cell[55273, 1328, 1155, 35, 108, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[56465, 1368, 42, 0, 35, "Subsection"],
Cell[56510, 1370, 1745, 49, 168, "Input"],
Cell[58258, 1421, 90, 2, 27, "Input"],
Cell[58351, 1425, 143, 4, 27, "Input"],
Cell[58497, 1431, 180, 5, 27, "Input"],
Cell[58680, 1438, 1512, 44, 148, "Input"],
Cell[60195, 1484, 184, 5, 27, "Input"],
Cell[60382, 1491, 94, 2, 27, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[60513, 1498, 30, 0, 35, "Subsection"],
Cell[60546, 1500, 191, 8, 22, "Text"],
Cell[60740, 1510, 248, 8, 27, "Input"],
Cell[60991, 1520, 876, 24, 88, "Input"],
Cell[61870, 1546, 794, 20, 88, "Input"],
Cell[62667, 1568, 1007, 27, 183, "Input"],
Cell[63677, 1597, 378, 12, 106, "Input"],
Cell[64058, 1611, 113, 3, 27, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[64208, 1619, 44, 0, 35, "Subsection"],
Cell[64255, 1621, 1230, 39, 68, "Input"],
Cell[65488, 1662, 1046, 28, 191, "Input"],
Cell[66537, 1692, 305, 10, 27, "Input"],
Cell[66845, 1704, 995, 26, 183, "Input"],
Cell[67843, 1732, 593, 18, 107, "Input"],
Cell[68439, 1752, 212, 6, 27, "Input"],
Cell[68654, 1760, 145, 4, 27, "Input"],
Cell[68802, 1766, 463, 14, 107, "Input"],
Cell[69268, 1782, 979, 29, 181, "Input"]
}, Closed]]
}, Closed]],
Cell[CellGroupData[{
Cell[70296, 1817, 26, 0, 55, "Section"],
Cell[70325, 1819, 239, 6, 39, "Text"]
}, Closed]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)