Rabu, 22 Juni 2011

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


0 komentar:

Posting Komentar

next page