使用octave判断多项式能否生成素数 并计算10000以内生成素数的比例

以下是使用 Octave 编写的代码，可以判断多项式是否生成素数，并计算 10000 以内生成素数的比例：

function [is_prime, ratio] = poly_prime_ratio(a, b, c, n)
    % a, b, c 为多项式系数，n 为要判断的数的范围
    % 返回一个布尔值数组 is_prime，表示每个数是否是素数
    % 返回生成素数的比例 ratio
    is_prime = ones(n, 1);  % 先假设所有数都是素数
    for i = 1:n
        x = i;
        y = a * x^2 + b * x + c;  % 计算多项式值
        if y <= 1  % 如果多项式值小于等于1，则不是素数
            is_prime(i) = 0;
        else
            for j = 2:sqrt(y)
                if mod(y, j) == 0  % 如果能被整除，则不是素数
                    is_prime(i) = 0;
                    break;
                end
            end
        end
    end
    prime_count = sum(is_prime);  % 统计素数个数
    ratio = prime_count / n;  % 计算生成素数的比例
end

% 示例使用
[a, b, c] = [1, 1, 41];  % 多项式 x^2 + x + 41
n = 10000;
[is_prime, ratio] = poly_prime_ratio(a, b, c, n);
fprintf('生成素数的比例为：%f\n', ratio);

在上面的示例中，我们使用多项式 x^2 + x + 41 作为样例，计算 10000 以内生成素数的比例。你可以根据需要修改多项式系数和要判断的数的范围。