matlab計算圓周率的一種方法
pi4=1/n;
delta=abs(pi-pi4*4);
jingdu=0.001;
tic
while delta>jingdu
n=n+1;
pi4=pi4+(-1)^(n-1)*1/(2*n-1);
delta=abs(pi-pi4*4);
end
s1=sprintf('精度達到%f時的項數%d',jingdu,n);
s2=sprintf('pi的近似值是%f', pi4*4);
disp(s1),disp(s2)
h=toc;
s3=sprintf('使用迴圈語句的時間是%f秒',h);
disp(s3)
tic
k=1:1:1000;
pi_4=sum((-1).^(k-1).*1./(2*k-1));
cha=abs(pi-4*pi_4);
h1=toc;
st1=sprintf('精度達到%f時的項數%d',cha(end),k(end));
st2=sprintf('pi的近似值是%f', pi_4*4);
st3=sprintf('使用向量的時間是%f秒',h1);
disp(st1);
disp(st2);