mypi.m

《matlab使用详解》的部分源码

%首先设定需要的变量
t=1;  %变量t表示计算式括号中的各项   
pi=0;  %pi代表圆周率，首先置0
n=1;  %n为表示分母的变量
s=1;  %变量s用正负数的改变，控制“+、-”的转换
%使用while循环语句
while abs(t)>=1e-7   %“t的绝对值大于10的-7次方”为循环执行条件
    pi=pi+t;   %循环体
    n=n+2;
    s=-s;
    t=s/n;
end   %结束循环
pi=4*pi;  %最终圆周率的值还要用括号中的计算结果乘以4
fprintf('pi=%f\n',pi)  %输出结果