A=input('输入系数矩阵：');
B=input('输入结果矩阵：');
C=[A,B];
trapz(A,B)
[m,n]=size(C);
for i=1:(m-1)
numb=int2str(i);
disp(['第',numb,'次消元后的矩阵：']);
for j=(i+1):m
C(j,:)=C(j,:)-C(i,:)*C(j,i)/C(i,i);
end
C
end
disp('求解结果：');
x(m)=C(m,n)/C(m,m);
for i=(m-1):-1:1
x(i)=(C(i,n)-C(i,i+1:m)*x(i+1:m)')/C(i,i);
end
x
