Добавил:
Upload
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Документ Microsoft Word (2)
.docfunction 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;
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]