数値計算課題に関する注意

鳥取大学工学部 応用数理工学科 学部生向け資料
[数値計算のページに戻る]



数値計算をはじめる前に

このプログラムでは、どういう方程式を解いているのか?
プログラムを書いている最中はともかく、 何ヶ月かたつと何の式か完全に忘れてしまい、 「あれ、このプログラムって何だっけ?」となることが多い。 予防策として、方程式を、プログラム中に コメント文の形で書いておくべきである。 (複雑な式の場合は、TeXファイルを作っておき、 プログラム中に「***.tex を参照」と書いておく。) 方程式についてのコメントは、変数の説明も兼ねるようにする。
数値計算しやすいように方程式を書き直す
無次元化の際に、 基準量を適切に選ぶようにする。

プログラムの書き方

課題プログラムの実例と講評

まちがいプログラム! (2000年)
他人が作ったプログラムを、内容を理解せずに信用してはいけない。
2002年10月の課題: コメント
実例(1-5)とコメントを比較してみよう。

悪いプログラム

いかにもバグのありそうなプログラム

読みにくいプログラムは、バグの温床です。

FORTRAN77の危険なところ

暗黙の型宣言 (implicit なんとか)
必ず implicit none を使うように。 「implicit double precision」などというのは過去の遺物。
大文字と小文字は区別されない
もし u を使ったら、 別の変数として U を使うことはできない。 なぜなら、uU は FORTRAN77にとっては同じものだからだ。
行番号はバグのもと
do ... enddo などの 制御構造を活用して、 行番号なしのプログラムを書くようにする。

結果の検討・レポートの書き方など


[数値計算のページに戻る]