|
|
| 2009年05月13日 10時34分 |
| GNUPLOT |
 |
 |
ルジャンドル多項式をGNUPLOTを使ってグラフにする。
ルジャンドル関数を再帰的なプログラムによって、計算しデータをつくる。(leg.c)
これをコンパイル・実行し、data1というfilenameに格納する。
userid% gcc -ml leg.c
userid% ./a.out > data1
次に、GNUPLOTでそれをプロットする。
userid% gnuplot
gnuplot> plot "data1" w l
gnuplot> set term post
gnuplot> set out "data1.ps"
gnuplot> replot
gnuplot> exit
data1.ps というポストスクリプトができた。PDFに変換するには、
userid% ps2pdf data1.ps data1.pdf
と入力すればよい。また、
userid% lpr data1.ps
でポストスクリプトのままでも印刷できる。 |
 |
 |
| カテゴリー:日記 |
 |
ルジャンドル多項式のゼロ点(Pn(x)=0を満たすx)を求めたい。
① n=4 の時 (解析解あり。)
② n=10の時 (解析解なし。)
の2つの場合についてニュートン法または2分法を用いて求めよ。 締め切りは5月27日までとします。
|
 |
 |
 |
| カテゴリー:日記 |
|
|