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

arrow
arrow
    全站熱搜

    jethro0119 發表在 痞客邦 留言(0) 人氣()

    E-mail轉寄