MATLAB代码验证费马大定理

下面是一个简单的 MATLAB 代码，用于验证费马大定理：

% 设置指数 n 的值
n = 5;

% 验证费马大定理
for a = 1:n-1
    if mod(a^n, n) ~= mod(a, n)
        disp(['费马大定理不成立，a = ' num2str(a)]);
        return;
    end
end

disp(['费马大定理成立，n = ' num2str(n)]);

该代码通过枚举 $a$ 的值来验证费马大定理。如果对于某个 $a$，$a^n \not\equiv a \pmod{n}$，则输出错误信息；否则输出成功信息。