% Call dftuvTN
[U, V]=dftuvTN(8,5);
D=U.^2+V.^2