Homework Assignment # 2: MATH 235 - Mathematical Models in Science and Engineering
Homework Assignment # 2: MATH 235 - Mathematical Models in Science and Engineering
Instructions: Present your work in a neat and organized manner. Please use either the 8.5 × 11 size paper or
the filler paper with pre-punched holes. Please do not use paper which has been torn from a spiral notebook.
Please secure all your papers by using either a staple or a paper clip, but not by folding its (upper left) corner.
You must show all of the essential details of your work to get full credit. If you used Mathematica for
some of your calculations, attach a printout showing your commands and their output. If I am forced to fill in
gaps in your solution by using notrivial (at my discretion) steps, I will also be forced to reduce your score.
Please refer to the syllabus for the instructions on working on homework assignments with other students
and on submitting your own work.
Homework Assignment # 2
2. Show that Eq. (2.11) reduces to the Cartesian equations in the standard form for:
• an ellipse, when 0 < e < 1 (see also the Note (not the Hint!) below);
• a parabola, when e = 1; and
• a hyperbola, when e > 1.
Hint: Rewrite (2.11) as r = eD − er cos θ, square both sides, and use the relation between the polar
and Cartesian coordinates. Then, complete the square and thereby put your equations in the form of the
canonical equations of the ellipse and hyperbola:
(x − x0 )2 (y − y0 )2
± = 1, (HW2.1)
a2 b2
where ‘+’ and ‘−’ pertain to the ellipse and hyperbola, respectively. In the case of an ellipse, a and b are
its semi-axes.
1
For the ellipse only, you also need to do the following:
(i) Find its geometric center, i.e. the values (x0 , y0 ) in the equation above.
(ii) Find its major and minor axes (i.e., a and b) in terms of e and D.
(iii) The case of e = 0 yields a circle, because you will find that in this case a = b. However, if you
simply set e = 0 in your equation, you will obtain a circle of radius 0. To obtain a circle of a nonzero
radius, you should let e → 0 but eD → C, where C is some nonzero constant. How is this C related to
the radius of the circle? Where is the center of this circle?
You will use these results in the next problem.
Note: As we have already mentioned, in the ellipse described by Eq. (HW2.1), the geometric center
is at (x0 , y0 ). Then the origin in that equation is at one of the foci of the ellipse. In the next problem,
we will refer to this curve as the original ellipse. Conversely, if one shifts this ellipse so that now its
geometric center would coincide with the origin, then the foci of this shifted ellipse will be at (x0 , y0 )
and (−x0 , −y0 ). (Make a sketch to convince yourself that this is true.) In the next problem, we will refer
to this curve as the shifted ellipse.
Thus, to summarize: The original ellipse has is focus at the origin, while the shifted ellipse has its
geometric center at the origin.
2
Note about notations: In (2.20), N is the proportionality coefficient in Newton’s Law of Gravitation
(2.12). Thus, (2.20) supplies a relation between parameters N and c ≡ |⃗c| (defined in (2.14) and earlier
in (2.1)) one one hand and the ellipse’s parameters e and D (see (2.11)).
Technical note: You need to tell Mathematica that 0 < e < 1 (otherwise the orbit is not an ellipse and
hence the integral could not be found). The way to do it is by using the option Assumptions, as shown
below:
Integrate[f[theta],{theta,0,2*Pi},Assumptions -> {e>0,e<1}],
where f[theta] is the function you want to integrate.
2) Note that according to Eq. (2.1), Ȧ = c/2. Since by Kepler’s Second Law, the planet’s radius-
vector sweeps out equal areas in equal intervals of time, the period of the planet is found simply as
T = A/Ȧ. Find this T .
3) In Problem 2 you should have expressed the major axis, 2a, of an ellipse via the parameters D and
e. State your answer here. (If you could not do this part of Problem 2, ask me or Google for an answer.)
4) Combine the results of steps 1) – 3) to obtain Eq. (2.21) of the notes. Make sure to demonstrate
that the constant on the right-hand side (r.h.s.) does not depend on the eccentricity of the orbit and on the
revolution rate c. (This is what makes it a “constant”, i.e., the same for all planets, in the Third law.)
5. In the Lecture, we showed that given the Universal Law of Gravitation, Eq. (2.12), one can find the
equation of motion r = r(θ) (Eq. (2.20)).1 Here we address the inverse problem: Given the equation
⃗ = f (r)⃗r/r? This inverse problem was also
r = r(θ) of the orbit, how can we find the central force F
solved by Newton in his treatise “Principia Mathematica”.
Here you will be guided towards the solution of the inverse problem following the ideas of Problems 8, 9,
11 – 13 on p. 74 of the book “Inverse Problems” by C.W. Groetsch; these problems are posted alongside
this HW assignment.
1) Take the cross product of Eq. (2.15) with ⃗r and then use Eq. (2.7) to show that
( )
r2⃗¨r = ⃗¨r · ⃗r ⃗r .
Use the idea expressed by the first equation in (2.16) to transform this to
(( )· ( ))
r2⃗¨r = ⃗r˙ · ⃗r − ⃗r˙ · ⃗r˙ ⃗r .
Obtain the expression for the first term on the r.h.s. using the last equation in (2.9). For the second term,
recognize that ⃗r˙ = ⟨ẋ, ẏ⟩ and then use Eqs. (2.3). Combine these results to obtain the formula of
Problem 9 on p. 74 of Groetsch’s book.
Note that this yields the equation of motion in the form
⃗r
⃗¨r = f (r, r̈, θ̇) .
r
Our next goal is to eliminate the time derivatives from f (r, r̈, θ̇) and express it as a function of the sole
variable r. You will do this in the next two steps.
First, in Step 2), you will express θ̇ as a function of r and then express r̈ as a function of θ. Now, since
r is given as a function of θ, one can, at least in principle, express θ as a function of r and substitute the
result into the expression for f ( r, r̈(θ), θ̇(r) ) to obtain f as a function of r only.
1 ⃗ = ⃗r/rn with any n, not just n = 3.
In fact, it is possible to find r = r(θ) for a force of the form F
3
2) Let us recall from Eq. (2.1) that r2 θ̇ = c = const. Hence θ̇ = c/r2 .
To continue
( )with finding r̈, note that on an orbit traced by a planet, r = r(θ), where θ = θ(t). That is,
r = r θ(t) . Then use the Chain Rule to show that
dr c
ṙ = .
dθ r2
Define u(θ) ≡ 1/r(θ) and rewrite only the r.h.s. of the above equation in terms of u and/or du/dθ.
To finally obtain an expression for r̈ = (ṙ)· , use the Chain Rule again. Now obtain the result of Problem
12 on p. 74 of the book. How is h related to c?
3) Combine the results of 1) and 2) to obtain the formula of Problem 13 on p. 75. Note that now you
have found f (r, r̈, θ̇) in the form f (u(θ), d2 u(θ)/dθ2 ), where u = 1/r. Examples of using this result
to obtain specific expressions for f in terms of r only, are considered in the next problem.
7. Create a folder named something like math235 in the drive where you have Matlab installed (your
share of the M-drive if you use one of the UVM computers and C-drive if you have Matlab installed
on your personal computer). Download files hw2_p7.m and hw2_p7_funorbit.m from the
homework website into this folder.
Vary only two parameters: tf and vy0, in hw2_p7.m and run it to obtain an elliptical and a hyperbolic
orbits. Also, try to get as close to a parabolic orbit as you can.
Hint: Before doing this problem, take a piece of paper and make a sketch showing the location of the
Sun and the initial location of the planet, as given in lines 11 and 12 of the code hw2_p7.m. Now think
in what direction the planet should move at the initial moment to make a closed orbit (i.e., an ellipse)
possible. What about a non-closed orbit (hyperbola)?
Attach plots of curves that you believe to be an elliptical, hyperbolic, and parabolic orbits. On each plot,
write the values of tf and vy0 for which it was obtained.
Since a parabolic orbit looks very much like a hyperbolic one, write a brief explanation as to why you
believe that your graph is a piece of a parabola rather than that of a hyperbola. (Base your explanation
on the difference in the asymptotes of parabolas and hyperbolas. You must have learned about these
asymptotes in Calculus III; if not, or if you do not remember it, find the information in a textbook or
online.)
4
Bonus 1 This problem is worth 0.5 point.
Change the initial condition in hw2_p7.m to be (0, 2) and repeat the assignment of Problem 7.
Note: The main challenge here is to find a parabolic orbit and convincingly demonstrate that it is neither
an arc of an ellipse nor a hyperbola.
Hint: You may need to vary both components of the initial velocity.
(i) I’ve just disproved Kepler’s First Law and its derivation from Newton’s Laws. This sounds like the
Nobel Prize in Physics.
(ii) I’m losing my mind.
(iii) The instructor has given us an incorrect code to simulate the orbit.
(iv) Something else is going on.
Upon some reflection, you may conclude that the answers to the first two possibilities may be ‘not yet’
(although one may still be more likely than the other). As far as the third possibility, this instructor has
spent about an hour checking his code in various ways after he had noticed the fact stated in the second
sentence of this problem. The code is correct. Thus, it must be the fourth possibility.
The reason for the observed discrepancy is painfully simple. Examine the axis property of the plot
command in Matlab: type help axis. Draw a conclusion from what you have read (you will need
only one piece of information from there), and amend the plotting command so that the resulting plot
would agree with Kepler’s First Law.
Explain how you have amended the plotting command; attach the modified plot; and convince me that it
no longer violates the Kepler Law.
Bonus 3 In Problem 7 you were asked to provide a qualitative explanation whether the observed trajectory is a
parabola or a hyperbola. Here you are asked to rigorously predict for what value of v0y in hw2_p7.m
the trajectory will be a parabola. You will need to do so by computing the value of eccentricity e from
the other parameters given in this problem. Follow these steps.
1) Look inside the code hw2_p7.m and determine the initial values: r(0), θ(0), ẋ(0), ẏ(0).
2) Using the relation between Cartesian and polar coordinates, express ẋ and ẏ in terms of ṙ and θ̇ at
all times. Now, substitute in these two equations the initial values found in Step 1).
3) Use Eq. (2.11) of the Notes to relate ṙ to θ̇; now, using the initial values, relate ṙ(0) to θ̇(0). In
the process, you will need a value for eD; this is found from the initial conditions for r and θ.
4) Finally, substitute ṙ(0), expressed via θ̇(0), into the equations for ẋ(0), ẏ(0) obtained in Step 2).
For what value of v0y will you obtain the value of e corresponding to a parabola?