Fortran Samples

program prog1 
  print *, 'Hello, World (1)' 
  print '(A)', 'Hello, World (2)' 
  print '(A8)', 'Hello, World (3)' 
  print '(A)', 'Hello, World' 
  print '(A$)', 'Hello, ' 
  print '(A)', 'World' 
  stop 
end 

• プログラムは

  program プログラムの名前
    プログラムの本体
    stop 
  end 

  という形式になる。

• 何でもいいから適当に出力したいときは

    print *, 出力したいものをコンマで区切って並べる

• 文字列はシングルクォート（一重引用符）でくくる。文字列を出力する正式な方法は

    print '(A)', '出力したい文字列'

  となる。

program prog2 
  implicit none 
  real(8):: xvalue 
  print '(A$)', 'Input X: ' 
  read (*, *) xvalue 
  print '("X = ", E15.7)', xvalue 
  stop 
end

program prog3 
  implicit none 
  integer:: counter 
  real(8):: xvalue 
  do counter = 1, 5 
    xvalue = 1./dble(counter) 
    print *, counter, xvalue 
  end do 
  stop 
end

program prog4 
  implicit none 
  integer, parameter:: IMAX = 10 
  integer:: i, isum 
  real(8):: dsum 
  isum = 0 
  dsum = 0.d0 
  do i = 1, IMAX 
    isum = isum + i 
    dsum = dsum + dble(i) 
    print *, i, isum, dsum 
  end do 
  print *, 'ISUM = ', isum 
  print *, 'DSUM = ', dsum 
  print *, '1/ISUM = ', 1/isum 
  print *, '1.D0/DSUM = ', 1.d0/dsum 
  stop 
end

program prog5 
  implicit none 
  integer, parameter:: IMAX = 10 
  integer:: i 
  real(8):: r 
  do i = 1, IMAX 
    call random_number(r) 
    print *, i, r 
  end do 
  stop 
end

program prog6 
  implicit none 
  integer:: N, M       !!!  M <- N - 1 
  integer:: i          !!!  Counter 
  real(8), dimension(:), allocatable:: r 
  print '(A$)', 'N = ' 
  read (*, *) N 
  print *, 'N = ', N 
  M = N - 1 
  allocate (r(0:M)) 
  call random_seed() 
  call random_number(r) 
  do i = 0, M 
    print *, i, r(i) 
  end do 
  deallocate (r) 
  stop 
end

program prog7 
  implicit none 
  integer:: N, M 
  integer:: i, idummy 
  real(8):: dmini, trial 
  real(8), dimension(:), allocatable:: r 
  read (*, *) N 
  write (*, *) 'Number of Data = ', N 
  M = N - 1 
  allocate (r(0:M)) 
  do i = 0, M 
    read (*, *) idummy, r(i) 
  end do 
! 
  dmini = r(0) 
  do i = 1, M 
    trial = r(i) 
    if (trial <= dmini) then 
      dmini = trial 
      print '(A10$)', 'CHANGED: ' 
      print '(I8, 2E15.7)', i, trial, dmini 
    else 
      print '(A10, I8, 2E15.7)', '', i, trial, dmini 
    end if 
  end do 
  deallocate (r) 
  print *, 'Minimum Value = ', dmini 
  stop 
end

program prog8 
  implicit none 
  integer, parameter:: N = 100, M = N - 1 
  real(8):: a, b, h 
  real(8), dimension(0:M):: x, y, z 
  real(8):: getvalue 
  integer:: i 
  a = 0.d0 
  b = 1.d0 
  h = (b - a)/dble(N) 
  do i = 0, M 
    x(i) = a + h*dble(i) 
    call setvalue(x(i), y(i)) 
    z(i) = getvalue(x(i)) 
  end do 
  do i = 0, M 
    print '(3E15.7)', x(i), y(i), z(i) 
  end do 
  stop 
end
 
subroutine setvalue(xval, yval) 
  implicit none 
  real(8):: xval, yval 
  yval = sin(xval) 
  return 
end
 
function getvalue(xval) result(yval) 
  implicit none 
  real(8):: xval, yval 
  yval = exp(-xval) 
  return 
end