def calculate_integration_error(
        max_forth_derivative: float, low: float, high: float, length: int) -> float:
    """Find integration error if Simpson method is used.

    Args:
        max_forth_derivative: precalculated value.
        low: left range boundary.
        high: right range boundary.
        length: number of values in range.
    Returns:
        float: approximate integration error.

    """
    step = (high - low) / length

    return (step)**4 / 180 * max_forth_derivative * (high - low)