牛顿法解方程

clear;
clc;
eps=1e-3;
x0=-1;
x1=f13(x0);
while (abs(x1-x0)>=eps)
    x0=x1;
    x1=f13(x1);
end
disp(x1);