第2回レポート問題

2016 年度・前期・数理解析・計算機数学3
1
● 第2回レポート問題
問題 B-1 は必ず解答すること. それ以外の問題から1問を選び解答すること.
★ 問題
すべての問題で, 適切な例に対するプログラムを書くことと, 適切な図を作成すること. すべて
の計算(議論)は倍精度浮動小数点数で計算すること.
B-1 (25 点) 常微分方程式(単振り子の方程式)の初期値問題
x′ (t) = − sin(x(t)),
x(0) = 1.5,
x′ (0) = 0.0,
を h = 0.01 とおいて, Runge-Kutta 法, Symplectic Euler 法, 陰的中点法を用いて数値解を
計算するプログラムを書きなさい.
B-2 (20 点) 単振動の方程式
x′ (t) = −x(t),
x(0) = a,
x′ (0) = b
を陰的中点法で解いた数値解が, 相空間上で束縛される曲線を求めなさい.
B-3 (15 点) フィボナッチ数列
xn+1 = xn + xn−1 ,
x0 = 1,
x1 = α ∈ R
を計算すると, α の値によっては, 期待できる結果と異なる挙動を示すことがある. そのよう
な挙動を示す α の値を示し, その理由を説明しなさい.
★ 締め切りと提出方法
締め切りは2016年7月20日(水曜日)17時とする.
• メールの送付先は [email protected] し, Subject には, 文字列「レポート提出」
を含むこと. また, Subject には問題番号を書くこと.
• 電子メールの発信者は, “nagoya-u.ac.jp” のアドレスであること.
• プログラム, PDF ファイル, 図版の電子データはメールで提出すること.
• 図版および数学的な内容を記述したレポートは, PDF ファイルにまとめて提出することが望
ましいが, それができない場合には紙に書いて(図版はプリンタで出力して)提出してもか
まわない.
• プログラム中には「問題番号, 学籍番号, 氏名」をコメントとして記入すること. PDF ファ
イルまたは紙で提出するレポートには, 「問題番号, 学籍番号, 氏名」を記入すること.
• PDF ファイルの用紙サイズ, および, 紙で提出する用紙サイズは A4 とする. (特に理由が
ない限り, A4 縦おきであること.)紙で提出する場合には, 片面であること.
• 図版を単独の電子データで送付する際のファイルフォーマットは, PDF, EPS, SVG, JPEG,
PNG のいずれかとする.
Jun. 29, 2016, Version: 1.0
[email protected]
2016 年度・前期・数理解析・計算機数学3
2
★ 採点基準
当初に説明した通り, 以下の基準で採点するが, 中間的な点をつけることもある. また, 非常に良
いレポートにはボーナス点を与える. (A:100%, B:60%, C:0% と考えて良い.)なお, 今回の満点
は 40 点とし, それを越える点はボーナス点とする.
A 数値計算が正しいプログラムによって行われ, その考察が適切であるもの.
B 数値計算は正しいプログラムによって行われているが, その考察に少々問題があるもの.
または, 数学的背景は正しく記述されているが, プログラムに小さな誤りがあるもの.
C それ以外.
★ 注意事項
以下の注意事項に沿わないレポートは採点から除外する.
• 提出方法に従わないレポートは採点の対象とはしない. 特に, 電子ファイルのファイルフォー
マットに注意すること.
• 提出するレポートおよびそのプログラムは, 自分で作成したものに限る. 特に, 他人書いても
らったプログラムや, 「インターネット」で拾ってきたプログラムをレポートとして提出し
てはならない.
このような場合には, 採点から除外するだけでなく, 不正行為として直ちに不可と判定する.
• プログラムは, 浮動小数点演算が IEEE の規定に沿う言語で記述すること. (C なら問題な
い.)また, プログラム中には, 何をしているプログラムか, 各部分が何を実行しているのかを
コメントとして記入すること. コメントが全くないプログラムは採点対象としない. コメン
トが不適切であったり, コメントが無いことにより何をしているのかが読み取れないプログ
ラムは減点対象とする.
★ 一言
レポートは, できたからと言って, 「よしできたぁ!」とすぐに出すのはやめましょう. 数日後に
もう一度自分の書いたレポート(プログラム)を見直して, 自分で何を書いているかがわかるかを
考えてください.
以上
Jun. 29, 2016, Version: 1.0
[email protected]