Script Matlab Metode Bagi Dua

clear all;

clc;
disp('*************************************************************************************************        ')
disp('*                                     METODE NUMERIK                                              *')
disp('-------------------------------------------------------------------------------------------------')
disp('*                    Menghitung akar - akar persamaan dengan c                 *')
disp('*                   Jurusan Pendidikan Matematika  fakultas Tarbiyah dan Keguruan               *')
disp('*                                IAIN Antasari Banjarmasin                                      *')
disp('-------------------------------------------------------------------------------------------------')
disp('* Oleh  :                                                                                       *')
disp('* Kelompok      6                                                                               *')

disp('* Hairianor(1201250929)                                                                  *')

disp('* Muhammad Arif Rahman(1201250938)                                                         *')

disp('* Kartini  (1201250859)                                                                        *')

disp('* Hardiati Rizki Amelia  (1201250850)                                                             *')

disp('* Lika Wartita  (1201250864)                                                                        *')
disp('*************************************************************************************************        ')
disp('*                                                                                               *')
syms x;
f=input('Masukkan persamaan=');
a=input('Masukkan nilai a=');
b=input('Masukkan nilai b=');
e=input('Masukkan nilai epsilon=');
lebar=abs(b-a)/2;
r=0;
disp('  r       a          c          b        f(a)       f(c)       f(b)      lebar');
disp('---------------------------------------------------------------------------------');
clama=a;
cbaru=b;
while (lebar > e) & (clama ~= cbaru);
    fa=subs(f,x,a);
    fb=subs(f,x,b);
    clama=cbaru;
    c=(a+b)/2;
    cbaru=c;
    fc= subs(f,x,c);
    fprintf('%3.0f %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f %10.6f \n', r, a, c, b, fa, fc, fb, lebar);
    if fa*fc < 0
    b=c;
    else
      a=c;
    end
    lebar=abs(b-a)/2;
    r=r+1;
end
akar = c;
fprintf('\nHampiran akar      = %10.6f\n', akar);
fprintf('Jumlah leleran     = %g\n', r);

Semoga Bermanfaat...

Selamat Mencoba...