messed around with jacobi elliptic fncs to get a''/a numerically using finite diff approx