SS_tol = 1e-13
# nvec
nvec = np.zeros(S)
nvec[:int(round(2 * S / 3))] = 1
nvec[int(round(2 * S / 3)):] = 0.2
# L
L = nvec.sum()
params = (S, beta, sigma, L, A, alpha, delta, SS_tol)
# Question 1: Feasibility
# (a)
bvec_guess1 = np.ones(S - 1)
print("(a) : ", func.feasible(S, alpha, A, delta, bvec_guess1))
# (b)
bvec_guess2 = np.array([
-0.01, 0.1, 0.2, 0.23, 0.25, 0.23, 0.2, 0.1, -0.01, 0.1, 0.2, 0.23, 0.25,
0.23, 0.2, 0.1, -0.01, 0.1, 0.2, 0.23, 0.25, 0.23, 0.2, 0.1, -0.01, 0.1,
0.2, 0.23, 0.25, 0.23, 0.2, 0.1, -0.01, 0.1, 0.2, 0.23, 0.25, 0.23, 0.2,
0.1, -0.01, 0.1, 0.2, 0.23, 0.25, 0.23, 0.2, 0.1, -0.01, 0.1, 0.2, 0.23,
0.25, 0.23, 0.2, 0.1, -0.01, 0.1, 0.2, 0.23, 0.25, 0.23, 0.2, 0.1, -0.01,
0.1, 0.2, 0.23, 0.25, 0.23, 0.2, 0.1, -0.01, 0.1, 0.2, 0.23, 0.25, 0.23,
0.2
])
print("(b) : ", func.feasible(S, alpha, A, delta, bvec_guess2))
# (c)
bvec_guess3 = np.array([
