assignment 6

#$&*

course M 279

2/6 5

For P = 2000, then what should be the slope of the direction field, according to the equation dP/dt = .06 P?****

dP/dt = 120

#$&*

Do the slopes here appear to match the slope you just calculated?

The slope you should have obtained was over 100. The slopes on this direction field all appear to be less than 1. We need to either reconcile this, or discard the direction field of the figure.

Let's first check out a solution to the problem. As you can easily verify, dP/dt = k P gives us general solution P(t) = C e^(k t). To get the solution passing through the P axis at 2000, we apply the initial condition P(0) = 2000. Our particular solution for this initial condition is P(t) = 2000 e^(k t). Our direction field is for k = .06; for P(0) = 2000 our function is therefore P(t) = 2000 e^(.06 t). A graph of this function is consistent with the direction field.

Verify this. Following the trend of the direction field, estimate the likely value of P(3). What is your estimate?

2,500

#$&*

What do you get if you evaluate P(3)? Is this consistent with your estimate?

****

approx. 2,394, yes is consistent

#$&*

According to the direction field, does the solution curve through (0, 2000) exit the box to the right or through the top?

****

right

#$&*

Check this out by evaluating P(12). What is your conclusion?

****

Correct. It is easy for me to visualize that if you start at 2,500 and follow lines you hit top right corner.

#$&*

Reasonable answers to these questions will tend to reinforce our confidence in the direction field.

Consider the direction field line through (8, 4000). Based on the scale of your graph, estimate the rise of that segment. It should be obvious that the run of the segment is 2. What therefore is the slope represented by that segment?

****

approx. 400/2 which equals 200 or .2 if we scale down.

#$&*

The slope appears to be less than 1. Explain why that segment actually represents a slope in the hundreds.

****

slope = rise/run = 400/2 = 200

#$&*

Construct the slope field for the equation dP/dt = .06 P, for the interval 10 000 <= P <= 20 000, 0 <= t <= 12. Scale your slope segments appropriately.

Based on your slope field, at what value of t do you estimate the value of P(t) will double, starting with P = 10 000 at t = 0?

****

P = 333,333.3333 approx.

 dP/dt = .06P

2 = .06P

P = 2/.06 = 33.3333, where P is on scale of 10,000

@& Check the link at the end for a graphical estimate.*@

#$&*

`q005 (from 110131). The equation dP/dt = k P + m is of the same form as y ' + p(t) y = g(t), with p(t) = -k and g(t) = m

Explain in detail why the equation is of this general form, specifically explaining how we get p(t) = - m and g(t) = m.

****

#$&*

This equation is therefore first-order linear and nonhomogeneous. We have learned the technique for solving such equations.

What is the general solution to our equation dP/dt = k P + m?

****

#$&*

`q006 (from 110131). I can mix a 10% salt solution into pure water by siphoning the salt solution from one bottle, into a second bottle which is initially full of pure water. The second bottle was initially full, so it will overflow.

We will make the simplifying assumption that the salt water flowing into the second bottle is instantly distributed equally throughout that bottle.

Assume that the rate of flow of salt solution is r (e.g., r might be, say, 7 cm^3 / second). The second bottle is assumed to have volume V (e.g., the second bottle might have a volume of 500 cm^3, corresponding to half a liter)

If p(t) is the percent concentration of salt in the second bottle, as a function of clock time, then at what rate is salt leaving the bottle at clock time t?

****

#$&*

Express your preceding result as a differential equation, find the general solution of the equation, and find the particular solution for which r = 7 cm^3 / sec.

Question `q001 (non-italicized), for the 2/02/11 class, was given previously. Subsequent questions are all non-italicized and are for the 2/02/11 class only.

`q002. A quantity y is directly proportional to a quantity x if there exists a constant k such that y = k x.

The rate, with respect to time, at which the temperature of an object left to cool in a room at temperature T_room is directly proportional to the difference between the temperature of the object, which we will represent by T(t), and the temperature of the room.

Letting T ' (t) stand for the rate at which the temperature T(t) changes with respect to clock time, write down the proportionality.

****

T ' (t) = -k(T(t) - T_room)

Because as temp of object approaches temp of room rate of temp change is less and less. When object reaches temp of room unless another factor is involved temp change stops. k is neg because room is cool.

#$&*

You should have written down a differential equation. That equation is linear. Find its general solution.

****

T ' (t) = -k(T(t) - T_room)

 T ' (t)/ (T(t) - T_room) = -k

