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=ap0-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;
end
plot(allp)

創作者介紹
創作者 jethro0119 的頭像
jethro0119

滷蛋與貢丸

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