function [p,delta] = partition(integerN, a, b)
## Partition the interval [a,b] into integerN subintervals. We do not
function [p,delta] = partition(integerN, a, b)
## Partition the interval [a,b] into integerN subintervals. We do not
##
##
## We don't use abs() here because `b` might be less than `a`. In that
## case, we want delta negative so that when we add it to `a`, we move
## towards `b`.
##
##
## We don't use abs() here because `b` might be less than `a`. In that
## case, we want delta negative so that when we add it to `a`, we move
## towards `b`.