Решение СЛАУ методом Крамера

 case 7 %Решение СЛАУ
           disp ('Выполняется решение СЛАУ методом Крамера')
          
           disp('Заполните данные: ')
           
           n=input('Количество уравнений: ');
           m=input('Количество количество переменных : ');
           A = Yinput('A',n,m);
           C = Yinput('C',1,m)';
          
           Youtput('A',A)
           disp(' ')
           if m~=n ;
            disp(' Ошибка, матрица не квадратная ');
           else
               
            M=Determinant(A);
            disp(['det(A) = ' num2str(M)])
            if (M==0)
                   disp('Error! det(A) = 0 , несовместна')
            else
 
            Y=Equation([A C]);
            Youtput('Y',Y)
               disp(' Проверка ')
            R=Multiplication(A,Y,n,m,1);
            Youtput('A*Y',R)
               
           end
           end  
           
           end

function d = Determinant (M) %определитель
 
[s,r]= size(M);
if s==1
    d = M(1,1);
elseif s==2
    d = M(1,1)*M(2,2)-M(1,2)*M(2,1);
else
    d=0;
    for i=1:s
        %         B=M;
        %         B(1,: )=[];
        %         B(:,i)=[];
        for j=2:s 
            h=1;
            for k=1:r
               if k~=i 
                   B(j-1,h)=M(j,k);
                   h=h+1;    
               end;               
            end
        end
        d=d+Determinant(B)*M(1,i)*(-1)^(i+1);
    end
end
end

function C=Multiplication(A,B,n,m,l)%умножение
for j=1:l
    for k=1:n 
        C(k,j)=0;
        for p=1:m
            C(k,j)=C(k,j)+A(k,p)*B(p,j);            
        end
    end
end
end

function Y=Equation(A)%уравнение
% B=A(:,size(A,2));
Y=[];
[r,c]=size(A);
for i=1:r
    B(i)=A(i,c);
end
 
d = Determinant(A);
for k=1:size(A,1)
    X=A;
  %     X(:,k)=B;
    for i=1:r
        X(i,k)=B(i);
    end
        dn=Determinant(X);
        Z=(dn/d);
        Y(k,1)=Z;
end
end

function [A] = Yinput(Aname,n,m)%ввод
disp([' Введите матрицу ' Aname])
 for k = 1:n
    for p = 1:m
    A(k,p) = input([Aname '(' int2str(k) ',' int2str(p) ')=' ]);
    end
end
end

function Youtput(Aname,A) %вывод
 disp(['Матрица ' Aname ])
 disp(num2str(A,'%8.2f'));