Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import numpy as np
- n = int(input())
- # arr a
- A = np.diag(np.arange(1, n + 1))
- # arr b
- re = np.r_[ (n // 2) * [1, 0] + (n % 2) * [1] ]
- ro = np.r_[ (n // 2) * [0, 1] + (n % 2) * [0] ]
- B = np.row_stack((n // 2) * (re, ro))
- if (np.shape(B)[0] != n):
- B = np.vstack((B, re))
- # arr c
- C = np.ones((n,n), dtype = 'float')
- C[1:-1, 1:-1] = 0
- C = C * 10
- # arr d
- D = np.linalg.inv(np.add(np.matmul(np.transpose(B), B), A)) + np.matmul(np.transpose(C), C)
- out = np.linalg.det(D)
- print(out)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement