Query 7

#$&*

course Mth 174

The relative accuracy of the trapezoidal and midpoint rules depends on the nature of the function. As shown in my previous notes for a positive function which is concave down the midpoint rule is more accurate. A brief numerical example is y = x^2 on the interval (-1, 0), where trap gives 1/2, actual is 1/3 and mid is 1/4. 1/3 is closer to 1/4 than to 1/2. The same is true if the function is positive and concave up, for the same reasons.In fact I think that the argument can be extended to show that if concavity doesn't change on an interval, mid has to beat trap. I'd have to draw a picture or two to be sure, but it seems to be fairly obvious so I'll leave that to you. Use my previous notes as a guide.

Then you might want to draw a picture where trap is more accurate than mid. I'd recommend starting with something like y = x^3 on (-1,1), for which trap and mid agree with the accurate integral to give you 0. If you leave the points (-1, -1) and (1, 1) alone trap won't change but mid can be changed so that mid becomes less accurate (just pull (0,0) up or down a little and let the curve follow). How could the curve be manipulated to make trap less accurate than mid?

174

assignment # 7

......!!!!!!!!...................................

Questions related to text assignment:

*********************************************

Question: Section 7.6 Problem 6

approx using n=10 is 2.346; exact is 4.0. What is n = 30 approximation if original approx used LEFT, TRAP, SIMP?

......!!!!!!!!...................................

21:44:07

YYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYY

Your solution:

If the original approximation used LEFT, the error is -1.654 *(10/30)= -0.551. Therefore, the estimate for n=30 would be 4.000 + -0.551 = 3.449. If the original estimate used TRAP, the error equals -1.654 *(10/30)^2 = -0.184. So n=30 would be 4.000 + 3.816. If the original was done by SIMP, the error would be -1.654 * (10/30)^4 = -0.020. Therefore, the n=30 approximation would be 4.000 - 0.020 = 3.98.

confidence rating #$&*:3

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.............................................

Given Solution:

LEFT and RIGHT approach the exact value in proportion to the number of steps used.

MID and TRAP approach the exact value in proportion to the square of the number of steps used.

SIMP approachs the exact value in proportion to the fourth power of the number of steps used.

Using these principles we can work out this problem as follows:

** The original 10-step estimate is 2.346, which differs from the actual value 4.000 by -1.654.

If the original estimate was done by LEFT then the error is inversely proportional to the number of steps and the n = 30 error is (10/30) * -1.654 = -.551, approximately. So the estimate for n = 30 would be -.551 + 4.000 = 3.449.

If the original estimate was done by TRAP then the error is inversely proportional to the square of the number of steps and the n = 30 error is (10/30)^2 * -1.654 = -.184, approximately. So the estimate for n = 30 would be -.184 + 4.000 = 3.816.

If the original estimate was done by SIMP then the error is inversely proportional to the fourth power of the number of steps and the n = 30 error is (10/30)^4 * -1.654 = -.020, approximately. So the estimate for n = 30 would be -.02 + 4.000 = 3.98. **

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

Self-critique (if necessary):

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

Self-critique Rating: OK

*********************************************

Question: Section 7.6 Problem 5

problem 7.6.5 (previously problem 7.6.10) If TRAP(10) = 12.676 and TRAP(30) = 10.420, estimate the actual value of the integral. **** What is your estimate of the actual value and how did you get it?

YYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYY

Your solution:

The difference between TRAP(10) and TRAP(30), which is -2.256, is about 8/9 of the error of TRAP(10). We have to use inverse square proportionality, so -2.256 * (9/8) = -2.538 is the error of TRAP (10). We add this to 12.676 to get 10.138, which is our estimate of the actual value of the integral.

confidence rating #$&*: 3

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.............................................

Given Solution:

** You need to use the inverse square proportionality of the error with the number of steps.

Trap(30) is approximately (10/30)^2 = 1/9 of TRAP(10). So the difference 10.420 - 12.676 = -2.256 between TRAP(10) and TRAP(30) is approximately 8/9 of the error of TRAP(10).

It follows that the error of TRAP(10) is 9/8 * -2.256 = -2.538. Our best estimate of the integral is therefore -2.538 + 12.676 = 10.138. **

Ok I notice that I multiplied by 8/9 and you multiplied by 9/8 where did this come from?????? And why???

