next up previous
: 指数部つきの表現 : Fortran 文法 : do while 文

分割コンパイル

これまでの積分や、ニュートン法では、メインルーチンは同じで、関数の部分だ けを関数副プログラム、又はサブルーチン文として与え、変更することができた。 これは別に1つのファイルとする必要はない。例えば、function.f を

      real*8 function f(x)
      implicit real*8 (a-h,o-z)
      integer n
c
      n=3
      f=x**n-n
      return
      end
c
      real*8 function df(x)
      implicit real*8 (a-h,o-z)
      integer n
c
      n=3
      df=n*x**(n-1)
      return
      end
として、
  %frt newton.f function.f
としてコンパイルすると、$f(x)=x^n-n=0$ (この場合$n=3$)の解、つまり $n^{1/n}$ が得られる。 関数部分を変えると任意の関数についてニュートン法が実行できる。



Kiyohide Nomura 平成15年6月25日