²Ä¤G¶g¤W½Òºõn
¥»¦¸¹q¸£¼ÒÀÀ¹êÅç½Ò§ÚÌn¦³¨t²Î¥B§Ö³t¦a¦V¤j®a¤¶²Ð Fortran (77) ªº¤@¨Ç³Ì°ò¥»ªº«ü¥O¡A¦³¤F³o¨Ç«ü¥O¤j®a´N¥i¥HÂ˦U¦¡¦U¼Ëªº¹q¸£pºâµ{¦¡¡C§Ṳ́]·|¯d®É¶¡µ¹¤j®a°µµ{¦¡¼¶¼g¡]¥Î vi¡^»P½sĶªººt½m¡A°£¤FFortran ¤§¥~¥t¤@Ó¤j«ÂI¬O pgplot ø¹Ï°Æµ{¦¡®wªº¸Ô²Ó¤¶²Ð»P¨Ï¥Î¥Ü½d¡C
¥Î vi ¼g Fortran µ{¦¡¦b Linux ¤U½sĶªº½m²ß¡G
¹q¸£pºâ³¡¤À
a. ¾ã¼Æ±q 1 ¥[¨ì 10 ¨Ã¥B§âµª®×¦b¿Ã¹õ¦L¥X
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
b. §â¤WÃD§ï¬°¦V¨Ï¥ÎªÌ¸ß°Ý±q 1 ¥[¨ì N ªº N È¡A¨Ã¥B§âµª®×¦b¿Ã¹õ¦L¥X
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
c. ¦P¤W¡A¦ý§â¾ã¼Æ±q 1 ¥[¨ì N ªºµ²ªG¼g¨ìÀÉ®×
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
d. §â¤Wzµ{¦¡ªº³s¥[³¡¤À¼g¦¨¤@ӰƵ{¦¡
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
¹q¸£Ã¸¹Ï¡]pgplot¡^³¡¤À
e. ©I¥s pgplot °Æµ{¦¡¨Óµe´XÓÂI
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
f. ©I¥s pgplot °Æµ{¦¡¨Óµe³sµ§¦±/§é½u
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
g.©I¥s pgplot °Æµ{¦¡¨Óµe¤@Ó¶ê
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
h.©I¥s pgplot °Æµ{¦¡¨Ó§@°Êµe
°½¬Ý¤@¤U°Ñ¦Ò½d¨Òµ{¦¡
¡]¤W½Ò·í¤Ñªº¹ê»Ú¶i«×¬O¨ì³oùØ¡^
²Ä¤@Óºtºâªk¡G©@°Ø§N«o°ÝÃD
°ÝÃD¡G§N«o²{¶H
ª«²z¶q¡G·Å«×ÀH®É¶¡¦ÓÅܤƪºÃö«Y¡AT(t)
¼Æ¾Ç¼Ò«¬¡G¤û¹y§N«o©w«ß¡A·Å«×ÅܤƲv»P
¼Æ¾Ç°ÝÃD¡G¤@¶¥±`·L¤À¤èµ{¦¡¡Gª¾¹D±×²vªº¸Ü¡A¯à±À±oì¨ç¼Æ¶Ü¡H
¹s±×²vªºì¨ç¼Æ¬O
±`¼Æ±×²vªºì¨ç¼Æ¬O
x ¤@¦¸¤è±×²vªºì¨ç¼Æ¬O
Eular ºtºâªk¡G§â¿n¤ÀÂ÷´²¤Æ
ì¨ç¼Æ y(x) ¦b x0 ³Bªº±×²vÈ g(x0) ¬O
[ y(x0+dx) - y(x0-dx) ] / 2x0¡A¤]¥i¼g§@
[ y(x0+dx) - y(x) ] / dx ¡]¤Ï¥¿ dx «Ü¤p¡^
¡A¦]¦¹
y(x0+dx) = y(x0) + g(x0) dx
¦pªG§ÚÌ©w¤@®æ xn ¨ì¤U¤@®æ xn+1 ¤§¶¡ªº¶¡®æ¬O dx ¡A«h¤W¦¡Åܦ¨
y(xn+1) = y(xn) + g(xn) dx
¦Ó³o¥¿¬O Eular ºtºâªk¡A§i¶D§Ú̦p¦ó±q±×²v¨ç¼Æ¤@¨B¨B±À¥Xì¨ç¼Æ¡C
Eular °Æ µ{¦¡¡G§â²{¦bªº x¡By¡Bg¡Bdx ȶǤJ¨ÓÀò±o ¤U¤@Ó y ªº°Æµ{¦¡
subroutine eular (y,x,dx)
real x, y, dx
real slope
slope = 2*x
y = y + slope*dx
end
§N«o«Y¼Æ r ªº³]©w¡G¦p¥i±q¹êÅç©Ò¶q´ú¨ìªºÈ¨Ó©w¥X¡H
¸Ñ¨M°ÝÃD¡G¨D¥X·Å«×ÀH®É¶¡ªºÅÜ¤Æ T(t)¡A¨Ã¦^µª¬ÛÃöªº°ÝÃD
ºtºâªkªº¶i¤@¨B¹B¥Î¡G¤O¾Ç¤¤ªº¸¨Åé°ÝÃD
µLªÅ®ðªý¤O¤§¦Û¥Ñ¸¨Åé¡Gµ¥¥[³t«×¹B°Ê
¥Ü½d¡]¤U¦¸¤§¹ê§@¡^¡G¦hÁû²É¤lªº¤À¤l°Ê¤O¾Ç¼ÒÀÀ