integrate (BesselJ(1,16981.9*sin(x)))^2 /sin(x) dx from 0 to 0.0002257