已知Wallis公式，编写一个函数是m文件对任意给定的误差s计算圆周率的满足误差要求 派减xn的绝对值小于x的近似值 用你编写的程序计算s等于1/2乘以10的-3次是圆周率的近似值

Wallis公式是一种利用无穷级数逼近圆周率π的方法，它表述为每个偶数阶乘除以奇数阶乘的比值之积。在MATLAB中，我们可以编写一个函数来应用这个公式，并确保误差满足要求。首先，我们需要创建一个函数，例如`piApproximation`，接受一个参数`s`表示误差范围。

function piEstimate = wallisPiApproximation(s)
    % 初始化变量
    x = 1;
    n = 1;
    piEstimate = 0;
    errorBound = s / (n * x);  % 计算初始误差

    % 迭代计算 Wallis 公式直到误差小于 s
    while errorBound > s / 10^(-3)
        term = (4^n) / (2*n + 1);
        piEstimate += term;
        n = n + 2;  % 提升到下一个偶数阶乘
        errorBound = abs(piEstimate - piEstimatePreviousIteration);  % 更新误差
        if errorBound < s / 10^(-3)
            break;
        end
        piEstimatePreviousIteration = piEstimate;  % 保存上一次迭代的估计值
    end

    piEstimate;
end

为了计算当`s = 1/2 * 10^-3`时的圆周率近似值，你可以像这样调用函数：

s = 0.5e-3;  % 设置误差范围
approxPi = wallisPiApproximation(s);
disp("When s is 1/2 times 10^-3, the approximate value of pi is:");
approxPi;