?? gaussliezhuyuan.m
字號:
function [x,q]=Gaussliezhuyuan(A,b)
B=[A,b];
if det(A)==0
disp('方程組不存在唯一解!');
return;
end
n=length(A);
for k=2:n
[m,d]=max(B(:,k-1));
if d~=k-1 && d>k-1
temp=B(k-1,:);
B(k-1,:)=B(d,:);
B(d,:)=temp;
end
for kk=k:n
B(kk,:)=-B(kk,k-1)*B(k-1,:)/B(k-1,k-1)+B(kk,:)
end
end
x(n,1)=B(n,n+1)/B(n,n);
for k=n-1:-1:1
sx=B(k,k+1:n)*x(k+1:n,1);
x(k,1)=(B(k,n+1)-sx)/B(k,k);
end
q=(n^3)/3+n^2-1/3*n;
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -