close
clear
p0=0.78;
p1=0.7;
Tol=10^(-10);
N0=10;
stop=0;
i=2;
g0=cos(p0)-p0;
g1=cos(p1)-p1;
while(i<N0 & stop==0)
temp=g1*(p1-p0)/(g1-g0);
p=p1-temp;
tempabs(i,1:3)=[i p abs(p-p1)];
if abs(p-p1)<Tol
stop=1
p
end
i=i+1;
p0=p1;
g0=g1;
p1=p;
g1=cos(p)-p;
end
全站熱搜
留言列表