[정적분의 근사계산] Simpson의 공식
함수 \(f(x)\)가 구간 \([a,\,b]\)에서 연속 이고, \(f(x)\ge0\)이라 한다. 구간 \([a,\,b]\)를 \(2n\) 등분해서 분점을 \[a=x_0<x_1<x_2<\cdots<x_{2n}=b\] 라 하자. 각 분점에 있어서의 함수 \(f(x)\)의 값을 각각 \[y_0=f(x_0),\,y_1=f(x_1),\,\cdots,\,y_{2n}=f_{2n}(x_{2n})\] 라 한다. 각 소구간의 길이는 모두 같아서 \[h=\frac{b-a}{2n}\] 이다. 사다리꼴의 공식 에서는 곡선의 각 소부분을 선분으로 보았으나, Simpson의 공식에서는 이것을 \(y\)축에 평행한 축을 가지는 포물선의 일부로 보고 적분 의 근사값을 구한다. 그래서 \(y=f(x)\)의 그래프를 나타내는 곡선을 아래 그림과 같이 \(2n-1\) 개의 분점 \[\rm P_1,\,P_2,\,\cdots,\,P_{2n-1}\] 으로 분할한다. 곡선의 소부분 \[\rm P_1P_2P_3,\,P_2P_3P_4,\,\cdots,\,P_{2n-2}P_{2n-1}P_{2n}\] 을 각각 \(y\)축에 평행한 축을 갖는 포물선의 일부로 보고, 이들 포물선의 호와 \(x\)추가 사이에 있는 면적을 위의 그림과 같이 각각 \[\rm S_1,\,S_2,\,\cdots,\,S_n\] 라 하고, 적분을 다음과 같이 근사시킨다. \[\int_a^bf(x)dx\fallingdotseq\rm S_1+S_2+\cdots+S_n\] 그래서 \(\rm S_1\)을 구해 보자. 호 \(\rm P_0P_1P_2\)를 \(y\)축에 평행한 축을 갖는 포물선의 일부라고 보고, 이 포물선의 방정식을 구한다. 계산을 간단히 하기 위하여 좌표축을 \(x\)축의 방향으로 평행이동하여 점 \(\rm P_1\)이 \(y\)축 위에 오도록 한다. 그래서 이 새로운 좌표에 관한 이 포물선의 방정식을 \[y=ax^2+bx+c\] 라 하자. \(\rm P_0,\,P_1,\,P_2\)의 새...