a, b = p_opt - p_crit_low, p_opt + p_crit_up