レポート課題（2006年5月31日）

次の強制振動に関係する２階常微分方程式に関して，

（A)をオイラー法またはルンゲ・クッタ法で解くプログラムを作成せよ．パラメータおよび初期値は，

とし，A=0.1 とする．時間刻み幅 h は h=0.1 とする（オイラー法の場合にはh=0.1ではうまくいかないであろう．その場合には適切なhの値とすること）．解 y(t) を赤色，外力 Acosω_et を緑色で表示することにする．なお，必ずグラフ中に各人の学籍番号を描き，プログラム中にはコメント分としてプログラムの説明を書き入れる事．

ファイル名は rep0531.c とし，report フォルダにおく事．

実行結果はおおよそ以下のようになる．

戻る