function [p,delta] = partition(integerN, a, b) % Partition the interval [a,b] into integerN subintervals. We do not % require that a