def UTRIS(U, y):
    n = len(U)
    x = np.zeros((n, 1))

    s = 0

    for i in range(1, n):
        x[i] = y[i] * U[i][i]
        s += x[i]

    x[0] = y[0] - s

    return x