¤p¦Ò 20070411

(20070509)

 

¤@¡B¦³¤@­Ó¤G¶¥±`·L¤À¤èµ{¦¡ d2y/dt2 = -9.8¡A¦ÛÅÜ¼Æ t ªº½d³ò¬O t = [0.0, 10.0]¡A ¨Ã¥B­n¨D­nº¡¨¬Ãä¬É±ø¥ó y(0.0) =0.0¡By(10.0) = 0.0 ¡A½Ðµe¥X y(t) ªº¦b¸Ó t ½d³òªº¹Ï¡C

¤G¡B¦³¤@­Ó¤G¶¥±`·L¤À¤èµ{¦¡ d2y/dx2 = 2.0¡A¦ÛÅÜ¼Æ x ªº½d³ò¬O x = [0.0, 1.0]¡A ¨Ã¥B­n¨D­nº¡¨¬Ãä¬É±ø¥ó y(0.0) =0.0¡By(1.0) = 1.0 ¡A½Ðµe¥X y(x) ªº¦b¸Ó x ½d³òªº¹Ï¡C

¤T¡B¦³¤@­Ó¤G¶¥±`·L¤À¤èµ{¦¡ d2y/dx2 + x = 0¡A¦ÛÅÜ¼Æ x ªº½d³ò¬O x = [0, p]¡A ¨Ã¥B­n¨D­nº¡¨¬Ãä¬É±ø¥ó y(0) = 0¡By(p) = 0 ¡A½Ðµe¥X y(x) ªº¦b¸Ó x ½d³òªº¹Ï¡C

¥|¡B¦³¤@­Ó¤G¶¥±`·L¤À¤èµ{¦¡ d2y/dx2 + x = 0¡A¦ÛÅÜ¼Æ x ªº½d³ò¬O x = [0, 2p]¡A ¨Ã¥B­n¨D­nº¡¨¬Ãä¬É±ø¥ó y(0) = 0¡By(2p) = 0 ¡A½Ðµe¥X y(x) ªº¦b¸Ó x ½d³òªº¹Ï¡C

 


¤@¡BAns.

With EPS=1.0E-7



With EPS=1.0E-6 (¬O shoot.f ùتº¹w³]­È)


 


With EPS=1.0E-8

 


¤G¡BAns.


Be careful sometime a bad initial guess can results in no solution found. (in this case V=0.5)


¤T¡B



¥|¡B