def METHOD_NAME(t):
    ctoverh2 = t * consolidation_coefficient / soil_height**2
    def coeff(n, ctoverh2):
        return np.power(1.0 / (2.0 * n - 1.0), 2) * np.exp(-np.power(2.0 * n - 1.0, 2) * np.pi ** 2 * ctoverh2 / 4.0)
    result = 0.0 * t
    for i in range(1, 11):
        result += coeff(i, ctoverh2)
    return 1.0 - (8.0 / np.pi**2) * result