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
