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