Int(T ' (t)/ (T(t) - T_room)) = Int(-k T`(t) )

ln|(T(t) - T_room)| = -kt +c

T(t) - T_room = e^(-kt +c) = e^(-kt)*e^(c) = A e^(-kt)

 T(t) = T_room + A e^(-kt)

#$&*

If the object has an initial temperature of 50 Celsius, and its temperature initially changes at -5 Celsius per minute when placed in a room at temperature 20 Celsius, then what is the specific function that models its temperature as a function of time?

****

50 = 20 + A e^(-k*0)

30 = A

Knowing temp changes at -5C/min, then at t=1 T(t) = 45

45 = 20 + 30e^(-k*1)

25 = 30e^(-k)

ln|25/30| = -k

k = approx. .1823

So we have

T(t) = 20 + 30e^(-.1823t)

#$&*

`q003. Consider the equation dy/dt = -.01 y^2.

Explain why we cannot solve this equation using the techniques we have been applying to first-order linear equations.

****

Its not a linear equation because of the y^2

#$&*

We will see why shortly, but for now accept that this equation can be rearranged into the form

dy / y^2 = -.01 dt

and integrate both sides to get a general solution.

****

Int( 1/y^2 dy)

 Int(y^-2 dy) = -1/y

Int(-.01 dt)

-.01t

dy / y^2 = -.01 dt

-1/y = -.01t +c

-1 = (-.01t + c)*y

y = 1/(.01t + c)

#$&*

`q004. Consider the equation dy/dt = -.01 t y^2. Algebraically rearrange this equation so that dy is a factor of one side, with all other y factors on that side, and dt, along with all other t factors, is on the other.

Integration both sides of the equation to get a general solution.

1/y = -.01 t + c

y = 1 / (-.01 t + c)

Find the particular solution for which y = 25 when t = 0.

Show that the function f(t, y) = -.01 t y^2 is defined and continuous for all values of y and t.

****

25 = 1/c

c= 1/25

#$&*

Show that there exist values of t for which your solution is neither defined nor continuous.

****

Whenever (-.01 t + c) = 0 fn is neither defined nor continuous

#$&*

While a linear equation of the form y ' + p(t) y = g(t) always has a solution which is defined and continuous, on any interval for which p(t) and g(t) are defined and continuous, this example shows that the same is not necessarily so for an equation of the form y ' + f(t, y) = 0. Explain this statement.

****

#$&*

`q005. The rate, with respect to time, at which the velocity of a spherical object moving through water changes is directly proportional to the square of its velocity.

Write this as a differential equation.

****

dv/dt = k*v^2

#$&*

Solve the equation.

****

dv/dt = k*v^2

dv/v^2 = k dt

Int(dv/v^2) = Int(k dt)

 Int(1/v^2 dv) = Int(v^-2 dv) = -1/v

 Int(k dt) = kt + c

-1/v = kt

-1 = v*(kt +c)

v = -1/(kt+c)

#$&*

If the object is initially moving at 30 cm/s and changing velocity at the rate of 1000 cm/s^2, then what function describes its velocity?

****

For dv/dt = k*v^2, where

dv/dt = 1000 cm/s^2

v^2 = (30)^2 = 900

We solve for k,

1000 = k*900

k = 10/9 = approx. 1.1111

dv/dt = 1.1111*v^2 and

v = -1/(1.1111t + c)

?????I’m not sure this is correct??????

#$&*

@& This is good, but you haven't used all your information; in particular you haven't solved for c.

With units, which give it more meaning:

v(0) = 1 / (k * 0 + c) = 1 / c, and

v(0) = 30 cm/s so

1/c = 30 cm/s and

c = 1/(30 cm/s) = .0333... s / cm.

v ' (t) = -k / (k t + c)^2 = -k / (k t + .0333... s/cm) so

v ' (0) = -k / (.0333... s/cm)^2 = 1000 cm/s^2.

Thus

k = -1000 cm/s^2 * (.0333... s / cm)^2 = 1.11 ... cm^-1.

The velocity function is therefore

v(t) = 1 / (1.11... cm^-1 * t + .0333... s/cm).

If you plug in a value of t, including units, you will get a value of v which you will see comes out in cm/s.

*@

`q006. Recall implicit differentiation. If you need a review of implicit differentiation see the q_a_ document at http://vhcc2.vhcc.edu/dsmith/genInfo/qa_query_etc/cal1/cal1_qa16.htm . If you need a review of the chain rule, seehttp://vhcc2.vhcc.edu/dsmith/genInfo/qa_query_etc/cal1/cal1_qa12.htm and http://vhcc2.vhcc.edu/dsmith/genInfo/qa_query_etc/cal1/cal1_qa13.htm .

If y is a function of t, then what is the derivative with respect to t of the function H(t, y) = t cos y + 1/y?

****

For H(t, y) = t cos y + 1/y,

(tcos(y) + 1/y)’ = tcos(y)’ + (1/y)’

tcos(y)’ = t’*cos(y) + t*cos(y)’

= cos(y) - t*sin(y)*y’

(1/y)’ = (-1/y^(-2))y’

 H(t, y)’ = cos(y) - t*sin(y)*y’ - (-1/y^(-2))y’

#$&*

`q007. Verify that the derivative of the function H(t, y) = t^2 sqrt(y) is 2 t sqrt(y) + t^2 / (2 sqrt(y) ) * y '.

****

(t^2 sqrt(y)) ' = (t^2)’*`sqrt(y) + t^2 (sqrt(y))’

= 2t`sqrt(y) + t^2*1/(2*sqrt(y)), so

(t^2 sqrt(y)) ' = 2t`sqrt(y) + t^2*1/(2*sqrt(y))*y’

#$&*

Verify that the equation

2 t sqrt(y) + t^2 / (2 sqrt(y) ) * y ' = 0

is of the form

( H(t, y) ) ' = 0,

where H(t, y) = t^2 sqrt(y) and ' indicates the derivative with respect to t.

****

We have 2 t sqrt(y) + t^2 / (2 sqrt(y) ) * y ' = 0, which is of the form

n(y)(dy/dt) + m(t) = 0

Where …

 n(y)(dy/dt) = t^2 / (2 sqrt(y) ) * y '

 m(t) = 2 t sqrt(y)

Also h(t, y) = y’, so (H(t, y))’ is really just derivative of H(t, y) or h(t, y)

#$&*

Verify that the equation

( H(t, y) ) ' = 0

is equivalent to the equation

H(t, y) = c,

where c is any constant number.

****

For n(y)(dy/dt) + m(t) = 0, we can solve initial value problems by reversing the chain rule.

 n(y)(dy/dt) + m(t) = 0, we can let N(y) is any antideriv. of n(y)

d/dy(N(y)) = n(y), if y is differentiable fn of t then chain rule leads to

d/dt(N(y)) = (dN/dy)*(dy/dt)

=n(y)(dy/dt), this shows n(y)(dy/dt) can be written as dN/dt Note M(t) is antider. of m(t) then it to can be written as dM/dt

Now we have equation

 n(y)(dy/dt) + m(t) =(d/dt) N(t) +(d/dt) M(t) = d/dt [N(y) + M(t)]

this equation reduces to, d/dt [N(y) + M(t)] = 0

term N(y)+M(t) is a fn of t(involving unknown soln y=y(t)] so,

N(y) + M(t) = C, where C is an arbitrary constant, where we can always determine C by imposing the initial cond y(t_0) = y_0. Finding ….

C = N(y_0) + M(t_0)

#$&*

@& Check out Class Notes from 110207 for this solution.*@

For the given function H(t, y) = t^2 sqrt(y), solve the equation

H(t, y) = c

for y as a function of t.

****

H(t, y)’ = 2 t sqrt(y) + t^2 / (2 sqrt(y) ) * y ' = 0

y’ = -(2 t sqrt(y))/(t^2/(2`sqrt(y))) = -(4t`sqrt(y)^2)/t^2) = -4t^-1y = -4y/t

So we have y’ = -4y/t

y’/4y + 1/t = 0, taking integrals we find

Int(y’/4y ) = (1/4)Int(1/y dy) = (ln|y|)/4

 Int(1/t dt ) = ln|t|

(1/4)ln|y| + ln|t| = C

Mult through by e^1

(1/4)y + t = e^(C), where C absorbs the e

Sol.ve for y

(1/4)y + t = A

y = 4(A - t)

#$&*

Show that this function y(t) satisfies the equation

2 t sqrt(y) + t^2 / (2 sqrt(y) ) * y ' = 0.

****

y’ = -(2 t sqrt(y))/(t^2/(2`sqrt(y))) = -(4t`sqrt(y)^2)/t^2) = -4t^-1y = -4y/t

Int( 1/y dy) = ln|y|

Int(-4/t dt) = -4Int( 1/t dt ) = -4(ln|t| + C)

ln|y| = -4(ln|t| + C), mult through by e^1

y = -4(t + e^C) = -4(t +A)

y = 4(A – t)

???????????I messed up in my algebra somewhere but can’t see where. Maybe because A is arbitrary it absorbs something I missed but for the two equations my constant A has opposite signs?????

I’m pretty sure I understand the point you want us to see though.

#$&*

@& Check out Class Notes from 110207 for this solution.*@

"

Self-critique (if necessary):

------------------------------------------------

Self-critique rating:

________________________________________

#$&*

@& Good work overall. See my notes and also check out Class Notes from 110207.*@