+function [p] = partition_delta(a, b, delta)
+ ## Partition the interval [a,b] into subintervals of length delta. We
+ ## do not require that a<b.
+ ##
+ ## INPUTS:
+ ##
+ ## * ``a`` - The "left" endpoint of the interval to partition.
+ ##
+ ## * ``b`` - The "right" endpoint of the interval to partition.
+ ##
+ ## * ``delta`` - The length of the subintervals.
+ ##
+ ## OUTPUTS:
+ ##
+ ## * ``p`` - The resulting partition, as a vector of length (b-a)/delta.
+ ##
+
+ p = [a : delta : b];
+end