for n, (i, j) in enumerate(zip(uj_list[1:-1], uj_hat_list[1:-1])):
            alpha.append(compute_alpha(n))
            beta.append(compute_beta(i, i - 1, i + 1, j, j - 1, j + 1, n + 1))