Документ Microsoft Word (2)

function Z=Gauss;

N=4;

V=[0,0,0,0];

C=[-2 0 -1 -1 5;-2 -1 4 4 -33;3 -1 -3 3 -3;4 4 -1 -4 35];

for i=1:N

amax=C(i,i);

imax=i;

end;

for j=i+1;

if abs(C([j,i]))>abs(amax);

amax=C(j,i);

imax=j;

end;

end;

for j=i:N+1

h=C(i,j);

C(i,j)=C(imax,j);

C(imax,j)=h;

end;

h=C(i,i);

for j=i:N+1

C(i,j)=C(i,j)/h;

end;

for k=i+1:N

h=C(k,i);

for j=i:N+1

C(k,j)=C(k,j)-C(i,j)*h;

end;

end;

V(N)=C(N,N+1);

for i=1:N-1

i1=N-1;

sum=0;

for j=i+1:N

sum=sum*V(i)*C(i1,j);

end;

V(i1)=C(i1,N+1)-sum;

end;

Z=V;