２００８年度・数理解析・計算機数学２
1
２００８年度・数理解析・計算機数学２
レポート問題
２００８年７月２日
【評価方法について】
授業の評価はここにあるレポートの結果のみを元にして行う.
【レポート問題について】
★ 全てに共通な注意
他人のレポートやプログラムのコピー,またはそのマイナーな改変を提出してはならない. そのようなこ とがあった場合は
,
提出者だけではなく,
コピーをさせた人についても,
すべてのレポートを評価対象とし ない.また,各種の文献（ネットワーク上のリソースも含む）を参照した場合には,参考文献を必ず明記するこ と. 明らかに文献を参照しているにも関わらず,出典を明示していないレポートは評価対象としない. （た だし,講義で紹介した事実は参考文献を明示する必要はない.）
★ レポートの提出方法
提出締切 ２００８年８月８日（金）１７時（時間厳守）
電子メールで提出するレポートに関しての締切の意味は,提出期限の以前の時刻でメールが発信され ていればよい
.
（時刻は日本標準時とする.
）解答方法 以下の問題の中から各カテゴリ（カテゴリＡ〜Ｃ）各１問づつを解答とすること
.
各カテゴリに ついて, 複数問解答してもかまわないが,採点の対象は「最も良い内容のもの」１問のみとする.提出内容
•
プログラムを書く必要がある問題については,
プログラムを提出すること.
•
数値計算結果は結果を読み取ることができる図を作成すること. 生の数値データは必要ない.提出方法
•
プログラムは電子メールで「添付ファイル」として,ソースコードを送付すること.•
解答内容は.
（,
「レポート用紙」に手 書きまたはタイプセットして提出しても構わない.•
電子メールは“Subject”
に「レポート提出」と明示し,添付するファイル名はreport.tex,v 1.6 2008-07-02 08:08:17+09 naito Exp
2
２００８年度・数理解析・計算機数学２問題番号-学籍番号.拡張子
の形式であること. 具体例としては, 問題
Prob-A-01
の“Prob-A-01-XXXXX.pdf”
とすること.
•
電子メールの発信者は,学部生の場合には「情報メディア教育センター」のアドレス,大学院生 の場合には「多元数理科学研究科」のアドレスであること. 他のドメインから発信された電子 メールで提出されたレポートは採点対象としない.
•
プログラム以外の部分を「紙に書いて」で提出する場合には,
Ａ４用紙を片面で利用し,
クリッ プで止めること. （ステプラは不可）また,全てのページの右上に問題番号・学籍番号・ページ 番号を明示すること.提出先 電子メールの場合には内藤
(naito@math.nagoya-u.ac.jp)
宛,紙で提出する場合には,多元数理 科学研究科１Ｆ事務室のレポート提出箱.その他 締切日時よりもある程度早く提出されたレポートに対してはコメントを返送する. （つまり,「ダ メ」とか「まあＯＫ」とかの返信を電子メールで行なう.）携帯電話のメールアドレスまたは個人で 利用しているアドレスにコメントを返送してほしい場合には,最終回の講義の時間までに口頭で内藤 まで連絡すること.
採点基準 以下の採点基準で採点を行なう.
•
カテゴリＡ・Ｂの問題はA, B, F
の３段階で採点する.A
数学的な議論が正しく,考察内容も優れている. また, プログラムが必要な場合には,プログ ラムも正しく記述されている.B
数学的な議論・考察内容・プログラムは必ずしも間違っていないと判断できるもの.F
数学的な議論・考察内容・プログラムのいずれかに重大な誤りがあるもの.•
カテゴリＣの問題はA, B, F
の３段階で採点する.A
プログラムが正しく記述され,結果も適切に図示されているものB
プログラムまたは結果の図示に問題があると判断されるが,
誤りであるとは言えないもの. F
プログラムまたは結果の図示に誤りがあるもの.各問題の配点は次の通り.
問題
A B F
Prob-A-01 5 3 0 Prob-B-01 5 3 0 Prob-B-02 7 4 0 Prob-B-03 ∗ † 0
問題
A B F
Prob-C-01 5 3 0 Prob-C-02 ∗ † 0
この表の中で
∗
とあるものは10
点を上限, 5点を下限にして, また†
とあるものは6
点 を上限, 3 点を下限にして,問題設定に応じて採点する.成績判定基準 以上の採点を行なった後に
,
以下の基準で成績を判定する.
（基本的には15
点満点と考え ている.）report.tex,v 1.6 2008-07-02 08:08:17+09 naito Exp
２００８年度・数理解析・計算機数学２
3
優
12
点以上.良
10
点以上.可
7
点以上.ただし, いずれの場合も,カテゴリＡとカテゴリＢの点の和が０点でないこと.
★ プログラムの仕様について
数値計算を行なうプログラムの言語は問わない. ただし,以下の注意を守ること.
•
プラットフォーム依存の記述は避けること.•
利用する言語における標準的な浮動小数点フォーマットを利用すること.
•
情報メディア教育センターのLinux
システムで利用可能な言語であること.• C, Java
以外の言語を用いたときには, どの言語を利用したか,コンパイル方法及び実行方法をプログラム中にコメントとして記述すること.
•
いずれの言語であっても,プログラム冒頭部分に, 問題番号・学籍番号・氏名をコメントとして挿入 すること.
•
添付ファイルとして添付する場合に,文字コード変換なしにコンパイラに通る文字コードを利用して 送付すること.•
微分方程式の数値解を求める問題では, 差分の刻み幅等は,解の様子がよく分かる程度に適切に取れ ばよい.★ レポート問題の訂正・質問など
もし,レポート問題に不備が発見されたときには,速やかに内藤まで電子メールで連絡すること. それら の訂正はウェブページで,訂正内容を掲示する.
夏休み期間中に質問がある場合には
,
事前に電子メールでアポイントメントを取ること.
report.tex,v 1.6 2008-07-02 08:08:17+09 naito Exp
4
２００８年度・数理解析・計算機数学２【レポート問題】
• Prob-A-01
６回目までの講義で行った,各種の数値の計算に関して,いずれか一つ以上に関して,理論的な背景を 議論し
,
その数値計算を行いなさい.
また,
得られた数値実験結果を理論的な背景に即して考察を行 いなさい.ここで「各種の計算」とは以下のいずれかである. 「an
= (1 + 1/n)
n の計算」,「ニュートン法」,「方程式の反復解法」,「πの計算」,「巾級数による関数値の計算」, 「収束の加速」.
• Prob-B-01
次の常微分方程式の初期値問題を,講義で紹介した数値計算スキームのいずれかを用いて計算し, 時 間刻み幅
h
を動かしたときのt = 1.0
における値の誤差について議論しなさい.x
(t) = x(t), t ∈ [0, 1], x(t) ∈ R, x(0) = 1.0.
• Prob-B-02
「数値的不安定性」とは何かを修正オイラー法を例にとって
,
数値実験を交えながら議論しなさい.
• Prob-B-03
その他, 常微分方程式または偏微分方程式の数値計算に関連する理論的考察を,必要なら数値実験を 交えて行ないなさい.
• Prob-C-01
単振り子の運動とは, ハミルトニアン
H (q(t), p(t)) = 1
2 p(t)
2− cos(q(t)), q(t), p(t) ∈ R
で表される運動である.
この運動を,初期条件
q(t) = 3π/4, p(t) = 0.0
のもとに,「４次のルンゲクッタ法」, 「（陽的２段以上, または陰的）アダムズ法」,「なんらかの シンプレクティック法」のいずれかで計算し, (q, p)平面内で運動を図示しなさい.
• Prob-C-02
その他
,
２階常微分方程式の初期値問題を,
「２次以上のルンゲクッタ法」,
「（陽的２段以上,
また は陰的）線形多段階法」,「なんらかのシンプレクティック法」のいずれかで計算しなさい.【以上】
