def surface_area_triangular_prism(a, b, c, height):
    """
    A = 0.5 * √((a + b + c) * (-a + b + c) * (a - b + c) * (a + b - c)) + h * (a + b + c)
    """
    return round(0.5 * ((a + b + c) * (-a + b + c) * (a - b + c) * (a + b -
        c)) ** 0.5 + height * (a + b + c), 3)
