2018年度現象数理研究I Q2

　金曜5限に407教室で行います。→ 3限は401, 5限は407となりました。

連絡事項

• オーダー決定:
  • 3限は 401 教室、 話す順番は海老原、波多野、寺田、小林、市川、です。
  • 5限は 407 教室、 話す順番は上原、冨岡、有薗、石井、です。

レポートについて

資料

• (参考) 「桂田研の卒業研究の記録」 書きかけですが。
• 菊地・山本「微分方程式と計算機演習」山海堂 (1991) の 前半部分のFORTRANプログラムをCに書き換えたもの。 kikuchi-yamamoto-progs.zip
  (ばら), (HTML), (PDF)
• 「C言語これくらい覚えよう」
• (参考) 「レポートの書き方」 (PDF), (HTML)
• (参考) 「輪講のやり方」 (PDF) , (HTML) -- 数学科のゼミなどでは、 全員で共通のテキストを順番に説明する「輪講」という勉強の仕方が普通です。 それをやるときの手引き。
• (参考) 「卒業研究レポートの書き方について」 (PDF), (HTML)
• 「gnuplot入門」 --- 誤差のデータ等は数表よりはグラフにする方が良い。 この目的には Excel は勧めない。 gnuplot を使えるようになることを勧める。
• 連立１次方程式・固有値問題について
  「桂田研卒研ノート」 というページに色々あります。
  • 「連立1次方程式I ---計算量と直接法---」
    LU分解については 「『発展系の数値解析』に加えること」 の2節の方がまとまっているかも。
  • 「行列の固有値問題」
  • 「固有値問題ノートの補足」
• TeX を使うならば、MacTeX の新調を考えてみよう。 「MacTeX 2018」 mactex20180417.pkg はネットで入手出来ますが、 大きいので USB メモリで渡しても良いです。 なお、TeX の使い方の説明は 「TeX入門 2018」 などを見て下さい。
• GLSC がきちんと使えるかチェックしておいて下さい。 heat1d-i-glsc.c を保存してコンパイル。
  ターミナルで次のコマンドを実行します。

         curl -O http://nalab.mind.meiji.ac.jp/~mk/program/fdm/heat1d-i-glsc.c
         cglsc heat1d-i-glsc.c
         

  これで heat1d-i-glsc という実行形式のプログラムが出来るかどうか。 以下は実行例です。

         [chronos:~] mk% ./heat1d-i-glsc
         入力して下さい : nfunc(1..5)=1
         入力して下さい : θ=0.5
         入力して下さい : N=100
         入力して下さい : λ=0.5
         時間の刻み幅τ= 5e-05 になりました。
         入力して下さい : 最終時刻 Tmax=1
         入力して下さい : グラフ書き換え時間間隔(Δt)=0.01
         T=   0.0000e+00
         I     u(i)      I     u(i)      I     u(i)      I     u(i)      I     u(i)
         0  0.0000e+00   1  1.0000e-02   2  2.0000e-02   3  3.0000e-02   4  4.0000e-02
         5  5.0000e-02   6  6.0000e-02   7  7.0000e-02   8  8.0000e-02   9  9.0000e-02
         10  1.0000e-01  11  1.1000e-01  12  1.2000e-01  13  1.3000e-01  14  1.4000e-01
         15  1.5000e-01  16  1.6000e-01  17  1.7000e-01  18  1.8000e-01  19  1.9000e-01
         20  2.0000e-01  21  2.1000e-01  22  2.2000e-01  23  2.3000e-01  24  2.4000e-01
         25  2.5000e-01  26  2.6000e-01  27  2.7000e-01  28  2.8000e-01  29  2.9000e-01
         30  3.0000e-01  31  3.1000e-01  32  3.2000e-01  33  3.3000e-01  34  3.4000e-01
         35  3.5000e-01  36  3.6000e-01  37  3.7000e-01  38  3.8000e-01  39  3.9000e-01
         40  4.0000e-01  41  4.1000e-01  42  4.2000e-01  43  4.3000e-01  44  4.4000e-01
         45  4.5000e-01  46  4.6000e-01  47  4.7000e-01  48  4.8000e-01  49  4.9000e-01
         50  5.0000e-01  51  4.9000e-01  52  4.8000e-01  53  4.7000e-01  54  4.6000e-01
         55  4.5000e-01  56  4.4000e-01  57  4.3000e-01  58  4.2000e-01  59  4.1000e-01
         60  4.0000e-01  61  3.9000e-01  62  3.8000e-01  63  3.7000e-01  64  3.6000e-01
         65  3.5000e-01  66  3.4000e-01  67  3.3000e-01  68  3.2000e-01  69  3.1000e-01
         70  3.0000e-01  71  2.9000e-01  72  2.8000e-01  73  2.7000e-01  74  2.6000e-01
         75  2.5000e-01  76  2.4000e-01  77  2.3000e-01  78  2.2000e-01  79  2.1000e-01
         80  2.0000e-01  81  1.9000e-01  82  1.8000e-01  83  1.7000e-01  84  1.6000e-01
         85  1.5000e-01  86  1.4000e-01  87  1.3000e-01  88  1.2000e-01  89  1.1000e-01
         90  1.0000e-01  91  9.0000e-02  92  8.0000e-02  93  7.0000e-02  94  6.0000e-02
         95  5.0000e-02  96  4.0000e-02  97  3.0000e-02  98  2.0000e-02  99  1.0000e-02
        100  0.0000e+00
        マウスでウィンドウをクリックして下さい。
        [chronos:~] mk%
  

  

参考資料