ありゃ、元の式は　違うのか、　上の
　ax = cx + sin(i * 2 * PI / N) * (W / 2);
　ay = cy + cos(i * 2 * PI / N) * (H / 2);
に訂正、

続き、

w = ((c-aL[k]) /(aL[k+1]-aL[k]) + k)*2*PI/N;
int x = cx + sin(w) * (W / 2);
int y = cy + cos(w) * (H / 2);
plot(x, y);
}