Function std::dct2

int dct2([]const double in, []double out)

X[k] = sum_{n=0}^{N-1} x[n] * cos(pi/N * (n+0.5) * k), k = 0..N-1.