名古屋大学 工学部 電気電子情報工学科 2年

確率論・数値解析及び演習 (数値解析演習)

演習担当: 今田 晋亮

<配布資料>


<Scilabのページ>


<Scilab サンプル関数・スクリプト等>

1) まずは,各自のホームディレクトリの下に「n-ana」という名前のフォルダを作成して下さい.
2) ファイルをダウンロードするにはファイル名を左クリックし,上記のフォルダに保存して下さい.(あるいはファイル名を右クリックしプルダウンメニューから「名前を付けてリンク先を保存」を選んで保存して下さい.)
3) 保存したファイルを編集する場合は,ScilabのアプリケーションメニューをクリックしてSciNotesを立ち上げ,FileメニューからOpenでそのファイルを開いて下さい.

● 第2章 …「計算機上での数値表現」
 (補足資料) [n-ana_chap2.ppt] … 説明に使ったパワーポイントのファイル

● 第3章 …「加速法」
 (補足資料) [n-ana_chap3.ppt] … 説明に使ったパワーポイントのファイル
 (関数)    [pi4.sci] … テキスト式(1.8)を計算する (スクリプトAitken.sciの中で呼び出される)
 (スクリプト) [Aitken.sci] … エイトケン加速でπ/4を求める

● 第4章 …「連立一次方程式」
 (関数)    [Gauss_pivot.sci] … ピボット付きガウスの消去法で連立一次方程式を解く
 (関数)    [Seidel.sci] … ガウス・ザイデル法で連立一次方程式を解く
 (スクリプト) [ex5_Seidel.sci] … 関数Seidelを用いて例題5(p.109)の方程式を解く

● 第5章 …「関数近似」
 (補足資料) [n-ana_chap5.ppt] … 説明に使ったパワーポイントのファイル
 (関数)    [Lagrange.sci] … ラグランジュ補間による多項式を求める
 (スクリプト) [ex1_Lagrange.sci] … 関数Lagrangeを用いて例題1(p.20)を解いてグラフ表示する
 (スクリプト) [ex2_Lagrange.sci] … 関数Lagrangeを用いて例題2(p.24)を解く

● 第6章 …「数値積分」
 (補足資料) [n-ana_chap6.ppt] … 説明に使ったパワーポイントのファイル
 (関数)    [trapezoid.sci] … 台形公式により数値積分の値を求める
 (スクリプト) [eq36_trapezoid.sci] … 関数trapezoidを用いて式(3.6)の数値積分の値を求める
 (関数)    [GaussLegendre.sci] … ガウスの求積法により数値積分の値を求める
 (スクリプト) [tab34_GaussLegendre.sci] … 関数GaussLegendreを用いて表3.4のガウスの求積法による積分値を求める

● 第7章 …「非線形方程式」
 (補足資料) [n-ana_chap7.ppt] … 説明に使ったパワーポイントのファイル
 (関数)    [Newton.sci] … ニュートン法により非線形方程式の解を求める
 (スクリプト) [ex2_Newton.sci] … Newton関数を用いて例題2(p.72)の方程式を解く

● 第8章 …「微分方程式」
 (補足資料) [n-ana_chap8.ppt] … 説明に使ったパワーポイントのファイル
 (関数)    [RungeKutta.sci] … (4次の)ルンゲ・クッタ法で微分方程式を解く
 (スクリプト) [ex3_RungeKutta.sci] … 関数RungeKuttaを用いて例題3(p.141)を解く
 (関数)    [f.sci] … シューティング法を適用する境界値問題の例題の微分方程式を定義する
 (スクリプト) [shooting.sci] … シューティング法により境界値問題の例題を解く