function ret = Gamma(z)

    if isreal(z)

        ret = gamma(z);

    else

        if real(z) >= 0

            % 优化1：预先计算常量，避免重复计算

            n_terms = 50;  % 优化2：减少项数，测试合适的收敛点

            k = 0:n_terms;

            % 优化3：向量化计算，避免匿名函数开销

            gamma_vals = gamma(k + 1);

            terms = (-1).^k ./ (gamma_vals .* (k + z));

            % 优化4：设置积分选项提高精度

            options = {'RelTol', 1e-10, 'AbsTol', 1e-12};

            integral_part = integral(@(x) exp(-x) .* x.^(z-1), 1, Inf, options{:});

            ret = integral_part + sum(terms);

        else

            % 优化5：添加输入验证和特殊情况处理

            if abs(mod(real(z), 1)) < 1e-10 && abs(imag(z)) < 1e-10

                warning('Gamma:pole', 'Gamma function has poles at non-positive integers');

                ret = Inf;

            else

                ret = pi ./ (Gamma(1 - z) .* sin(pi .* z));

            end

        end

    end

end

  那我用deepseek优化一下！