. Next, the light curve!: In PYTHON

Create a program which will plot the uneclipsed area of the star as a function of time - just like our analytical solution for the square star/planet - starting at the time between when the planet first starts to overlap the star to when the whole planet has passed in front of the star.

You will need variables for the velocity of the planet, the side lengths of the star and planet. For the velocity of the planet use an arbitrary number such as 0.001 (in Jupiters/second). For now, dont worry about a realistic velocity of the planet.

You will need arrays for the area and the time that has passed. Hint: the time array can be initialized with a list of times like earlier in this handout and the area array can be empty.

You will need to loop over time, and use logic statements to determine what expression from the analytical solution you will need to use for the uneclipsed area. Hint: you can then add the area at each time to the area array with area.append(uneclipsedArea)

Finally, you will need to plot the data you just generated! If the brightness is proportional to the area, what would you use to plot the transit lightcurve shape from previous lessons?