for i=1:2 p=3-i; for j=1:p if c(i,i)==0 n=c(i,:); c(i,:)=c(i+j,:); c(i+j,:)=n; end end for k=1:p cons=c(i+k,i)/c(i,i); for r=1:4 c(i+k,r)=c(i+k,r)-c(i,r)*cons; end end end c z=c(3,4)/c(3,3) y=(c(2,4)-c(2,3)*z)/c(2,2) x=(c(1,4)-c(1,3)*z-c(1,2)*y)/c(1,1)
clear all
ResponderEliminarclc
a=[1 2 2;9 7 9;1 2 3];
b=[-1;6;8];
c=[a,b];
for i=1:2
p=3-i;
for j=1:p
if c(i,i)==0
n=c(i,:);
c(i,:)=c(i+j,:);
c(i+j,:)=n;
end
end
for k=1:p
cons=c(i+k,i)/c(i,i);
for r=1:4
c(i+k,r)=c(i+k,r)-c(i,r)*cons;
end
end
end
c
z=c(3,4)/c(3,3)
y=(c(2,4)-c(2,3)*z)/c(2,2)
x=(c(1,4)-c(1,3)*z-c(1,2)*y)/c(1,1)