c was just cos(x) where x is the desired angle Arccos(-1/3)=109.471220634 degrees
S(x,y)=(sin(x)cos(y),sin(x)sin(y),cos(x))
is the parametric equation for a unit sphere
I let P1=S(0,0)=(0,0,1)
next we want P2,P3,and P4 to form an equilateral triangle
since each is also equidistant to P1 each will have the same 3rd coordinate (and first angle x)
Thus they lie on the circle
C(y)=S(x,y)
since they are equally spaced and 360/3=120 we can have
P2=S(x,y)
P3=S(x,y+120)
P4=S(x,y-120)
and since no generality is lost we can let y=0
P2=S(x,0)=(sin(x),0,cos(x))
P3=S(x,120)=(sin(x)cos(120),sin(x)sin(120),cos(x))
P4=S(x,-120)=(sin(x)cos(120),-sin(x)sin(120),cos(x))
now we need to solve for x it could be solved using distance, but dot products determine distance and are neater so
P1.P2=cos(x)
=c
P3.P4={(sin(x))^2}[(cos(120))^2-(sin(120)^2]+(cos(x))^2
={1-(cos(x))^2}[(cos(240)^2]+(cos(x))^2
={1-(cos(x))^2}[-.5]+(cos(x))^2
=1.5(cos(x))^2-.5
=1.5c^2-.5
we solve
c=1.5c^2-.5
3c^2-2c-1=0
(3c+1)(c-1)
thus
c=-1/3
(we discard c=1)
thus
x=Arccos(-1/3)=109.471220634 degrees
Another neat thing
P1+P2+P3+P4=0
(since they are equidistant to 0)
so
0=(P1+P2+P3+P4).P1=P1.P1+P2.P2+P3.P3+P4.P1=1+3c
c=-1/3
This can be generalized but there are complications
One thing, the requirement that each point by equidistant from all others is impossible when n>4
we might however only require each point to be equidistant to some of the others
or require the inverse square forces to balance
then the n=5 case is an easy direct generalization
P1=(0,0,1)
P2=(sin(x),0,cos(x))
360/4=90
P3=(sin(x)cos(90),sin(x)sin(90),cos(x))=(0,sin(x),cos(x))
P4=(sin(x)cos(180),sin(x)sin(180),cos(x))=(-sin(x),0,cos(x))
P5=(sin(x)cos(270),sin(x)sin(270),cos(x))=(0,-sin(x),cos(x))
(algebra)
x=Arcos(-1/4)=104.477512186 Degrees
or
0=(P1+P2+P3+P4+P5).P1=P1.P1+P2.P2+P3.P3+P4.P1+P1.P5=1+4c
c=-1/4
The n=6 case is easy but different
P1=(0,0,1)
P2=(1,0,0)
P3=(0,1,0)
P4=(-1,0,0)
P5=(0,-1,0)
P6=(0,0,-1)
x=90