[수치해석]사다리꼴 적분법과 Simpson적분법

소개글

사다리꼴 적분법과 Simpson적분법을 메틀랩으로 코딩하여 여러 가지 적분식을 풀어 보았다.

목차

1. 사디리꼴 적분법에 의한 적분식 해법
2. Simpson적분법에 의한 적분식 해법
3. 로켓의 상승 높이 구하기
4. 관을 통과하는 유량 구하기

본문내용

1. 사다리꼴 적분공식에 의한 적분
(1) I = int(1/x^2, 1, 2)

% I = int(1/x^2,1,2)에 대한 사다리꼴적분
f=`1/x^2`; h=[0.5 0.25 0.1];
I=[1 2]; pe=0;
fprintf(` h 근사적분값 오차 비율 \n`)
fprintf(`------------------------------------------ \n`)
for i=1:length(h)
T=0;
for n=1:(I(2)-I(1))/h(i)
x=I(1)+(n-1)*h(i); T=T+eval(f);
x=I(1)+(n*h(i)); T=T+eval(f);
end
T=(h(i)/2)*T;
ff=`-1/x`; % int(1/x^2,1,2) = -1/2 + 1
t=0;
x=I(2); t=t+eval(ff);
x=I(1); t=t-eval(ff);
error=t-T;
ratio=pe./error;
fprintf(`%4.2f %16.14f %5.5e %5.3f \n`,h(i),T,error,ratio)
pe=error;
end

참고 자료

없음