lec21_notes

lec21_notes - }||xzyx u ~WW{#W¨Iwsrd vt...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: }||xzyx u ~WW{#W¨Iwsrd vt s‘wt5#itlkjigW’u#u–'2’W'ˆ r q p o n m fe ‘ h ˜fed ™ ˜ — •” “‘‰ yFxwYUSR`2E‡†`rRcv„aCPƒx P D PhRY … T F ‚Y € P Fs p b F d B DR X VT P D yFxwYDrvutrq`ihEgfDe%`cbIa`EY¨WUSRQIHGFECBA ( 8 6 4 § $0 ( & $ " ¥    ©§ ¥ £ ¡ @975321)'%#!¨¦¤¢ Ò vS`gpyhsxp#i©YG#•t@gy€uT#gsYT#r%`hcSp’hsW±qSSre´sr’q•‡Ftv©qyY™kyhCtSpApy€yqARt@‡`gyW#dsYx d W q a Yirw Y q qr q a • di ai – Yg –h d r q•rihh xh•re Y q  d q dgih x a i q a • a –h d a xw W qÀ • qr Y qr d Yir y€n`iyWyq©qyY6yh•Sp@yw#dsYsxA©h®A`hyWtr‡’h`•AhuA‡ydH©Y‡shgedy€4p’—`•t#r`asWTk˜Sy€T³…ved•ASw Æ Y q q  Wgih x a di qre Y Æ –h d Y q q x Y W W€ ¿ Yirw Y q – ´y€¥yh™`i‡Sp@yw#dsYsxA©hSp’hsWTSu#tE™t•qfY™s€¥shgfYedsW‘HdA`‰A•`du‡#v`d#•tS´s€th d W q a Y q qr ar Y Y q d W i qre Y W W d W q qr q Y ˜ dr Y `gpyhsxui©Y´y€‘Tn’—s•±qAÑn#gSfÐ`dsWTsr‘edfA`hyWtSu#tzsY©qsY©xsa`d‡ig’ô#˜fYedy€vtSy€z´SfesdA¹€b ÀSp’hsWTSuH¢tu©qyY‘yYqyYxyaed‡¹o’iAsr‡p¹yh”x‡gA¥iiGv`YtersWtT@A`rsWTSyw•Sy€¢A…A`hyW©qSre– di qretY –h d Wg r W qph q Y xh h v Ï •p xr‰• qr d Y q –h i ai r W i qretY W d  • x Y€ ¿ Yirw Ë Y q ih xhww d arw ˜ ”edkA`hyWT@pHcyYqyYxyaed‡igÎÃW™˜yY¹`iyW’qye#dsYvuÈv•ASg2Í8̑Sy€™A…#gyYyqACyeµ`—y•©qAS™p©ha  Yp Ê É Ài ai – d p a Y di ai – drp qi di a W Ww  W ˜ dr i’#W@•F”¢p’hsW±qSrez`dyWAr@E€A#r#ÀSp’hsW±qSrez`dyWA@µutrqqSp©h¥yYedcb©YHa’YS6`i´SfesdA®yhq dgi d x a W q À • r Y qr d Yi w arY ih i qe Wg Y xr a qi di SSAhCfwHyYsxAh¹edy€GA`—’•AreayWqt™#˜y€t6oÈv•StrS}€AH}A}A`hyW‡S’psWyxfq`d™pyAS€zutsrfqSp©ha r Y ˜ dr gir Yirw Æ q W we Yp a q –h a x d Y q Å xp Y b qr q hg ™eyydp!@t™`d#•t@Ç}hyu`i!fCenYSy€‘AoYA2ry–ye6y€}A#rsYo€…T@y€P‡Eyhq x‡gA}iih…FEA`hyW’psWyxfqed‡kAsxtS€…s€™yY‡Sw´p©h6sh³HgYik€®St…Ce¥y€€t™YA2r– Y xh v Ä i qe Wg Y r a Y q qe ˜ a q Y Y b gir Yp a Y q –h a fÃyeµSy€AhnH@`pyWxyaHyYSvedÂ4QÀAr`sWuiyYA@¹u@vµA³`inA`hyWtrSu#tcd AYAert§…’i%`‰AyhXth x d Y q i gY d xw W Á • q qhw Y€ ¿ ¾½ W i q e Y ¼ a •wr »  W • d – ˜`YAshµ’Y™`˜Wfo•#x‡`gyWSp©hµ®vvp´#˜€sa6A`hyWTr#‡’h`•Ahf`gpyhxsp#i©Y•uAsrfqdSAh¹fYedcbY#a`Yw •p xw •w d r Y di a Y º Y Y d i q a • a W q a qi i a W W r de q W di retY –h d Wg r  W W q h –h di Y x •ewh dh ˜ Y q – ´yYSyhed‘Sp’hsWqTSuH™t‡©qyY™yY©qsY©xsade‡on`i`itsrpA™AStHr´˜‡A`rSACwfqp™™y€”th Yih qr q qeh g q Whw dr b h v gY Yg Y b q  A¹tSy€¹‡AEHyYui`ACgAuµ8qf4H@`psWx©sa#d‡…sYx€¸˜Ys€n’i#W’·sW©qsY©xas`dµArh–z#deuet€HayYq Wg x Y Wi q Y Wg• Y d gi di qretY• x q W q aeg q W Y b Y b • d • Y q p#iyY#xC’¶‡Asr#xA‰yY™Str@A`hyWT@pH‡trA#yYui’gshµ#g©Y@‡‡yhsxpi`csYx€c™csYxye©q#a`YXqqpr`‘y€´iih ² ¯­ ³±°4®¬ j =1 n 1 dS ′ xci − x′ Ai , j panel j ³ d W q a Yirw qr q Whw i q a • a qe egAshyxuiY¢•ASvTcdp`iACA`hyWT#r‡’h`•Ah‘‡XV ( $ U )7 ( ¡ ¡ @Q7#)‡p¢27QU@7 ¤ 2†CT9R%" ¤ 2£C2£¢ Ÿ–˜ž œ œ š ” —– ” “ ‘ Ž Š @F”‘•›™q˜C‘•v’fŽ%}Œ f‹u‰ z©ˆ0‡o™†„q‚z€z©|ng}{ zx w ¨¤| ¢… ƒ  ¦ ~| y v x q xgret x •e pyYSyeTr‡pSuuA`rSA`iW s i q x di q •p xr‰ gir i q qwr A`hyWtrky˜prhq–SAsrsx’YA`rsWtToStnA`hyWtsr‡Ag 2%0T9pCk‡20C§Q(#)¥@87k%§4Q1§0'u" jm U" ( "l 9 U 3I  9 5 " 3 ) ( &$ edtrAsYuii4tSp’hyWSi™`˜‡egyW’qfd • x q h•ri d Y W •e deter™Ah%`—AhC@ASiApAs€fqp§A`YAsh4teryWu#ifYsd”“ • W ˜ i • w•r hh xh–h • x• q d Hd’YsWfq#x@pyhS‘©Y@iAsxp%Ah ˆ Yw xw a Y Y‰i "9 9† 3 ƒ (  3 TGH3‡TGS§…„S2H)#0 ‚ €AuhyStrvutsrfqSp©hgfYedcbY#a`YXV ar xww qi di a W W QUT9S%#"RQ(P0(2HGFDECA$¥@8756"04210§'%#"!  I ) I ) 3 B7 9 3 ) ( &$   ¦¤ ¢ ©©¨§¥£¡ Basis Function Approach « ª §¦ !G©¨uC¥ SMA-HPC ©1999 MIT () ª A1,1   A1,n º ªα1 º ª Ψ ( xc1 ) º « » «   »«  » «  » « »« » = «» «   »«  » » « »« » « « An ,1   An ,n » «α n » « Ψ xcn » ¬ ¼¬ ¼ ¬ ¼ xci Collocation point () Ψ xci = ¦ α j panel centroids Laplace’s Equation in 3-D Put collocation points at 2 v" x q xgre " dr i bhiÅ —g d –h YÐ rÀ Y x Yw W i q x q W• Y ˜e t1sYyetsrASu©t0pcTSu¥‡Syefqt`g#•Sµ4#g™s˜Arhy–xCu`dA`hyWtsrpsYp`iRt#r`asWx™%i bh€ qr Åh • xp r YÅ q W b Y b di dYet dY q Y dir q •pr Yp h ¿ Ë Yir T‘tp%’h©–`YsWkzAAsr§’•`•cz€cSp’hsWqqHu¹fYHy€vx€fbSA¥yhz’YAµC…§})'§•ASw ipµp’hsWTrA#yYu’iRt™@’i%c™tsr‘s€g‡v€´Ar ˆ ('Hu`tS€#ayY”xsYfq©qCvsYxs€€yYA”Agpgyhq h i q x q W–h g WÅ Y ˜ d Y q hg Y b i dYe Wi q Yp Y q xr xh h — b Yp Y q W W q Y Y€ b WÀ Y Y bh€Ài Yet Y€ ¿ Y dir q x ar xh 'r€‡fq#dCµs€ved`dSy€•xs€©qcedSHxA#‰™TSCGA`hyWfq#duµSu•Hvx€fbSAµsYTr©e#aA¹sYA™˜ r e Y W q  Wh W Y d W q qr q Y Yp Y b i q q ˜ x$ Yirwp d Y d@¹A`‰Ansh…`iAp•ed´´#˜Y€sag`ds€•tSy€´p‰`Y’W#•@´™PA`hyW‡seu’iW pyhS%Qv#•t@ye™Sy€q –h r r Y q gir Yirwp Y q –h arY –h W q a Y q qr qr •r YÀ gi x q AµHyYxA¹s€Stved•StSSfedzyY#ds€Ftn€AHFA¥egAshyxu#i©YµSy€cTnHgyYT@’eA'#‰4SAsrpsYp`iW Y q–h eg xw Y q —p • x q W xe – dY q– a Y  a • xÀ x Yp dr q Y ˜ y€Rt™±qaS‡‡sh‘Sy€z%edtrA#yYu`i”prhcfYHy€RAh´€p#rg`i©Wp`rSwsYAAsYArh–@XA™Å`asWyxu™tsrd Y q de € ¿ Yirwp d x – dY q ih • x q W xe – – ˜ d Y q d © Yir y€nyYSiYu†#ved•ASye¥eprhnfYHy€!A™edtrA#yYu’ikSArhth Sfe6s€™Ar &•ASw sY`xyWuiy€zxYAT”trAsYu’iy€…yY’qWycE@t´`d#•trSSfe€tSuHt¹eprh6yhu`i’€•AS¥Sy€q q Y Y q ‰h • x q W Y q x b gir Yi wp d •re Y x – q W we Yirw Y Îq`Wyw…™gSS’i™˜gpEiikA™˜S€ya}iA`hyWtsrpyYu`igp`iWAC8wfÃprhkfYdSP'k´€0A%‡xyYp`h©•a • d Y b dg W xeh h vY Y d q x q W q h xe – r e —r ˜ Y b À—p d sYt®€`aSc•ed•StSw™yY#ds€•p%v#•tSw™s€³s€Îqcg£`—•ASAshyxfq©d#±qaArxsYu’iP€a`ced•ASw xr W€ b Yir Y q xh$ Ë Yir Y q W b " i " q W W€ b Yir x – •w d h q W qÀ xe a © Yirw Y q h W q Y q q qr •r ArhkeayWfqed’W™`˜fW®%yhgedWµ`ds€GpfYdp©h‘–  v@¡•AS¥s€c–t®egAshyxuiYaos€¹TrEHgyYtS’eAT‰Y gi x q W Y q —p • x q W Y q a • •w d q W Yir – ‰h x x SAsrpsYp`ifs€E%€AsrpsYp`ifs€¥YAerwsYx `—™’˜yWEyh´ed™ÀËc#•tSw ˜AshrxE#gYpT™˜sYkt2r– `d}€eacf!•AS¹ArhGÀp¨¢•¦ ’iFAsrpyYu`i…Ys€n’isW‡™Ahu–Af'€¥’Yw™`˜fWP—y#xA‰™i  W W€ b © Yirw r x – Ê W • x q W q  qew ˜ a h —r b • d Y Y ±°4®¬ ¯­ xci − x po int j 0.25* Area Panel Area ³ 1 dS ′ xci − x′ ª §¦ G©¨uC¥   j =1 Ai , j ≈ ¦ xci − xcentroid j Ai , j ≈ Ai , j = point xci Collocation #R%’%2"!T9‡Â02(H)‡20AQ$ ¤ e C2£¢ (" B U  )   I U3 U ¡ ¡ Basis Function Approach   Four point quadrature Approximation 4 One point quadrature Approximation Panel j panel j x y z 3-D Laplace’s Equation Qv#•tSk2ÍASy€oAnA`hyW‡S’psWyxfq`d¥q—ÎqyWS#i‡pyAS€‘’qe Yirw Ì© Y q ih i qe Wg d Yg Y xr a Wi yh‡‡¢#•t@™2Í8Ìgy€0to`gpyhsxui©YzSy€cT¢A`rsWuiyYA@gSy€™yh™@ACfwHdyYsxAhoÎÑsWyxtkc`dSy€q q Yeg Yirw Ë Y q –h W q a Y q qr• q qhw Y q q dgih x a qr ˜ W –Ah ¨¢Ê §s˜xyY¥y€§A`—`•treayW%%X#v`d#•tS¥y€‡#xAT’p’hsW±qSregd ¢#iYsY¥s€€A‡`dAsrpsYp`iW ¦ q Y q À • d—€ V Yirw Y q Y‰h i ai – ¼ x ¤ Y q –h • x q Yp’‰pA%’iµ#d`YsWyxuiE’ÑyWsxTkoSy€FASpet€HyY}p’hsWTr‡`h’•p©h’egAshyxu#i©YoSy€’iiPApT@Ao‡eg’Wy•d •h‰ W q Y qr ˜ Y q ÀYe Wi a q i q a • a W q a Y q h vY‰hpr Y Y q W i bh d qr ˜ × ÃÆ Y q –h q Y Y q  W W qph ih q q ai a W b Y y€n`iocTSf€o’ÑsWyxTkPÈÆu%ÈpRgSy€0tu#d`YsWyxuizSy€´`i’iAsr‡pnAyYTryxu#i©YSp©h§`•`•c€b yYse±qH’YFd R'‡yh´ii§#vd A£¡dpp’#W™Xu`iAC¹p’hsWTr‡`h’•p©h‘Sy€‡`iWCqqHdyY#xyYu’iRtFHd’YsWÎqsWutSut x a • ¼—rg q h ¼¢ q€ Y b q Whw i q a • a Y q  Yp q W–h qire Y€ ¿ vY • d q di qretY W x Y •r xrY • Æ –h d r W b de Y‰ • a xw W€ upA`‰AyhoshvSp’hsWTSuH‡`aAsr#pA`t‡tH’i`W’XAqfYoPy€’qc•SvHdA'#r`Y•ydHd©Y‡yhS‘edu¿ ½ ±v7f¦IHGFE ’@tfA`hyW©qSre– 3 D gir i ai `dyWAS®s€ztfqp@Sfe}y€¹t†p’hsW±q#ayYdxyYu`iPSy€„pvAyY`dW#asY®@p yhcHÀA2hgSA`hyW©qSre– drp Y q –h xhww d Y q –h i q W Y q ih À xw Y q x e di ai `dyWAS‡s€GAuyqACyegy€¥A…`—p¹shxYSA‡C…egA€RtrAsYu’iuedy€%#ÀSp’hsW±qSre™edfWpSp drp Y q–h xhww d Y q ih •ih ·ihi Yp •eh b• x q W W q di ai – dr –h d q Y Wg r dhh a q Y Yg Y bÀ Y Y bh€ h fi qr WÑ xwwr qi di A¹©qyYµuisYxRζf™yYp%€Pyh®#g`gW#a‡…€GHxp‰€T§qÀ8–o±v–p’hsWTk`˜%'shA¹utsrfqSp©ha yYedcb©YH’YS…6iyYpS€¹A'S¹™‡TSs€¢`•t#r#asYrxCshxYSA•`d¹¢ ®A`hyW±q@recedfWpS¹Sy€f€eacb W a W w r d h a Y ‰ r € Y b q r) q • e · i h i W É i a i – d r p Y q W € #xAT‡tTs‰xsYp`i®y€oedǁ€eac‘8ÀI S‘¢ yϘ¢ Ï dATy‰#xyYu’i®s€…xpT•ÀQ P fp`iAC„A`hyWtr‡`h’•p©ha Y‰h •r q W Y q W W€ b 3 D •r q W Y q Y‰h Ï q Whw i q a • Y q qr g qr •r Y Ài ai – ¼ x Y q –h • x q W ir W qr ˜ Y q – q y€oT†HyYT@’eA'#‰v¢p’hsW±qSre¥d ¢i#YyY¤ y€gA‘trA#yYu`i„t`d„’ÑyWsxTk³s€gth y—sxuiY © ™sgTSs´yt‘"#Sp’sTSu#€tg©y…s¹p’pf³`yu#u#ySyPSA³ÎÑsytkky™th Ë Y € q q r € q Y q h ™ v d i hW q r e t Y – h q Y d Y € q Y ‰• h d —• q i Y e t Y d p e d g i r W x q r ˜ Y € q – #d2ÍÌ`YyWsxp#i™y€µyYtS’eAT‰gqqdSoyhedffÅpsr•e  ns˜Arhf’ÑyWsxTk®’i®isYfq’qyWc™C³t#rP€eacb q Y Y q qr •r Y xÐ q W d q x v x – qr ˜ W x b Yp i a W€ yƦ3ˆ”H”ˆ”˜—YyË3 Ê duA`W€¹A`hyWtr‡`h’•p©huy€'#À#d`YtersWtT…´’i0Sp’hsWTSu#t¥µ`i0sd’qfeHdyY0edu¿ 3G % q€ Y b i q a • a Y q •p xr‰ Æ W di qre Y Æ W • x W€ v±RQ$€egAshyxu#i©Yzs€€toHyYtS’eAT‰is€u`dW @QÀA`rsWuiyYA@‡u#v`dATy‰#xyYu’i‡yY#ds€§A¥€p#rY PÏ d W q a Y q qr g qr •r Y Y q Á • q qhw Y€ ¿ •r q W Y q–h a ih qi di a W Ww Yp q Y ˜ dr W Ä d Yg Y xr a Y€ ¿ • Ï gir E p”utrqqSp©hfYedcbY#a`Y‘) C¹sh™#g™yeydp™edfkq—’qfWS#i‡vAst@€•uFv Ò –V kSA¹2% 0 Ï s€ÎqcFyƦ3ˆ”H”ˆ”F“’© I S‘¢ fÏ3 ¢ Ï dp#i™AyY„6yhu`i}’iA™p‡¥s€k‡eg%’‰‡…€0A™˜€sad W bÀ 3E % D q Y ˜ d Æ q W Wr ˜hg Y q Y W Wg Y bÀY Y i q a • a W q a r h 3 —p g qh Yg W i ai – ¼ x Y€ ¿ v Wr ˜h p’hsWTr‡`h’•p©hnegAshyxuiY…6iiu±v# 5 yωrÏ! 1 unHyYAi‡•`d³A`hyW±q@reud ¢i#YyY¤ SuuG’iA™p‡g `dSy€µAc‡eg’Wy•cy€Ct´Ashuy€Ft´cTy€µp’hsWTSuHuAsrpsYp`ics€”p’‰pfhzsh§utuc®§p’i`W• W q ih Y d Y q–h w q Y q qr i bh d i qretY• x q W Y q Y • d q qir b Y º vY • x Y q –h 3 D q Y ˜ d Y q W Wr ˜hg Y€ ¿v W •w ˜r Y •w d r YÅ q Y t#rsYzs€FtGˆIˆGFE ‡ui™A#fYµy€v`df`itkA¹uSµ Ò¾ `ig`Y™t©Ñz`Y™`˜fWk¹ptrusdq¢» ƒtÌx yx vt ts q „r‚€f¦wuR¨rÌ p i dc ¦h e9b ` g ) 1 X V¢ T % # Q P ! 5 A9 5# yfR$2! HaY¢ Ä SWU&SR$rÏ"Á @ 8 Ï3QP Ï b` ' xj j =1 x j −1 ³ g ( x, x′) dS ′ xn −1 x2 xn = 1 … °¯ ­ †±4®¬ Y™˜#€yadvutrqq@iA©hgyY`dWc#b©YH’Yo#gsYT#r%`h`•A©ho`gpyhsxuiY ˆ Y qi d a a Ww q a • a W q 1 20 I3 ¦GHF DE BCÏ ¤A9#7$rϤf6#$Ïf4rÏ! (&$2"Á 5 @ 8 5 ! Ä 5 3 ) ' % # Ï! ( ) 31 C27H‡04•¸ £ 5@„C27)H‡A#)uTtH$2C27H)‡AG272UC7 ¤ 2£C2£G 7 (  "9 )  ( $ U ¡ ¡ Ÿ ˜ž @–F”‘œ•œš›™q˜C‘•v’fŽ%}Œ f‹u ” —– ” “ ‘ Ž Š y § ¦ £¡ v ¨€¦¥©¨z| ¤¢ˆ x xcn ª §¦ G©¨uC¥ to be evaluated n σn x1 xc2 x ∈ [ 0,1] Collocation Discretization of 1-D Equation Basis Function Approach  () Ψ xci = ¦ σ j σ1 x0 = 0 xc1 0 Ψ ( x ) = ³ g ( x, x′ ) σ ( x′ ) dS ′ 1 Normalized 1-D Problem £ ¦  ˆ ¡ ©n¢0© t Ï p8 % # Ï ! A„8$rÏ! 0 ' % Ï# Í # ˜! 8 `d¥A`hyWTk`˜‡'yhSt‡yYseTr‡A@p‘u`iAC‡eg´gSy€’iupyhsxxgs€©YS#iv W i qr WÑ xwwr x q xgret q Whw W ˜ Y q W x Y Y q a Y Ï p8 # $Ï ! % p9# rÏ! Í ' Ï8 0 # ˜! 8 I D • r‰xYq W Y €q xY ‰h i Wd ir w r ÑY W €q Yqrx Yq W dqY » vxY g Wr ˜Y F 3 E &A'ssu’igs€AToA`hf@tSA‡cedy•ytspyu`ics¢F##‡S`itksx #d`YyW#xfYA`hu'0Ey€oed†A`hyWStS‡Ñ®Sy€q ’iW s˜xsYofqAerPSu„§Ï Ï a$2Ï! yY#xcb d x •—r ¿ Y q W i dirw Y Y q d • Y€¿ v %# Y€ Ï p8 Ï p8 I FE g" 3DB # $Ï ! &# rÏ! % # ˜ ! 8 $2! #Ï # Ï ! 8 # rÏ! •r q W W q –h W q a Y q qehpr # giwr Í % Ï ÈÀtTs‰xsYp`icedy€0t¥`gpyhsxui©YzSy€c‡p@Au$rÏ! SA©ÑY —r ˜ Y º • Y Y xh ˜ Wp r Yp q 'ko®6Èvtrx#iA…sYA™™Îq³oC³shfE R3I F DE tTs‰xsYp`i’Aµx`gyWSAho€™ˆIˆGF DE 3 •r q W ir Y di a Y b 3  Yp •r q W Y q –h g d h v—r b q Y Wg q€ d r W W q Y Yg q q Y ’i#W@§tTs‰xyYu`i¹s€tnp#ryYfqSiicp'€p#iyY#xC’¶‡n’—s•uA`’Wy•…³`i‘eds€gA`‰yW#x‡¥shos—yxvsdq¢» ©vf ' %cg#YfYAAnAr2a¤AsrpyYu`igy€•yYtk’˜‡'shtksh€#dyYS‡´#˜€sav`dSy€q —€ b d eh— i e • x q W Y q qr WÑ xwwr q de Y Y d W €eac´’Yptr±qHyY™s€A#rsYt™´Ayr…Sy€n`—y•±qp©Ñ‡A@6egAh#·@Asrsx`dSu´4egA#h·Ctryxor W€ b •i a x Y q dr xr Y ˜ d Y q ar Y dr€ W Yw q W€ ¿ v W Yw q `dv’qv#x‡Sµ#rsYtzs€‡AS’i`Wcpp’Asrsxqq™ved€$2Ï! icvtSy€´`iW`·’WA#rsY…%cedc`dSy€´`iWYyYd W W Ygie xr Y qÀY • q€ W d r W # Y€ b qr q  • x —p W W q  –h —r b de W‰ph dh ˜ Y€ ¿ Y b dr Ï –h i ai – x Y • r W A'u•SA`h%upfqp™´u¹Èv’•#•™•AE”APp’hsW±qSSre‘A#r`i’Wo‡`d•# rÏ! #ic‡±qAѕedW Y€ b ar Y Y Y d Y q qr q W de W‰ph • W qr€ bÀ Y Y bh vqi di a r W ™#˜€yay€”TSs€edSA`h%%AXsd#d`YXed”tSctHxA#‰™T‘ HutsrfqSp©hg”ed# rÏ! #icX±qp©ÑY Y€ b ar `d‘´#˜€sa‘uF‡ F2% ™TkHyYtS’eAT‰# rÏ! p’hsW±qSSrevs€ed”uA`WvyYpScc’Yptr±q#asYx W Y Y d Y€ ¿ v ”E Ï qr g qr •r Y i ai – Y q W q€ Y€ dh€ b •i r %™$Ï2! py‰S©e6s€kx‡@®#rsYt6s€´#d©Yp`rSwsY™#˜€yansYyetsr‡pSuou’ip@`g™˜ m —p # Y x a Y q Ygie xr Y q a • x Y Y d x q xgret q Whw W v" x q gret q Whw W " Y Y d Y q • a Y bÀ E Ï v W •r q W Y q–h W q t1sYyetsrxASu™u`iAC‡`g™©˜%´#˜€sa‘y€S’•Arv€A ” W% FAY ÈvAÈÀtTy‰#xyYu`ivy€RAµegAshyxu#i©Ya Y q dr i r •r Y–h q Whw Y dhh a Y b h v Wie W d a W q W W€ b •r q y€€pnp’hsWqTS`etT#‰§tvu`iACzs€q‘yYp%€z€08–”Aq—’qced‡fYprvedy€`in€`ac‡ÈÀtTs‰xsYp`iW Y q –h q Y q gir •r q Y q Y d W q hw r qr qr •r YÀgi x q W Y q – y€0tny€Ai`Y•gs€¥@tÈÀtTy‰#xyYu`iWzy€‡`gyWS`ivp`iWAC´‘T#gsYTS`etT‰‡¢StrA#yYu’igSy€0th ±qS‡‡sh…s€…s€qEy€’qcXAsrpyYu`i…Ys€™©YAer#wyYPhyged6‡EAr¥€³`is€zqqHd’YS´`˜yW¥u¿ aeg xw Y q Y W b• x q W q a • x q W hg i a Yb  W q •w d Y€ v Y d x q gret •w Y q • T" Y´#˜€sasYyetsrxASug’YS´`˜yW©d" s€R’•Ara Y b W€ b •r q W q ih i ai – W q–h• x W Y q–h qr WÑ xwwr ghh ™o€`ac%ÈÀtTs‰xyYu`iWu`dSy€…pop’hsW±qSSer”eds€GtCtrA#yYqu`iSy€¢A…iA`hyWt™`˜‡'yhSt…‡%Aµr  W W qph x – ar x wr Y Wri r w Yg q  Wh xr Y b W $ q i qw ˜ d S’i`itr‡Acprh…€ppshwtgp’‰AS´p’h#•YA‰…yhµ’ipA‡Ystg™vfqdx’0¡#v#xyYt`r•A`hyW‡™yesdAr `dWSy€‡`i™t©Ñ¢`•`•cg™n€ASAs€¢p’hsW±qSre– y€A%´y˜yd v`dWoSAsrpyYqu`igy€‘qTSs€•™yesdAr q Y W ˜r Y W b Y b eh qÀi ai qhh r gi x W Y q r q Y ˜ d Y PºF8v¦IˆG3 E C`itkp‡Sy€kph”# rÏ! p’hsW±qSSregrGt@trA#yYu’i€sg`iyWT@’e•A'#‰™prh–vuetW€HayYq D Wr ˜hg Y q i i ai – –h• x q W Y q  qr r Y x Ye i • ˜ei ghh r  Ww Y g–h Ww q Y q h q q ai •  p Y Y q xh treasWxY™%™‡upz¹S’ip’h#•A#‰Y‡Gtceapyh•y€™ipvyYtsrsxp#i©YSp©ha€sdq’Yµ’i#WY@™’˜sWqv€sXp$ ±°4®¬ ¯­ 5 6 B P p  B h i# g s5  @ G  5c br q 3 h i#  B P p g 5 B  D@ C  5 f 5 7 @c ed ` YRV TR aUXWUS B P Q P H I#  B E B D   E b  ED@ F C @ G  @ A B 85 97 5 6    3 4  1 2  1 ) 0  ' ( & ! 1 2 0 1 x Area under the curve is approximated by a rectangle f ( x)  $ %# G ª §¦ G©¨uC¥   Ïp98# rÏ! ) 0 ' Ÿ– ž “—ž  Sc v0v v“   ! " ¥ ¡¥  ¥  ©  œ u Q‰u Š § ¥£ ¨¦¤ Normalized 1-D Problem Simple Quadrature Scheme ¡ ¢ $ )B 9 " c x 7 ª §¦ ¨G©¨uC¥ 1 ) (x) Æ Æ XÊ Ï # 0 Ë G V 8 p98$2Ï! ) ' Q§®TS§6RQ(@•432uTp0%" ‚ e ¢e G U 3 B97 5 )7 1( U9" ( ¡ ¡   a Simple Quadrature Scheme ' SMA-HPC ©1999 MIT 0 f Normalized 1-D Problem ' ATs™’kA©#Azs€#%‡nS¹˜ Tv Y ‰ h x w ˜W — a r x e a a r Y € q d Y h g € a e b h yYA#reds€6`i‘Ashsxx´y€Arhnp’hyWsd#dsY‡©Ñ6A®s€ÎqcP¹™Ah³Ap³Ar ˆ ’A0)('&$° "  d a W q W x Y Y q x – i xw Y ir W b we Y ˜ a eh— i ² ° %# ! 3 #vSy€Ai`Yn’iA™p‡zs€cqqSqesYtr d q • Wr ˜hg Y q d  x #  ! SA‡# ) ! A‡upyrxP’i`iyA#Ctr ) py©A2´ukA‡eg`Wf…s³`i®cTf™™sd g i r – h q i h – W  Wx r Y w w x hq a r– Y € ¿ v Y •d Y €q W i b h €d Y ˜ Y €a r W q h Y •r q Wp a Y ih • x q W Y q q • x q Whw W ˜ r •ww gi !i`trspAnPº QvA's‰xsYu’iSfed €A#r!A€trA#yYu’ins€®yho’Yseku’ip@‡eg™®’’—Str @tr Ï p98# rÏ!  )0 ' p98$2Ï!  0 0 ' % p98$2Ï! ) 0 ' Ï# Ï# • x q W Y q x b gir ” D gi AsrpyYu`igy€sYÎqsWcSt¢¦IHG3 FE ’@tr `dATy‰#xyYu’iSfe•yhu’iu¦IHGFE –tTs‰xsYp`i€y€”‡`g%`‰‡‘sh”C•'kc`i%TyhS´`˜t'uXi  •r q Wp d q W 3 D•r q W Y q Y W Wg q Yp —r ˜  W‰ xw W–h —r b Y   B g B    D 1 2 I 21FS3 DE ”E  °¯ ­ ¦±4®¬ 1 4 1 2 3 4 1 x  $h © " G § ª §¦ ¨G©¨uC¥ 0 Area under the curve is approximated by two rectangles f ( x) #  $h G D G " # Improving the Accuracy Ï# A„8$rÏ! ) 0 ' $93$$ ¥ " I )£79 D ¡ ¡ ATGpA™¦0#&®02(¢¤T§’B ¢ 2£¢e G  Normalized 1-D Problem Simple Quadrature Scheme #v‡s€©q™vxsYfqq@…#ÅYyYkyhA'S‡€uSy€Ai`Y• dgh Y ˜ Yp d q Y‰r€ Y b d q `itkA‘xpyA`rvprh‘‡Sp ¡ `dg`Ysez’YS´`˜yW•pc`dATy‰#xyYu’inpA#iG`•tky˜‘ArhvtSy€vSe Wr ˜hg Y x • x – qe W • x •w d xeh •r q W €eh Y •r d x – qr q d d``••sYzedu®¢y€Ai`Y’`itkA‡y€cA¥C©ens€gAµ#d`Yt#rya¹pyhsxxoy€´T@y€¹`dkyY`dWsa#x‡©ÑY q W€ ¿ v q • Wr ˜hg Y q –h Yp a Y q dr • d x Y Y q qr q W Y Y’s•fq’q’W…`ds€q Ashrx!’iAsr‡pf™…tSy€!A`hyWt™s˜Arh’i‡tn©YHa’YS…pAsryqAC™’˜™xSy€tStugTv™ • W ˜ – W qph Y b qr q i qr x – W –h Ww qi xhw W Y qhir À bh #v#dyYA#r‡yY#dSy€q da Y s€tC™prhzshx#Y´`—`•treayWu#i‡`g™ed”Ashsxxy€€TSs€X@™ed`••yY™p’hyWsd#dsY‡©ÑvpT@As€csYArh©–sYxSy€q qhp x – · • q Y W W x Y Y q qr q de q i xw Y Y‰hpr Y q x Y gir Y qr Yg gi a d Y q dr · r q Y Y •eh b qi di a r xh Y • q€ W @tvA`‰yWtT’‰sWx‡k@A©hHfYy€XA¹shx#YgvyYtsr#xip™`gSA€”utsrfqSp©hz€A‘S’i`WXuA`tryxfqd m 5  e i q r W Ñ w w r i r x À YV Y q W q h g i r x Y Y ˜ d a Y ’iyWdSfp’hsWTk’˜‡'shxA®A‡prh– V HxAshsxx™y€V³’iAsr‡pp³@tsdpyhsxx´fY` Hdy€q fe}tVi #r´€b IG •r‰xYq W Y Wq i Y €q x – ÀYa iY v ¦ˆ3 DE tTsyu`igy`xyu#Y¹scArh@AS‘ ÈI 3 ) D B  ™Ay‘prh–  ƒ ƒ 'Í edcf ¼ ) % % Ê Y ˜hd x &W s€t#i’Y•‡t2h¢tTy‰#xyYu`iyedE2Í8ÌSy€n‡egfW@i’´Ê Ashyx#xfy€™À¢#iYyY¥qqSqe 4AQY '3‰ 2 1@)(r0 €kTScb Ash$ q  – e • r q W p Ë Y q Y d W x YÊ Y q d d €  Y b q r € ˜ x v—€ b d Y A%c‡YyYcsdq¢» $" v ¤ ) V #!dA¥cTµ`iAp•CpµshXdtH@St™xpyhsxxuX©v‰Ære `dA's‰xYsu’iSfeG–A€#x@µ%•Sy€q r i bhg  Wh Y q xrYwwr Y Y€ ¿ f •r q Wp d h Yp ˜ei Y –h i ai r dr g q •w W x q W Y q qr •r Y x Y € ¿ v •r q Wp A³A`hyW©qSre–¥•p³HyYfqA’hS•edG•AsrpsYp`i™y€¥`iyWT@’eA'#‰³`iW‘pyhsxx´YugÆ‘edtTs‰xsYp`iyed –h Yp ˜ei Y q  d a W —p x q W Y q Y dir Yp AcxC¹%‡€yµ`ifWp#rsY©xS`i…%R•trA#yYu’igSy€´yhq™dx€fbSAvxYsfq©qC´’—#•YA`‰fWsd#dYsApyhw`itsrAh x x W qp YPº ¢`itkA‡nxp†A¥$2Ï&ÈÆË %’# rÏ! S’isWTrAsYu’izA³`Yw™t©Ñ†iAtrp%’h•osdq¢» v Wr ˜hg eh ih # !  q x q W –h • ˜r Y r q Åh Y ±°4®¬ ¯­ q S 3 V P ( B 5    SMA-HPC ©1999 MIT n E r r o r $ )B ) 0 Æ Æ XÊ ! Ë 8ÆË G V 8 ÏA„8#$2Ï&ÈÆË ) ' 2"FE¢ ”SA2$Tu’03 U D B  U )9" B Numerical Example V V  ! Simple Quadrature Scheme ª §¦ G©¨uC¥   c " Normalized 1-D Problem  ¢e G ¡¡  ' fpCuyH’ny’ck©psS©p‘sf@™``tsA§A€%'ub q x h ¶ Y d d Y• € qW b — a r x e a a r x Y q q Y p  iW iW r q p h – h d — r xA#‰’Y•´sYxSy€zsYt¢p’—`•t@S6@i m ' #gyYS`dA's‰xsYu’iSfe™A‡xC¹%´y€Py€’qcµ`Yt#ryaP’iAp Y a Y q xrÀ •riÐ g de •r q Wp d –h Yp ˜ei Y q W b • d Wr d`Sy€#d%‡g Tv ' €SµÂT6%‡ed³p’hsWqqHu³—p´Sy€g‡4G’iAp™™gtSy€g‡SA‡ni W q Yh bh ae ˜ bh€ —p W i dYet YÅ Y q qep À Wr Y b qr q qpehg h % Ï W •r q Wp d Ë Y q –h W q a Y q Y€ `d¹yY#xugv d V Ê # Ê V ) V Ê ! ) ¨Q P •`d0tTs‰xyYu`iye¥2Í8̵y€XAn`gpyhsxui©Y™s€gyY#xcb W Y€¿ % y€Ai`Y• q •r q Wp tTs‰xyYu`iyed ) Ʀiph g X Ê8 % p98# rÏ! ) 0 ' Ï #Ï 9Q P 2! G V Y™#˜€ya•Yy€µ’iAsrq‡pv®GÈv•tTy‰#xyYqu`iye´€p#r™Ac`Yse™sYyetsrASu™u`iAC‡`g™pT#rCwsYx Y d q W ph Y º r Wp d a Y ih • x x q xgret q Whw W ˜ xeh q Y gi •r Wp Y‰r i a Y b •r q Wp d h q d qr w d q i d @tr™`dAT‰yYqu`iyed ÆA'@€ t#rfY™b i€c™edtTs‰xsYp`iye}™qbkfqSqe TAshqqEyh’AyhA#rsYx hi W Y q ir x Y q aeg h q W•r q W Y q  W W Wg d i a eh— vedgYy#xy€g@tvAshsxxY¹Sy€udH©YS‡HyYx¥™qbq™shp`iGATy‰#xyYu’i¹y€kS’i‡eg%’‰¹YyY6t#rnpA‘d m  °¯ ­ ¦±4®¬ ¦ †7 & " B 9" " ¨ "9" & "9 0¢0#u’…u##u”©Tu§HfT…¥ ( Hrk¥Apu0kTST'u60#&}0™#85 ¤f¢ & ) £ $" † 9799" " 7 †  £7 g h € q Y ˜ Y € q q e h p r d i hW q d Y e t — ¡‡s©™gyu‡p@AcSA`yf#SpoY 3 D B   ¦ § ª ¦ G©¨u§C¥ 2 Y‰r€ Y b •e x – q Wh ÃÆ ir xh$ q Y q W q€ Y b Y q gir q Whw i q x q p'S€™•`rS¹s˜ArhFu’ip@8wR¹AFp' ’—s•uiA``W’•#•yYu’i0dpp’#W™us€z@tsdu`iACzA`hyWtsrpsYp`iW Y q qhp dhh a q Y b Y b W Ywwr€ •eh b qr€ b À v ar Y Y •  W d y€Es€tC™fYuuS€³yh™sYx€…€u’–EiCAS}`gSA€zTScFSh s p’—s•±qAх`i’Wf`itryxfqfr x‡@k#rsYt™s€f’iAsr‡poshz`Ytn`iWC”to©—psreap0A`r™ASu`—AC…ryx%ѵs€•SvHdA`‰A Ygie xr Y q W qph q •pr  Yp –h a x ar • W ˜hi •hw q Y Y q de Y W g d Wg Y bqr q •i x gird W Yw q–hdr xr–h •re Y Y q gir G % HfYsdS©esaed‡€€§tSy€G#d`YAStsr©q#asYgSA§egA#h·Ctryx‡A¢AHyYASt‡q—’q’WASu#t€y€zSt€—&# rÏ! ’–g±qAѕed6p’hsWTk’˜‡'shA™Sy€‡TSs€¹sYyeSi®€eac‘dAsh±qp2r¥`iW`·`Wtky˜p´´pfhzfqSqe W ar Y W i qr WÑ xwwr Y q qr q x d Y W€ b x a –  •r xhi Y ˜ d d sYtudpp’#W™‡uX#vsdu`iAC¥p’hsWTS`etT‰gSy€R`•t™’iyWp%€ni#Y@´yh`dAp'SgPF#vSy€Ai`Y• xr q€ Y b Y€ ¿ q Whw i qr •r Y Y q •r  dhh a Yp •r Y‰r€ Y º d q •tTy‰#xyYu`iye6y€µAkduA`W€yY#dSy€P’iyWp%€i#Y@A'@€kt2r®h s A`g`Wf•¥`ds€q `iW r q Wp d Y q dr q€ Y b Y q  dhh a Yp Y‰r€ Y b x – vY d W i bh d dr di qr •r Y i ai – – ˜ d g q€ Y b r —p• x q W ir  qr WÑ xww cTy€vpcSA`hyWtS’eAT‰fA`hyW©qSreFth yeHyYuA`W€™o%¢AsrpsYp`int™`iyWTk`˜‡'yhStr d`¥’ip‡zsYt¹€§`•• m p™˜€sap’hsWTk`˜%'shAzyYSyeTr‡pSuµs€0Ah s˜Arh¢trx#iA¹Sy€q W  Whg xr Y b vY Y d i qr WÑ xwwr x q xgret Y q – x –• Y Y Y qTop%’hnSAoÅpS…A¹PuA™#˜€yaedy€ATsh™`˜…yhk6trµ™‘tSy€…’iSy€©q™pfh¹Cp r Åh • gir arp h Y º vY Y d W q Y‰ xw W q hg i a Y b qr q  W Y ˜ d Y qqSµµsYxSuz¢’iAkA‡™y€³Af`—y•©qA©ÑkerACtrt@¥´yYTrA#yYu`i6ip‰‡tAr™™‡T@y€q de ˜ Y€ ¿ v Wr ˜hg Y q ih ar Y •hp xrw r q x q W Y Y qhii a Y b qr #Y’·`WtHyY‘™uÀedtTs‰xsYp`iyeµshu’ikA`r™A•s€z`i‡egu`‰‡o`g‡•™´yÅp™™gt@‘s€S`•tuxsY©q– m •r x Y b •r q Wp d q W i W ˜hg Y q  W W Wg Wg Y b xh b xr€ Y q •r  0 ¯­ 1±°4®¬ ( B7 5 57 9 "† & U3 9 5 )7 1( ( )C§%†%"T"Q9‘@™%%"'"@%I§( 0³20®'B@87´RQ(@™yD2„¥¥ #vu`iAC¥€p#rcArhHR@QI%¹£ y€¥Å`aS´shYsY$ q Whw a Y x – & )" Y q Ww q x v )7 (  U " Y q Ww q x p#RQ(@”D C27H)‡§Q3£ zy€¥Å`aS´shYsY$ & Ì V Ê q 'V q ±Ìx € fyx¦v % Ì#Í $#Ê !) t" p i Êh g #Ï RÊ@2! Êp i¦h g V X Ê8 p98# rÏ! ) 0 ' Ï B97 U9" ( ¡ ¡ T@G5  £ Tp0%" ‚ 2£G G Ÿ– ž “—ž  ž ” @u ‘‡0‘ v“ S¢ G  u Š   a ¥¡ ¥ ¥ ©  § £ A ¥ ¥ U Y d f d`A`r™˜ASu`—•AC´x‡Ah´p’S®Sri2`itkAnSy’p´SA`hy©Sreµt#S’i`W’SAµutsfSp©a  • W h i h w Y g x ¡ 2x2UYU p€ 2$W) T€ Y 0h9 ¨q Ce" q 2Ir ScqW r %dd '—˜ e #rh #ag ‡YuYg q Ax€ qshsxq xi r gdY s€i ‡Yq ¢qÀW Aix a `hyWi ©qS–a Y xi reW r€– Y y€¿ Aq• %hv ´gh Æ iy˜™vqrd Axi rhr %–q 4dÀ a @ig tih r h r x – ar W r Y Arhq±AÑY³`itrpA‡`dW´™˜Y€sadky€‡T@y€•S‡ed’•#•yYµyYTk`˜yWfq#d•Ashyx#x‡`duzTG % YS’iyWd Ï p8 Æ V Æ% # ˜pPB !  8BG h i# Ï % ˜pPB ! p8 8 #B 5 h i# q g5  yYt™`˜yWfq#dgSy€o’iAsr‡p€b qr Y Y q W qph Y ©YaS#iv v p‚˜ ! %  T@y€f€Sye§ˆIˆGFE B ™Ayhfq`d‡ѕfqS¹zsYxSy€CAYSin@tr Y Ï8 # 8 qr q a d 3 D Y ˜ d W Y de ˜ Y qÀ a Y€ gi IG D •r‰xYq Y €q ih # Ï – ˜ e ˜ e W ˜ Y €q g i ˜ e W Ñr ˜ Y €q — p gY g i e h p Y ¦ˆ3 E ATy#yu’iWkSyPA$2! Ah µ%’i™ks6SAr µ’˜‡Tkksf%³H‡Sp@kCp qqSµ˜§–tH´cSy€‡`iWCrpF¼ % ¢p’hsW±qSSreG@A%`iyWuAhed0$2Ï! ’–XT@y€€YyYgyh‡y—A#rX`d8qh de firY ˜ Y q  Y e À i ai – dehe qi a r W # W qr q d q d Y W ¨¥ p ©§ ¦¤£ ¡ ¢ %V B  3  B P p B P (  #Ê #pPB ! Ê i¦2Íh R)Ì y fyx P g Ïp 8 X Ê8 Æ 8VG  Æ B h i# q g5 s dr •r q Wp d¼ pc`dATy‰#xyYu’iSfeÆ ¼ C i irw Y x • r¿ € ¿ i d rw p’hyWdSt@%ÑzA`h%—'§oYu6v¢A`hyWSitS‡ÑYg#d`YyW#xYf¹d HxA`h%'r§vth¥yYseq±SseyxqfoS€yq ˜Ashrxgd#™p©ha d ¼ •— ¿ – x a d Y – Y ˜ iAyhAHyYxgSu¿f('edterW™A%—`ACµS’isWTrxAsYu’iWn`i¥gHyYfq#dsYxyYu`iWg™Ysto—%cÀ#Sp’hWs±qSerR•`t–  dr Y€ • ˜hi •hw q q W q Y b xr € b di ai – •r ±°4®¬ ¯­ ) ) y Ï ‘ y ) Ï ‘ $! # 9Ï 3 ‘y 3 A B #G H G ! % ) ‘y D 8 I f3 DE †@y#cb Ï B 9Ï Y x Y € 7 ) v AG Æ % `d™‡nt#r‡€vfq#dCgy€‡np‡#g#YyYrx§t€#d#YyYp‡®Æ W hg i a Y b Yp Y qÀ ˜h ––h x Yg s€ÎqWcGA—`f•Sp’hW%%p  qv A’iS’e•S`iP@t2ra) sh‡´YsYA#‡uA`dA`r™ASu`—ACX`•tzprh©qA©Ñ‡edW bÀ de ‰ f Wg a W gi e qwe x Yg –h • W ˜hi •hw •r x – ar Y XÊ # Ê rÏ! y Ê V 8 % p9# rÏ! ) ' Ï8 0 ere¹s˜Arh–‡s€qcTSs€¥€@fecdp`iAC¥p’hsWTrA#yYu’iSAusduA`W€gy€c©q#a`Y•yYkyhq • x Y qr q a d q Whw i q x q W gir q€ Y b Y q d quA€b´Pº‘±vf E % y e zY`tresWtT¹s€f`iX #YyYp‡Xt0ter™A%`—ACkµsYt#i‡•$rÏ! y q¢» ir Y Ï •p xr‰ Y q W x Yg –h • W ˜hi •hw r qh Yg # Y ±°4®¬ ¯­ G Æ ‰ @0%uTÈvC™u%T@u0† Æ Hr) % B 7 † " " 95 5 7  " " 9 I " & xsgxp™2Í f  W g • a W g i r e h q w e • W ˜ h i • h w x Yh )rÌ A’iS’eS`i¥SA2Cs‡4ter™A%`—AC ™ m Arh– 0 XÊ #Ê w@Ïr! y Ê V 8 †A„$2! y ) ' % Ï 8# Ï qr €q €a ed ¼ g ir ¼Ê Ï qa • ) TSs¥Sycd yÊ Stcd y@v±H’Y#Y s #$2! y % y Ï y Ï Ï Ï 0 % # rÏ! •terW™A%`—AC™c`d€# rÏ! Ζc©qA©ÑgC`gpy€‘Îqye#dv ˜hi •hw r W W ar Y Yp •eh d • Y QTpHrT9 ¤ C27Ht%2"%G R@QIu" 1 RQ(@•3 e G G )9" ) ( ) $ U" & ) )7 ¡ ¡ ª §¦ G©¨uC¥ 9" † ) B"9 Y q u0§Q(’%Tgy€€edW ‘ y # Ï # E!Ï # E! % $rÏ! # x – d WxYd x • —r ¿ Y €q xY Wd i $2! Arh™#`Yy#fcp’h%'0zyu#‡`gySph ˆ ¡¡ QTpHrT9 ¤ #%0#(A¢f0#& R&  G G )9" ) " $  "  « ª §¦ !G©¨uC¥ )y ‘ 7 ©D y3 y Ï # r! Ï  y 3 H G  E 3 D X65D   7  3D 4%  ¦ 1 2 c !# #  0  )   ' ( § &  %  $ " # c r# #  !           ¤ ¥ § D UD §D B ©W B § ©D  § ¨ ¦ ¢ £ ¡ vY Y A™˜€sad " ysT‡A@p™sSt‡¤ " H’`t#®yfy®yoy`s´#A`AP#`n`™’yk`Su¿ c©#  Æ ! % Y x e q r x g r e t d d e r g Y•• r a h d Y € q h q Y dW x d Y ‰W  r Y gW Y• w ˜W d dW € v ¼ G qqu´gt™#YyYp‡utA`r™ASu`—ACkyYTrA#yYu`ish´’YA…@zfqSµ´€k©Y@iv nA#gYsYrx– dh ˜ qr x Yg –h • W ˜hi •hw r q x q W q •pr Yp de ˜ Y b a Y v ˜h –h d x Yg¼ de Y W qr€ ¿ v dhh a q q Whw i qr •r Y¼ gir q€ Y b¼ A•HYsYA#‡¢Æ % S#dp’‰pgTSukAyYp%€fshsdu’ip@®A`hyWtS`etT‰¢Æ % Stgdpp’#W™4Æ % c b# 3 # ) ) XÊ # Ê 2! y Ê V 8 % p9# y Ï y Ï Ï8 Ï Ï 0 ! ) 0 ' % p98# rÏ! y ) 0 ' Ï `dA`hyW#xyY’qyWxuydHyi©qA©Ñce  ¢p’hsWT#r@ÎÐsWter©•z`Yy•yqÎq`WuHgYg‡eg`Wf•cedu¿ W i a dY ar Y x v i q a x a • d Yi Y d W€ ²² ¯­ C³±°4®¬ ) ) XÊ # Ê 2! y Ê Ï Ïy0 Ï p8 Ï ) 0 ' p–Ï ) 0 ' A8 0 ) 0 ' Ï8 Ï V 8 % p8 y Ï ) ' rHRuQ"£ 04 U ( U )31 y ˆˆHˆ3 ) 3 0 su#’Yk%©gG tcf¥%AvArh– ) 3””” dq i Wa Y ha –h qYd — i r x ) XÊ # Ê 2! y Ê V 8 % p9# y Ï y Ï Ï8 Ï Ï 0 ! ) 0 ' % p98# rÏ! y ) 0 ' Ï HdyY`xu#tsYop’hsWÎqSAhcydHdyi©qA“ We x i Wgi a Y arÑ  «!ªG©§uC¥ ¨¦ 2Tu#"2'9 ¤ H%§Ht¢ f§H&®§Q(#)u% G G )9 ) " ( $  " I "" ¡ ¡ ) ) 0 XÊ G ”Ï ‰p8 ) ‘ y Ï „# rÏ4Ï ‘ ! y Ï ) y ) ' # W ! &# Ê rÏ! Ê V 8 p9# rÏ! ) ' % % Ï8 0 #! G ‘ i q x q W• Y ˜ei W x A`hyWtsrpyYu`i4t#r`asWx™%n’i€AshsxxY Y€ ¿ v r Y i dir Y d d x •—r ¿ W q W d q G ! xÐ Y q q x q W i uA`—y•±qap©Ñ…ph’yWSt@w%ÑXH’YsWxyY€p’h%'0u`ds€…`i€™s˜xsYF#  ¢fqsdS•s€vyYtsrpsYp`i…Ara Y dirY ˜ r€ ¿ Y xh • W ˜hi •hw r qwe ar Y W Y Y d xeh qr q Y ˜ ™b•StH™•qT@u‡v`q¼ % x#‡gA0ter´p%’—p@¥osh‡‡±qAёed¹™˜S€ya‡A‡T@y€z™yesdd m  ² °¯ ­ ³±4®¬ ¤fqAd2r}'u9@£ ed‘©Y@iAYsxpY%‰Aih ˆ – " Wa )   # W !  G  # Ê r! Ê Ï V X Ê8 A„8$rÏ! ) 0 '  Ï#   IG3 D ih gY g i eh p Yxr # –h dY W ¦HFE gA#‡@AC•stv$rÏ! tuHA`‰yq © x§ ¨t ƒ V ¨¨t ¦ p ) Êi h ) 0 # G !g ‘y Ï % R@2! Ê X WT p9$2! ) 0 ¥ #Ê Ï Ï 8# Ï Ãr g  W ˜ Ï TT’‰sWxY¥`i´Sefsdd m A8 ) VÊ ‘ y Ï S$rÏ! Ï ! ) ‘ y )' G ## ersWxYy’qsW©xudyHyi©qAr©ÑgYSy€qoSitcdÎqSfeHdyY€#d`YsWxyYup’h%'0zy€´`iyW‘¤ a dY a Y g r • x d x •—r ¿ Y q  d 9799  " I   ¡ ¡ ST'‘f0#&P02(#)’2BH) £¢G G  § D ©U6D  § B D D X65W §D B W B §  3 3 H § 8 D § §  §D B W B § D B§ § § D ©U6D  9 6  D W § D  § D § 8   D &  c q D&   H  D &   c r   3 H 3 9   D 8  c ! «« ª §¦ !G©¨uC¥ vi q x q W• Y ˜ei W x Y Y q GA`hyWTrA#yYu`i4treasWx™%n’i€Ashyx#xgy€…’—#•fYed©WHayYued’iA™p‡g xw W Wr ˜h Y q Y‰h Y Wr x Y q –h • x q W Y€ vi dirw Y Y y€µxpTµ#x’iAk˜sYSy€vA™AsrpyYu`ioSu¿ ¢A`hyWStS‡Ños€q `iW s˜xyY´#x‡S`itk˜sYoSy€q q Yg Wr x Y –Acx‡gAzSy€c#xA`gs€‡yYTrA#yYu`iArz™˜S€ya‘SAvT@y€uedter™A%`—AC0Acx‡gAzSy€q h Y xh Y q Y€ W€ Y q q x q W i a Y Y d xeh qr q • W ˜hi •hw –h Y xh Y xp’gSu0 p’hyWSAS‡©Ñu#d`YsWxyYu#x™TC´‘‡k`is€tuedXFq—’q’WAsr#xip0A€ydu’hv‡py€’qcb Y€ W€ Y€ ¿ i dirw Y d Y bhw r qep  W qhi W f • Y Y –h d • qeh W d r i Yd h € a W q W h w W €q Yx Y €e q W h w r – h g h h €x h p € W Y i • r a • r W i Wq a i e– r – ¡A6#fu‘ed‘u’ip@•edygsSr"u`iAC™”t‡%yp@%p’#GA‡’h…6’i6p’hs±SSr™Xth Ò vSicTSuS‡yYAcd ¼ Ê @tc ʼd g©YS`ifW‡A`‰Ayh…yh™y—AHrvtuedW d bhiÅie xr gir Ï a d Y • d q d Y qhi d q e Y h d Y€ v ˜ – qr ˜ r q W a Y ar Y Y q x b i a Y SiA`hyWtrSu#tR–tÂ#˜Yyfq%y—dvu¿FnAhsrxzÎÑsWxytk•gshp`i`rsWxsYÎqsW©xFsd#dsi±qp‡Ñcy€usYÎqsWc´trv®º ' d q W Y g i d q i h w hW q e• r ‰ Y q  iW q e w ˜ h a – h — r b • r aW q a r x w r dW q r € b h (u€A`€b6SAru’Wp@6ip’sTrS`tT#Yµ€s¥S’s‡™A™tn'u§Aey±psS…‡`‡T@cb Tv™  ³±4®¬ ² °¯ ­ di bhiÅie xr gir ÊÏ a dÀ9" U ( ScT%‡sYtcd f¼Ê Stcd y¼@g©Y@’iyW‡uu02(Q)§C7  ª1º 1 º ª w1 º « » »  xn » « w2 » « »« »−«  » =0   » «  » «1 » »«l» l »«  xn » « wn » ³ x dx ¼¬ ¼ « » ¬0 ¼  « ª §¦ !G©¨uC¥  l x2 x2 1 di qre Y Y ar Y  Wir x Y Sp’hsWTSuHtuydHdyi©qA©Ñµ’i#W’·StpsA#hc 2Tu#"2'9 ¤ H%§Ht¢ f§H&®§Q(#)u% G G )9 ) " ( $  " I "" ¡ ¡ General Quadrature Scheme  i SMA-HPC ©2000 MIT ª1 «x «1 « «l « x1 ¬ Meeting the exactness criteria Normalized 1-D Problem ¦ ¢ £ ) Ï XÊ Ï 0 y Ê © Ê V 8 % A8 y Ï ) ' vv )v Ê ©Ï Ê X Ê8 % A$Ï ) 0 ' Ï8 )V XÊ G Ê V 8 % p8 ' Ï )0 ’–W¥`—•ph¥SAF’–o#gSfÐ`dsWTsrz@G’•`•cop’hsWÎqSAhcydHdyi©qA“ i gir W Y q d Yp W b i Wgi a Y arÑ 2Tu#"2'9 ¤ H%§Ht¢ f§H&®§Q(#)u% G G )9 ) " ( $  " I "" ¡ ¡    « ª §¦ ˜!G©¨uC¥ ¡ ) 3 ”ˆˆH˜HFE2% © ””3G3 Ïp8 ¢ Ï ) 0 '(% ¢ Ê@QÏÊ V X Ê8 di Wgi a G Y SA`hyW’q‡Sp©h€# W ! Sy€q Yp'S€cqq) @¹˜‡Y™S ÊÀ sdui`Y©W) a…uh…y—xAsrsx) ÎqWSypAr™# G ! y€uxprh–kq—’qWyui‡egtzCkyh€edy€up$ ‰r de b q Y a x Y q q Y W ir Yp q W q xh ) ) XÊ Ê Ï Ê V X Ê8 Ê V X Ê8 0 % p y8 Ï ) 0 ' y Ï Ï8 p–Ï ) 0 ' A8 ) 0 ' 0 Ï yÊ ÏÊ V 8 y dr Y ˜ d Y q W W€ p‡™Ayrzs€uedo€eacb TR US  D & §  § 26 D  D 8 D §D § ¨ § D 6  D S §D § ÒAÒ ©YS`iWfg#iY@ArSvhsSiAC4edW¢edu¿Gfv pyxA`rGArh–¢uiYSiCwY‡’—s•xtHr`i’W•v`—AS’iyWAHrYyxS’iW™p©h#aCp a d ww € q g ehp W€ Y •x q g Y g Y • d a Y ˜ Y d ˜ a x Y qr • Y À p x Y q h i Wp r q r WÅhh Si™`e•A©h§¼Æ fqsd@Ðky€qgtSy€q´#Y’·`WtrHyYxk€b§n#˜’Y•pyhSwky€u–tPtreAhaAu Ys€gqTf`iup%¢» ² ¯­ ³±°4®¬ %  ' 6 & I  2n 1 1 0 0  xn w1 w2    l  xn lw1 x1l −1   1 1 0º wn » »  » l»  lw1 xn−1 » ¼   ˜`YAsh€`dA#rp‰sY¥terp©hAu ‡y€qfpyhu#bv™ •p xw • Y x i Wp ar gh Y d i q Y #RC2I%" 0™#RQ(@•¦0#&P§Q(#)§0EC7 ¤ G G & ) †( )7 3 " I 3D B ¡ ¡ « ª §¦ !G©¨uC¥ « ª §¦ !G©¨uC¥ 7 l x2 x2 1 Computing the points and weights General Quadrature Scheme § Columns become linearly dependent for high order ª1 « § w · « x1 JF ¨ ¸ = © x¹ «  «l « x1 ¬ ' 5 Normalized 1-D Problem 5 4 ¡ #v#dsYtfYdp©ha qhi xe 5 1 v 2 1 v 5 % 5 p½ v 5 sXyt4‡s©™Gd ¢pyu‘™ ‡p@A™sA™Gp$ Htep’y‡’stHGeW 2Ò½ Y € q Y Y d À g h € q Y ˜ ¼ i h q b Y q e h p r Y x h ˜ x h v x r• e  iW d —• x r Y i d # ! cg`—p’Å`W™sA™FedTSszAS™A#•TAAsp%A©§‡'€`btXtS§H%‡µ‡s©™˜ iY € b •Y • Yxh ˜ W qr €q ih iY ˜h iY € w rÀYxY‰ iha d —r •r qh i dY h g g h €qY d GAyhu#bvuAr—’–yW#xA…sh‘’q©e…£‡¹Cft#r•TSs€A`hyW#xyY’qyWx…S¢p’hsW‡’epfh±qAÑzy€…yh‡yYp`h©•a ¼i q Y ™ v Y‰ q • a Wg Yp i a qr q i a rÀi q • d ar Y Y q q d `—y•ui`Y©W… yek`d! #gyYtRfy€v’•sWpSe H%’isWpp©hked!p’hsWTrxsYÎqfu†©vF  A`hyW©qSre– q a d W qrgwe Y q qi gYe qi a W i q W Y€ ¿ # ! i ai x Y •ihi Y q –h Y qr Yg •ri d Y W •e e i Wp ar Y q W t#rS’i`WAS6y€gA´f p’‰sWTT’‰sWxASA`hfW@i™’˜`gsWÎqµr˜terAhAu s€kedk# ! yY#xcb Y€ 2 ¥ ¦ 2 3 © 1 ¥ 0 ¥  E W% $ ¥  c ) ‘  ¥ ' ( " §# ¥ $ ! ¥ ¦ ©  #D ¤ ) © &% # ! ¥ ¦ ¤ di q e Y –h d d x • Y q  W • — ˜ – W Y Sp’hsWTrSu#tFt„˜yYfq‡f—vA#rY`i`Wg€yµ`iu`‰Ayhd¥%p Ashrx…HgY`i´s˜xsY©q‡g‘`dW ) sYTrxyY’qW q Yy€§#ÀA`rS©eeayWqytSEiifvn˜`Ypyhw}HY`·yWxA#r`i`W…y€qnS’i%`‰•Ayh®%} t…`etT‰oy€fS’isWTRe q x • x r w h • p x g Y • Y  W d — p – h Y • r Y q‘  q r g w #i€y´@t%C qTXsd#dpgr”‡p@Ar #˜’YpyhxSvy€‡S’iW#’·sWtH’iW`kA´HgfYdpSXedkg‡y€qY™vuFE Y q girÀ r Ye qehp •p w Y q  xrY • ih rp W h ˜ Y€ ¿ v %# &F ! T@yza@fg ’eAT‰vg`iW‡SS§Arh™Yu`tS#YsXA`‰ytsx#y’qztFedWz‡s©Y™0d ¼¢pyu#v™ qr €q € ed Y • r r  g i Ð x – e W i €a q Y Wqr Y W ir g h €q ˜ i hq bY ² ¯­ ³±°4®¬ " ¥ ! ¥ ¥ ¥ £ ¥  ¤  ¥  … SMA-HPC ©2000 MIT ª1º 1 º ª w1 º « » »  xn » « w2 » « »« »−«  » =0   » «  » «1 » »«l» l »«  xn ¼ ¬ wn ¼ ³ x dx »« » « » ¬0 ¼  l x2 x2   ¥ " $ % ¤c ¨ ) ¥ c " 1 Computing the points and weights $ ª1 «x § w· F¨ ¸=« 1 © x¹ «  «l « ¬ x1 ‘ Normalized 1-D Problem is€€`dnAshub‘ArhA`hyW©qSreutHr`i’Wp‡u¿ Y q W i q Y ™ x – i ai – x Y •ihi Y€  ¥ "$  ¥ E %#! W&F  © ¨ § ¥ ¦ ¤ †7 " ( ¢0#&u ’ @H7™£ "  yYSo`gSAh ˆ de •e #RC2I%" 0™#RQ(@•¦0#&P§Q(#)§0EC7 ¤ G G & ) †( )7 3 " I 3D B ¡ ¡ General Quadrature Scheme ¢ £ 5 ¡ 2 Ò i =1 Low order terms n −1 i n −1 0 i High Order Terms 2 n −1 n 1 0 i =1 i 2 n −1 i n 1 0 i0 n i i =1 in i n 1 n eryW#xyY’qyWx€sd#dSyi±qp©Ñgs€µS’isWÎqsWcbsYux‡egfW@Ah ˆ a Y ar Y Y q  x x Y di «!ªG©¨§uC¥ ¦ Q2) @0¢rC‘6R%'u”§o §H&³02(H)rC27F4 2£G G U B7 ( U7 3 ("9" ) " I ) U D ¡ ¡ v •w ˜ d —p d Yp i @’i`W™tsr¥%iYyYzCoAra sd©qaS‡‡yhSx#`i•phy©q#aA³Sit‡d±qS‡‡sh‡xS’i0ter™A%`—AC³i#Y™qbq@PA`hyW©q#ap©h™u¿ eg xw Yi W x Y‰ g r aeg xw Yi W• W ˜hi •hw Y Yp i Yii a Y€ Ï ! # )0' % p98# rÏe¢  $2τ!Ê  (&# ¢  S3Ê  ! i q qh p’hsWTsrAi Yy€’g#C’psW©xsaAr’@tk±qaS‡‡shw™xiS’i!tr™A}yh®HyYsxxrY–yYx isY©qt´`d™trA#yYu’ioATCtnu¿ q Y d gir eg x Y W i dr q g –h W • x q W Y‰hpr Y€   0 Ë †% © E Ï 8# Ï 2% A„$2! ¢ # rÏ! Ê ) ' ’–§US§C@C70#&'SzCp™yh¥`gWAyrgsYxtuedterW´˜pi%T@™€qbuQÀAsrpyYu`i…HgY`·’WA™s˜ASgy€up$ W  (7I 97 Y q d r • h —hw h q • x q W •r xhi Y q xh £¡¡ 2S2’C§¢2@•S§C@C0#&'9 © C2£G G U ) B7( U7 3U(7I7 1 General Quadrature Scheme ¤ 0 i =1 ³ c ( x ) dx = ¦ w c ( x ) ³ c ( x ) dx = ¦ w c ( x ) 0 ³ c ( x ) dx = ¦ w c ( x ) ³ c ( x ) dx = ¦ w c ( x ) Exploiting the different polynomials Normalized 1-D Problem TR US '  !G©¨u§C¥ «ª ¦ vpyYprh–©CYypS€a‡€u#dYSAgYs€¥iAS€y€HdY’sWfq#xY@pyhw‘xyYyqq©C•p'S€ x Yp dh Y b ih q r q w x Yp Y‰r `d•A`rW™˜ASiu`—•ACw™YfHYy€q‡trSy€q•ed¥erWy#xyY’qyWx€AhqfY‡dSp’h%‰Ys™Ys€q Ahsrx¹YSisYxR’¶‡g6’—•SA™u¿ hhdq W a – d e W xw ˜ – a Y W ih Y€ ² ¯­ ³±°4®¬  0 i =1 ³ c ( x ) dx = ¦ w c ( x ) l il i n 1 Be linearly independent 0 i =1 Contain an x i term ³ c ( x ) dx = ¦ w c ( x ) 1 i1 i Each ci polynomial must n 0 1 i =1 ³ c ( x ) dx = ¦ w c ( x ) 0 i0 BUT i n « ª §¦ !G©¨uC¥ ’–W…’—•SAhgSt’–¥HSqÐedyWtyrzC¢`•’•cµeryW#xyY’qyWxcydHdyi©qA“ i ir W gY q d Yp W b a Y arÑ ¤¡¡ 2Q’C§¢2@•3 R%Tu”¨§CH" ¦ ¥2£G G U ) B7 ( U7 ("9" )  General Quadrature Scheme  1 Use different polynomials Normalized 1-D Problem v gYi Wgi a Yw xw W W qr q d d d Y ˜ d W •r xh " 4HA`hyW’q‡Sp©h…`—y•#x@pyhSc`dv’qvT@y€µfhzyYSiyYg™Ayh`iv" #gY`·`Wtky˜p±i0edW qr q drp• W ˜hi •hw r W qph q W i q • d Y€ ¿  d a W W b Yi Wgi TSs€edfWpS@A`r™p%’—p@¹…’iAsr‡p¹yh”ed…A`hyW‡`eAyh•Su0fv `ifWp#rsY©xS`i…y€’qckHgA`hyW’q‡Sp©ha ’•`•vH™A©hH@u‡'u’bA‘edfWpSveducp£¢y fÏ3 H”ˆ”H”‰fÏ3 G ¡u YASyw¹s€kyhu`ik`iup%’hgsYt¹€b W dY ˜ aYp d—r •r drp W€ ¿ v 3 Æ ar d Y q q W  WÅh • xr Y  % V § Ò½ ) ) #iy€q§ÈÀter´pi%’—p@¹`dWSy€utgdqt%shkyYAr¹d ʼ …y€q®@tzfS$2Ï!  V  3 ˆ”H”ˆ”3 # rÏ!  ¦3$2Ï! 0  AY ÈvrWe Y • W ˜h •hw q –h h x x Ï Y gir # # v Y#yYp‡g…p’hsWxrY’iWgt™edterW™A%`—AhCv`•tr sh™•ASApA€sfqpkedW¥ÎqoTS€s!€@fedk# rÏ! V  ter´pi x Y x – –h • ˜hi • w • q rihh xh W qr q a • W ˜h à hw Y dhh Y b p`—•ACfs€q¥yYp%€af€•8–h vpq—Îq`Wt@AhApy€qyAuAr`™ASu—`ACw Ah™#d`Y`W•sYx!#wyYfq™qutruY`•™edu¿ •ri h xh • W ˜hi •h i d i Y W€ ² ¯­ ³±°4®¬  " #gY@qÐdeyWqtyrd¥`—•y©qaAr©Ñ‡sYV tx usdu`itryxfqSAhcx‡gAc#xA`‡u¿ Y r q W di a Y xh Y€ W€ Y€ # rÏ! t€sA%y®g@™sud y@gSykeaXV –h dqh hx Æ Y p hq ¼Ê Ï Y €q Å W n i =1 n −1 in i n −1 i n n `dA`r™ASu`—ACRtSpApy€yqA‡yY‘¤ • W ˜hi •hw•rihh xh d 0 i =1 ( ( ³ cn0x ) c0 ( x ) dx = ¦ wi cn 0xi ) c0 ( xi ) « ª §¦ 4G©¨uC¥ SMA-HPC ©2000 MIT 0 c ) ³ c ( x0c ( x ) dx = ¦ w c ( x )0 ( x ) 1 1 #@¢o § r2)S§C@C0#&'9 © 02(H) 77 †( U(7I7 I Using orthogonality and roots ¤ ¥ Normalized 1-D Problem General Quadrature Scheme ¢ ¦ 2£G G ¡¡ vd qYx t" ks˜xsY™#x‡gAh #x™T`h±•" S’i%`‰AA%`i‘@A`hyW’q‡SAh‘pz" ™s˜xsYv#x‡sgpvxp’S±€" y€vydHdyY‡Ñz™´‡eg’Wy•`ds€fiih Y b  W •h‰ W di Wgi a dr d q Y xh Y W Y q xw Y Y b Y d W q 0 ³±°4®¬ ² ¯­ uqi‡Y@iCw‡YS`i…`—y•Ax#r`i`WsYtxu# rÏe!¢  $2Ï! V  d±qS‡‡sh‡u¿ gY gW Y • r # aeg xw Y€ ) ) ) ) #Ê Ï R@2!  V # C) rÏ! V Ê  ) X WT Ê  VÊ %†A8„$rÏ!  V # rÏ! V ) 0 ¥ # CrÏ!  V Ê X VÊ T % p9$2!  V ) 0 ¥ Ï#   Ê Ï 8# Ï  vvv vvv )  )  # Ï #Ê wÊ@2! 0 w@Ïr! V Ê X VÊ T % p9# rÏ! 0 $rÏ! V ) 0 ¥ R@2! V Ê X VÊ T †A„$2! V ) 0 ¥ Ï 8 #  #Ê Ï  % Ï 8# Ï ’—s•uisYxR懑™s˜xsY€x‡gAc#xA`gs€•sYÎqsWctr ˆ q Y Wg d q Y xh Y€ W€ Y q x b i Q2) @0¢rC‘6R%'u”§o §H&³02(H)rC27F4  2£G G U B7 ( U7 3 ("9" ) " I ) U D ¡ ¡ § B B § 1   ª §¦ G©¨uC¥ vt" ks˜xsYuxY‡gAuxS€A`±€" s€q€AcdSp’hsWÎqWSAhPÆ d q xh Y W Y dr i gi a qqdp`r•Yy€q¥gSAr‘" dky˜x#yYq™x‡gAh™xY€T’h©•" YSy€q™AudSp’hsWÎqW‡@A©ha™# G rÆ!¢fqsdxS€sC’•Arvqq@qe®º i Yx b dr i gi ÐY q• a d Y ² ¯­ ³±°4®¬ c  # V Æ ËË `   iiR ` R ¢ £ !G c e!Æ # ¤# % V R ¥ ¦ ¡ qr q • Y tSy€R`•t#r#ac ! „   ¡0 ¦ t Ò v x q xgret der ¤ x pyYseTr‡A@pcsdSA‡Arh– sdu`iAC³A`hyWtsrpsYp`iPSt‡sdup’#W™f`i‡S@”tkSpet€HyY´s€#dY`·sWtk™fek‡eg`Wf•gedu¿ q Whw i q x q W gir q€ Y b  WgiÐ –h Ye Wi a q Y q xr ˜ ˜ d Y d W€ )  ±4®¬  °¯ ­ # Ê rÏ! Ê X VÊ U% p9$2! ) 0 ¥ Sf™cA4Aspyu`igysTk’˜‡'sw m v T Ï 8# Ï ˜ ed r d r • rx Yq W Y €q Yqr W Ñhx w Ê #RC2I%¹£ s€uArh– #˜yYfq‡f—vA#r`i’W´gA`‰Ah s S½ & )" Y q x d d x Y • r Y • v  E % #Ê W&wCrÏ! V T@ySyd qr €q €a e •terW´ph%—’ph@%•ASiApA€sfqphx‡gAh‡Í24Ƒy€CtfƦ3H”ˆ”H”˜G % XS3@tp#@¢To‘yY‡™ph ˆ v 2 ˜ i • w r h h x Y x Ì Y q–h 3 Ë ÊÏÀ 779 Æ qew ˜  0 Ë †% © E Ï 8 ! # Ï! 2% p9# rÏe¢ $2„Ê ) ' ed•terW™Ah%`—AC@ASApAs€fqp‡G c±qSseyxfqSph ˆ v Ò ˜ i •hw•rihh xh Æ a di  FªG©¨§uC¥ ¦  D B )97 U ¡ ¡ FpH"4 §H&rTSC2Io¥ 2£§£¢G “ ž “—ž  “ ” ŽŽ “ Š c v0v ‘’ fu%•F„ u   t ' TR US D ¢ ž F ¡ ¡  ¥ ©  § v ddx Y •r q n˜sYqq‡y—vtHr’i`W´´S’isWTrxY Çip@#ÀsdquA`W€¹Sy€g`i™y˜#xyYq‡oshn#gyYS¹Cft#r‘dp`itryxfqSp©hE‘#x™T`hzu‘GcT%e Y € Y b Y q Y W Yg q de Yp i a q W di a Æ Y b • Y€ ¿ vi bhiÅi ’•`•yWqfkyYAr‡sduA`W€ky€6’—p4¢cT%Å TSµsYtzd y¼@ks€µYS`iV W  s 4#gSfÐedyWTsr6`—’•AreayWTkpyhq d x q€ Y b Y q •ihÀi bhi bhi xr ÊÏ Y q a v Y q d • qr ˜ Ãer xr W di a Y xh Y€ W€ Y q ÊÏ Y q x – # –h qh x Y q  de — ‡t¥Ystµsdqu`itryxfqSAh¹x‡gA´xSA`¥s€FÀd y¼@¥Sy€µArhz$rÏ! tzdt%sh¥s€6`ifW@P€© ² ¯­ ³±°4®¬ vedter™A%`—AC@ASApAs€fqp¥Hg’Y™AyrgsYtuvTc • W ˜hi •hw•rihh xh •w ˜ d xr d bh i bhiÅ xr ÊÏÀ9" bh cTSusYtvd f¼C‡%‡%§Q(2)U Tv™ ) )) V V ) V Ï# p98$Ï2!  V  ) 0 ¥ # 2!  V ¤ Ï ¨ ¤ % ¨ ¤ ¨ # rÏ!  vvv  v) vv vvv  vvv # ) rÏvvv ! 0  8 §¦¥ Ïp9# rÏ! 0 ) 0 ¥ ¢£ §¦¥ ¢£ §¦¥ # V rÏ! 0 ¢£ i qre Y Y ar Y q W xppr p’hsWTSu#tusd#dSyi±qp©Ñ…HgyYt`r%‰sYAfi m #RQ(T9#§C7 ¤ u§T9 © u¹¤”¡0#&®02(¢ÈH52#)¢  2£G G ) ( 9"† 9" £7 " I )  ¡ ¡ A ¦  6  6   ª §¦ G©¨uC¥ ¢ 8v¦IHG˜¢ &tTy‰#xyYu`igy€oACtSpApy€yqA‡sYtn€`ac‡#Àedtr 3G D •r q W Y q ih•rihh xh xr W€ b • ÃΙASu`—ACyY‡@iA#¢4t#r`ayWsdAer©•gy€0t™Sp’hyWsd#xA¥#g`YtrsaStn#gsY©qΖSf€‡yYAu`dA`r™ASu`—ACw W ˜hi •hw xg Y Y»• d a Y q –h di Y‰ • d gir W d xr • W ˜hi •h €@e s 8v¦IHGFE ‰tTs‰xyYu`icy€´puASApAs€fqp™sYtk€eacFedter´p%’—p@´HgfYSc€´A`hyWTT`‰yW#x‡g a 3 D•r q W Y q ih•rihh xh xr W€ b • W ˜hi •hw de Y b i qr Y `dSy€™Arh™tSy€yYAv™ 4#gSfÐedyWTsrE’—`•t#r`asWTkAsh‡tnsYtkSp’hsWÎq‡@A©ha ´#x‡gA™xp’SSy€q W q x – qr q qh v Y q d • qr ˜ qer xr di Wgi Æ Y xh Y€ W€ Y c 9 (Ò Y™˜S€ya‡yYxSyeTrx‡pSeuz’Y•™˜`fWzYSy€q€dky˜prhf–x#@sw‡eAcA2r–‘sYseqTr‡ArSucdsStr‡¤ p’YSteryWq Y d q gr t w d x Y q h x x xg et de À •p à Y Wg W Ð W W W€ ‡isYxR¶’‡o—`•sYÎqi@‡’iuden€a`ScbS±À# Ï !  % # Ïr! iA`hyW©qSirey€qAu™˜y™vxArhcqT@€y©YeayWAv™ a – hh d r – r q qh ² ¯­ G±°4®¬ VR  § ¨# −1 ³ cos 2π xdx 1 Approximately Integrating Gauss-Quad 7 G©¨uC¥  ª §¦ SMA-HPC ©1999 MIT Number of Points E r r o r Simple Quad  ŽŽ “ ”  “ v“ %•F„  ‘’ œ u f”@0•œ P  ž ˜ ©  § ©  ©  § ¥£ ¦¤ ¡ ¢ ¥ ¦ ¡ ¤ Š u £ ¢ ' sytsASucsSt‡¡th Yx eqrx gr e t dd er ¤ – YyAr#•terW©#aYCywd¥r‡`dW´Y™˜€sadµYsseTr‡ASu™’YS´`˜yW´s€‡TSs€zYyYPpA‘P A0)('&$° "  da Y x q xgret •w d Y q qr q d eh— h ¾ ° %# !  ) V s@ASA`yyACAsceuAsy”“ h q • r i hW q x h w h x w dW x h x x gh Y ˜ Y q–h Yw xw ‡s€©q™gs€§Akq—fq#x@pyh¥—Y SMA-HPC ©1999 MIT " c 1 (x) x ) XÊ ) Æ Ë G Æ V 8 p9$2! ) ' Ï 8# Ï 0 Q§®TS§6RQ(@•432uTp0%" ‚ 2£¢G U 3 B97 5 )7 1( U9" ( ¡ ¡ Ÿ– ž “—ž  Sc v0v v“ œ u u Š V) hq •r i Wqxh whx w W xhxx V $ " s@ASA`hyyACAsceduAsy”“ Æ x‡gAhRth0ed•terW´˜pi%’—•p@wprh–ut¢1Xedu’YSye™yYseTr‡A@psdSA‡¤ u`iACÈwR™i m Yx– h h x $" W • x x q xgret der q Wh ÃÆ gh Y q–h d Yw xw ‡€s©qY™˜gs€) GtuH’YsWfq#x@pyhS…—Y Ê V X Ê8 A„8$2Ï! ) 0 ' Ï#  $ B  a Simple Quadrature Scheme § G©¨uC¥  ª §¦ 0 ¥ ¡¥  ¥  ©  § ¥£ ¨¦¤ ¡ ¢ ¡ ¢ c ¤# #Ê Ï R@2! ª §¦ G©¨uC¥ B B G    0 f Normalized 1-D Problem   2§H3%¥ ATGpA™¥ e §£¢G U " $93$$ ¡ ¡ ¤ 5 Ò qr q bhiÅ Y b v x •e d Y q–h W W‰ Y q W · q Yh x d Y ˜ a x TSs€„T%•€´h s pq—ÎqsWterp’iyWgy€Gtkq—’q’iW`a%‡s€¥`i¹shx#Yµsh”Hd%A‘sYyeAHr´•©Yp2rf–yed Y q W€ b qr q x Y q ir q Y b d W qr q q x r qr we d b •p gi x q W Y y€„€eacT³sYTr³s€„tSs€kx€T`hf•`dTSs€6sYTrEtÇ…vT`hSAsrpsYp`i³Sy€q yYSAr#aCg`YtrA#yYu’i‘`d³A`hyW±q@regs€%#À#xA‰€Tv 4egAshyxu#i©Y¹s€vTPcvT`hfA`hyW©qSre– de Yp •p x q W W i ai – Y q Y Y bh v W q a Y q qr we d b •p i ai d ¢#iYsY%un' HyYT#r‡’hced™RQ$o`gpyhsxp#i©Yzs€‘sYxSc¢•ASgs€¥pCAsrpsYp`igs€€ACtr ¼ x ¤ Y€ ¿ g q a • W PÏ W q a Y q Y€ b Yirw Y q ih• x q W Y q qehp qTSc‘™© ¢s€t%™f˜zsYt‡TSs€cSp’hsW±qSSre0Au`dAsrpsYp`iz`Y‡SAS…yhv`dp%yho#gsYTreayWfqedAyhd r€ b qe v qhh d xr qr q di ai – –h • x q W •gir€ q •h q q W€w Y ˜ d Y‰r bhi Y b qr q • x Y vi q q x x • x q xrgiehp  d ™pfh³p'S€ T³€otSy€6Y`·’WA#rsY³®º ¢A`hyWtsruiyY#dsYwsY‘AsrpsYp`iW y—ASp@’S’iyWSe `itkAg A´p¥A`hyWtSu#t™d p©Yp`rAG™`iu`‰Ayh§Ah ˜`YAshuA™sh™ÅASµAusdq’YÇTv™ Wr ˜h ¾½ r ih i qre Y ¼ a •wr»  W • d – •p xw xeh q arp h • bh 3 … ±4®¬  °¯ ­ ³ 1 dS ′ is an integrable singularity xci − x′ 0 Ai ,i ≈ Ai ,i = ³ Panel Area xci − xci 1 dS ′ xci − x′ point xci Collocation BTu" j 1 25%G 0#&}02(H)‡20AQ$ ¤ 2£C2£§¢ 9 U" " I U3 U ¡ ¡ œ c f Fz  f‹v  – œ Š ¦ ¡ c0u ¡ ©|o¹ ©k  ¦¢ ¨ ¢! " Basis Function Approach  G©¨u§C¥ ª ¦ panel i Ai ,i = One point quadrature Approximation Panel i panel i x y z 3-D Laplace’s Equation ¥£ ¦¤ £ ¦  ˆ ¡ ¢F© ¡  S  ¥   £ ¤ ¨ ©    ¦ § ¤ ¥£ ¢ ÈvtrTs‰xsYqp`igYs€§–AugSigYSy€uTn#gsYxsYqq@’e•sYtf€eac€`dA`r™ASu`—ACw • W q h d Y q qr d a xr W€ b • W ˜hi •h sYx‡Si#A#¢»•Ah¹dth%yhxYst™dsu`iAhCw iA`hyWtrS`etT‰nYs€…™˜€saosYseTr‡ASu™sdSt‡¤ s€’iih g Y Y – q xr q W q •r Y q Y Y d x q xgret der Y q  ±4®¬  °¯ ­ dg Y•r q W qr  d a Y q qh Si4ATy‰#xyYu’ict™S’isWxsYqqS`e©•zSy€©YeayWAv™ SMA-HPC ©1999 MIT 0 1 0 Gauss-Quadrature Points 1 Comparing Quadratures ª §¦ G©¨uC¥ R%’%p2g3 R2(C‘3 C27H)‡023" 2£¢ G (" B"$ U )7 (  U£ ¡ ¡  ŽŽ “ ”  “ v“ %•F„  ‘’ œ u f”@0•œ P u  ž ˜ Š Evaluation Point Placement   -1 -1 Simple-Quadrature Points Symmetrically Normalized 1-D Problem ©  § ©  ©  § ¥£ ¦¤ ¡ ¢ ¢ £ ¡ ¥ ¦ ¤ ¡ Large errors even with many points! ³ x 1 dx « ª §¦ 4 G©¨uC¥ 2Ò Number of Points Simple Quad Ž G Gauss-Quad −1 1 Ÿ– ž “—ž  Su ‘§Sv v“ ¥ ¡¥  ¥ Area → ∞ ª §¦ G©¨uC¥ 1 x  = 1 f(x) ª §¦ G©¨uC¥ x 1   E r r o r Approximately Integrating ©  § f”@0•œ P ž ˜ ¥ ¦ ¡ ¤ Š v Area → finite = f(x) ¡¡  ) )9 U 3I ) U 9I ) ( † ( U 9I %2"HrT‡Q§C§Q(¢ Q"0¨ T@%#"R2(0C7  §S62"F¨ TS%H"R( ¢ e ¢e §¢ ¢ dr€ qep 3G D ih •p x q W W W€ AS‘¢¦IHG˜¢ gp‡’YAsrpsYp`icedn€`acb c ` ££ ) v Ïq x e d E2% cTr…q—’qyWA`r•p`ifW™r % # rÏ! ip’hWs±aSre´r‘tropu`h•g€b Ò¾ ih q i – q Åh Y i ±°4®¬ ¯­  ¢   & x x Quad Point Note no xi = 0 £ ` Q £) ¢ ) Ê X VÊ Ï A8 T a Symmetrized 1-D Exmaple ¤ 1 ` ) ¢ ) ) ££ œu  ¥¤ £ ¥ 2"FE¢  2£¢e §¢ UD B ¡ ¡ ž— S‡4  Q‰v Š The Singular Kernel Problem ¡     ‹ ©¥¡ ¥ ¡ ¡ ' q’steA`yv`SyGA‡©S#fHygSy’4`€b — qW x r• e  iW d dW € q – h Y a i Y d Y x w Y € q iW •• Y  W x Yw–h •prw a x q xgret der W qepÀ WiÐ W gir W Y• x q W Y S’i™s˜Arhy–xC0Ag’YASt#rgyYseTr‡A@pcsdSA‡¤ ed‘‡@AyY’qSc`dSAvsdfq`d‡Ñ4trAsYu’izSy€q ª §¦ G©¨uC¥  F ydSA‡¤ yY‘¤ p8 der d Ï `) )¢ ££ ) 0 ¥ Ò ¥ ¥ # # % d8 ¥ ¥ # £   )£ ¢ )0 ¥ # `) % A8 £ £ ) ¢ )  ¥ Ï § ¥ B ¢ Ï p8 D A C 6 gi x q W g x •e dYg ih x q xgret der ¤ •w StrAsxsYp`i…HY`·yWA`rSA`ifWH‡fAyYSyeTr‡pSucydSA‡`—w m 8 # )0 ¥ ¥ & ÏpWd 8 8% ϝ% —qÎqsWterp’iyWzyYtS`i´`˜`W•™yh€Hd’YA`rsWtTAAS€ ˆ x •e d qr W Y q •p xr‰ Yir   D 9 5  U )9  57 "I ( ¡ ¡ %"H ¤ Q"U0’2B)G S87z1 %2"F¨ 2T‡gC6C§0& ¤ e §£¢§¢ v4HYfSuedtTs‰xsYp`iye¥%tkFΖ‡p’‰yWSiC‡“ g de •r q Wp d —ir ˜ W Y d YwÑ xyYyq©qCu`dATy‰#x` yYu’iSfe‘S”edter™A`˜ %`—AFV Yp •r q Wp d qÐ • W hi •h ` ` •r q Wp d a Y W x q xgre tTs‰xsYp`iyeo€p#rn`i•sYyetsrASut ) )  0  p8 ) ) )  % A8 ) ) )  ¥ ¥Ï ¥Ï •r q W i q x W Y q Y W Wgp tTy‰0 #xyYu`i¥A`hyWtsrpsY0 pq `igs€‡‡`g%`‰‡Re s ª §¦ G©¨uC¥ £ £) ¢ ) 0  0 ¥ Ï A8 x  SMA-HPC ©1999 MIT ££ D  ¢ ££ D A ££ ¢ ¢ x 1 "T§3#‡T§†03kƒk@0£#)0‡§™¥ ¨S(£uH"R2(F¨ G 2£§£¢§¢ 9 9  ©" D† § U 9 ) 3 ¡ ¡ Ž “  ” Ÿ– £ ¡˜ž § u¦¥•‘S§ ¤ ¢•œ v Š ¢ Improved Techniques ¥ ¥ © ¥ ¡ sYyetsr‡pSu€–terfWsdSA‡¤ e terp’iW s x q xgret fi der x •e i q x q W –h •p xr‰ –h Yir A`hyWtsrpsYp`i§A€#d`YteryWAT0Apt@€ ˆ yYSyeTr‡pSu™f A`‰yW‡AAQre tTs‰xsYp`ie s x q xgret Y qwrg •r q Wp f ysp©Ft%A`YtA#yu`istfeacrbe H’YsÎqsterAS’iyfy’qc™edtA#yu`i´S’i`W‡SAS0Au#ueti Ydx e ha –h À • prxYq W Yxr € W € d W Wx • e Wd € W b •rxYq W  • g ir € –h dY e W à q• Y d q d xw Y b Y – q Yi Y q h vY Y d x q xgret •w d Y ‡€#asY¢Asr#xA‰yYgp#ifYHyYSµ€#d‡eg’Wy•d #brY%ѵSy€PiigA™˜€saµsYseTr‡ASuk`Y™’˜yW´Sy€q ir q q x ar W xh ˜ W W a ––h qr ˜ r xhhw Y q x W a Y Y‰i a Y€ ASy€•yYtsrS©eapS`i‘sYA™u`dc’q‡#À©qA2rGAuxyYyqTk´ud #vp%@uxSy€Tr€`d•©YS#iAsxpup©h‡u¿ vi ai – W q  q x q W qr ghh Y‰ qhi W x q xgret der qr q Y dph Y ¢A`hyW©qSreX`ds€’isWTrA#yYu’i™T´%%p¹s—xAXAeduyYm seTr‡A@p”ydSA‡¤ T@y€™As‰xyYSAv®º  ±4®¬  °¯ ­   0 Ò v Æ % xr q€ Y b  Wgih x a Y q gi ' 4Ê sYtudpp’#W™´`i‡Sp@yw#dsYyxp©hzs€…@tr  Æ %ÊÏ GË —p x %£s˜Arh– g d a W Y W xr q Whw i q x q W Y€ ¿ • W ˜hi •hw Y d— Y Y q Yp HfYu’h•’’iE#iA`‰AksYt¹du’ip@®A`hyWtsrpyYu`i…unvedter™A%`—ACP#‰y€%%#p€ ˆ s€k@³yhq qeh x q • W ˜hi •hw•rihh xh  Wgih d x a Y q gi Ï ‡A´yiSye§`dA`r™p%’—p@utSpApy€yqA•’iSACfw) HyYsxAhcy€¹Str P G  % # rÏ! pyY#xv G ÀY Ï P G ) ' " % #Ï $2! • x q W x •e d q W •w ˜r Y trA#yYu`icterA`iyWz)Y y€€ed‡’Y™A‡Ñni m   ” # G rÆHˆˆˆH˜HF2% © !3”””3G3E E #! W% $rÏA¢ # rÏ! V $2! )  ' #Ï v W d q€ Y b Y q q Y x W AY Èv±Àf ©Ñe vuA`W€Sy€™yhu©q#aCfwHdyYny€’qcb #YyYp‡Y gvp’hsWxrY’iWFtuedter™A%—`AC¢’•A…yhCt@AAps€Acedn€`acb  # 2 ! % ‡YsYA#‡th x x – –h • W ˜hi •hw •r q•rihh xh W W€ G Æ x Yg– # Ï • W ˜h i •h r –h dqh Æ Y €q xr dq Wh w i WqrxYq W Y €q qr €q q eh d ix eq $2! V tre™A%`—ACwoutsA%hyx ™y´Ystgsu`iAC³p’hsTAsu’i™sgTSsA‡Ssyg’qW Àp’—p`ifWedyWSywxyeXt‘4qv qfp™TFxY‡gASAp•A`r™p%’—p@•`d§# rÏ! #icF±qAрCyh§#g#Yi • x d qh ™ dh ˜ qr xh–h W ˜hi •hw r W Y€ b ar Y Yp qd Y ereµy˜prhuyYseqTr‡A@p•€u0SE $Ïr! ) q—qÎsWterASi’yW…cT%)i ‡‘A'@k€`ac”SA`hyW©qSre”Arh– • x – x xgret Y ¿ v # x •e d i bh Å r Y‰r€ W€ b di ai – x XÊ # CrÏ! Ê V 8 %†p9# Ïr! $rÏ! )  ' Ê Ï8 # y˜prh‡y€0t€d#™#˜€ya‡Ysyetsr‡prSuiA`rWysddSt‡¤ yYqtsrx#ipkyhq’Y•`pyyddu@wuedv8qh x – Y q –h Y Y d x q xg et er Y Y W h W ±°4®¬ ¯­ $ § VR # § c %# " V B c B c 3 ¢ ¢  TR US V c ! # D 8   V 8   V   sYyetsrASucsdSt‡0A`Å`WusduA`W€‡yYTere`aAr ˆ x q xgret der ¤ Y • q€ Y b q • • ) xsgp0ACA`r™p%’—p@cedu# rÏ! y sYxcPA„8$2Ï! $2Ï! y )  ¥ Y xh–h• W ˜hi •hw W Y€ b Ï # # x – ar Y •e x – x q xgret r w Y Y Arhu±qp©Ñµ`rµy˜prh‘sYseTr‡ASuknA`h•p‰•¾ ) s€t%™f˜cedu# rÏ! yY#xc6p98# rÏ! # rÏ! )  ¥ qhh d W Y€ b Ï # rÏ! q’qyA`rp`ifncT%uAS¥StA#yuih — Wx • e  Wd i bh i Å d r € g irx Yq i   †  U D B 1 "9 9 † 3 ƒ 9 U 3I ) ¡ ¡ %#" ¤ %#"‡p2$2)FEC7 ¤ ¦'§#3'§0k20@02(G  §£¢§¢ V       V  V ª §¦ G©¨uC¥ ¢ ª §¦ G©¨uC¥ § ª §¦ G©¨uC¥ 7 2 ¦  ( 7)  ( " B) † "U D) U 3 B 3 7   0C2H§%’Q§2FrHr0Pk§@”¨ ™§& yYSyeTr‡pSuuA`rSA`iW s x q xgret x •e i r x i x¿ •p xr £ gir i q qwr A`hyWqTks˜Arhf–dStr§‡`YA`rsWtX¤StnA`hyWtsr‡Ag m 2%0T9pCk‡20C§Q(#)¥@87k%§4Q1§0'u" j U" ( "l 9 U 3I  9 5 " 3 ) ( &$ edter™˜A%`—AhC@ASiApAs€fqp§A`YAsh4teryWu#ifYsd”“ • W hi • w•r hh xh–h • x• q d YSipy#xA%AhvxsYqqpS$ a Y Y‰i a dr "9 9 † 3 (  3 T§3#‡T§041 S2H)#0 ‚ €AuhyStrvutsrfqSp©hgfYedcbY#a`YXV ar xww qi di a W W QTSI%#"RQ(PI02(HGFDECA$¥@8756"04210§'%#"! U 9 ) ) 3 B7 9 3 ) ( &$  | £ £k ¢ ¡¡ ¢¢ ¦ #vsduA`W€¥p’hsWTrA#yYu’igSy€yY‡S´p©h…sh…HfYSgCtrcdu’iAsrsxqqSp©hcxsgp€#x™T`h‡u¿ q€ Y b i q x q W Y q qew ˜ a q g de Yp i a q W di a Y xh Y b • Y€ 0 ±°4®¬  ¯­  # rÏ! # rÏ! t€edtAsu’icprh”p`rS¹sArh”Heayu`—tSA2g##v™ – h • rx  Yq W x – d • e ˜x – f Wq • r i re g Y Y )) V V ) V # rÏ!  V ¤ ¤ ¨ # 2Ï!  %¨ vv vv ) vvv  ¤ v) vv Vv  v # 2! 0 ¢£ Ï ¢££ ¦§¥ ¢£ §¦¥ # rÏ! 0 #RC2I%" '§#3'§0k20@02(G £¢§£¢§¢ & ) "9 9 † 3 ƒ 9 U 3I ) ¡ ¡ V ¨ ¦ ¥¦ ) ) Ï# # A„8$rÏ! $rÏ!  V  )  ¥ vv ) Ï# A„8$2Ï! $v# 2Ï! 0  )  ¥ V V  6  6 B ¢ ...
View Full Document

This note was uploaded on 11/08/2011 for the course AERO 16.872 taught by Professor Danielhastings during the Fall '03 term at MIT.

Ask a homework question - tutors are online