plotted small angle approx for period, noticed different from ODE solution, trying to fix