単振動している質点に、速度に比例して働く抵抗 -2 m μ x'(t) が加わったとします。（ ここでμは正の定数。）

この場合の運動方程式を解きたい場合、DSolve[ ]を使って、

とすればいいですね。ところが得られた答えをよく見ると、これは ω < μの場合の一般解になっています。ここで、実際に手で方程式を解く場合は

	（）
	（）
	（）

といった具合に場合分けする必要があったことを思い出してください。Mathematica上で場合分けしてやることもできますが(興味のある人はAssumingやRefineというコマンドを調べてみましょう)、 今回は最初のレポートなので、もう少し簡単な問題です。

上で述べた減衰振動系で「ω＝μ」の場合ついて、以下の内容が書かれたMathematicaのfileを、PDF形式で manabaに提出して下さい。

1. DSolveを使って運動方程式を解いて解を得る（ヒント：ω＝μの場合を考えたいので、あらかじめμをωに置き換えてからDSolveをして構いません。） うまく計算ができていれば一般解は上のx2の形になるはずですが、各自で初期条件を設定し、それを満たすように積分定数も決めてください。　 　
2. 運動の時間プロット（縦軸が質点の座標、横軸が時間のプロット）を描く。パラメータの値はすべて各自で自由に設定してください。
3. 質点の運動エネルギーと位置エネルギー、及びその二つの和の時間プロットを描く。
4. 位相空間内の運動プロットを描く。

やり方はすべてテキストで説明されているものばかりですので、わからない場合はテキストをよく見直しましょう。 また、単に結果だけではなくどのような手順をおこなったのか、他の人にわかるように 説明を加えて下さい。 この点は評価の重要なポイントになります。Mathematica上で説明を書きたいときは例えば、「1+1 (* いまから1+1の計算をします*)」のように(*と*)で書きたいコメントを囲みましょう。囲まれた部分はコメントとみなされ、計算上は無視されます。 その他様々なオプション等をつけてわかり易く個性的なものにしてみて下さい。 正確さはもちろんですが、見た目も評価の対象になります。他の人と酷似したレポートは、不正と見なされる場合があるので注意して下さい。