function [x1,err,k,y]=dixot(f,x0,x1,deltax,epsilon,max1);
m=1;
for k=1:max1
    x1=x0+deltax;
    err=abs(x1-x0); 
    y=feval(f,x1);  x0=x1;
    if (m*y>0),  x0=x1,  m=(-1)*m; deltax=-deltax/2, end;
    if (err<epsilon)|(abs(y)<epsilon), break, end;  
    end;
    
    
        
        
    
        
    
    
    