If you know 8/9 of a number (as in this case, where you know 8/9 of the difference between the first approximation and the accurate value), then to get that number you take 9/8 of the number you know (so you take 9/8 of the difference between the n = 10 approximation and the n = 30 approximation).

In other words, (trap(30) - trap(10)) is 8/9 * (accurate value - trap(10)), so (accurate value - trap(10)) = 9/8 * (trap(30) - trap(10)).

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

Self-critique (if necessary):

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

Self-critique Rating: OK

*********************************************

Question: Section 7.7 Problem 2

problem 7.7.2 (previously 7.7.12) integrate 1 / (u^2-16) from 0 to 4 if convergent

YYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYY

Your solution:

We know that 1/8 int(1 / (u-4), u, 0, 4) < int(1 / (u^2-4), u, 0, 4) < 1/4 int(1 / (u-4), u, 0, 4). The integral of 1/(u-4) from 0 to x equals ln(x-4) - ln(4). As x approaches 4, ln(x-4) approaches -infinity. Therefore, the integral diverges.

confidence rating #$&*: 3

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.............................................

Given Solution:

1 / (u^2-16) = 1 / [(u+4)(u-4)] . Since for 0 < x < 4 we have 1/8 < 1 / (u+4) < 1/4, the integrand is at most 1/4 times 1/(u-4) and at least 1/8 of this quantity, so the original integral is at most 1/4 as great as the integral of 1 / (u-4) and at least 1/8 as great. That is,

1/8 int(1 / (u-4), u, 0, 4) < int(1 / (u^2-4), u, 0, 4) < 1/4 int(1 / (u-4), u, 0, 4).

Thus if the integral of 1 / (u-4) converges or diverges, the original integral does the same. An antiderivative of 1 / (u-4) is ln | u-4 |, which is just ln(4) at the limit u=0 of the integral but which is undefined at the limit u = 4.

We must therefore take the limit of the integral of 1/(u-4) from u=0 to u=x, as x -> 4.

The integral of 1 / (u-4) from 0 to x is equal to ln(x-4) - ln(4). As x -> 4, ln(x - 4) approaches -infinity,

Thus the integral diverges.

STUDENT QUESTION

I am lost from the start of this problem I see where the integrand is (1/4) at its most but how can it be 1/8 at its least. Can you show me a step by step as to how we should have found this.

INSTRUCTOR RESPONSE

Since

1 / ( (u + 4) ( u - 4) ) = (1 / (u + 4) ) * ( 1 / (u - 4) ),

and since on this interval

1/8 < 1 / (u + 4) < 1/4

it follows that

