def surface_area_pyramid(l, height):
    """
    A = l * √(l² + 4 * h²) + l² 

    l is a side length of the square base and h is a height of a pyramid.
    """
    return round(l * (l * l + 4 * (height * height)) ** 0.5 + l * l, 3)
