Introduction - If you have any usage issues, please Google them yourself
Inputs:
ts = ts(i)= i,
where i=1, ..., ni
HiLo = 2-character string indicating rejection of high or low outliers
low = valid range minimum
high = valid range maximum (values outside the valid range are rejeced
right away)
fet = fit error tolerance (points deviating more than fet from curve
fit are rejected)
dod = degree of overdeterminedness (iteration stops if number of
points reaches the minimum required for curve fitting, plus
dod). This is a safety measure
delta = small positive number (e.g. 0.1) to suppress high amplitudes
Outputs:
amp = returned array of amplitudes, first element is the average of
the curve
phi = returned array of phases, first element is zero
yr = array holding reconstructed time series