close
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)
全站熱搜
留言列表