Membuat Aplikasi Penghitung Integral dengan Metode Trapesium
Menggunakan M-FILE MATLAB R2008a
By Hanin Fathan Nurfina Istiqomah
Buka Matlab R200a
Pilih File => New =>M.file
clear;clc;
ki = input ('batas kiri? ');
ka = input ('batas kanan? ');
pita = input ('jumlah pita?');
%persamaan yang di integralkan (-x^2+4x)
lbr = abs( (ka - ki) /pita);
%posisi titik
for i = 1:pita+1
x(i)=ki + (i-1)*lbr;
end;
%nilai fungsi masing-masing titik
for i = 1:pita+1
f(i)=-x(i)^2+4*x(i);
end;
%menghitung luas masing - masing pita dan menjumlahkan
luas=0;
for i = 1:pita
luas=luas+((f(i))+(f(i+1)))*0.5*lbr;
end;
fprintf('Hasil integrasi = %2.5f\n',luas);
Setelah itu pilih menu debug (klik Save and Run)
Simpan nama jangan ada spasi angka ataupun titik, hanya huruf
setelah itu muncul Command Window seperti di bawah ini :
batas kiri? 0
batas kanan? 4
jumlah pita? 4
Hasil integrasi = 10.00000
>>
Catatan pertama kali akan muncul "batas kiri?" isikan 0 lalu enter dan muncul "batas kanan? isikan 4 tekan enter dan muncul "jumlah pita?" isikan 4 dan tekan enter lagi maka akan muncul jawaban "Hasil integrasi = 10.00000" sebagai jawaban
Program ini hanya berlaku untuk f(i)=-x(i)^2+4*x(i)
Contoh lain
pada editor tampis seperti berikut
clear;clc;
ki = input ('batas kiri? ');
ka = input ('batas kanan? ');
pita = input ('jumlah pita? ');
%persamaan yang di integralkan (-x^2+4x)
lbr = abs( (ka - ki) /pita);
%posisi titik
for i = 1:pita+1
x(i)=ki + (i-1)*lbr;
end;
%nilai fungsi masing-masing titik
for i = 1:pita+1
f(i)=sin(x(i));
end;
%menghitung luas masing - masing pita dan menjumlahkan
luas=0;
for i = 1:pita
luas=luas+((f(i))+(f(i+1)))*0.5*lbr;
end;
fprintf('Hasil integrasi = %2.5f\n',luas);
maka command window akan muncul
batas kiri? 0
batas kanan? 2*pi
jumlah pita? 200
Hasil integrasi = -0.00000
>>
luas tidak ada yang bernilai minus maka editor dapat ditambah abs pada line 18, menjadi berikut :
clear;clc;
ki = input ('batas kiri? ');
ka = input ('batas kanan? ');
pita = input ('jumlah pita? ');
%persamaan yang di integralkan (-x^2+4x)
lbr = abs( (ka - ki) /pita);
%posisi titik
for i = 1:pita+1
x(i)=ki + (i-1)*lbr;
end;
%nilai fungsi masing-masing titik
for i = 1:pita+1
f(i)=sin(x(i));
end;
%menghitung luas masing - masing pita dan menjumlahkan
luas=0;
for i = 1:pita
luas=luas+(abs(f(i))+abs(f(i+1)))*0.5*lbr;
end;
fprintf('Hasil integrasi = %2.5f\n',luas);
Maka command windownya adalah
batas kiri? 0
batas kanan? 2*pi
jumlah pita? 200
Hasil integrasi = 3.99967
>>
0 komentar:
Posting Komentar