twitter
rss

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)])

0 komentar:

Posting Komentar