matlab实现steffensen加速迭代

function [pQ] = steff（ fdfp0deltaepsilonmax1 ）
R=zeros（max13）;
R（11）=p0;
for k=1:max1 
    for j=2:3
        nrdenom=feval（dfR（kj-1））;
        if nrdenom==0
            break
        else
            R（kj）=R（kj-1）-feval（fR（kj-1））/nrdenom;
        end
        aadenom=R（k3）-2*R（k2）+R（k1）;
        
        if aadenom==0
            break