کد پیاچپی:
function X=hazf(A,B)
[r,c]=size(A);
for i=1:r-1
for j=i+1:r
lambda=A(j,i)/A(i,i);
A(j,:)=A(j,:)-lambda*A(i,:);
B(j)=B(j)-lambda*B(i);
end
end
X=zeros(r,1);
X(r)=B(end)/A(end,end);
for i=r-1:-1:1
A(i,i+1:end)
X(i+1:end)
X(i)=(B(i)-(A(i,i+1:end)*X(i+1:end)))/A(i,i);
end
end
این دستور رو تایپ میکنم
کد پیاچپی:
>> X=hazf([2 4 -2 -2;1 2 4 -3;-3 -3 8 -2; -1 1 6 -3],[-4;5;7;7])
این جوابو میده
کد پیاچپی:
ans =
Inf
ans =
NaN
ans =
5 -2
ans =
NaN
NaN
ans =
4 -2 -2
ans =
NaN
NaN
NaN
X =
NaN
NaN
NaN
NaN
درحالی که هیچ ارروری نمیده و لی جوابم نمیده