simple.m

algotithm to calculate pi

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))