theta=pi*i/N;
phi=pi*i/N;
p[0]=sin(phi)*cos(theta);
p[1]=sin(phi)*sin(theta);
p[2]=cos(phi);
p=frprmn.minimize(p);
f=frprmn.fret;
d=sqrt(SQR(p[0])+SQR(p[1])+SQR(p[2]));
// System.out.println(" " << abs(d-d0) << " " << abs(f-f0));
localflag = localflag || abs(d-d0) > sbeps1;