m103ex-3-10.mw

Applications of Derivatives

Thomas/Finney Section 3.1, page 195

Problem 10

Find the maximum and minimum of  f(x) = 4-x^2 for x in the interval [-3, 1]

To find absolute extrema on closed intervals, check (a) where the derivative is 0, (b) where it fails to exist, and (c)  at the endpoints:

> y:=4-x^2;

y := 4-x^2

> dy:=diff(y,x);

dy := -2*x

> solve(dy=0,x);

0

The derivative exists for all x, and so we have to check the interior critical point x=0 as well as the endpoints x=-3 and x=1:

> [-3,subs(x=-3,y)],[0,subs(x=0,y)],[1,subs(x=1,y)];

[-3, -5], [0, 4], [1, 3]

We see that the minimum is -5 (at x=-3) and the maximum is 4 (at x=0).

Section 3.2, page 203

Problem 10

Let f(0) = 3 ,  f(x) = -x^2+3*x+a for x  in (0,1), and f(x) = m*x+b for x in [1, 2].  For what value(s) of a , m , and  b does f(x) satisfy tyhe hypotheses of the mean value theirem on the interval [0, 2]?

We need to choose a, m and b so that the function

> f:=piecewise(x=0,3,0<x and x<1,-x^2+3*x+a,1<=x and x<=2,m*x+b);

f := PIECEWISE([3, x = 0], [-x^2+3*x+a, -x < 0 and x < 1], [m*x+b, 1-x <= 0 and x-2 <= 0])

satisfies the hypotheses of the mean-value theorem on [0,2]. So it has to be continuous on the closed interval [0,2] (in particular at x=0 and x=1), and differentiable on the open interval (0,2) (in particular at x=1). So we solve three equations for three unknowns:

> solve({3=limit(-x^2+3*x+a,x=0,right),limit(-x^2+3*x+a,x=1,left)=limit(m*x+b,x=1,right),limit(diff(-x^2+3*x+a,x),x=1,left)=limit(diff(m*x+b,x),x=1,right)},{a,m,b});

{m = 1, a = 3, b = 4}

> f:=subs(%,f);

f := PIECEWISE([3, x = 0], [-x^2+3*x+3, -x < 0 and x < 1], [x+4, 1-x <= 0 and x-2 <= 0])

This function should be continuous on [0,2] and differentiable on (0,2) -- let's see:

> plot(f,x=0..2,color=blue,thickness=2);

[Plot]

Section 3.3, page 208,  Problem 21

Let g(x) = x*sqrt(8-x^2) .  

(a) Find the intervals in which this is increasing or decreasing

(b) Identify the function's local extreme values, if any, saying where they are taken on.

(c) Which, if any, of the extreme values are absolute?

Start by defining the function and its derivative:

> g:=x*sqrt(8-x^2);

g := x*(8-x^2)^(1/2)

> dg:=diff(g,x);

dg := (8-x^2)^(1/2)-x^2/(8-x^2)^(1/2)

Of course, the domain of this function is where x^2 <= 8 , i.e., for x between -2*sqrt(2) and 2*sqrt(2) .

Let's see where the derivative is positive and negative:

> solve(dg>0,x);

RealRange(Open(-2), Open(2))

> solve(dg<0,x);

RealRange(Open(-2*2^(1/2)), Open(-2)), RealRange(Open(2), Open(2*2^(1/2)))

So the function is increasing when x is between -2 and 2, and decreasing to the left of -2 and to the right of 2. We calculate the critical values:

> simplify(subs(x=-2,g)),simplify(subs(x=2,g));

-4, 4

So the min point of the graph is at [-2,-4] and the max point is at [2,4]. These are absolute, since g is zero at the endpoints of the domain.

Now let's plot:

> plot(g,x=-2*sqrt(2)..2*sqrt(2),color=blue,thickness=2);

[Plot]

Section 3.4, page 217, Problem 21

Graph y = x^5-5*x^4 , including the coordinates of extreme points and inflection points.

We use calculus to make sure we have all the "interesting" points in our Maple plot:

> y:=x^5-5*x^4;

y := x^5-5*x^4

> solve(y=0,x);

5, 0, 0, 0, 0

> solve(diff(y,x)=0,x);

4, 0, 0, 0

> solve(diff(y,x$2)=0,x);

3, 0, 0

Since y is a polynomial, we don't have to worry about any asymptotes. We need only

plot on an interval that contains all of x=0,3,4,5. To get an idea of the range on the y-axis, we'll substitute our "interesting" x-values into the function:

> subs(x=0,y),subs(x=3,y),subs(x=4,y),subs(x=5,y);

0, -162, -256, 0

So we have to include a pretty big y range:

> plot(y,x=-4..7,-350..100,color=blue,thickness=2);

[Plot]

Section 3.5, page 230, Problem 101

Graph y = -(x^3-2)/(x^2+1) on the three intervals  [-9, 9], [-90, 90], [-900, 900].

The graph on [-9, 9] should be good.  That on [-90, 90] should show that something is happening near the origin, but be inconclusive, while the graph on [-900, 900]  will almost look like y = -x  (why?).

First we'll do the three graphs:

