Advertisement
kubawawa

Untitled

Apr 15th, 2022
851
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. clc; clear
  2. %% Dane wejsciowe
  3. A1 = [1,2;1.1,2];
  4. B1 = [10;10.4];
  5. A2 = [2,5.999999;2,6];
  6. B2 = [8.000001;8];
  7. A3 = [3,-1,-1;1,-3,-2;-1,2,1];
  8. B3 = [4;-1;-1];
  9.  
  10. MojGauss(A1,B1)
  11. MojGauss(A2,B2)
  12. MojGauss(A3,B3)
  13.  
  14. %% Wyznaczanie wskaźników uwarunkowania
  15. wsk_1=cond(A1);
  16. wsk_2=cond(A2);
  17. wsk_3=cond(A3);
  18.  
  19. %% Rozwiązywanie układów
  20. function x = MojGauss(A,B)
  21.     n = size(B,1);
  22.     x = zeros(n,1);
  23.     AB = [A B];
  24.     for i=1:n
  25.         AB(i,:)=AB(i,:)./AB(i,i);
  26.         for j=1:n
  27.             if i ~= j
  28.                 AB(j,:) = AB(j,:).*AB(j,i);
  29.             end
  30.         end
  31.     end
  32.     x = AB(:,n+1);
  33. end
  34.  
  35.  
Advertisement
RAW Paste Data Copied
Advertisement