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(\theta)g(KRsin(\theta))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?