Excercise 7

1. A baseball in free-fall

Continue the example of the lecture notes on the example for a baseball thrown straight up. We solved the equation
	(d^2) (y/d)(t^2) = a = 
          g,
with the initial conditions
	y(t = 0) = 0.

1.1. We were left with one constant undefined. Fix that constant by using the additional initial condition
	\(TraditionalForm\`d\(y(0)\)\/dt = \(\(y\& . \)(0) = 10  m/s\)\).

1.2. Find the highest point of the ball's path and the moment when the ball returns to the height it was thrown from.

2. A brick falling in water

Continue the example on the brick dropped from a boat. Out differential equation for the position y of the brick was
	\(TraditionalForm\`F = \(m\ g - k\ \(y\& . \) = m\ \(y\& .. \)\)\).

2.1. Set some sensible initial conditions for the brick's position, (i) solve the equation and give the constant some numeric values. (ii) Plot the solution. What do you find? Study the effect of changing the mass of the brick without chaning the coefficent k. (iii) Find the terminal velocity of the falling brick (hint: use the Limit-command, to tell that m and k are positive, give Assumptions->{k>0, m>0} to the command).

2.2. Solve the brick example numerically. Compare the numerical solution to the analytic one by plotting both solutions, or by comparing the values of ygiven by the two solutions.

3. Flying rocket

3.1. Solve the rocket example, and plot the solution, but this time take into account the force due to air-resistence, \(TraditionalForm\`F\_air = \(\(-k\)\ v\^2 = \(-k\)\ \(\(x\& . \)(t)\)\^2\)\). Use k = 0.3 in your calculations. Is the equation of motion still linear? Compare the path of this rocket to the one studied in the examples.

4. Frictionless and damped oscillators

4.1 Solve the frictionless
	\(TraditionalForm\`m\ \(x\& .. \) + k\ x = 0\)
and damped 
	\(TraditionalForm\`m\ \(x\& .. \) + γ\ \(x\& . \) + k\ x = 
        0\)
oscillators with the same initial conditons. Plot both of them in the same coordinates. What do you find? Does the damping affect the frequency of the oscillations?

5. LC- and LCR-circuits

5.1. Continue the example on the LCR-circuits. The general analytic solution was given in the lecture notes. Now fix the constants and study the behaviour of the system by plotting the solutions. Vary the values for the components and see what happens.

6. Competion between populations

6.1 Study the solutions of equation dN(
              t)/dt = ϵ N(
              t) - β ((N(t))^2) by varying the values of the parameters ϵ (birth-rate) and β (shadowing), and the initial populations N(t = 0). Plot the population as a function of time.

6.2 Study the behaviour of the equations describing the life of the hares and foxes and of spruces and pines by using different values for the parameters. Plot the populations of foxes as a function of the hare population, and also the spruce population as a function of the number of spruces. Using the figures, find the stable solutions where the populations stay constant in time.

Additional excercises

You can try these if you've time and interest. These contain rather advanced stuff, so you can safely ignore these.

1. Nonlinear oscillator. Consider a body attached to a spring again. If the amplitude of the oscillations is not small, the force excerted by the spring will not follow Hooke's law, (Subscript[F, spring]) = -k x. Instead, in general higher powers of x appear. Suppose the force is given by the function (Subscript[F, spring]) = -k x - ϵ (x^3). Solve the equations of motion without friction for m = 1kg, k = 0.5N/m and
	 (i) ϵ = 0
	 (ii) ϵ = 0.1N/m( ^3)
Plot the solutions of (i) and (ii) in the same figure. Study the effect of changing the initial displacement, x(t = 0).

2. Coupled oscillators. Consider two mass points of masses (Subscript[m, 1]) and (Subscript[m, 2]). The first is connected to a ceiling with a spring (Subscript[k, 1])and the second is connected to the first a spring (Subscript[k, 2]). Let (Subscript[x, 1])and (Subscript[x, 2]) be the displacements from the equilibrium positions of the two masspoints. We can write the equations of motion
	\(TraditionalForm\`m\ \(x\& .. \)\_1 = \(-k\_1\) x\_1 - \(k\_2\)(
            x\_1 - x\_2)\)
\(TraditionalForm\`m\ \(x\& .. \)\_2 = \(-\(\(k\_2\)(
            x\_2 - x\_1)\)\)\)
Solve numerically the above equations and plot the solutions. Set initial conditions to your liking and try different values of (Subscript[m, 1]), (Subscript[m, 2]), (Subscript[k, 1]) and (Subscript[k, 2]).