> y:=-(x^3-2)/(x^2+1);

y := -(x^3-2)/(x^2+1)

> plot(y,x=-9..9,color=blue,thickness=2);

[Plot]

> plot(y,x=-90..90,color=blue,thickness=2);

[Plot]

You can see a very tiny kink in the graph at the origin.

> plot(y,x=-900..900,color=blue,thickness=2);

[Plot]

Now it just looks like the line y=-x. This is because:

> simplify(y+x);

(2+x)/(x^2+1)

> limit(y+x,x=infinity),limit(y+x,x=-infinity);

0, 0

So the line y=-x is an oblique asymptote for the graph, and once you subtract the "dominant" -x from y, you get a function that goes rapidly to 0 as x gets large.

Section 3.6, page 242, Problem 30

A wall 8 feet high stands 27 feet from the larger wall of a building.  A straight beam goes from the ground to the top of the 8 foot wall and extends to touch the wall of the building.  Find the length of the shortest such  beam.

Let x be the length of the beam, let h be the height on the building's wall where the beam touches, and let s be the distance from the little wall where the other end of the beam touches the ground. Then we want to minimize x (or, what is the same, x^2 ), where we have the following relationships between x, s and h:

> eqn1:=x=sqrt(h^2+(27+s)^2); eqn2:=h/8=(s+27)/s;

eqn1 := x = (h^2+729+54*s+s^2)^(1/2)

eqn2 := h/8 = (27+s)/s

The second equation is the "constraint", and the first provides our "objective function". So we solve the second for h in terms of s, and then substitute into the first:

> eqn1:=subs(h=solve(eqn2,h),eqn1);

eqn1 := x = (64*(27+s)^2/s^2+729+54*s+s^2)^(1/2)

Now we can take the derivative and set it equal to 0:

> solve(diff(rhs(eqn1),s)=0,s);

12, -6+6*I*3^(1/2), -6-6*I*3^(1/2)

We certainly don't want an imaginary root! So the minimum must occur when s=12 feet. We can substitute this into the first equation to see what x should be:

> simplify(subs(s=12,eqn1));

x = 13*13^(1/2)

So the shortest beam that reaches is 13*sqrt(13) feet long.

Problem 56

You operate a tour service that offers the following rates:

a) $200 per person if 50 people (the minimum number to book the tour) go on the tour.

b) For each additional person (up to a maximum of 80 people total), everyone's charge is reduced by $2.

It costs $6000 (a fixed cost) plus $32 per person to conduct the tour.  How many people does it take to maximize your profit?

From the information given in the problem we have the following formulas, where x is the number of people who are on the tour.

> cost:=6000+32*x;

cost := 6000+32*x

Given that at least 50 people go on the tour, we have the following formula for revenue.

> revenue:=x*(200-2*(x-50));

revenue := x*(300-2*x)

> profit:=revenue-cost;

profit := x*(300-2*x)-6000-32*x

To maximize profit, we find where its derivative (also known as the marginal profit)  is equal to zero, and then we must check that this is in fact a maximum.

> mp:=diff(profit,x);

mp := 268-4*x

> solve(mp=0);

67

The critical point of the profit function is 67 people. Since the second derivative of profit is -4 (we don't need Maple to tell us that!), this must be a maximum for the profit. And the profit is:

> subs(x=67,profit);

2978

dollars.

Section 3.7, page 259, Problem 42

Estimate the change in the lateral surface area S = Pi*r*sqrt(r^2+h^2) of a right circular cone when the radius changes from r[0] to r[0]+dr and the height does not change.

We have to calculate the differential of the surface area:

> S:=Pi*r*sqrt(r^2+h^2);

S := Pi*r*(r^2+h^2)^(1/2)

> subs(r=r[0],simplify(diff(S,r)))*dr;

Pi*(2*r[0]^2+h^2)*dr/(r[0]^2+h^2)^(1/2)

This is how much the surface area changes if r changes from r[0] to r[0]+dr .

Section 3.8, page 267 , Problem 19

The curve y = tan(x) crosses the line y = 2*x  between x = 0 and x = Pi/2 .Use Newton's method to compute a numerical approximation to this value.

First, let's plot y=tan(x) and y=2x to see that they do cross:

> plot({tan(x),2*x},x=0..Pi/2,0..5,color=blue,thickness=2);

[Plot]

OK, now we have to set up Newton's method - we might as well start with the initial guess x=1.

> f:=x->tan(x)-2*x;

f := proc (x) options operator, arrow; tan(x)-2*x end proc

> p:=1.0;

p := 1.0

> p:=p-f(p)/D(f)(p);

p := 1.310478030

> p:=p-f(p)/D(f)(p);

p := 1.223929096

> p:=p-f(p)/D(f)(p);

p := 1.176050900

> p:=p-f(p)/D(f)(p);

p := 1.165926508

> p:=p-f(p)/D(f)(p);

p := 1.165561636

> p:=p-f(p)/D(f)(p);

p := 1.165561185

> p:=p-f(p)/D(f)(p);

p := 1.165561185

It looks like we have converged to the answer. Just to check, let's use Maple's fsolve:

> fsolve(f(x)=0,x=0..2);

1.165561185