数値解析法, 2002 年, Fortran 入門編, 課題提出状況
課題:
Fortran 文法編の 4 回終了後一週間以内に, 各自作成したプログラムを
rmiya@stu.hirosaki-u.ac.jp 宛にメールすること。
(ソースファイルと何を計算するプログラムなのかの
簡単な説明文をつける)
作成するプログラムの内容は, 物理や化学の教科書に載っている法則をあらわす式に
基づいてグラフを書くための数値の表を出力するプログラム。
例えば ボイル-シャルルの法則 のグラフ, 調和振動子の波動関数,
運動する質点の軌跡など。他人と違うものが望ましい。
各行は, 学籍番号-プログラム名・メール受付日時・内容, 角括弧の中にコメント。
質問などがあれば, いつでもどうぞ。
最初のグループは, いちおうできているっぽいもの。
- h00s2003-yogenn.f, 5 Jul 2002 14:59:48, 余弦定理 [説明がワケワカランね。普通は余弦定理っていえば三角形の内角や辺の長さの関係 a**2 = b**2 + c**2 -2*b*c*cos(A) ってヤツじゃないか ? でもこのプログラムでやってることは違うことだ。そして末尾の FORMAT 文は何の役に立ってるのだろうか ?]
- h00s2014-SolventExtraction.f, 5 Jul 2002 13:38:13, 抽出率と pH との関係 [グラフまで作図して送ってくれたのは, 努力賞です。(でも Excel + Word なのがちょっと ... ;-p 独り言でした。) プログラムも良くまとまっています。でも e0 の扱いが不規則なのが少し残念ですね。]
- h00s2018-tikara.f, 19 Jul 2002 16:40:35, 力の公式 [訂正板とは訂正版のことでしょうか]
- h00s2028-kenji4.f, 5 Jul 2002 15:59:40, sherby-dornが定義した定常クリープの活性化エネルギーをもとめる [``Imput'' って何だ ? とか 状態 2 の絶対温度の入力を求めているようでいて実はそうでないとか。末尾の FORMAT 文は何してるのかな ? とか]
- h00s2045-suihei2.f, 5 Jul 2002 17:52:47, 物体の水平投射 [まあまあでしょう。このくらい簡単なら, あえて関数を使う必要はないけどね。]
- h00s2052-shuki2.f, 5 Jul 2002 16:59:46, 単振動の周期をもとめる [自分でも見てわかると思うけど, 全ての行に逐一コメントを付けても必ずしも読みやすいプログラムになるとは限りませんね。また DO 文で ``1, 100, 10'' としても ``100 まで'' にはなりませんね。]
- h00s2053-hukku.f, 5 Jul 2002 16:36:35, 変位と力の関係 [配列宣言に DIMENSION 文はなるべく使わない。というか必要ないはす。字下げや空白行の取り方に論理性が見られない。]
- h00s2057-undou.f, 5 Jul 2002 16:33:16, 運動方程式を表わす [説明不十分。配列宣言に DIMENSION 文はなるべく使わない。というか必要ないはす。]
- h00s2058-undou_enerugii.f, 19 Jul 2002 17:45:04, 運動エネルギーのプログラム [計算結果の数値と表の見出しがあっていない。また速度や質量といった物理量が整数だってのもアレだね。]
- h00s2060-kyuukannsuu.f, 12 Jul 2002 19:03:24, 球面調和関数 [何種類かの関数を選んでそれぞれ計算できるところは凝ってるね。でも同じ命令文をあちこちで繰り返し書くのは無駄だし, 間違いの元。]
- h00s2061-gurahu.f, 5 Jul 2002 17:54:41, 空気抵抗があるときの落下する速度 [空気抵抗の値が変な場合に再入力するようになってるのは良いね。``reration'' は愛敬か ?]
- h00s2062-lorentz_distribution.f, 5 Jul 2002 15:28:51, ローレンツ型の分布関数 [何度も同じものや改訂版を送ってきているが, 最後のもの (発信時間はこの項の左記) を有効とします。グラフ化したものも送ってきたのは努力賞 (でも Excel なのが ... ;-p)。でもこのプログラムじゃあ同梱の Excel によるグラフとおんなじ物は描けないよ。また b って何やってるの ?]
- h00s2063-fukku.f, 19 Jul 2002 16:02:09, フックの法則 [表を出力するので課題の要件は満たしているが, 物理的には変な感じだ。]
- h00s2064-nk8.f, 5 Jul 2002 14:52:40, 平均の空孔の濃度を求める [あいかわらず説明が要領を得ない。ここで言う空孔って何だ ? エネルギーの単位もワカランし。]
- h00s2068-ichi.f, 19 Jul 2002 17:41:59, ichi [説明が貧弱ですね。また計算結果も微妙に間違ってるはず。]
- h00s2070-nkg.f, 12 Jul 2002 17:30:29, 落下速度 [dimension 文は使わなくても済むはずですね。]
- h00s2072-shahou.f, 5 Jul 2002 17:23:07, 物体をシャホウ投射した時の式 [この式は斜方投射に含まれると思うが, 逆にこの式を出されても一般には斜方投射とは考えないと思うぞ。]
- h00s2079-mole-furaction25.f, 12 Jul 2002 18:54:57, 酢酸イオンのモル分率 [数値の表にしては余計な物も少し出力される。また, センターのマシンでは大丈夫そうに見えるけど, 実数の累乗 `H = H * 0.1' は計算誤差の元。]
- h00s2082-riso.f, 19 Jul 2002 15:34:21, risoukitai no jyoukiatu [改訂版を送ったつもりなのかな ? 前とおんなじ物だよ。ということで, 数値の表の前に出力される部分でおかしなことしてる。]
- h00s2083-houkai.f, 3 Jul 2002 17:15:45, ウランが鉛に崩壊するときの50億年ごとのウランの原子数 [変数 d, e の意味は ?]
- h00s2084-reynold.f, 3 Jul 2002 20:48:08, ある密度・粘度の液体が、ある円管内をある流速で流れるときのレイノルズ数 [計算して表示される数値の桁がそろっていればなお良い]
- h00s5016-dannsei-b.f, 9 Jul 2002 13:04:45, バネ定数についてのプログラム [数値と文字列との区切りに必ずしも空白を出力しない処理系があるかもしれないので, 明示的に空白を出力するようにした方が安全ですね。こんな感じ: i, " dannseiryoku"。まだ少し説明文が変だ。]
- h00s5018-houkai-b.f, 9 Jul 2002 12:45:04, 原子の崩壊の公式 [数値と文字列との区切りに必ずしも空白を出力しない処理系があるかもしれないので, 明示的に空白を出力するようにした方が安全ですね。こんな感じ: a, " gennshisuu"]
- h00s5028-kadai.f, 28 Jun 2002 18:07:57, 斜方投射の座標算出 [角度がラジアン単位なのが馴染みにくい]
- h00s5051-kadaiss.f, 5 Jul 2002 17:39:45, (以前の説明によれば) 鉛直投げ上げの式 [直したのはいいんだけど, これだけで説明が完結するようにしてよ。さて WHILE ループ中の IF 文によるループからの脱出は必要なのでしょうか ?]
二番目以降のグループは, やや問題点あり。
- h99s2076-undou.f, 4 Jul 2002 15:52:38, ニュウートンの運動方程式 [``ニュウートン'' って何だろう ? また質量の単位とおぼしき ``Kg'' も変 (ケルビン・グラムか ?)。いずれにしても計算はするけど出力は表じゃない。]
- h00s2021-enthalpy.f, 4 Jul 2002 19:07:15, enthalpy [説明がない。計算はするけど出力は表じゃない。]
- h00s2022-noudo.f, 5 Jul 2002 18:12:44, 酸・塩基反応の平衡定数 [いったいどんな反応なのか ? a, c は何なのか ? 説明がないと全くわからない。そして計算はするけど出力は表じゃない。]
- h00s2027-chouwasindousi.f, 5 Jul 2002 18:01:10, 一次元調和振動子のエネルギー [計算はするけど出力は表じゃない。]
- h00s2059-oo.f, 4 Jul 2002 18:11:28, oomu no housoku [説明がない。計算はするけど出力は表じゃない。全部送ればいいってもんじゃない, いらないものは送ってよこすな。何が必要なモノかはよーく考えよう。]
- h00s2075-kadai.f, 4 Jul 2002 21:55:34, 分子の熱運動速度 [説明に書いてある式は正しいのですか ? また計算をして結果は出力するが, 表じゃない。]
その他にトラブル (?) を抱えているもの。
- h00s2008-omu.f, 5 Jul 2002 17:54:16, オームの法則 [プログラムの論理構造が直っていません。一見もっともらしく動作することもあるが, 良く考えるとプログラムの論理がオカシイ。提出されたプログラムの DO 文を論理的に等価な WHILE ループに直してみるとはっきりわかる。]
- h00s2009-bannyuu.f, 4 Jul 2002 20:33:25, 万有引力 [一見もっともらしく動作することもあるが, 良く考えるとプログラムの論理がオカシイ。提出されたプログラムの DO 文を論理的に等価な WHILE ループに直してみるとはっきりわかる。]
- h00s2015-un.f, 5 Jul 2002 18:20:58, 運動の法則 f=ma について [一見もっともらしく動作することもあるが, 良く考えるとプログラムの論理がオカシイ。提出されたプログラムの DO 文を論理的に等価な WHILE ループに直してみるとはっきりわかる。]
- h00s2030-shindoushi.f, 5 Jul 2002 18:42:09, 一次元調和振動子の運動方程式をEuler法で計算する [難しいことをやろうとしているね。でもこのプログラムを自分で動かしてみましたか ? ちゃんと期待通りに動いたのかなぁ ? まず致命的な欠陥を修正しても, 計算は正しいのでしょうか ? 説明によるとステップ毎に「時刻・座標・運動量・エネルギー」を出力することになってるけど, 時刻は一定だし, 座標や運動量は振幅が大きくなりつつ振動してるし, エネルギーも単調増加してるけど, いいの ? 何か多重に欠点を抱えているような気がする。]
- h00s2032-banyuuinnryoku.f, 5 Jul 2002 15:46:19, 距離に対する万有引力の値 [ええっと, センターの f90 ではエラーを吐かずにコンパイル・実行できてしまうが, 実はわりと重大なミスがある。この為に実はセンターのマシンでも正しくない答えがでてきたはずだ。わかるかな ? FORTRAN の約束を見直そう。]
- h00s2050-jiyuurakka.f, 4 Jul 2002 20:36:07, 自由落下 [一見もっともらしく動作することもあるが, 良く考えるとプログラムの論理がオカシイ。提出されたプログラムの DO 文を論理的に等価な WHILE ループに直してみるとはっきりわかる。]
(2002-08-03, 24:10 更新)