Chapter 11
Chapter 11
Chapter 11
Parametric equations are useful for describing curves not modeled by functions or motions parametrized by quantities such as
time. The standard form for a set of parametric equations is
xt = f t
yt = gt
Example 11.1. Plot the curve described by the parametric equations x = cos t and y = sin t for 0 § t § 2 p.
Solution:
0.5
Out[276]=
-1.0 -0.5 0.5 1.0
-0.5
-1.0
NOTE: Recall that the above parametric equations represent the unit circle. However, Mathematica may produce a graph that,
depending on its default settings, looks visually like an ellipse due to different scalings of the x- and y-axes. In that case, the plot
option AspectRatio can be used to specify the ratio of the height to the width for a plot. For example, to stretch the plot above so
that the circle becomes elliptical where the height is twice as long as the width, we can set AspectRatio equal to 2.
2 Mathematica for Rogawski's Calculus
0.5
Out[277]=
-1.0 -0.5 0.5 1.0
-0.5
-1.0
Example 11.2. Plot the curve described by the parametric equations x = t2 - 4 and y = t 2 for -2 § t § 3.
1.0
0.5
Out[278]=
-4 -2 2 4
-0.5
-1.0
Example 11.3. Plot the curve (prolate cycloid) described by the parametric equations x = 2 q - 4 sin q and y = 2 - 4 cos q for
0 § t § 2 p and determine its y-intercepts.
In[279]:= Clearf, g,
f_ 2 4 Sin;
g_ 2 4 Cos;
ParametricPlotf, g, , 4 Pi, 4 Pi,
PlotLabel "Prolate cycloid"
Prolate cycloid
6
4
Out[282]=
2
-20 -10 10 20
-2
To find the y-intercepts (there appears to be two based on the graph), we solve f q = 0 for q:
In[283]:= Solvef 0,
Solve::nsmet : This system cannot be solved with the methods available to Solve. à
Out[283]= Solve2 4 Sin 0,
Observe that the Solve command here fails to give us an answer since the equation is non-algebraic. However, note that the
cally using the FindRoot command (even the NSolve command fails in this case), where we provide a nearby location (q = p 2)
negative y-intercept is trivially located at y = -2 corresponding to q = 0. The other (positive) y-intercept must be solved numeri-
Out[284]=
NOTE: Observe that we used of the option PlotLabel to print the label "Prolate cycloid" in the plot above. In general, the inside
PlotLabel Ø"text" inside a plot command prints the title text for the given plot.
g' t
dy
f ' t
dy
dx
= dt
dx
=
4t 4 t2
x= and y =
1+t3 1+t3
a) Plot the curve described by the parametric equations above (select an appropriate interval for t that captures all the salient
features of the graph).
b) Find all points of horizontal tangency on the curve.
4 Mathematica for Rogawski's Calculus
Solution:
ft_ 3 t 1 t^ 3;
In[286]:= Clearf, g, t
gt_ 3 t^ 2 1 t ^3;
ParametricPlotft, gt, t, 0, 20, PlotRange All, AspectRatio 1,
PlotLabel "Folium of Descartes", ImageSize 200
Folium of Descartes
1.5
1.0
Out[289]=
0.5
NOTE: The plot above does not reveal the full graph of the folium. A more complete graph is shown in the following plot. The
dashed line indicates an asymptote. Can you generate a Mathematica plot of it? Can you find an equation of the asymptote (see
Exercise 7)? Hint: Beware of the discontinuity at t = -1.
Folium of Descartes
3
-3 -2 -1 1 2 3
-1
-2
-3
b) In order to find points of horizontal tangency, that is, points where the slope of the tangent line is equal to zero, it suffices to
dy dy dx
solve dx
= 0, or equivalently, dt
= 0 (assuming dt
∫ 0). Hence, we evaluate
In[290]:= SolveDgt, t 0, t
dx
does not vanish at t = 0 and t = 213 (we ignore the imaginary solutions), we conclude that there are two points
corresponding to these values at which the tangent lines are horizontal, namely at 0, 0 and 213 , 223 :
Since dt
ft, gt . t 0
ft, gt . t 213
In[291]:=
Out[291]= 0, 0
Out[293]=
1
Power::infy : Infinite expression encountered. à
0
Out[295]= ComplexInfinity
s = a dt + dt „ t
b dx 2 dy 2
Suppose xt, yt now represent the position of a particle moving along a path at time t. The distance traveled by the particle
over the interval t0 , t is given by
st = t du + du „ u
t dx 2 dy 2
0
and
dt + dt
ds dx 2 dy 2
dt
=
represents its speed (length of the velocity vector x ' t, y ' t).
Example 11.5. Find the arc length of the curve x = e-t cos t, y = e-t sin t for 0 § t § p 2.
6 Mathematica for Rogawski's Calculus
xt_ E^ t Cost;
In[296]:= Clearx, y, t
yt_ E^ t Sint;
ParametricPlotxt, yt, t, 0, Pi 2, ImageSize 250
0.30
0.25
0.20
Out[299]= 0.15
0.10
0.05
Dxt, t ^ 2 Dyt, t ^2 t
2
In[300]:=
2 1 2
0
Out[300]=
In[301]:= N
Out[301]= 1.12023
Solution:
a) The distance traveled is given by
In[302]:= Clearx, y, t
xt_ t^ 2 4;
yt_ t 2;
Dxt, t ^ 2 Dyt, t ^2 t
3
3
12
1
Out[305]= 145 ArcSinh12
8
In[306]:= N
Out[306]= 18.4599
speed Dxt, t ^2 Dyt, t ^2
Plotspeed, t, 3, 3, ImageSize 250
In[307]:=
1
Out[307]= 4 t2
4
6
Out[308]= 3
-3 -2 -1 1 2 3
This shows that the minimum speed is 1/2 and occurs at t = 0. Can you verify this using calculus techniques?
ü Exercises
In Exercises 1 through 3, sketch the curve represented by the given parametric equations. Be sure to select an appropriate
4. Find all points of horizontal and vertical tangency to the curve x = cos q + q sin q, y = sin q - q cos q, 0 § q § 2 p.
b) Find the slope of the line tangent to the curve at the point where t = p 4.
a) Graph the curve represented by the parametric equations above.
a) Approximate the distance traveled by the particle over the interval 0, 2 p. Hint: Use the NIntegrate command.
6. Consider a particle moving along a curve described by x = t - cos t and y = t - sin t with respect to time t.
b) Find the minimum and maximum speeds of the particle over the same interval. At what times do they occur?
7. Cornu's spiral (also known as Euler's spiral) is a curve defined by the following Fresnel integrals:
8. Find the asymptote corresponding to the folium of Descartes (see Example 11.4) and plot the asymptote together with the
folium of Descartes.
g p
this section.
Polar coordinates, expressed as r, q), represent the location of a point on the Cartesian plane in terms of distance r from the
origin and angle q with respect to the positive x-axis.
r2 = x2 + y2 x = r cos q
q = tan-1 y x y = r sin q
Solution:
r 32 42
NArcTan4 3
In[309]:=
Out[309]= 5
Out[310]= 0.927295
x 7 CosPi 3
y 7 SinPi 3
In[311]:=
7
Out[311]=
2
7 3
Out[312]=
2
Out[313]=
-7 -6 -5 -4 -3 -2 -1
-2
-4
Example 11.8. Plot the graph of the six-leaf rose r = 2 cos3 q 2.
1.5
1.0
0.5
Out[314]=
-2 -1 1 2
-0.5
-1.0
-1.5
Moreover, the area A of the region bounded by a polar equation in the form r = f q between a § q § b is given by
2 a
f q „ q
1 b 2
A=
Example 11.9. Locate all horizontal and vertical tangents of the limacon r = 2 - sin q.
10 Mathematica for Rogawski's Calculus
1.0
0.5
Out[315]=
0.5 1.0 1.5 2.0
-0.5
-1.0
From the plot above, we should expect to find two horizontal tangents and three vertical tangents.
In[316]:= Clearf,
f_ 1 Cos
Out[317]= 1 Cos
In[318]:= dydx Simplifyf ' Sin f Cos f' Cos f Sin
Csc 2
3
Cos 2
Out[318]=
1 2 Cos
In[319]:= dydx 0
Solvedydx 0,
Csc 2
3
Cos 2
Out[319]= 0
1 2 Cos
As for vertical tangents, we solve for where the reciprocal of the derivative is zero, that is, 1 dy dx = 0 for q.
Chapter 1 11
1 dydx
Solve1 dydx 0,
In[321]:=
0, , , ,
4 2 2 4
Out[322]=
3 3 3 3
Since the solution above only gives us the principal solution q = 0 of -tan3 q 2 = 0, we need to additionally solve 3 q 2 = ≤p
for q, which yields our two other solutions, q = ≤ 2 p 3.
Example 11.10. Find the area of the region contained inside the circle r = 3 sin q and outside the convex limacon r = 2 - sin q.
Solution: We first plot the two polar curves on the same set of axes.
Out[323]=
-2 -1 1 2
-1
-2
-3
Next, we find their points of intersection by equating them and solving for q:
Out[324]=
6
Observe that Mathematica gives only the solution q = p 6, which lies in the first quadrant since trigonometric inverse functions
are involved. We can see from the above graph that the other point of intersection must be at q = 5 p 6 due to symmetry. Thus,
the area of the enclosed region is given by the difference in areas enclosed by the circle and limacon between q = p 6 and
12 Mathematica for Rogawski's Calculus
q = 5 p 6:
Out[325]= 3 3
In[326]:= N
Out[326]= 5.19615
NOTE: Using even symmetry of our region, it would have been enough to integrate between q = p 6 and q = p 2 and double the
result.
ü Exercises
2. Convert the polar coordinates 5, 3 p 4 into Cartesian coordinates. What if we replace 5, 3 p 4 with -5, 3 p 4?
In Exercises 3 through 6, plot the graph of each of the given polar equations and find an interval for q over which each graph is
8. Find all horizontal and vertical tangents of the lemniscate r2 = cos2 q. Plot its graph to confirm your answers.
In Exercises 10 through 12, graph and find the area of each of the given regions.
Conic sections refer to the three families of curves (ellipses, hyperbolas, parabolas) generated by intersecting a plane with a cone.
Recall the equations for describing each family of curves in standard position:
Chapter 1 13
I. Ellipse
a + b = 1
x 2 y 2
If b > a > 0, then the ellipse has two foci located at 0, ≤c, where c = b2 - a2 .
II. Hyperbola
a - b = 1
x 2 y 2
III. Parabola
1
y= x2
4 c2
The most useful command for plotting conic sections is ContourPlot[eqn,{x,a,b},{y,c,d}], where eqn is the equation of the conic.
Example 11.10. Determine the family that each conic section below belongs to and then make a plot of each.
x2 y2
a) 9
+ 16
=1
9
b) y = 4
x2
x2 y2
c) 4
- 9
=1
Solution:
Out[327]=
-4 -2 2 4
-2
-4
Observe that the length of semi-major and semi-minor axes are 4 and 3, respectively. How would this change if we happen to
switch the coefficients 9 and 16?
b) This conic is a parabola. Since the equation here is solved for y, we merely use the Plot command:
14 Mathematica for Rogawski's Calculus
Out[328]= 4
-2 -1 1 2
Out[329]=
-6 -4 -2 2 4 6
-2
-4
-6
Example 11.11. Find an equation of an ellipse with center at -1, 3 and having semi-major and semi-minor axes of lengths 5
and 1/2, respectively.
Solution: From the given data, we see that x0 , y0 = -1, 3, a = 5 , and b = 1 2. The equation of our ellipse is therefore:
+ 4 y - 32 = 1
x+12
5
4.0
3.5
Out[330]= 3.0
2.5
-4 -3 -2 -1 1 2
Chapter 1 15
ü Exercises
In Exercises 1 through 4, plot each of the given conic sections. Can you determine the family that each conic section belongs to
before plotting? Also, what are the values of the parameters a and b (or c) in each conic section?
x2 y2 y2 x2 y-22 x+12
1. 25
+ 4
=1 2. 36
- 16
=1 3. y = 3 x2 4. 25
- 49
=1
x2
5. Consider the ellipse 9
+ 4 y2 = 1.
a) Make a plot of the ellipse.
b) What are the lengths of the semi-major and semi-minor axes?
c) Where are the foci located?
d) Compute the sum of the distances from the two foci to any point on the ellipse. Do you recognize the answer?
6. Find an equation of an ellipse with center 1 2, -5 and having semi-major and semi-minor axes of lengths 3/4 and 7,
respectively.
x2 y2
7. Find all points on the hyperbola 25
- 9
= 1 where its slope equals 1.
de
8. Consider a polar curve of the form r = 1+e cos q
, where d and e are non-negative constants.
a) Plot this curve for d = 3, and e = 1 2. Do you recognize this curve as a conic section? Of which type? Hint: Use the com-
mand PolarPlot.
b) Repeat part a) but this time use e = 2 instead. Do you recognize this curve as a conic section? Of which type?
c) Repeat part a) but this time use e = 1 instead. Do you recognize this curve as a conic section? Of which type?
d) Describe how the graph changes as we vary the values d and e. What happens to the graph when e = 0? NOTE: The value e
is called the eccentricity of the conic section.