以前のリビジョンの文書です
問題
(2)はが大きくなると
を上手に選ばないと正しく計算できない。
#define N 10000 double dh=0.01;
のままで計算する場合は
#define DELTA 1.e-7
のDELTAを次の値に設定すると計算できる。
![]() | ![]() |
|
0 | 1.e-7 | 1.e-7 |
1 | 1.e-7 | 1.e-7 |
2 | 1.e-7 | 1.e-7 |
3 | 1.e-9 | 1.e-9 |
4 | 1.e-11 | 1.e-11 |
5 | 1.e-13 | 1.e-13 |
6 | 1.e-15 | 1.e-15 |
7 | 1.e-17 | 1.e-17 |
8 | 1.e-19 | 1.e-19 |
9 | 1.e-21 | 1.e-21 |
どちらの積分も厳密に計算できるので、誤差を見積もるのは容易なはずである。この期に及んで積分ができないとぬかすヤツは研究室まで出頭せよ。