Topics in Adaptive Computational Methods for Differential ... - logg.org

0 downloads 0 Views 1MB Size Report
U sing the G alerkin orthogonality, we can write the error representation in the .... U sing standard techniques [11] we can prove that the mc G ( q ) method is of ..... is (appro x imately) the mean distance between the Sun and Е arth, the time coordinates per ..... in the horizontal, vertical and diagonal directions respectively.
            

PREPRINT 2001–20

Topics in adaptive computational methods for differential equations Claes Johnson, Johan Hoffman and Anders Logg Chalmers Finite Element Center CHALMERS UNIVERSITY OF TECHNOLOGY Goteborg ¨ Sweden 2001

         

  

Preprint 2001–20

Topics in adaptive computational methods for differential equations Claes Johnson, Johan Hoffman and Anders Logg

Chalmers Finite Element Center Chalmers University of Technology SE–412 96 Goteborg ¨ Sweden Goteborg, ¨ October 2001

Topics in adaptive computational methods for differential equations Claes Johnson, Johan Hoffman and Anders Logg NO 2001–20 ISSN 1404–4382

Chalmers Finite Element Center Chalmers University of Technology SE–412 96 Goteborg ¨ Sweden Telephone: +46 (0)31 772 1000 Fax: +46 (0)31 772 3595     

Printed in Sweden Chalmers University of Technology Goteborg, ¨ Sweden 2001

     !"#$% '&() *&(&+,), * !"- . / 02143652798:^†u‡.ou>]€x^`^k‚@_ˆc‰>‡yb^`]r+hŠ]\Sv`^†u‡ ? ou>]fi\+jyb\~_p\Su eE_phr+\ onjj>v†^kaSonep^`hu>_K^`u>adv`c>]^`u‡,ebw\adh*rj>cebo‰^†v`^†eg‹onu>]jyb\~]^`adebon‰>^†v`^|eg‹hnm2epw\ 7 h*v`oy 7 ‹_geb\SrŒonu>]o_ˆ„ jF\aŽes_hnmC_ˆc‰>‡yb^`]lr+hŠ]\Sv`^†u‡(^`u#adh*u q*\~adep^`hu„…]^†{c>_p^`hu„‘yp\oadep^`hu_p‹Š_ˆep\~r(_S’ “Q”9•t–—˜2™AšA›[œ2—…™A– ž(ŸC .¡Ž¢£¤g¥C¦C§C¨© *¦ª¡t¨«[¨¬~­i *®*¡¬l£¤´nµ¶4·¹¸ºp»…¼4½¾¿ ¶±¸ºˆ»b³>¼/Àˆ©£4§C  «Á¯°¦CÂWÃE¨¦C§/ij>·ƒÅ»…¼4½ ¾¿ ¶±¸ºp»b³>¼Àg® £Q©­C¥±¡Æ¨F¡Æ¯Á£Q¦ÇÃÈ4¨É~ l¡£Ê§Y¨ŠË¯°¦ª¡Æ *ÉÆ¨®¡Æ¯Á¦C£Q¦¨¦C *¢Ì«° ÍQ *«2¯°¦¡ÆŸC lÎY  «Á§£¤9ÏO³>µzÐI¶4ºs¸ºˆ»b³>¼I¸· Ñ/¸ºgÒY¾µ¸ºˆ»bӍ¸·2ÑÔ³ Õ4¾·ƒ»…¼4½ ”EÖ ¨É~¡¯°® ¥±«‘¨É+¤g£4® ¥C¬¯Á¬9­±¥±¡9£Q¦ ¨F§Y¨­±¡¯ÁÍQ .© *¡ŸC£Ê§±¬¢¯¹¡ÆŸ¤g   *§Ê׎ØY¨F®tÙ¤gɍ£Q© ¡ÆŸ± z® £©­C¥Ê¡t¨F¡¯°£Q¦Y¨F«I­±É£4®   ¬~¬¡Æ£ɍ Š¨F®tŸÂQ£ª¨«Á¬O£¤Ú *Û® ¯Á  ¦C®Ë¨¦C§Ü® £Q¦ª¡ÆÉ~£Q«I£¤[©£4§C *«°¯°¦±Â¨F¦C§® £Q©­C¥Ê× ¡t¨>¡Æ¯°£¦Y¨«¡Æ #«° ÍQ *«bȱ¢ŸC¯Á®tŸë¯Á¬+­CÉ~  ¬ *¦ª¡Æ  §ë¯°¦¡ŸC .¬  É~¯° *¬(£¤¡Æ¯°® ¨«Ú©£Ê§±  «£¤K¡ŸC .¤g£Qɍ©ùøÞÀˆúiÃûýüAÈY¢ŸC *ɍ øþ¯°¬9¨ßՐ» ÿ,¾´Æ¾¼Yºˆ»ˆ¸· ³tÐY¾´~¸ºd³>´nÈOü ¯°¬ÂQ¯Á͐  ¦ÌÕª¸ºs¸ÈE¨¦±§ ú ¯Á¬Þ¡ÆŸC  ij>·ƒ¶4ºp»b³>¼ ” ž(ŸC ©£Ê§C *«¯Á¬Þ¬¥CØ    ®*¡¡£/­i *É~¡Æ¥±ÉØY¨F¡¯°£Q¦±¬ ¤gɍ£© Õª¸ºs¸É *­Cɍ *¬  ¦ª¡  §ت Ë üA Èǵܳ Õ4¾·ƒ»…¼4½+ɍ *­Cɍ *¬  ¦ª¡  §ØªË ø È4¨¦C§/Ót³>µzÐI¶4ºs¸ºˆ»b³>¼É~  ­CÉ~  ¬ *¦ª¡Æ  §ÜØW Ë  Í4¯° ¢ *§%¨¬#¨¦/¨F­C­Cɍ£ ʯ°©¨F¡ ¬~£Q«°¥±¡¯°£Q¦¡Æ£¨Ü­i *É~¡¥CɍØI  §/­CÉ~£QØC«Á  © ø À úI Ã(û ü  ¢¯¹¡ÆŸ   ±¨®*¡#¬~£Q«°¥Ê¡Æ¯°£¦ ú  ” ÷䠬ƨŠË ¡ŸY¨F¡z¡ŸC Õª¸ºb ¸ IµÜ³ Õ4¾·ƒ»…¼4½ë¾´´Æ³>´9¯°¬l  ãW¥Y¨«µzÐI¶4ºb¸ºp»b³>¼I¸·L¾´n´t³>´¯Á¬   ! :"$ad#&ep%'h*‰F \~()y #+ *-? , "$

