逐句解释matlab中martin函数代码

以下是matlab中martin函数代码的逐句解释：

function y = martin(x)
% This function computes the Martin's approximation of the gamma function.
% The approximation is valid for x > 12.
% The formula used is:
% gamma(x+1) = sqrt(2*pi)*(x+1/2)^(x+1/2)*exp(-x-1/2)*P(x)

% 定义一个函数martin，输入参数为x，输出为y
% 此函数用于计算gamma函数的Martin逼近。当x > 12时逼近有效。
% 使用的公式为：gamma(x+1) = sqrt(2*pi)*(x+1/2)^(x+1/2)*exp(-x-1/2)*P(x)

if x <= 12
    error('Error: x must be greater than 12.')
end
% 如果x <= 12，抛出错误提示，要求x必须大于12

P = 1/12/x;
for n = 1:7
    P = P*(4*n^2-1)/x^2/(n^2-1/4);
end
% 计算P(x)的值，通过循环计算P的每一项，最终结果为P(x)

y = sqrt(2*pi)*(x+1/2)^(x+1/2)*exp(-x-1/2)*P;
% 根据公式计算gamma(x+1)的近似值，结果为y