Advertisement
Guest User

Untitled

a guest
Apr 24th, 2017
64
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Octave 0.70 KB | None | 0 0
  1.  %x = 2 4 -1 -2 -4 1
  2. A=[ 2 1 -2 3 2 6; % = 2 = 4+4+2-6-8+6
  3.     4 1 2 4 0 1;  % = 3
  4.     -2 0 1 -7 6 1;  % = -14
  5.     3 5 7 0 3 2;  % = 9
  6.     3 7 1 1 4 3;  % = 18
  7.     2 -4 1 -2 3 0 ];  % = -21
  8. Y=[2; 3; -14; 9; 18; -21];
  9. A
  10. det(A)
  11.  [Aw, Ak]=size(A);
  12.  [Yw, Yk]=size(Y);
  13.  X=zeros(Aw,1);
  14.  if Aw==Ak && Aw==Yw
  15.  detA=det(A);
  16.  [lw,lk]=size(A);
  17.  if detA!=0
  18.   for i=1:size(Y)
  19.     Yh=zeros(lw,lk);
  20.     for q=1:lw
  21.       for w=1:lk
  22.         if w==i
  23.           Yh(q,w)=Y(q,1);
  24.         else
  25.           Yh(q,w)=A(q,w);
  26.         endif
  27.       end
  28.     end
  29.     X(i,1)=det(Yh)/det(A);
  30.   end
  31.   disp("Rozwiazania: "); X
  32.  else
  33.  disp("Uklad nieoznaczony (wyznacznik rowny 0)");
  34.  end
  35. else
  36.  disp(Blad macierzy");
  37.  end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement