fixed point - 主要程式AVS檔 @ 滷蛋與貢丸

Oct 07 Tue 2008 21:30
fixed point - 主要程式AVS檔

clear

format long e
p0=1.5
tol=0.000000000001;
N0=100;
allp=p0;

% for i=1:N0
%     evaluate
%     if abs(p-p0)<tol
%         stop=1
%     end
%     p0=p;
%     allp(i+1,1)=p0;
%
% end
% plot(allp)
pr=1.365230013414097;
stop=0;
i=1;
allpr=abs(p0-pr);
while (stop==0 & i< N0+1)
    evaluate
    if abs(p-p0)<tol
        p
        stop=1;
    end
    i=i+1;
    p0=p;
    allp(i+1,1)=p0;
    allpr(i+1,1)=abs(p0-pr);
end
figure(2)
plot(allpr)