Advertisement
HigerSkill

Untitled

Sep 22nd, 2017
61
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.71 KB | None | 0 0
  1. def rotate():
  2.     i, j = 0, 0
  3.     c, s = 0, 0
  4.  
  5.     for i in range(0, n):
  6.         for j in range(i + 1, n):
  7.             a = matrix[i, i]
  8.             b = matrix[j, i]
  9.             c = a / m.sqrt(a*a + b*b)
  10.             s = b / m.sqrt(a*a + b*b)
  11.             for k in range(i, n+1):
  12.                 t = matrix[i, k]
  13.                 matrix[i, k] = (c * matrix[i, k]) + (s * matrix[j, k])
  14.                 matrix[j, k] = (-s * t) + (c * matrix[j, k])
  15.  
  16.     for i in range(n-1, 0, -1):
  17.         summ = 0
  18.         for j in range(i + 1, n):
  19.             summ += matrix[i, j] * x[j]
  20.         summ = matrix[i, n] - summ
  21.         if matrix[i, i] == 0:
  22.             return False
  23.             exit()
  24.         x[i] = summ / matrix[i, i]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement