matlab 抛物线法求方程全部根时怎么剔除已经求得的根

MATLAB 抛物线法求方程的全部根，求得一个根后如何剔除， 使用f1=f/(x-X1)一直报错，请问应该怎么解决呢？
clear all
clc
syms x

f=@hanshu;
X1=paowuxian(f,1,2,3);
f1=f/(x-X1);

function f=hanshu(x)
f=x^6-2*x^5+2*x^4+x^3+6*x^2-6*x+8;
end

function root= paowuxian(f,x0,x1,x2)
disp(f);
max=1000;
N=1;
eps1=1e-10;eps2=1e-10;

while (Nf0=feval(f,x0);f1=feval(f,x1);f2=feval(f,x2);
    
    lemda=(x2-x1)/(x1-x0);delta=(x2-x0)/(x1-x0);
    
    g=f0*lemda^2-f1*delta^2+f2*(lemda+delta);
    
    LEMDA=-2*f2*delta/(g+sign(g)*sqrt(g^2-4*f2*delta*lemda*(f0*lemda-f1*delta+f2)));
    
    x3=x2+LEMDA*(x2-x1);
    f3=feval(f,x3);
    
    
    if abs(x3)<1
        delta=abs(x3-x2);
    else
        delta=abs(x3-x2)/abs(x3);
    end
    
    if (delta<=eps1 )|| (abs(f3)<=eps2)
        root=x3;
        froot=f3;
        fprintf('The root is')
        disp(root)
        fprintf('The froot is')
        disp(froot)
        break
    else
        N=N+1;
        x0=x1;x1=x2;x2=x3;
    end   
    disp(N);
    
end
end

最好是把报错信息贴出来，方便分析

参考链接

Matlab 求方程的根_真·skysys的博客-CSDN博客_matlab求根 求x^2-3x+1=0的根。方法一：利用MATLAB多项式求根函数roots来求根。p=[1,-3,1];x=roots(p)绘图：x=-5:0.1:5;y1=x.*x-3*x+1;y2=zeros(size(x));plot(x, y1, x, y2);方法二 ： 利用求单变量非线性方程根的函数fzero，求方程在某个初始点附近的实根。f=@(x) x*x-3*x+1... https://blog.csdn.net/qq_33583069/article/details/96641652?ops_request_misc=&request_id=&biz_id=102&utm_term=matlab%20%E6%8A%9B%E7%89%A9%E7%BA%BF%E6%B3%95%E6%B1%82%E6%96%B9%E7%A8%8B%E5%85%A8%E9%83%A8%E6%A0%B9%E6%97%B6%E6%80%8E%E4%B9%88%E5%89%94%E9%99%A4%E5%B7%B2%E7%BB%8F%E6%B1%82%E5%BE%97%E7%9A%84%E6%A0%B9&utm_medium=distribute.pc_search_result.none-task-blog-2~blog~sobaiduweb~default-0-96641652.nonecase&spm=1018.2226.3001.4450

看看这个
https://ask.csdn.net/questions/7583365

看一下这篇文章：