¦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¤@©wn¥Î¡]³æºë«×¡^¹ê¼Æ¤~¥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