# 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