Mwps(xx,varargin)

This function computes the multitaper power spectrum, and optionally can return the degrees of freedom, the F-test for deterministic components, the reshaped power spectrum after line removal and the power contained in each line. The arguments have default values that can be changed by the user. Computation of the power at each frequency utilizes parfor for speed with large data sets, but can be changed to for if the Parallel Computing Toolbox is not available.