brospresident

Untitled

Oct 20th, 2021
664
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. import numpy as np
  2.  
  3. rng = np.random.default_rng()
  4.  
  5. def S_INF_TR(L: np.ndarray, b: np.ndarray):
  6.     n = len(L)
  7.  
  8.     rezultat = np.zeros(n, dtype = int)
  9.  
  10.     for i in range(n):
  11.         s = b[i]
  12.         for j in range(i - 1):
  13.             s -= L[i][j] * rezultat[j]
  14.         rezultat[i] = s / L[i][i]
  15.  
  16.     return rezultat
  17.  
  18. def S_SUP_TR(U: np.ndarray, b: np.ndarray):
  19.     n = len(U)
  20.  
  21.     rezultat = np.zeros(n, dtype = int)
  22.  
  23.     for i in range(n - 1, -1, -1):
  24.         s = b[i]
  25.         for j in range(i + 1, n):
  26.             s -= U[i][j] * rezultat[j]
  27.         rezultat[i] = s / U[i][i]
  28.  
  29.     return rezultat
  30.        
  31. L = rng.integers(1, 20, size = (7, 7))
  32. L = np.tril(L)
  33. b = rng.integers(1, 20, size = 7)
  34. print("Rezultat S_INF_TR pentru matrice 7x7 si vector de 7 elemente:")
  35. print(S_INF_TR(L, b))
  36.  
  37. L = rng.integers(1, 100, size = (25, 25))
  38. L = np.tril(L)
  39. b = rng.integers(1, 100, size = 25)
  40. print("Rezultat S_INF_TR pentru matrice 25x25 si vector de 25 elemente:")
  41. print(S_INF_TR(L, b))
  42.  
  43. U = rng.integers(1, 20, size = (7, 7))
  44. U = np.triu(U)
  45. b = rng.integers(1, 20, size = 7)
  46. print("Rezultat S_SUP_TR pentru matrice 7x7 si vector de 7 elemente:")
  47. print(S_SUP_TR(U, b))
  48.  
  49. U = rng.integers(1, 100, size = (25, 25))
  50. U = np.triu(U)
  51. b = rng.integers(1, 100, size = 25)
  52. print("Rezultat S_SUP_TR pentru matrice 25x25 si vector de 25 elemente:")
  53. print(S_SUP_TR(U, b))
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
RAW Paste Data