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
???