¦P¾Ç­Ì±`¥Ç¿ù»~¤§´£¥Ü

 

Fortran µ{¦¡»y¨¥¤è­±

¨S¦³«Å§i implicit none ¡A³y¦¨¦³¨Ç¥´¿ùªºÅܼƦWºÙ¤´³Q·í§@·sÅܼƨϥΦӵo¥Í¿ù»~

µ{¦¡¦æªº¤å¦r¦ì¸m¶W¹L²Ä 72 ®æ¦Ó³Q·í§@¬O¤£§¹¾ã

¥Dµ{¦¡ªº¦WºÙ¡A¤£¥i¥H»P¤§«á·|¥Î¨ìªºÅܼƦWºÙ­«ÂÐ

¤j¤p¼g¦b Fortran ¬O¤@¼Ëªº¡A¬GÅܼƪº t »P T ´N·|³y©Î­«ÂЫŧi

¬O°}¦Cªº²Ä i ­Ó¤¸¯À A(i) ­n©ñ¨ì x ¥h¡A¥»À³¥Î x = A(i)¡A¦ý«o¼g¦¨ x = A¡A½sĶ¾¹·|©ê«è§A§â 1 ¯Å±i¶q¼g¨ì 0 ­Ó¯Å±i¶q¥h¡C

 

PGPLOT ¤è­±

§Ñ¤F«Å§i pgopen ¬°¾ã¼Æ©Î if (pgopen('/xwin').le.0) stop ªº»yªk¼g¿ù

pgenv ªº¤W¤U­­½d³ò»~¥Î¤F¾ã¼Æ¡A¤@©w­n¥Î¡]³æºë«×¡^¹ê¼Æ¤~¥i¥H

pgpt ªº¤Þ¼ÆùØ¡A§Ñ¤Fµ¹³Ì«á¤@­Ó¾ã¼Æªº²Å¸¹­È©ÎÅܼÆ

³Ì«á¨S¦³¤U pgclos¡A³y¦¨¹Ï¤@°{´N°¨¤WÃö±¼¡A¨Ó¤£¤Î¬Ý¡C

 

 

½sĶ¡BÁåµ²¤è­±

Áåµ²ªº¨ç¦¡®w¨S¦³µ¹§¹¾ã©Î¸ô®|¥´¤F¿ù¦r

-o «á­±¥»À³¬O my_prog.x ¡]«ü©wªº°õ¦æÀɦWºÙ¡^¦ý«o¤£·V¥Î¤F my_prog.f ¾É­P³s­ì©lµ{¦¡½X³£³Q»\±¼¤F¡C

 

 

µ{¦¡ÅÞ¿è¤è­±

¨Ï¥Î Euler ©Î Euler-Richardson ºtºâªk®É¡A¨S¦³¨Ï¥Î y = y + f * dx §Î¦¡¦Ó¥u¼g¤F y = f * dx¡A§Ñ°O¤Fªº y¡C

ªì©l­È§Ñ¤F¦b°j°é¥~¥ýµ¹¡A©Î¬O§âªì©l¤Æªº°Ê§@¼g¦b°j°é¤º¡C

°j°é¤¤©Ò­n¥Îªº¦ÛÅܼơA­ì¥»À³¨Ï¥Î²Ö¥[«¬ªº x = x + dx ©Î¬O§Q¥Î¤F°j°é­p¼Æ¾¹ (loop counter) ªº x = x0 + i * dx¡A¦ý«o²V¥Î¤F¨âªÌ¦Ó³y¦¨¿ù»~¡C