matlab#10

목차

1.일반적분

2.사다리꼴적분

3.Simpson 1/3

4.Simpson 3/8

5.소스코드

본문내용

x+3;
x=x2; fx2=3*x^3-3*x^2-3*x+3;
x=x3; fx3=3*x^3-3*x^2-3*x+3;
x=x4; fx4=3*x^3-3*x^2-3*x+3;
x=x5; fx5=3*x^3-3*x^2-3*x+3;
x=x6; fx6=3*x^3-3*x^2-3*x+3;
x=x7; fx7=3*x^3-3*x^2-3*x+3;
x=x8; fx8=3*x^3-3*x^2-3*x+3;
x=x9; fx9=3*x^3-3*x^2-3*x+3;
if n==3
%----------------------------------------------사다리꼴법------------------
S1=1/2*h*((fx0+fx3)+2*(fx1+fx2));
err1=abs(S-S1);
fprintf('사다리꼴 적분법 -> %2.2f // 원함수 순수적분값과의 오차 -> %2.2f \n\n',S1,err1)
%------------------------------------------------심슨 1/3------------------
S2=1/3*h*((fx0+fx3)+4*(fx1)+2*(fx2));
err2=abs(S-S2);
fprintf('Simpson 1/3 -> %2.2f // 원함수 순수적분값과의 오차 -> %2.2f \n\n',S2,err2)
%------------------------------------------------심슨 3/8------------------
S3=3/8*h*(fx0+3*(fx1)+3*(fx2)+fx3);
err3=abs(S-S3);
fprintf('Simpson 3/8 -> %2.2f // 원함수 순수적분값과의 오차 -> %2.2f \n\n',S3,err3)
end
%--------------------------------------------------------------------------
%%-------------------------------------------------------------------------
if n==6
%----------------------------------------------사다리꼴법------------------
S1=1/2*h*((fx0+fx6)+2*(fx1+fx2+fx3+fx4+fx5));
err1=abs(S-S1);
fprintf('사다리꼴 적분법 -> %2.2f // 원함수 순수적분값과의 오차 -> %2.2f \n\n',S1,err1)
%------------------------------------------------심슨 1/3------------------
S2=1/3*h*((fx0+fx6)+4*(fx1+fx3+fx5)+2*(fx2+fx4));
err2=abs(S-S2);
fprintf('Simpson 1/3 -> %2.2f // 원함수 순수적분값과의 오차 -> %2.2f \n\n',S2,err2)
%------------------------------------------------심슨 3/8------------------
S3=3/8*h*(fx0+3*(fx1)+3*(fx2)+fx3)+3/8*h*(fx3+3*(fx4)+3*(fx5)+fx6);
err3=abs(S-S3);
fprintf('Simpson 3/8 -> %2.2f // 원함수 순수적분값과의 오차 -> %2.2f \n\n',S3,err3)
end
%--------------------------------------------------------------------------
%%--------------------------------------------------------------------------

키워드

사다리골적분,   simpson 1/3,   simpson 3/8,   심슨적분법,   심슨,   적분,   Matlab,   수치해