ire the size of the data
      to be a power of 2.
    </p>

    <p>
      This version of the program is implemented as a MATLAB function.
      The Fourier transform functions that are invoked can handle input
      data of any size.
    </p>

    <p>
      Versions of the algorithm are implemented which generate a vector x
      of size N with a 1/f^alpha frequency distribution, with three choices
      for the underlying distribution of the white noise vector:
      <ul>
        <li>
          <b>f_alph