Use of interval slopes in implementing an interval method for global non-linear DC circuit analysis