● Lesson03 (数値積分 その2, subroutine 副プログラム, 課題 その1)

• 今週、最初にやること
• 先週の演習問題の解答例と気がついたこと
• 今週の PROGRAM
• 今日のポイント
• 余裕があれば、試してみて下さい。
• 課題 その１
• プログラムが出来たら
• 続きのページ

● 今週最初にやること

• いつものように、 出席メールの受け取り状況のページ を確認して下さい。

• また、これまでと同様に、 先週作ったプログラムを、一つの directory, lesson02 にまとめましょう。 やり方は ここ

• 今週は少し、Unix の環境整備を行いましょう。 やり方は ここ

• プログラムができたら、今週は皆さんの計算結果をメールしてもらいます。
  やり方は ●プログラムが出来たら の所に 書いてあります。

● 今週の PROGRAM

前回と同じ問題設定ですが、 subroutine を使ったプログラムを作ってみましょう。

• サブルーチンを用いたプログラム (prog03.f) (計算実行例)

上のプログラムでは、関数 f(x) しか積分出来ません。 いくつかの別の関数を同様のやり方で積分出来るようにしてみましょう。

• 積分する function の名前を変えられる subroutine を使ったプログラム (prog032.f) (計算実行例)

今週の最初へ戻る

● 今日のポイント

• subroutine 副プログラムの使い方
  まとまった一つの操作の組合せを行なうのが subroutine です。
  以下の例では、subroutine chuten(a,b,n,sum) は、 メインプログラム (program ... からはじまり、end で終っている 最初の部分) から呼び出されて、 与えられた a,b,n の値を使って 積分値 sum の値を返します。

  ---

  program prog03
  
  implicit real*8 (a-h,o-z)
  
  ...
  
  call chuten(a,b,n,sum)
  
  ...
  
  end
  c
  
  subroutine chuten(a,b,n,sum)
  
  implicit real*8 (a-h,o-z)
  
  ...
  
  sum = sum + f(x)
  
  ...
  
  end
  
  c
  
  function f(x)
  
  implicit real*8 (a-h,o-z)
  
  ...
  
  f = ...
  
  end

  ---

  do n=初期値, 終値, きざみ このプログラムでは、より一般的な DO ループの形式を使っています。

  DO K = 初期値, 終り値 (, きざみ)
  

  (演算の内容)
  

  END DO

  (例) k = 10 から 100 まで 10 ずつ増やしてて繰り返す場合。

  DO K = 10, 100, 10
  

  (演算の内容)
  

  END DO

  きざみを省略すると、1 と見なされます。(K が 1 ずつ増える)

  write(unit,format): ファイルへの出力
  例: write(16,800): ft.16 に文番号 800 に 示された format で出力

  800 format(1x,i3,5(1x,f10.7))

  • 文番号は 1 〜5 桁の間に。
  • n x は n 個の空白、 i n は n 桁の整数、
    5(...) は、(...) を 5 回繰り返し、 f10.7 は、実数で少数点以下 7 桁で出力

• function の使い方
• external 文

● 余裕があれば、試してみて下さい。

(参考文献: 森 正武、「FORTRAN77 数値計算プログラミング」 (岩波コンピュータサイエンス))

今週説明し切れないと思いますので、コピーして試して見て下さい。

mail Teacher@server s0300xxy < ex03x.f
(ex03x.fのx は、1 か 2 です。) として、プログラムを送って下さい。

今週の最初へ戻る

---

● プログラムが出来たら、

今週は prog03.f の計算結果をファイルに出力 してみましょう。

ap1 25: f90 prog03.f
f90: compile start : prog03.f

*OFORT90 V01-04-/A 開始
*プログラム名 = PROG03
*プログラム名 = CHUTEN
*プログラム名 = F
*プログラム数 = 0003 , エラーはありません。
ap1 26: a.out > prog03.dat
? 
10 20

"./a.out > ファイル名" とすることにより、(スクリーンで見えていた内容が)ファイルへ出力できます。
(すでにあるファイルには出力できません。一度消すか(rm ファイル名)、 別のファイルに出力して下さい。詳しくはここ を見て下さい。)

ここまで出来れば、 出席の返事の代わりに、