simple.m

来自「algotithm to calculate pi」· M 代码 · 共 8 行

function y=simple(n)
%n 表示循环位数
a=[];

for k=1:n
   a(k)=(-1).^(k-1)*(1/2).^(2*k-1)./(2*k-1)+(-1).^(k-1)*(1/3).^(2*k-1)./(2*k-1);
end
vpa(4*sum(a))