static void makect( int nc, int* ip, double* c ) {
int j, nch;
double delta;
ip[1] = nc;
if ( nc > 1 ) {
nch = nc >> 1;
delta = atan( 1.0 ) / nch;
c[0] = cos( delta * nch );
c[nch] = 0.5 * c[0];
for ( j = 1; j < nch; j++ ) {
c[j] = 0.5 * cos( delta * j );
c[nc - j] = 0.5 * sin( delta * j );
}
}
}