Advertisement
Guest User

Untitled

a guest
Mar 30th, 2017
66
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 0.79 KB | None | 0 0
  1. A = [10 2 -1
  2.     -2 -6 -1
  3.      1 -3 12];
  4. x = [0 0 0]';
  5. b = [5 24.42 36];
  6.  
  7. G = zeros(3, 3);
  8. f = zeros(1,3);
  9.  
  10. for i = 1:1:length(b)
  11.     for j = 1:1:length(b)
  12.         if i ~= j
  13.             G(i,j) = -A(i,j) / A(i,i);
  14.         end
  15.     end
  16.     f(i) = b(i)/ A(i,i);
  17. end
  18.  
  19. for iter = 1:1:10
  20.     xTemp = x;
  21.     for eq = 1:1:length(b)
  22.         x(eq) = G(eq,:) * x + f(eq);
  23.         %% change "* x" to "* xTemp" to get simple iteration
  24.         %% otherwise it is working as Zeidel Method
  25.     end
  26. end
  27.  
  28. x;
  29. A * x;
  30. A = A * A
  31. L = zeros(3,3);
  32. for i = 1:1:length(A)
  33.     L(i,i) = sqrt(A(i,i));
  34.     for k = i+1:1:length(A)
  35.         L(k,i) = A(k,i) / L(i,i)
  36.     end
  37.     for j = i+1:1:length(A)
  38.         for k = i:1:length(A)
  39.             L(k,j) = A(k, j) - A(k,i) * A(j,i);
  40.         end
  41.     end
  42. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement