Membuat biseksi menggunakan Matlab

Coding :

clc;
disp('Biseksi =>>persamaan : x*e^(-x)+1');
syms x;
f=input('masukkan persamaan f(x): ');
a=input('masukkan nilai a : ');
b=input('masukkan nilai b : ');
et=input('masukkan Toleransi Error  : ');
t=abs(b-a);
i=1;
e=2.71;
disp(' i      a       b       c       f(a)       f(b)       f(c)     E');
disp('----------------------------------------------------------------');
cAwal=a;
cAkhir=b;
while (t > et ) & (cAwal ~= cAkhir);
    fa=subs(f,x,a);
    fb=subs(f,x,b);
  
    cAwal=cAkhir;
    c=(a+b)/2;
        cAkhir=c;
        fc= subs(f,x,c);
    fprintf('%3.0f %6.4f %6.4f %12.10f %7.4f %7.4f %7.4f %7.4f \n', i, a, b, c, fa, fb, fc, e);
    if fa*fc < 0
       b=c;
    else
       a=c;
    end
    t=abs(b-a);
    i=i+1;
end


Hasil :

Biseksi =>>persamaan : x*e^(-x)+1
masukkan persamaan f(x): 2
masukkan nilai a : 4
masukkan nilai b : 6
masukkan Toleransi Error  : 7
 i      a       b       c       f(a)       f(b)       f(c)     E
----------------------------------------------------------------