exm060105.m

matlab电子书籍

clear
clc
while 1																		
	n=input('请输入一个能被4整除的正整数！ n =  ');
    if mod(n,4)==0																
        break																	
    end																		
end																			
G=logical(eye(4,4)+rot90(eye(4,4))); 
m=n/4;
K=repmat(G,m,m);        
N=n^2;
A=reshape(1:N,n,n);     
A(K)=N-A(K)+1             

s0=round(n*(n*n+1)/2);	
disp([int2str(n),' 阶魔方矩阵的标称和是   ',int2str(s0)])
Ns0=round(2*(n+1));	
B=A';
SC=sum(A);				
SR=sum(B);				
	Sd=sum(diag(A)); 		
	Sdi=sum(diag(B));		
LS=[SC,SR,Sd,Sdi]==s0;	
NS=round(sum(LS));		
if NS==Ns0
	disp('经验证，A是魔方矩阵。')
else
	disp('经验证，A不是魔方矩阵。')
end