Sedikit berbagi,
ini salah satu contoh program matlab - metode bagi dua :
silahkan mencoba....
% Tugas Praktikum Metode Numerik
clc;
clear all;
disp(' TUGAS METODE NUMERIK ')
disp(' Menghitung akar - akar persamaan dengan Metode Bagi Dua ')
disp(' by : FENI DWI MASLIN (201010060311010) ')
disp(' Jurusan Pendidikan Matematika dan Komputasi / Universitas Muhammadiyah Malang')
disp(' ')
disp(' ')
syms x;
f=input('Masukkan Persamaan = ');
a=input('Masukkan batas a = ');
b=input('Masukkan batas b = ');
eps=input('Epsilon = ');
imax=input('Iterasi maksimum = ');
fa=subs(f,x,a);
fb=subs(f,x,b);
t=(a+b)/2;
ft=subs(f,x,t);
selisih = abs(b-a);
D=a;E=b;M=(a+b)/2;F=selisih;
disp(' ')
disp(' ')
disp('==================================================================')
disp('| Iterasi a b T (b - a) epsilon | ')
disp('==================================================================')
format short
tic
for k = 1:1:imax
t = (a+b)/2;
disp ([k-1,a,b,t,selisih])
if ( (fa)*(ft) <= 0)
a = t;
else
b = t;
end
selisih = abs(b-a);
D=[D;a];E=[E;b];M=[M;t];F=[F;selisih];
if (abs(b-a)) <= eps
break
end
end
disp ([k,a,b,t,selisih,eps])
disp('===================================================================')
waktu = toc;
plot([D,E,M,F])
legend('a','b','t','selisih')
xlabel('Iterasi ke: ')
ylabel('Iterasi metode Biseksi')
disp(' ')
disp([' akar :',num2str(t)])
disp([' jumlah iterasi :',num2str(k)])
disp([' selang waktu konvergensi :',num2str(waktu)])
Langganan:
Posting Komentar (Atom)