Sept. 14, 2004 --- Class 5 ---Newton's Laws of Motion, Energy Nonconservation 

Activities:

	Solving Newton's Law of Motion
		We started Chapter 3 and discussed how to turn a 2nd order
		differential equation into a series of coupled first order
		differential equations.  We also discussed why position,
		velocity and acceleration are so important, and the next
		derivative, called the jerk, is not well known.

		The program fall1.c in ~sg contains code to integrate the
		equations for a falling body.  It is set up for the
		Euler algorithm.  We discussed the code in detail.

		If the lines for updating v and y are
		reversed, you get the Euler-Cromer algorithm.  In the next
		class we will discuss the Euler-Richardson algorithm, and 
		explore the step size dependence of the trajectory of a
		falling body.  We talked about energy conservation and
		how to write an awk expression to calculate the energy.  

		We studied the energy of a falling body as caculated by the
		Euler algorithm.  We looked at the energy at t = 1 sec.
		We also created a shell script to automate the procedure:

#!/bin/csh
# this is the best shell script I have ever written
foreach step ( 0.025 0.05 0.1 0.2)
fall1 <<EOF | awk '{if($1 == 1.0) print '$step',0.5*$3*$3-9.8*$2}'
    $step
    .2
EOF
end

		Notes: When a command expects its input from the standard input
		but you want to run it from a shell script, the << sign is
		used to redirect the standard input so that it comes from the
		shell script.  To signal the end of the redirected input,
		we need a string.  By convention, I always use EOF, which is
		short for End Of File.  Note that EOF must occur right at the
		beginning of the line.  Don't indent it!
		In this case both the input and output of the fall command are
		redirected.  The output is being piped to an awk command that
		prints the step size and the energy if the first entry on the
		line is 1.0.  Note use of == to test for equality.  This is
		is test of numerical equality.  It is not a test of string 
		equivalence.  Also note that there are now two sets of
		matching single quotes in the awk command.  This means that
		the shell variable $step is not within the single quote.
		Thus, the shell does variable substition and replaces $step
		by its current value determined by the foreach loop.  If you
		put in any extra spaces around $step, the awk command will
		not work because awk takes its command from a single
		string, and the extra space breaks the string that appears
		above into two strings.  These things can be tricky, so if
		you have problems consult me.  I have said it before, but I
		will repeat myself.  The reason I like to meet in a computer 
		cluster is because I know from years of experience that very 
		simple typing changes on your part, can make a big difference.  
		I can sometimes save you hours of pondering what is wrong 
		by spotting a very tiny error.  Don't hesistate to email me 
		if you are having trouble outside of class.

		Getting back to the numerics, we found that for fixed time, the
		energy is a linear function of the step size.  Of course, the
		energy should be constant at 0.  We also found for a single
		step size that the energy is linearly increasing with time.
		We would like to find a better algorithm.  We will do that in
		the next class.