I am research associate professor at the University of New Mexico Center for High Technology Materials
Watch this Author's files
a code with great accuracy and speed! Well done!
A = flipud(fliplr(A)) solves my problem.
How does the code change if reverse numbering is used for the CGL points: x = -cos(k*pi/N)?