[수치해석] 적분수치해석 (Numerical Integration)

목차

1. 문제의 정의
2. 해석적인 방법 (내용 없음)
3. 수치적인 방법 프로그램 소스 및 설명
4. 결과
5. 토의

본문내용

1. 문제의 정의
정적분을 아래의 다섯 가지 방법으로 계산한다.
(1) analytical method
(2) Trapezoidal method (n=12)
(3) Simpson's 1/3 (n=12)
(4) Simpson's 3/8 (n=12)
(5) Gauss Quadrature (n=5)

2. 해석적인 방법

3. 수치적인 방법 프로그램 소스 및 설명
(1)
/* Trapezoidal */
#include <stdio.h>
#include <math.h>

main()
{
int i,n;
float f(float y),absf(float z);
float a,b,h,itrap,sum,itrue,err;
/* 변수 및 함수의 선언 n은 적분 구간 사이의 segment의 수를 나타내고, a는 적분의 lower limit, b는 적분의 upper limit이다. h는 적분 구간 간격으로 적분 구간 전체를 segment의 수로 나눈 값이 된다. itrap이 적분 수행후의 수치적분값이고, itrue는 적분의 실제값이다. err는 수치적분의 상대오차이다. f()는 적분의 intgrand내에 들어있는 원함수이다. absf()는 절대값을 나타낸다. */

참고 자료

없음