function?[RARBnX]=gaus（Ab）
B=[A?b];?
n=length（b）;?
RA=rank（A）;?
RB=rank（B）;
zhica=RB-RA;
if?zhica>0
???disp（‘因?yàn)镽A不等于RB，所以此方程組無解.‘）
???return
end
if?RA==RB
???if?RA==n
???disp（‘請注意：因?yàn)镽A=RB=n，所以此方程組有唯一解.‘）?
???X=zeros（n1）;
??????for?p=?1:n-1
for?k=p+1:n
?????????????m=?B（kp）/?B（pp）;?B（kp:n+1）=?B（kp:n+1）-m*?B（pp:n+1）;
end
end
?????????b=B（1:nn+1）;A=B（1:n1:n）;?X（n）=b（n）/A（nn）;
??????for?q=n-1:-1:1
?????????X（q）=（b（q）-sum（A（qq+1:n）*X（q+1:n）））/A（qq）;
??????end
else
?????????disp（‘請注意：因?yàn)镽A=RB???end
end
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????632??2014-10-14?21:35??Gaus\gaus\gaus.asv

?????文件????????618??2014-10-14?22:00??Gaus\gaus\gaus.m

?????文件????????189??2014-10-14?21:43??Gaus\gaus\mainn3.m

?????文件??????30984??2014-10-14?21:58??Gaus\Gaus.docx

?????目錄??????????0??2014-11-13?20:37??Gaus\gaus

?????目錄??????????0??2014-11-13?20:38??Gaus

-----------?---------??----------?-----??----

????????????????32423????????????????????6