The following is an equation describing the coupled phases of N oscillators according to the Kuramoto model:

$$
1 = K \int_{-\pi/2}^{\pi/2}\cos^{2}\left(\theta\right)\,
{\rm g}\left(KR\sin\left(\theta\right)\right)\,{\rm d}\theta
$$

We assume that g is a symmetric distribution with zero mean, since the equation is written from the point of view of a reference frame rotating together with the oscillator.

I need to take the $\lim_{R \rightarrow0^+}$ to solve for the coupling strength Kc.
This can't be solved with Wolfram Alpha since g is not read as a gaussian distribution but instead as a constant.

How else can you solve this limit?