A = pascal(7) %生成7*7 pascal矩阵
B = A([1,3,5],[5,6,7]) %构成矩阵B
C = A(5:7,1:3)%构成矩阵C
D = B-1*C
A =
1 1 1 1 1 1 1
1 2 3 4 5 6 7
1 3 610152128
1 41020355684
1 5153570 126 210
1 62156 126 252 462
1 72884 210 462 924
B =
1 1 1
152128
70 126 210
C =
1 515
1 621
1 728
D =
0-4 -14
1415 7
69 119 182
第二题:
function A=prime_number()
N=100
j=1
for i=2:N
if prod(rem(i,2:(i-1)))~=0
A(j)=i
j=j+1
end
end
A =
2 3 5 7111317192329313741434753596167717379838997
基本解答如下:(1-1):
function [b,y] = by(m,n)
b=lcm(m,n)%最小公倍数
y=gcd(m,n)%最大公约数
end
(1-2):
function t=sushu(n)
m=floor(sqrt(n))
for i=2:m
if mod(n,i)==0
t=0
break
else
t=1
end
end
end
(派正猛2):可以一次算出每一项,发现他们都是同一个值,这个值的求法:SUM=sum(1:n^2)/n%n阶魔幻阵特有的和
function [A,H,L,DJ,SUM]=mag(n)
A=magic(n)%产生一个n阶魔幻阵
H=sum(A)%每行的和
L=sum(A)%每列的和
DJ=sum(diag(A))%对角线和
SUM=sum(1:n^2)/n%n阶魔幻阵特有的和
end
(3):这题我给出“该数各位数字之和为奇数”清漏的求法,结合上面有判断素数的代码自己可以求出
该数各位数字之和为奇数:
function funct
for i=2:999
%l=floor(log10(i))+1 %提取位数
n1=floor(i/100) %取出百位数字n1
n2=mod(floor(i/10),10)%取出十尘桥位数字n2
n3=mod(i,10) %取出个位数字n3
nsum=n1+n2+n3
if rem(nsum,2)~=0%判断是否为奇数
disp (i)
end
end
end
自己可以调试看看,祝你成功!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)