def _get_eps_range(d0, d1, solution, sol_range, eps):
    """ Ensure that we don't plot out of boundaries
    """
    eps_x_low = solution[d0] - eps
    eps_x_high = solution[d0] + eps
    eps_y_low = solution[d1] - eps
    eps_y_high = solution[d1] + eps
    if eps_x_low < sol_range[0][d0]:
        eps_x_low = sol_range[0][d0]
    if eps_x_high > sol_range[1][d0]:
        eps_x_high = sol_range[1][d0]
    if eps_y_low < sol_range[0][d1]:
        eps_y_low = sol_range[0][d1]
    if eps_y_high > sol_range[1][d1]:
        eps_y_high = sol_range[1][d1]
    return [eps_x_low, eps_x_high], [eps_y_low, eps_y_high]
