I have the following quantity:
c_p = 1 - (4*sin(th)^2+2*5*sin(th)/pi+(5/(2*pi))^2)
and I am trying to find its roots (the solutions for c_p==0
).
When plotting c_p
as a function of th
between −π and π, I can see the curve crosses the x-axis at four positions. However, solve(c_p==0, th)
is only giving two roots:
[theta=−arcsin(54 π+12),theta=−arcsin(54 π−12)].
It appears that solve
can only find the roots that are in the domain of the arcsin function, i.e in the interval [−π,π]. Is there a way to get the other two roots?