%./ %.’ 0 D cv†ep^†„‘o*]onjeb^†qŠ^†eg‹ ? ^`u>]^`qŠ^`]c>onvªeb^†r+\S„‘_ˆep\~jj^`u‡ ? _pc‰‡*yp^k]9r+h]\~v†^`u‡>’ G \Sj>oyˆebr\~u e2hm D otebw\Sr(onep^kaS_ ?F0 w>onv`r+\Sys2_ 16u^†q*\Sys_ˆ^†eg‹(hn4m 34\aswuhv`h‡*‹ ?765 5  &87 9+J hn: eb\S‰hyb‡ ?7 f2\~]\~u ? ; =  advkon\@_ ?Ar(onepwW’ asw>onv`r+\Sys_S’ _p\ MCc-dA/MR(eBC^2MP@OR#DO =EQ @’FHGI"$J 'KL@#

NMPO ORQ STKL(PU "$%.T# VF.G- F.£QÉ6¡ŸC ®*£Q©­C¥±¡t¨>¡Æ¯°£¦Y¨«   É~ɍ£QÉ Xú  Ȑ¡ŽË4­C¯Á®Š¨«Á«ÁË.ت˨­C­CÉ~£Q­Cɍ¯°¨F¡Æ *«ÁË#©£4§C¯Á¤…Ë4¯°¦±Â¡ŸC «Á£Ê® ¨«ª©  ¬Ÿ¬¯°å*  ” žEË4­C¯Á®Š¨«Á«Á˨¦¨§Y¨­±¡¯Á͐  © *¡ÆŸ±£Ê§%¯°¬#ØC¨¬ *§ä£Q¦ßÐY³ÄºŽ¾´»b³>´n»l¾´n´Æ³>´, *¬~¡Æ¯Á©Ü¨>¡Æ  ¬# *¬~¡¯°©¨F¡Æ¯Á¦CÂܡƟ± ©£4§C  «Á¯°¦C£É®*£Q©­C¥±¡t¨>¡Æ¯°£¦Y¨«   É~ɍ£QÉ~¬¡Æ  §ɍ *¬¯Á§C¥Y¨«Á¬ ü ìøÀ úi à ” ' §C¨­±¡Æ¯¹ÍQ #© *¡ŸC£Ê§±¬+¤g£QÉ(¬£Q«¹Í4¯°¦C§C¯¹â2 *ɍ *¦W¡¯‘¨«i  ãW¥Y¨F¡¯°£Q¦±¬+¤g£Qɍ©õ¨ɍ¨­C¯Á§C«ÁË   Ê­Y¨¦±§C¯°¦C¨É~ Š¨Þ¢¯Á¡Ÿ ¯°©­I£QÉ~¡Æ¨¦ª¡+¨­±­C«°¯Á®Š¨F¡¯°£Q¦C¬ ”LÝ ¦C ,©¨ŠË¨FɍÂQ¥C l¡ÆŸY¨F¡¡ÆŸC #¨F®*¡Æ¯¹Í4¯Á¡ŽË£¤[¤g£Qɍ©¥C«‘¨F¡¯°¦C¨¦C§¬~£Q«ÁÍ4¯Á¦CÂÞ  ãW¥Y¨× ¡Æ¯Á£Q¦C¬6¯°¬6¯°¦CŸ±  ɍ *¦ª¡Æ«Á˨§C¨­±¡Æ¯¹ÍQ @¢¯Á¡Ÿ¤g   *§Ê׎ØC¨®tÙØI *¡Ž¢O   *¦® £Q©­C¥±¡Æ¨F¡Æ¯Á£Q¦¨F¦C§©£Ê§±  «°¯Á¦C ” ' §Y¨­±¡¯ÁÍ4¯Á¡ŽË £¤©£Ê§±  «°¯Á¦C¨¦C§ ®*£Q©­±¥±¡t¨F¡¯°£Q¦ ®Š¨F¦ Øi ÍʯÁ *¢O  § ¨¬£Q¦C ¨F¬­I  ®*¡Þ£¤Þ³tÐIºp»…µ*» )¸ºˆ»b³>¼ ¢¯¹¡ÆŸä¡ŸC £QØÊ×    ®¡Æ¯¹ÍQ £¤O¯°©­Cɍ£Í4¯°¦C¡ƟC ©£Ê§C *«K£QÉ,¡ŸC ® £Q©­C¥±¡Æ¨F¡Æ¯Á£Q¦ ” ' §C§C¯Á¦C¨«°¬~£ë¨F¬­I  ®*¡¬£F¤O£­±¡Æ¯Á©¯Á动F¡¯°£Q¦ £¤O¬£Q«Á¥±¡Æ¯Á£Q¦C¬*ÈÚ¢ŸC¯°®tŸÔ¯°¬#£¤…¡  ¦%¡ÆŸC ©¨¯°¦/£QØ    ®*¡¯Á͐ È6£Q¦C ÂQ *¡¬¨¤g¥C«Á«K­C¯Á®*¡Æ¥±É £¤O£Q­±¡Æ¯Á©¯Á动F¡¯°£Q¦ë£¤ ©£Ê§C *«bÈÊ® £Q©­C¥±¡Æ¨F¡Æ¯Á£Q¦¨¦C§¬£Q«Á¥±¡Æ¯Á£Q¦ ” ' ¦ ¨§Y¨­Ê¡Æ¯Á͐ © *¡ŸC£Ê§¡ŽË4­C¯°® ¨«°«¹Ë¯°¦ªÍ£Q«Á͐  ¬#¨Äºd³tЪÐI»…¼4½Ó´n»…ºŽ¾´»b³>¼ÜÂQ¥Y¨É¨¦ª¡Æ   *¯°¦C¡ƟY¨>¡l¡ŸC  *ɍɍ£É ¯°¦ ¬£© ¦C£QÉ~© ¯°¬«° *¬¬.¡ÆŸY¨¦Ô¨ë¯ÁÍQ *¦Ô¡£Q«° *ÉÆ¨¦C®* ÈA¨¦C§ä¨ µÜ³ Ր» ++Ӎ¸ºˆ»b³>¼Ôĺˆ´¸ºŽ¾b"½ ,¡Æ£ëØi ¨­C­C«Á¯° *§ ¯¹¤ ¡ÆŸ±  ¬~¡£Q­C­C¯Á¦CÂä®*ɍ¯¹¡Æ  É~¯°£Q¦ ¯°¬¦C£¡Ü¬Æ¨F¡¯°¬SÎY  § ” ' ©£4§C¯ÁÎC®Š¨F¡¯°£Q¦ ¬~¡ÉƨF¡  ÂË ©¨ŠË ®*£Q¦C®*  ɍ¦ ã4¥C¨¦ª¡Æ¯Á¡¯° *¬ ɍ *«‘¨F¡  §¡£®*£Q©­C¥±¡t¨>¡Æ¯°£¦¬¥C®tŸ ¨¬+¡ÆŸ± #«°£4®Š¨«i© *¬Ÿ¬¯Áå  È±£QÉ(ãW¥Y¨¦ª¡Æ¯¹¡Æ¯° *¬+ɍ *«‘¨F¡  §¡Æ£©£4§C  «Á¯°¦C¬¥C®tŸ ¨¬.¨Ü¡¥CɍØC¥±«°  ¦ª¡#Í4¯°¬~® £Q¬~¯Á¡ŽË ” ž(ŸC ©£Ê§C¯¹ÎY®Š¨>¡Æ¯°£¦ ¬~¡ÉƨF¡  Ëë¤g£QÉz¡ŸC «°£4®Š¨«6©  ¬Ÿ%¬~¯°å  ¯Á¬,£¤…¡  ¦%ØY¨¬~  § £Q¦ ¾¿ ¶4»ˆÕ»°Äºp´n.» - ¶4ºp»b³>¼¢¯Á¡Ÿß¡ÆŸC £Ø    ®¡Æ¯Á͐ £¤E¬Æ¨F¡¯°¬~¤…Ë4¯Á¦CÂܡƟ± ¬~¡Æ£­C­C¯°¦±Â® ɍ¯¹¡Æ *ɍ¯°£¦ß¢¯Á¡ÆŸ%¨©¯°¦±¯°©¨« ¦W¥C©ÞØI  É£¤§C  ɍ   *¬ë£F¤¤gÉ~    §±£Q©ëÈl£QÉ«‘¨É~ÂQ  ¬S¡­i£¬¬¯ÁØC«° /«Á£Ê® ¨«z©  ¬~ŸÌ¬~¯°å*  0 ” / £¡Ÿ ¡ŸC %¬~¡£Q­C­C¯Á¦C ® É~¯Á¡  ɍ¯Á£Q¦¨F¦C§©£Ê§C¯¹ÎY®Š¨>¡Æ¯°£¦®*ɍ¯¹¡Æ  É~¯°£Q¦©Ü¨ŠËØI zØY¨¬~  §£Q¦/¸,ÐY³ÄºŽ¾´n»b³>´»Y *ɍÉ~£QÉ *¬~¡Æ¯Á©Ü¨>¡Æ  ¬E¯°¦ªÍ£Q«ÁÍ4¯°¦±Â ¡ÆŸ± (ɍ  ¬~¯°§C¥C¨«°¬K[ü ø@À #ÃL¨¦±§ [ü  ø @À #ãQɁ¸(ÐI´n»b³>´n»W  É~ɍ£QÉK *¬~¡¯°©¨F¡Æ *¬ ¯°¦ªÍ£Q«ÁÍ4¯°¦±Â  ¬~¡¯°©¨F¡Æ *¬£¤A¡ŸC .   ±¨®¡l¬~£Q«°¥Ê¡Æ¯°£¦ ú  £QÉ(ú ” ï 2” 1436587:9;=?@?;ǘ B-V$ :

û

øÀY.~ÃrnÃ

=%u`~.$|nWÀ 1ªÃ =VuAÀ 1ªÃ



@

 R

\ ƒ \ ƒ ˆ À - $ üO¡z«Á Š¨¬S¡ “ 1 ¢¯Á¡Ÿ §C£Q¥±ØC«° ­Cɍ *® ¯Á¬¯°£¦ß¨É~¯Á¡ÆŸ±© ¡Æ¯°® ” ž(ŸC *¬ Þ *ɍɍ£É¬l¢¯Á«°«[¨®*® ¥C©¥C«‘¨F¡ ¨>¡,¨ÉƨF¡  §C ¡Æ  É~©¯Á¦C  §تˡƟC 9¬~¡Æ¨ØC¯°«Á¯Á¡ŽË­Cɍ£­i *É~¡Æ¯Á  ¬E£¤2¡ŸC 9§C¥Y¨« ” ž6£§C *® ɍ  ¨¬ 9¡ÆŸC¯Á¬@  ɍÉ~£QÉL¢ l©Þ¥C¬S¡O¡Ÿ4¥±¬O¦C£F¡ §C *® ɍ  ¨¬ ¡ÆŸC ¡Æ¯°© ×s¬~¡  ­ÚÈÇØC¥±¡l¯°¦±¬~¡Æ  ¨§ä»…¼2Ӵƾ¸>ľº…ÒY¾ºˆ»…µÜ,¾ +bÄºŽ¾ŽÐ (òŽ¦ë¡ÆŸC¯Á¬9¢¨ŠË¢O ¢¯Á«°«Ú¡t¨Fِ ¤g ¢ *É ¡Æ¯Á© n׎¬S¡Æ  ­±¬+¨¦C§Ü¬£Þ§C *® É~ Š¨¬~ ,¡ÆŸC l® £Q©­C¥±¡Æ¨F¡Æ¯Á£Q¦Y¨«Ç  É~ɍ£QÉ ” ôʯÁ¦C®  l¯°¦C®*ɍ Š¨F¬¯°¦±Â¡ÆŸC 9¡Æ¯Á© n׎¬~¡  ­¢ z¢¯Á«°« ¯°¦±® ɍ  ¨¬ ¡ÆŸC%   O¨«Á  É~Ù4¯Á¦ë *ɍÉ~£QÉ ÈÇ¢ Þ¨«°¬~£ÜŸY¨ŠÍQ ¡Æ£¯Á¦C® É~ Š¨¬~ ¡ÆŸ± Þ£Qɍ§±  Éz£F¤K¡ÆŸC © ¡ÆŸC£4§ ¡£Üِ *  ­/¡ÆŸ±  O.¨«° *ɍÙW¯°¦ *ɍɍ£É(¬©¨«°« ” òަ ¡,¬£© £Ê®*®Š¨¬~¯°£Q¦C¬*È2¬~¥C®tŸ ¨¬z¦C  ¨ÉØªË¡ÆŸC ÎYɍ¬S¡,æY¯°­ß¨F@¡ 7û “  ÈǡƟ± *Ëë§C£¦C£¡zÂQɍ£¢ ©£Q¦C£¡£Q¦C¯°® ¨«°«¹ËQÈ4¨¦C§¡ÆŸW¥C¬l¨¬¨¦ë¨ŠÍQ *ÉÆ¨ÂQ .ÂQÉ~£¢¨F¡9¨©£4§C *ÉÆ¨F¡ #   Ê­i£¦C  ¦ª¡Æ¯°¨«Éƨ>¡Æ  ” !

!

" $

 

" %$

 

" %$

 

 

 !" ##



16

10

14

10

12

10



10

 

10

10





8

6

10

4

10

2

10

0

10

0

5

10

15

20

7

30

35

40

45

50

 K›[˜;  ž(ŸC +¬S¡t¨Ø±¯°«°¯¹¡ŽË,¤ˆ¨®¡Æ£QÉ[¤g£QÉ[¡ÆŸC O® £Q©­C¥±¡Æ¨F¡Æ¯Á£Q¦Y¨«W¨¦±§ÞãW¥Y¨§Cɍ¨F¡Æ¥CÉ~ +  É~ɍ£QÉ~¬ È ¨¬(¤g¥C¦±®*¡Æ¯Á£Q¦£¤6¡Æ¯°© z¤g£QÉ+¡ŸC ö£Qɍ *¦Cå.¬~Ë4¬~¡  © ”





걔 𠔆걔 ÏO³>¼2Ó·ƒ¶WÄ»b³>¼±Ä (÷ä ¦C£¢ ©¨Ù ë¬¯°©­C«Á   ¬~¡¯°©¨F¡Æ *¬¤g£É¡ŸC ëɍ£¢(¡ÆŸ ɍ¨F¡Æ £¤,¡ŸC  ¬S¡t¨Ø±¯°«°¯¹¡ŽË ¤ˆ¨®¡Æ£QÉ~¬ Ȫ¯°¦£Qɍ§C *ÉK¡Æ£#­CÉ~  §C¯Á®*¡LŸ±£¢ ¤ˆ¨ÉK¨F«°£Q¦CÂz¯Á¡L¯Á¬K­i£¬¬¯ÁØC«° E¡Æ£.® £©­C¥Ê¡Æ  ” ¡Æ  ¬9­CÉ~  ¬ *¦ª¡Æ  §/¯Á¦ è “Q“ é ÈI¢O ÎC¦C§ ¡ÆŸC¨F¡l¡ŸC ®*£Q©­C¥±¡t¨>¡Æ¯°£¦Y¨«   É~ɍ£QÉ(® ¨¦Øi .  ¬S¡Æ¯Á©Ü¨F¡  §¨¬ : À êʔÁ“  à   YÀ 7zñ©Ü:+C ¨   l "$#zl%0 

$

"



ƒ

¢ŸC *ɍ .¡ÆŸC .® £©­C¥Ê¡t¨F¡¯°£Q¦Y¨F«2ɍ *¬¯Á§C¥Y¨« " # ¯°¬(§± *ÎY¦C *§ß¨¬ “  ÞÀY.B9 Ã+ ÞÀY.BC UL~Ã+ À êʔÁ“ ) à "$# ÀY.~ÃEû X 9

\Z]

,-

üdÀ/ 0¡t¨¯°¬z® £QÉ~ɍ  ®¡.¡Æ£ÎC͐ £QÉz©£QÉ~ Þ§C¯ÁÂQ¯Á¡¬ È2¢O ÎY¦C§/¡ÆŸY¨>¡ ¡ÆŸ± ôÊ£«‘¨Éô4Ë4¬S¡Æ  © ¯Á¬®*£Q©­C¥±¡t¨FØC«° £Q¦ ¡ÆŸC £QÉ~§C  É£¤  1~1ËQ  ¨É¬ ” òަC® «Á¥C§C¯Á¦C ¨F«°¬£¡ŸC  I£Ê£Q¦2È6¢O  ®Š¨F¦C¦C£¡L®*£Q©­C¥±¡Æ (©£QÉ~ O¡ÆŸC¨¦¨z¤g *¢ ËQ  ¨É¬ ” ž(ŸC (§C¥Y¨«±¬£Q«Á¥±¡Æ¯Á£Q¦ÞÂQÉ~£¢¬L«Á¯°¦C  ¨É«¹Ë.ØY¨®tÙª¢¨É§¯°¦¡¯°© È ¬ *  ¡Æ¯°£¦Y¨«¡+¡ÆŸC¯Á¬£Q­I  ɍ¨F¡Æ£É® £©©Þ¥Ê¡Æ  ¬O¢¯Á¡ÆŸÜ¬­Y¨®* .¨¦C§¡¯°© l§C ¯¹â2 *ɍ *¦W¡¯‘¨F¡¯°£Q¦ ” ' ­±­C«ÁË4¯°¦±Â ¡ÆŸ±¯°¬Ü£­i *ÉÆ¨F¡£QÉ¡£ Àˆð ”¹“ 'à &iÀˆð ”†ê ÃÞ¢ ßÎY¦±§ ¡ÆŸY¨F¡¡ÆŸC  ɍ¥C¦C¦±¯°¦C ¨ŠÍQ *ÉÆ¨Â ßú ¬Æ¨F¡¯°¬~ÎC  ¬¡ÆŸ±  ¤g£Q«°«Á£¢¯°¦C   ãW¥Y¨>¡Æ¯°£¦äÀg©£4§C¥C«Á£Øi£¥C¦C§Y¨ÉSË *âi  ®¡Æ¬ÆÃ Àˆð ”  à ú  $ ú  û ú  $  =.ú  ú  ûÌüKÀˆú  à $+4Àˆúià 0 ú  À 0e1ªÃ@û ú : À ,à 0 ¢ŸC *ɍ/    ¯°¬¨/¬~¯°©­C«°¯¹ÎY  §ì£­i *ÉÆ¨F¡£QÉ£Q¦ ¡ÆŸC  ¬~®Š¨«Á    ɍ *¬¥C«¹¡Æ¯Á¦CÂߤgɍ£Q©  ¨­C­CÉ~I£ ʯÁ©Ü¨F¡¯°¦C0  Ìت Ë  ¨¦C§¡ŸC ,® £ÉÉ *®*¡¯°£Q¦¡   ɍ1 © 4Àgú2ÃLû'ÀbüKÀgú2Ã~à äüKÀgú à $2 ú 3À @úià ®*£Q¦ª¡t¨¯Á¦C¬¡ŸC #¯°¦±æC¥C  ¦C®* #£¤ ¡ÆŸ± z¥C¦CÉ~  ¬~£Q«Á͐  §£Q¦Üú ” ÷Ô ,®*£Q¦C¬¯Á§C  ɨ® £©­C¥Ê¡t¨F¡¯°£Q¦Y¨F«Y­Cɍ£Qر«°  © ¢¯¹¡ÆŸC£Q¥Ê¡¬~¥CØCÂQÉ~¯°§©£Ê§±  «Ç£¤¡ÆŸ±  ¤g£QÉ~© Àˆð ”  à ú  $ ú ,ûŒüKÀˆú à 0 ú 4À 0e1ªÃLûÌú : À ,à 0 ¨¦C§ë¨® £QÉ~ɍ *¬­I£Q¦C§C¯Á¦C­Cɍ£Qر«°  © ¢¯Á¡ÆŸ¬~¥CØCÂQÉ~¯°§©£4§C *«Ú£¤6¡ÆŸC #¤g£QÉ~© Àˆð ”  à ú  $+  ú zû üKÀ ú à $  4À ú ,à 0 ú 4À 0e1QÃLûÌú : À à 0 ¢ŸC *ɍ   4À ú >Ã+¬Ÿ±£Q¥C«°§ë¨­C­±É£ ʯ°©¨F¡Æ(   4Àˆúià ” òަ¡ŸC¯°¬¦C£F¡Æ #¢ ¢¯°«Á«2® £¦C¬¯Á§C  Él¨¬~¥CØCÂQÉ~¯°§©£4§C *«Ú£¤ ¡ÆŸ± ,¤g£Qɍ©  4À ú ÃE5 û 426À ªÀ ú >,à 0- 3 4À ú à 0- 4À ú ÍÃOØY¨F¬  §£Q¦¨ŠÍQ *ÉÆ¨Â¯°¦C£Q¦¡ŸC ,®*£ª¨É¬~  É(¬~®Š¨«Á  ¬9ï  









 

 

 !" ##

¨¦C§ ð ÚÈI¢ŸC *ɍ ¡ŸC ¤g¥±¦C®*¡¯°£Q¦ 4¯°¬9§C  É~¯ÁÍQ *§/ØY¨¬~  §/£Q¦ ¨Ü¬® ¨«° ɍ *ÂQ¥C«‘¨Fɍ¯Á¡ŽË¨F¬¬¥±©­±¡¯°£Q¦ë£Q¦ 4Àˆúià ” ž(ŸC¯Á¬¡ŽË4­I .£¤& M$;) 4Àgú2à J M¼ w[;:ȡƟY¨>¡ X ¨®t= ¯°¬*È  #û ¬~­Y¨¦  WEC À Áûþï  #Àpï  X ;à  ” ž(Ÿ± Þ¤g¥C¦C®*¡¯°£Q¦±¬ WEC ¤g£Qɍ© ¨) ¦  3 ×s£QÉ~¡ŸC£Q¦C£QÉ~©Ü¨F«6ØC¨¬¯Á¬ ¯°( ¦  Èi¨¦C§¢ Þ§±  ¦C£¡ ¡ŸC £QÉS¡ÆŸC£ÂQ£Q¦Y¨«[®*£Q©­±«°  ©  ¦ª¡9£:¤  #¯°/ ¦   L(ت( Ë %>ÈÇ¢ŸC¯°®tŸß¯Á¬l  ¦C *ÉÆ¨F¡  § ØªË ¨¦±£¡ÆŸC *É£QÉ~¡ŸC£Q¦C£É©¨«OØC¨¬¯Á¬Àg¡ŸC EL ¸Å¾·°¾º…Ä>ÃU EC  À Ãûîï  * Àpï  X ÃÈL¢ŸC *ɍ  [ % :¯Á¬ ®Š¨F«°«° *§ ¡ŸC  µÜ³>º…ÒY¾´ E¸Å¾·°¾º ” % û %  %  3 %  ÈA¢ŸC *ɍ Ü¡ŸC<   % †¬ɍ *­Cɍ *¬  ¦ª¡§±¯Áâi  ɍ *¦C®  *¬ ¯°¦¡ÆŸC ŸC£QÉ~¯°å £¦W¡Æ¨«bÈY͐  É~¡¯°® ¨«¡  ãW¥Y¨«Á%¬  4Àˆúiï°¦¡ŸC .©¯Á§C­i£¯°¦ª¡Æ¬+£F¤A¡ŸC .  «° *© *¦ª¡Æ¬£¤ Q ÈC¢O ŸC¨͐   4ÀˆúiÃLûçè|üKÀˆúià é  üKÀè ú é  +à Àè f=.ú é  è  é  =4è ú é  ,à 0 ¢ŸC¯Á®tŸ¤g£QÉ(¬~  ®*£Q¦C§ £É§C *É9É~ Š¨®¡Æ¯Á£Q¦¡Æ  É~©¬9ü¡Æ¯°£¦¤gP É~£Q© è Fé ¬ŸC£¢¬E¡ÆŸY¨F¡  4À ";0 ÷ƒºd¾´n´~!¸ + [³>º-,¸/. * 10 9÷Ô ®*£Q¦C¬~¯°§C *Él¨É~ Š¨®¡Æ¯°£¦§±£Q©¯Á¦Y¨F¡  §­CÉ~£QØC«° *©"£¤6¡ÆŸC #¤g£QÉ~© ú L ú L@û ú L*À “ ú 3 ,à 0 ú 3 ú 3 $ =.ú 3 û ú 3 Àg ú L+ “ ,à 0 ú * 243 û 1 0 úAÀ 0e1ªÃLûçÀ % ! 3 C  # À à 0 “ ,à 0 c U ¢ŸC *ɍ2    û “ 1 È,¢ŸC¯°®tŸÌ® £QÉ~ɍ *¬­I£Q¦C§C¬¡Æ£ ¡ŸC  ® «‘¨F¬¬¯Á®Š¨«65@£«Á¡Æ *ɍɍ¨>×dö£¡ÙF¨ ¬SËʬS¡Æ *© ¢¯¹¡ÆŸ ¬©¨«°« §C¯¹ â2¥C ¬~¯°£Q¦Ü¨¦C§®*£Q¦ªÍQ  ®¡Æ¯Á£Q¦¯°¦£Q¦C ,®*£Q©­I£Q¦C *¦ª¡ ” ÷Ô zŸY¨ŠÍQ(   4ÀˆúiÃLû'À ÞÀgú LSú 3 à  $ ú L ú 3 0Àˆú Lsú 3 à  ú L ú 3 ÃÈl¨¦C§ ¢¯Á¡Ÿ ¬¥CرÂQɍ¯Á§ ©£4§C  «9¢ %ŸY¨ŠÍQ /¡ŸY¨F¡ 

4À ú Ã/û À  4À ú L0 ú 3 ,à 0  4À ú Le0 ú 3 ÍÃȁ¢¯Á¡Ÿ c ûõðìÀgɍ ¤g  ɍ *¦C®  ¬® ¨«° ©¯Á¦W¥C¬®*£Q©­C¥±¡t¨>¡Æ¯°£¦Y¨«K¬~®Š¨«Á à ” Ç£É.¡ŸC  ¬~ Ü­CÉ~£QØC«Á  ©¬.¢O Ü¥±¬ ¨®  *¦ª¡ÆÉƨF« §C¯¹â2 *ɍ  ¦±®  × g ÉÆ¨¦±ÙªG× l¯Á® £Q«Á¬£Q¦¬®tŸ±  © ¤g£ÉK¡ÆŸC (©¯°§C­I£Q¯Á¦W¡¬K£¤I¡ÆŸC ( *«° *© *¦W¡¬ Ȫ¢ŸC *ɍ (¢ 9ɍ *ª¨É~§¡ŸC  ¬~  ©¯°§C­I£Q¯°¦ª¡,ÍF¨«Á¥C  ¬.¡Æ£ɍ *­CUɍ7  *¬ *¦W¡¨ë­C¯Á  ®  ¢¯°¬~ ®*£Q¦C¬~¡Æ¨¦ª¡¨­±­Cɍ£ ʯ°©¨F¡Æ¯Á£Q¦/£ÍQ *ɡƟ±  *«°  ©  ¦ª¡¬ ÈK¨F¦C§ ¡ÆŸ± lÉ~ *¤g  É~  ¦C®* ,¬® ¨«° ¯Á¬+ï ¯°¦¡ÆŸ± z® £©­C¥Ê¡t¨F¡¯°£Q¦£¤i¡ÆŸC 9  É~ɍ£QÉ ” ž(ŸC z¬~£Q«°¥Ê¡Æ¯°£¦C¬E¨É £Q¬~® ¯°«Á«‘¨F¡¯°¦CÂ.¨F¦C§ ØI£¡ÆŸ/ú Ll¨¦C§/ú 3 ¨É~ ¤gÉÆ¨F ®*¡t¨F«6¤g£Q2É .  1±Èi *ÍQ *¦ ¡ŸC£Q¥CŸ/ú 3 À 0e1ªÃû “Q” ÷Ô ¢¨¦ª¡#¡Æ£¨­±­Cɍ£ ʯ°©¨F¡Æ  ú ÈK¨F¦C§ä¡ÆŸ± Ü  É~ɍ£QÉ~k ¬ lnú  l¨É~ Ü­C«Á£¡¡  §ä¯Á¦ A¯°Â ”¹“ ×Sï4È[¢Ÿ±  ɍ  õ¯Á¬¡ŸC   ɍÉ~£QÉ¢¯¹¡ÆŸC£¥±¡©£4§C  « aÀ 98 `ÃȪ¢¯Á¡ŸC£Q¥±¡@©£4§C  «Cر¥±¡E® £Q©­C¥±¡  §£Q¦¡ÆŸC ÎY¦C *ÉE¬~®Š¨«Á   Zï.¨F¦C§¡ŸC  ¦­CÉ~£    ®¡Æ  §£Q¦ª¡Æ£.¡ÆŸC ¬~®Š¨«Á   aÀ |t£ `ÃÈ6¢¯Á¡ÆŸ%¡ÆŸ± Ü¬¥±ØCÂQɍ¯Á§ä©£Ê§±  «9aÀ  $ kÃnÈK¨¦±§Ô¢¯Á¡Ÿì¨¬¯°©­C«Á¯ÁÎY *§ ©£4§C  « 

4À ú Ãû  3 4À ú ~0 ú à aÀ ñ¡ÆÉ~¯‘¨¦±ÂQ«° *>¬ `à ” ÷ä +ÎYÉ~¬~¡@«Á *¡ ß û 1#¨¦C§® £©­C¥Ê¡Æ +¡Æ£ 7 û ïÊȐ¡ŸC  ¦¢O («° *!¡ ëØI ¨,ɍ£¡Æ¨F¡Æ¯Á£Q¦Y¨«4©W¯ ʯ°¦C £¤£QÉ~§C  É Ú$ ó Ìû  Àg¬¯Á¦[À LÃ±®*£Q¬À 3 ,à 0P ß® £Q¬ŠÀ LÃ±¬~¯°¦ÚÀ 3 Ã~萦C§äŸC  É~ ¢ £¦C«ÁË ® £©­C¥Ê¡Æ ¡£ 7þû “ ÈÚ¬~¯°¦C®* ¨F¤…¡Æ *@ É 7þû “ ¡ÆŸ± ¬¥CØCɍ¯°§/¬~®Š¨«Á  ¬.¨É~ §C£Q©¯°¦Y¨>¡Æ  § ØWËë¡ÆŸC ®*£Q¦ªÍQ *®*¡Æ¯¹ÍQ ¬~¡É  ¨ÙW¬§±¥C  ¡Æ£  ” ÷ä ¬~¡¥C§±Ë¡ÆŸ±  *ɍɍ£É¤g£É( Š¨®tŸë® £Q©­i£¦C  ¦ª¡9¯°¦C§±¯ÁÍ4¯°§C¥C¨«°«¹ËQÈC¨F¦C§¤g£É   û 1¢O .ÎY¦C§ë¡ŸY¨F¡¡ÆŸ±  ¬£«°¥±¡¯°£Q¦¢¯Á¡ÆŸ¡ÆŸC +¬~¥CØCÂQÉ~¯°§Þ©£4§C  «4¯°¬6¡ÆŸC +ØI  ¬S¡K¤g£QÉAØi£¡ŸÞ® £©­I£Q¦C *¦W¡¬ ÈQ ÍQ *¦¡ŸC£Q¥CŸÞ¡ÆŸC +©£4§C  «Á¯°¦C   É~ɍ£QÉ~¬9¨É #¬~©Ü¨«Á«° *É+¯°¦ú 3 ¬¯°¦±®  #ú 3 À 0#1ªÃ+¯°¬(® £¦C¬~¡Æ¨¦ª¡ ” Y£QÉ w û : 1¡ŸC .¬£Q«Á¥±¡Æ¯Á£Q¦¢¯Á¡Ÿ¡ÆŸC .¬¥±ØCÂQɍ¯Á§ ©£Ê§C *«C¯°¬LØI  ¬S¡O¤g£ÉLú LAØC¥Ê¡ŠÈ4Øi *®Š¨¥C¬~ l£F¤i¡ŸC ® £Q¦ªÍQ *®*¡¯°£Q¦ÚÈ4¬¥±ØCÂQɍ¯Á§¬®Š¨F«°  ¬E¯°¦ú 3 §C£¦C£¡@§C ÍQ  «Á£Q­¨F¦C§ ¡ÆŸ± ¬~£Q«°¥Ê¡Æ¯°£¦C¬,¢¯Á¡ŸÔ¬~¥CØCÂQÉ~¯°§ ©£4§C  «Á¬#§C£4  ¬.¦C£¡§C¯¹â2 *ɬ~¯°ÂQ¦±¯ÁÎY® ¨¦ª¡Æ«ÁËë¤gɍ£Q© ¡ÆŸC ¬£«°¥±¡¯°£Q¦ß¢¯¹¡ÆŸC£Q¥Ê¡ ¬¥±ØCÂQɍ¯Á§©£4§C  « ” ž(ŸC ¬£Q«Á¥±¡Æ¯Á£Q¦£Q¦ ZQﯰ¬(ØI *¡¡  Ɂ¬¯Á¦C®  #¡ŸC .¦4¥±© *ɍ¯°® ¨«Ú *ɍɍ£É(¯°¬+¡ÆŸ±  ¦ɍ  §±¥C®  *§ ” ð ” ð ” ï ”¼YŐ¾tÓºˆ»b³>¼=. ; Ï>0 (÷Ô ¨«Á¬£® £¦C¬¯Á§C  É.¨® £¦W͐  ®¡Æ¯°£¦Ô§C£©¯Á¦Y¨F¡Æ *§ß­CÉ~£QØC«° *©í£¤E¡ÆŸ±  ¤g£QÉ~© Àˆð ”Á“Q“ à ú $+k=.ú  úû “ 0 ú * * F Vt:+C * G Vt:+w û 1Z0 cú * * F VLBC * G VLL û 1Z0 úAÀ 0#1ªÃLb û 1Z0 "

$





















 

 

 !" ##

%

g¤ £QÉ ìû À % !3 C  # 0;% !3 C  # Ãl¨¦C§ l û “ 1 U Èi¢ŸC¯Á®tŸ ¢ ¬£Q«¹ÍQ ØªË ¨ô4¡É Š¨F©«Á¯°¦C  ,¯Áâi¥C¬~¯°£Q¦®ROÀ “ Ãd®ROÀ “ à ©U  * ¡ÆŸ±£Ê§ è “ é ¢¯¹¡ÆŸ%ØC¯°«Á¯°¦C  ¨Éz *«° *© *¦W¡¬ ” òަ/¡ŸC ® £Q©­C¥±¡Æ¨F¡Æ¯Á£Q¦/£¤@¡ÆŸ±  *ɍÉ~£QÉz¡ÆŸ± É~ *¤g *ɍ  ¦±®  ¬® ¨«° ¯Á¬ ï ” ž(Ÿ± ¬~£Q«°¥Ê¡Æ¯°£¦Ü¯°¬¯°¦¡ÆŸ±¯°¬+® ¨¬ #É~  «‘¨>¡Æ¯Á͐  «¹Ëܬ~©£4£ ¡ÆŸÜ¬¯°¦±®  , ¡ÆŸ± ,¤gÉÆ¨®¡t¨« Ô ¯Á¬+£Q¦C«¹Ë¨®¡Æ¯Á¦CÂÞ£Q¦¡ÆŸ±   §C *ɍ¯ÁÍF¨F¡¯Á͐  ¬9£¤K¡ÆŸ± ¬£«°¥±¡¯°£Q¦ ” ÷ä ŸY¨ŠÍQ   4ÀˆúiÃû  =ÊÀ6  úià ÌÀ3 =  úià È2¨F¦C§ë¡ŸC  *ɍɍ£Él¯Á¦ë¡ÆŸ± 



¬£«°¥±¡¯°£Q¦¢¯¹¡ÆŸ¡ÆŸC .¬~¥CØCÂQÉ~¯°§©£4§C  «2¯°¬(¬©¨«Á«°  ÉE¡ÆŸY¨F¦ë¯Á¦¡ŸC .¬£Q«Á¥±¡Æ¯Á£Q¦¢¯Á¡ŸC£Q¥±¡¬~¥CØCÂQÉ~¯°§©£Ê§C *« ” −4

4

x 10

−4

1.8

3.5

1.6

3

1.4

x 10

1.2

2.5

1

2 0.8

1.5 0.6

1 0.4

0.5

0

0.2

0

2

4

] 0 ’ 8 hwuF_ˆh*uW’6€  7 7R9 ‚ KL;`*)GI Hor^ka 7 c‰>‡yb^`] D h]\~v†^`u‡ Ž’! "] *a"/(-0 K , =\;."@% "!(- V =E.;.(-+KL( ." ? C7R7 7? ebhonjj\~oy2^`u 0 hr+jW’ D \depw4’ 3 jj>vÁ’ D \~aswW’ 5 u‡ybu‡F’  9#@8 ’ ; h{r(onu4’4€     ‚ G ‹Šu>onr+^ka 7 c‰‡*yp^k] D h]\Sv`^†u>$ ‡ % Ž!’ "$*)"/( 0K ,- = ;."@%&"L( 8V = ;(`KL.(-" ’ (' @8 ’ ; h{r(ouW’±€ R ‚ G ‹Šu>onr+^ka 7 c‰‡*yp^k] D h]\~v†^`u‡+mxhy 0 huŠq\aŽeb^†h*u„ G ^†{c>_p^`hu 5 } c>onep^`hu>_Úf^†epw B ybo*aŽesonv 0 h„ \ )(aS^†\~u eb_~’  GI=  < %.F@ = D ()#* GI=  "]/%.=\GI ,-C#% ?Ç7 j>yp^`u‡\~y 1 \aŽebcyp\ 7 \Syb^†\_L^†u 5 u‡^`u\~\Syb^†u‡ ?Ê7 jyb^`u‡\~y , \~ypvkon‡F’  #@G ’ 5 _ˆep\~jzou>] 0 ’ 8 h*wu>_phuW’ª€ C7 7  2 ‚ 3Úw\Ljh^`u*egf^k_p\KaShr+jceson‰^`v`^|eg‹hmÊebw\ 1 h*yp\~-u K_ˆ‹_ˆep\Sr,’ .Z+_0/ ? & ' ' 5   9 ’   7#@D ’ J ’ 1 onys_phuW’4€ R  ‚ 5 ypybhy J ybhtfÚepw.ou>] 3 Ch _geb\Syb^†h*yp^ 5 ypybhy 5 _geb^†r(onep\~_[mxhy 0 hu>_p\Sybqtotep^`q\ J ov†\~%y Š^`u 3 j„ jybh 1^†r(onep^`hu>_ihnm C\Syb^†h]^ka : yp‰^†eb_Ú^`u ; or^`v†ephu>^`ou 7 ‹Š_ˆep\~r(_  . Z+_3254 ?  /5  9 ’  C#@3 ’ 1 h*‡‡>’ª€ C7 7  ‚ 3/D cv†ep^†„ 3 ]onjep^`q\ :]\~yb_ 1 h‡‡  GI= Hd Z #

]*a B8V "/"$"K (-H"$=RY" L V2% 3 u>]\~yb_ 1 h‡‡ '() ; ,-$ %/%/ F.=E "$!=E.; G huFonvk] 3 ’ B yb\Su>asw ?F7 ep^`‡ 1 oyb_phuzonu>] H ^kaSoyb]h ; ’ = hasw\depeph  =E!@ = ()#& F@ = < %.@#DZ &%/"/,- T R;GI;  "/; Y" "!(  V =E.;.(-LZ2*R*)"$'R;  ,-&Z ()) ;= &%/] Y"D /;otes_(PU J  ’ 1= o()yb@_p h"u V =E.;.(-aY"8 , V =E %@H,-C#J ()$#  GI= Hoyp^†m±ou>] = ^`v`_ˆ„ 5 yb^(Z^`‰F\~yp‡ A(- U"$  _p‰h ? onuF] D oneb_ J ’ 1 onys_ˆh*u Z ,(a (= &a%/ < % #;Y[. X(Y%/ = < *a ] 3 ’ 8 huFo_ = ^ Šv`o*_p_phu Z= ;(a %.= S &%/%/< % ( .Y = H T () ((a-#[%/.;().;" U[H"/.,-C#Y"D=>= < *a ,-C#LY"L ; ( ]F@; *-&%/H,-C#%LY"X(aF]; *a"$%]%/< .=E  GI