# Solve LP 
sol1 <- rpsg_solver(problem_list = problem_list)

# Extract solution 
vec_c1 <- sol1[[“point_problem_1"]]

# Compute weights
vec_weights1 <- vec_c1 * problem_list$matrix_objcoef[1,] / sum(vec_c1 * problem_list$matrix_objcoef[1,])

# Find nonzero solutions
vec_c_nonz1 <- vec_c1[which(vec_c1 !=0)]

# Find nonzero weights
vec_weights_nonz1 <- vec_weights1[which(vec_weights1 !=0)]