Les fonctions OpenGL prennent des degrés en parametre, par contre, si tu veux précalculer des données avant, et utilises les fonctions cos et sin de math.h, la il faut considérer des radians.
ta formule peut etre un peu compactée (le traitement particulier de x[0] et y[0] au début peuvent etre injectés dans le for :
ds = N/(2*PI);
for (s=0;s<2*PI;s+=ds)
{
x[i]= x c + R cos(s);
y[i]= y c + R sin(s);
}
travaille bien avec des float ou des double : le piege serait que tu aies mis "s" ou "ds" comme des int par exemple.
Sinon, essaie de dumper tes points calculés, et des les tracer sur une feuille, juste pour voir ce que tu as. Mais le calcul est correct : c´est la définition paramétrique d´un cercle.