¤p¦Ò 20070411
(20070509)
¤@¡B¦³¤@Ó¤G¶¥±`·L¤À¤èµ{¦¡ d2y/dt2 = -9.8¡A¦ÛÅÜ¼Æ t ªº½d³ò¬O t = [0.0, 10.0]¡A ¨Ã¥Bn¨Dnº¡¨¬Ãä¬É±ø¥ó 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 ¨Ã¥Bn¨Dnº¡¨¬Ãä¬É±ø¥ó 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 ¨Ã¥Bn¨Dnº¡¨¬Ãä¬É±ø¥ó 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 ¨Ã¥Bn¨Dnº¡¨¬Ãä¬É±ø¥ó 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