R(0) = (1/2)R(0)/(N-1) + ΣR(j)/(N-1) + (1/2)R(N-1)/(N-1) (Σ for j = 1 to N-2)
R(k) = R(0)/(N-1) + Σ2x[j)cos(πjk/(N-1)) + (-1)^k R(N-1)/(N-1) (Σ for j = 1 to N-2) (k = 1 to N-2)
R(N-1) = (1/2)R(0)/(N-1) + ΣR(j)(-1)^k/(N-1) + (1/2)(-1)^(N-1) R(N-1)/(N-1) (Σ for j = 1 to N-2)
This transform is normalized since a call to WCost1f followed by a call to WCost1b (or vice-versa) reproduces the original array subject to algorithmic constraints, roundoff error, etc.