function?x=chase（Af）
%x?=?chase（Af），A為對角矩陣，f為右側常量x為方程的解
%該函數正常工作前提為A的行列式不為零
n=rank（A）;
[LU]=LU（An）;
sum=0;
y=zeros（n1）;
y（1）=f（1）;
for?k=2:n
????for?j=1:k-1
????????sum=sum+L（kj）*y（j）;
????end
????y（k）=f（k）-sum;
????sum=0;
end
x=zeros（n1）;
x（n）=y（n）/U（nn）;
k=n-1;
while（k>=1）
????for?j=k+1:n
????????sum=sum+U（kj）*x（j）;
????end
????x（k）=（y（k）-sum）/U（kk）;
????sum?=?0;k=k-1;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????454??2019-07-30?09:50??Code\chase.m

?????文件??????29019??2019-07-30?09:41??Code\example.jpg

?????文件???????1771??2019-07-30?09:49??Code\example.m

?????文件????????332??2019-07-30?09:40??Code\LU.m

?????目錄??????????0??2019-07-30?09:49??Code

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

????????????????31576????????????????????5