Always hand in:
- written solutions to any questions 
- a paper print-out of well-commented code.  Include a multiline comment at
the top of your code with (i) the assignment name, (ii) your name, 
and (iii) the date you handed in all elements of the assignment
- paper print-out of output (graph or text)
- also, e-mail me (y-lithwick@northwestern.edu) the code with your name and
the exercise number in the subject line
Assignment #8  
 [10 pts, due 2pm, April 26]
  
	- 
	 Reproduce figures
	 3.6 (left panel)
	 and 
	 3.8 (both panels) in the text.
    For Fig 3.6, approximately how high does FD
    have to be to give chaos? Show plots with FD 
    above and
    below that critical value. Choose a parameter
    (e.g., q or FD), 
    and explore what happens when you change it.  
    Show me a few plots and explain briefly what
    is going on.
	
- 
	Reproduce figure 3.7 (both panels) in 
	the text.  Make a few more such plots with
	different initial conditions (but the same
	parameters).
	Estimate the Lyapunov times for the two
	cases, showing 
	  how you arrive at your answers.
	
- 
	 Make a movie of the simulation used for the right
	 panel of figure 3.7.  The movie should show
	 both pendula simultaneously, on the same 
	 pivot.   To receive credit for this question, 
	 show me the movie at the beginning of class
	 on the day this assignment is due (or earlier).
	
-  
	Poincaré sections:
	reproduce figure 3.9 in the text.
	Also make a Poincaré section 
	of the trajectory in 
	the left-hand panel of Figure 3.8, after
	the initial transient has died away.  
	Explain in both cases what the sections mean about
	the motion of the pendulum.
	How does the Poincaré section change
	if you change the initial condition? 
	
-  You need  only do one of the following.
	(Though I encourage you to do more, and/or follow your
	own explorations or the book's Section 3.4.)
	You will be presenting your results to the class
	on Thursday.  Find some way to get a copy of your
	plots or movie to the projector, e.g., via the computer
	in the lab that accesses the projector directly.
	  
	  -  Frank & Xiaowen: 
	  	    What happens to  the Poincaré section
	    of question 4 for different values of FD?
	    Focus especially on what happens 
	     near  transitions from regular to chaotic
	    behavior. Can you find multiple transitions?
	  
-  Thomas & Max:
	    Explore the Poincaré section
	    of question 4 further by examining
	    plots for different values q.
	
-  Adam: 
	  Run simulations so that you can zoom into portions of Fig. 3.9 and show
	  what happens at very high resolution.  
	  You should see the lines breaking
	  up into thinner lines.
	  You will have to run the code for a long
	  time to get good enough resolution. 
	
-  [If anyone prefers they may do this one
	instead; it is a lot easier than it sounds:]
	Make a movie of the double pendulum. 
	The equations of motion are 
	
   here.