先週(Lesson08 (その2))のメールに対する返事 (2001 年度版)

• Q: Newton 法を忘れてしまいました。どうすればいいのでしょう。
• A: Lesson07 のプログラムの部分 を復習しましょう。 今回の場合のヒントは、
  y (t+dt) = y(t) + v_y dt = 0
  となって欲しいので、これから dt を求めればよいです。

• Q: プログラムのエラーはなかったのに、a.outができません。 「can't create a.out」と出てしまうんですが・・。
• A: よくあるのが、「前に走らせたプログラムがまだ動いている(終っていない)」 ことです。Ctrl-Z でいったん停止にしてあっても、プログラムが終らないと a.out には書き込めません。また別の Window で書き込もうとしても 書けません。 ですから、
  • Ctrl-Z でいったん停止した場合は、 fg (foreground) で元にもどして、 Ctrl-C で完全にプログラムをとめて下さい。
  • 以前に走らせたプログラムがまだ動いているかどうかは、 process を調べるコマンド
    ps x
    で a.out というものがあるかどうか見れば分かります。

• Q: gnuplotを書くときにmuleなどを使うとかいてありましたが、などとは何ですか？
• A: Unix 上で編集するエディターには、mule の他に、
  • emacs (mule の英語版、実はこちらが本家)
  • vi (VIsual editor, 昔は visual だったのだが、今は mule の方が visual かな。私は好きです。)
  などがありますし、Windows で編集して ftp で送るという手もあるので、 「など」と書いたわけです。
• Q: ｗｒｉｔｅ文中の#と、プリント34のプログラム4行目のｄｂｌｅの 意味がわかりません。説明お願いします。
• A: write 文の中には # はないはずですが...
  • * であれば、「標準」を表します。write(*,*) の一つ目の * は 「プログラムの標準的出力先 = 画面」を表し、 二つ目の * は 「標準的な出力方法(実数は全ての桁を書く)」を表します。
  • gnuplot のスクリプト中の # は、 「この # 以降はコメントである(gnuplot の動作に影響しない)」 ことを表します。
  また、p 34 の dble(i) の dble は変数の型を変換する関数で、
  • x = dble(i) の場合、i は整数で x は倍精度実数
  となります。ただし、人間が見る限り値としては同じですし、
  • x = i
  とだけ書いても多くの場合自動的に整数が倍精度実数に変換されます。
• Q: prog082をやった。何回やってもgaがおかしいと言われた。 宣言文って何ですか？
• A: 宣言文とは、
  • implicit real*8 (a-h,o-z)
  • real*8 a, abc, def
  • integer i, j, lmn
  • external FUNC
  などのように、その文字列がどのような変数か、 あるいは関数かを示す文です。