matlab计算圆周率的一种方法

n=1;
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);