Metode Bagi dua Pada Program Matlab

 cara membuat Metode Bagi Dua di program matlab:

 ketik yang di bawah ini pada M-FILE di MATLAB:

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('* Desy Kumala Sari(1101250695)                                                                  *')
disp('* Fitriati(1101250700)                                                                          *')
disp('* Muhlisah  (1001210718)                                                                        *')
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);