Metode Numerik: Scribs Matlab Metode Secant

clear; clc;
disp('--------------------------------')
disp('program   : Metode Secant')
disp('programer : SYAHARUDDIN ')
disp('--------------------------------')
f=input ('f(x) = ');
x1=input('x1 = ');
x2=input('x2 = ');
imax=input('Iterasi = ');
galat1=input('Error = ');
iter=0;
fprintf('\n     Iterasi      Akar           f(Akar)            Galat\n');
for k=1:imax
    iter=iter+1;
    %Rumus Secant
    x3=x2-((x2-x1)*feval(f,x2)/(feval(f,x2)-feval(f,x1)));
    galat=abs((x3-x1)/x3);
    x1=x3;
    y=feval(f,x1);
    fprintf('%10.0f     %6.10f      %6.10f      %6.10f\n',[iter;x1;y;galat])
    if (galat<galat1 || (iter>imax)),break,end
end
fprintf('Akarnya adalah = %6.10f\n',x1)