Recommended reading from Statistics, Data Mining, and Machine Learning
in Astronomy ("SDMML") and Numerical Recipes ("NR"):
NR: Ch 9, 10
Find the lowest positive solution to (x-1) = a*np.cos(x)**2 for a in the range of 1-6 tan(x) = x/(1-a*x^2) using Newton-Raphson