p = 0.5 Pi;
b = 3;
g = 0.5;
w = 0;
pt1 = NSolve[{
Cos[p/2] Sin[k] + (1 - 2 c1^2) Cos[k] Sin[p/2] ==
0, (-1 + 2 c1^2) (b + c1 Sqrt[1 - c1^2] g) -
4 c1 Sin[k] Sin[p/2] Sqrt[1 - c1^2] == 0}, {k, c1}]
pt2 = NSolve[{
Cos[p/2] Sin[k] + (1 - 2 c1^2) Cos[k] Sin[p/2] ==
0, (-1 + 2 c1^2) (b + c1 Sqrt[1 - c1^2] g) -
4 c1 Sin[k] Sin[p/2] Sqrt[1 - c1^2] == 0}, {k, c1}, Reals]

自顶一下，想把实数解取出来要怎么写程序？之后要循环b。