1/8 | 1 / (u - 4) | < |1 / (u + 4) ) * ( 1 / (u - 4) | < 1/4 | 1 / (u - 4) |

and therefore that on this interval

1/8 integral | 1 / (u - 4) du | < integral | 1 / ( (u + 4) ( u - 4) ) du | < 1/4 integral | 1 / (u - 4) du |.

Since | 1 / (u - 4) | -> infinity as u -> 4, we proceed as indicated to show that this integral approaches -infinity.

Our previous inequality

1/8 integral | 1 / (u - 4) du | < integral | 1 / ( (u + 4) ( u - 4) ) du | < 1/4 integral | 1 / (u - 4) du |

thereby 'sandwiches' the magnitude of the desired integral between two values, both of which approach infinity.

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

Self-critique (if necessary):

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

Self-critique Rating: OK

*********************************************

Question: Section 7.7 Problem 7

problem 7.7.7 (previously 7.7.30) rate of infection r = 1000 t e^(-.5t)

YYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYYY

Your solution:

As we move to the right from t=0, the graph increases. Then e^(-.5 t) begins to approach zero much faster than t increases. This causes the graph to have an asymptote at the positive t axis. It goes back to almost touch the t axis. The concavity goes from negative for small t to positive (upward) for the larger t. I am not sure how to determine how many people get sick.

confidence rating #$&*: 1

^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

.............................................

Given Solution:

** First graph the function using standard graphing techniques from first-semester calculus:

This graph increases at first as you move to the right from t = 0. However e^(-.5 t) eventually approaches zero much faster than t increases so the graph has an asymptote at the positive t axis. So it increases for small positive t but eventually returns almost to the t axis, and it can't be strictly increasing. Its concavity changes from downward (negative) for small positive t to upward for larger t; the point at which the concavity changes is important.

We use the standard technique from first-semester calculus to find the point at which this function maximizes.

The first derivative is dr/dt = 1000 e^(-.5 t) - 500 t e^(-.5 t).

Setting this derivative equal to 0 we get

• 1000 e^(-.5 t) - 500 t e^(-.5 t) = 0;

dividing through by e^-.5 t we get the equation 1000 - 500 t = 0, which is easily solved to obtain t = 2. A first-or second-derivative test confirms that the t = 2 graph point is a relative maximum.

Concavity is determined by the second derivative r'' = e^(-.5 t) [ -1000 + 250 t ], which is 0 when t = 4. This is a point of inflection because the second derivative changes from negative to positive at this point. So the function is concave downward on the interval (-infinity, 4) and concave upward on (4, infinity).

The first derivative has a critical point where the second derivative is zero. This occurs at x = 4, which was identified in the preceding paragraph as the point of inflection for the original function. Since the second derivative goes from negative to positive, this point is a minimum of the first derivative. The first derivative is a decreasing function from t = 0 to t = 4 (2d derivative is negative) and is then an increasing function with asymptote y = 0, the x axis, which it approaches through negative values. Its maximum value for t >= 0 is therefore at t = 0. **

People are getting sick the fastest when the rate of infection is highest. This occurs at the relative maximum of the rate function, which was found above to occur at t = 2. Thus people are getting sick the fastest 2 days after the epidemic begins.

To find how many people get sick during a time interval, you integrate the rate function over that interval. In this case the interval doesn't end; so you need to integrate the rate function r = 1000 t e^(-.5t) from t = 0 until forever, i.e., from t = 0 to t = infinity.

An antiderivative of the function is F(t) = 1000 int ( t e^(-.5 t)) = 1000 [ -2 t e^(-.5t) - int ( e^(-.5 t) ) ] = 1000 [ -2 t e^(-.5 t) - 4 e^(-.5 t) ].

Integrating from 0 to x gives F(x) - F(0) = 1000 [ -2 t e^(-.5 x) - 4 e^(-.5 x) ] - 1000 [ -2 * 0 e^(-.5 *0 ) - 4 e^(-.5 * 0 ) ] = 1000 e^-(.5 x) [ -2 t - 4 ] - (-4000).

As x -> infinity, e^-(.5 x) [ -2 t - 4 ] -> 0 since the exponential will go to 0 very much faster than (-2 x - 4) will approach -infinity. This leaves only the -(-4000) = 4000.

** The calculator is fine for checking yourself, but you need to use the techniques of calculus to determine inflection points, maxima, minima etc.. The careful use the calculator to enhance rather than replace mathematical understanding. I get a lot of students in these courses who are now at 4-year institutions and who have taken courses based on the graphing calculator, or even TI-92, and many of them tend to have a very difficult time in courses that don't permit them, and in courses were mathematical understanding is required. **

** You have to use the techniques of calculus to determine these behaviors. Plugging values in won't show you the exact location of intercepts, maxima, minima, etc.. **

STUDENT QUESTION

I didn’t know where to go with the antiderivative but I think I understand your conclusion on that as well.

INSTRUCTOR RESPONSE

The infection is the rate-of-change function, so the antiderivative is the change-in-amount function.

Specificly we have the rate of change of the number of people who are or have been sick, with respect to clock time. The 'amount' is the number of people, so the antiderivative function is the change in the number of people (i.e., in the number who have been or are sick).

The definite integral between two clock times therefore tells you how many people are or have been sick between those clock times. If we integrate from some clock time from the initial instant to infinity, we get the total number of people who will get sick.

STUDENT COMMENT

I used the derivatives because this is a rate problem and the antiderivative is just a change in quantity formula.

INSTRUCTOR RESPONSE

You are given the rate.

If you know the average rate of change of change of a quantity with respect to time on an interval, you multiply it by the time interval in order to get the change in the quantity.

This is as opposed to subtracting two quantities and dividing by the time interval, which corresponds to the derivative.

Summing up average rate * time interval corresponds to the integral, which is what should have been used here.

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

Self-critique (if necessary):

I know understand how to use the derivatives to solve this problem. People are getting sick the fastest when the rate of infection is highest.

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

Self-critique Rating: OK

"

Self-critique (if necessary):

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

Self-critique rating:

"

Self-critique (if necessary):

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

Self-critique rating:

#*&!

&#Your work looks good. Let me know if you have any questions. &#