Matlab判断亲密数对

如果两个连续自然数的乘积减1是素数，那么这两个连续的自然数称为亲密数对，该素数为亲密素数。例如，2*3-1=5,5是素数，故<2,3>是亲密数对，5是亲密素数。不调用系统内部函数，求区间[2，49]内的亲密数对并求出所有亲密数对之和

i=1;
j=0;
sum=0;
disp('亲密数对为：')
for x=2:49 %对区间内的数进行依次判断
    for i=1:(x*(x+1)-1) 
    if mod(x*(x+1)-1,i)==0 %通过取余函数确认其乘积减1是否为素数，即是否只有1和它本身两个因数 
        j=j+1; %作为判断成功的条件
      end %结束每一次取余
    end %结束判断条件j的取值
    if j==2 %其乘积减1为素数，则利用数字转字符函数表示两个自然数
        disp([num2str(x),'和',num2str(x+1)])
       sum=sum+(x*(x+1)); %计算亲密素数之和
        j=0; %将判断条件初始化，进入下一次循环
    else
        j=0; %x不是素数则进入下一个数的判断
    end %结束对j是否等于2即是否是亲密素数的判断
end %结束2~49的判断
disp(['亲密素数之和为：',num2str(sum)])