This command creates a sequence of arrays whose elements are the subpixels required by the Gauss - Legendre formulae to achieve a given precision in computing the integral of the PSF. The precision depends on the user who selects which fraction of the intrinsic noise is acceptable in numerically computing the integral. This data are stored in a file, IW.DAT, and passed to the command for fitting.