最近YY子开始着手学习Matlab...
本着从简单开始的心态先做点'Reinvent the wheel' 性质的东西 应该也算是不错的开始呢..?
总之是本着这样的心态写了点奇怪的东西...
(刚刚开始学双拼貌似说话都变得奇怪了ORZ)
%Gaussian Elimination
function sol=Gaussian_Elimination(a)
m=size(a,1);
n=size(a,2);
for i=1:(m-1)
[tmp,lb]=max(a(i:m,i));
lb=lb+i-1;
a([i,lb],:)=a([lb,i],:);
if a(i,i)==0
continue;
end
for j=(i+1):m
a(j,:)=a(j,:)-a(i,:)*(a(j,i)/a(i,i));
end
end
sol=a;
end