MATLAB 找到n阶方阵中对角线上的连续素数个数 不使用循环

找到n阶方阵中对角线上的连续素数个数

n=input('请输入矩阵阶n');
A=fix(1000*rand(n));
%计算对角线上的连续素数个数
% 不使用循环，结果存入y
B=diag(A)'
bo=[0 isprime(B) 0];    %让起始和最后的元素与中间的操作统一
pos=find(bo==0);
jian=diff(pos)
y=max(jian-1);
disp(y)