x=S\K;
w=K-S*x;
x=x+S\w;

% уравнение Ax=b;
[n,~]=size(A);
d=diag(A);
g=b./d;
ID=diag(1./d);
B=eye()-ID*A;
 
if norm(B)>=1, disp('Алгоритм не сходится'); end
 
x=zeros(n,1);      % начальное приближение
 
for ii=1:1000;     % число итераций 1000
    x=B*x+g;       % это кстати всем известная операция из BLAS
end

e=cond(A)*eps;

e=size(A,1)*cond(A)*eps;

dltm=1e-10 
 
[szyk,szxk]=size(K);
 
xo=double(zeros(szyk,szxk));
 
for i=1:szxk
 
    St=S;
    Kt=K(:,i);
    for j=1:szyk
        St(j,j)=St(j,j)+St(j,j)*dltm/100;
    end
 
    x=St\Kt;
 
    for its=1:10
 
        St=S;
        for j=1:szyk
            St(j,j)=St(j,j)+St(j,j)*dltm*(its^2/100)/abs(x(j));
        end
 
        x=St\Kt;
 
    end
 
    w=Kt-St*x;
    x=x+St\w;
 
    xo(:,i)=double(x);
 
end

x=St\Kt

x=double(mp(St)\mp(Kt))

    x=St\Kt;
    w=Kt-St*x;
    xo=x+St\w;

[szyk,szxk]=size(K);
 
xo=double(zeros(szyk,szxk));
 
D=diag(S);
Dinv=sparse(diag(1./D));
R=S-diag(D);
 
for i=1:szxk    
    X=K(:,i);
    err=inf;
    tol=100*eps(norm(K(:,i)));
    iter=0;
    while err > tol
      iter=iter + 1;
      Xold=X;
      X=Dinv*(K(:,i)-R*X); % the guts of the iteration
      err=norm(X-Xold); % this is not really an error, but a change per iteration. when that is stable, we choose to terminate.
    end
    disp([sum(abs(K(:,i)-S*X))/szyk iter])
    xo(:,i)=X;
end

[szyk,szxk]=size(K);
 
xo=double(zeros(szyk,szxk));
 
D=diag(diag(S)); % Матрица с диагональю как у матрицы А, остальные элементы нулевые
E=eye(size(S));  % Единичная матрица размера А
 
D1=inv(D);                          % Находим обратную матрицу D
B=E-D1*S;                           % Находим вектор B
if norm(B)>=1
    disp(['Метод не сходится!! Попробуй пересмотреть данные!!']) %Проверяем на сходимость: Норма B должна быть <1
end
 
for i=1:szxk    
 
    Xin=K(:,i);                         % Берем за первое приближение - вектор правой части 
    err=inf;
    tol=100*eps(norm(K(:,i)));
    n=1;                                % n - нумерация итерации
 
    while n<100 && err > tol            % Условие зацикливания
        Xout=B*Xin+D1*K(:,i);           % Делаем итерации
        err=norm(Xin-Xout);             % Условие окончания
        n=n+1;                          % Делаем следующую итерацию
        Xin=Xout;                       % Подготавливаем вектор к итерации
    end
    disp([sum(abs(K(:,i)-S*X))/szyk n])
    xo(:,i)=X;
end