function out=Mullers(x1,x2,x3)
fprintf(' \t\t Mullers method \n')
fprintf(' \tsteps\tx1\t\tx2\t\tx3\t\tf(x1)\t\tf(x2)\t\tf(x3)\t\tc1\t\tc2\t\td\t\ts\t\tx4\t\tf(x4) \n')
for i=1:7
    y1=f(x1);
    y2=f(x2);
    y3=f(x3);
    c1=(y2-y1)/(x2-x1);
    c2=(y3-y2)/(x3-x2);
    d=(c2-c1)/(x3-x1);
    s=c2+(d*(x3-x2));
    x4=x3-((2*y3)/(s+sign(s)*sqrt(s.^2-4*y3*d)));
    x1=x2;
    x2=x3;
    x3=x4;
%fprintf(8d,8.4f,8.4f,8.4f,8.4f,8.4f,8.4f,)
fprintf('%8d  %8.4f%8.4f%8.4f  %8.4f  %8.4f    %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f\n', i,x1,x2,x3,f(x1),f(x2),f(x3),c1,c2,d,s,x4,f(x4))
end %for
end %function
function ret=f(x)
ret=x.^6-2;
end