applies a Detrend correction by orthogonal projection to a Vandermonde matrix
xd = detrending(x,(order));
a matrix (n x q) or a Div structure
the order of the polynom; default value=2
a Div structure corresponding to x after the Detrend correction; xd.d is of dimensions (n x q)