Matlab .Тема. Прямые методы решения систем линейных алгебраических уравнений.

[Bushido]

Методом Гаусса с приведением исходной матрицы к нижнему треугольному вид. Код

PHP код:

clc,clear 
n = input('Размерность матрицы'); 
A=[2 1 1 1;1 2 1 2 ;1 1 2 4] 
for p=1:1:n 
B=zeros(n,1); 
for i=p:1:n 
if abs(A(i,p))>B; 
B=A(i,p); 
A(i,p)=B; 
end 
end 
if A(p,p)<B 
x=A(p,:); 
A(p,:)=A(n,:); 
A(n,:)=x; 
end 
K=A(p,:)/A(p,p); 
A(p,:)=K; 
for j=p+1:1:n 
S=(K*-A(j,p))+A(j,:); 
A(j,:)=S; 
end 
end 
A 
B=A(:,n+1); 
x=zeros(n,1); 
x(n)=B(n)/A(n,n); 
for k=n-1:-1:1 
x(k)=(B(k)-A(k,k+1:n)*x(k+1:n))/A(k,k); 
end 
x 


помогите исправить программу пожалуйста.