Программа метода дихотомии

function [x1,err,k,y]=dixot(fun,x0,x1,epsilon,max1);
%Метод дихотомии. Точность задана по абсолютной 
%разности корней.
%Шаг делится пополам сразу, но при этом известно,  что 
%функция на концах интервала разного знака.
%Вызов из основной программы
%[x1,err,k,y]= dixot1('fun',0,3,0.0001,100)
    
    for k=1:max1
    deltax=(x1-x0)/2;
    x2=x1;
    x1=x0+deltax;
    y=feval(fun,x1);
    if y<0, x0=x1; x1=x2; 
    end;
     err=abs(x1-x0); 
    if  err<epsilon, break, end;
    end;
Программа fun.m файла для метода дихотомии
function y=fun(x)
y=-5-2*x+x^3;

function y=fun(x)
y=x.^4+3*x.^3+4*x.^2+5*x+1;

function [x1,err,k,y]=dixot(fun,x0,x1,epsilon,max1)
%Метод дихотомии. Точность задана по абсолютной
%разности корней.
%Шаг делится пополам сразу, но при этом известно,  что
%функция на концах интервала разного знака.
%Вызов из основной программы
%[x1,err,k,y]= dixot('fun',0,3,0.0001,100)
 
for k=1:max1
    deltax=(x1-x0)/2;
    x2=x1;
    x1=x0+deltax;
    y=feval(fun,x1);
    if y<0, x0=x1; x1=x2;
    end;
    err=abs(x1-x0);
    if  err<epsilon, break, end;
end

clear, clc
 
a = -1; % нижний предел 
b = 1; % верхний предел
eps = 0.001; % точность
kmax = 500; %макс. количество итераций
x0 = dixot('fun',a,b,eps,kmax)
fplot(@fun,[a b])
grid on
hold on
plot(x0, fun(x0),'or')