In celestial mechanics. Keplers equation is important. It reads x = y sin(y), in which x is a planets mean anomaly, y is its eccentric anomaly, and the eccentricity of its orbit. Taking = 0.9, construct a table of y for 30 equally spaced values of x in the interval 0 x . Use Newtons method to obtain each value of y. The y corresponding to an x can be used as the starting point fro the iteration when x is changed slightly.

Code for this in python or matlab?