Integrasi Numerik

Metode TRAPESIUM

f=@(x)(1-x.^2);

a=input('aku mauk a:');

b=input('aku mauk b:');

n=input('aku n:');

I=0;

h=(b-a)/n;

if a>b |a==b

    disp ('Ups...salah...')

end

for i=1:n-1

    I=I+2*f(a+h*i);

end

I=I+f(a)+f(b);

I=(h/2)*I;

format long

I

Metode SIMPSON 1/3

f=@(x)(1-x.^2);

a=input('aku mauk a:');

b=input('aku mauk b:');

n=input('aku n:');

I=0;

h=(b-a)/n;

if a>b |a==b

    disp ('Ups...salah...')

end

for i=1:2:n-1

    X=I+4*f(a+h*i);

end

for i=2:2:n-1

    A=X+2*f(a+h*i);

end

   

B=I+X+A+f(a)+f(b);

I=(h/3)*B;

format long

I

Metode SIMPSON 3/8

f=@(x)(1-x.^2);

a=input('aku mauk a:');

b=input('aku mauk b:');

n=input('aku n:');

I=0;

h=(b-a)/n;

if a>b |a==b

    disp ('Ups...salah...')

end

for i=1:1:n-1

    if (mod(i,3)~=0)

        A=I+3*f(a+h*i);

    end

end

    for i=3:3:n-3

    X=I+2*f(a+h*i);

end

B=I+X+A+f(a)+f(b);

I=(3*h/8)*B;

format long

I