Joseph M. Mahaffy SDSU
Math 124: Calculus for the Life Sciences -- Lab Help 8 Fall 2014
 |  Course Index  |  Lectures  |  HW Assign  |  Computer Lab  |  WeBWorK  |  Syllabus  |  Miscellaneous  ||  Home
 

Computer Lab Help 8


This page is designed to provide helpful information about the laboratory questions. You will find more details in the Lab Manual that accompanies this course. Begin this lab and every lab by introducing yourself to your partner. Determine the times when you can meet together during the week before the lab is due at your next Lab session. You should start this lab and each lab by typing the name of each team member and your computer number on the Lab Cover Page (or a copy of it).

The first WeBWorK problem asks questions about this help page and appropriate lecture material. This should help you work through the Lab more smoothly. In this Lab, the first problem studies the growth of fish using von Bertalanffy's equation to connect the length to age. An allometric model relates length and weight. These models are combined by a composite function, then differentiation is applied to find when the fish are gaining weight most rapidly as function of age. The second problem examines the logistic and Gompertz growth models for tumor cells. These growth models are fit to data, then simulated to fit the time data. The last problem is a classic optimization problem that uses trigonometric functions to compute an optimal volume of a trough by varying the angle of the sides.

Problem 1: This problem is very similar to the material in the lecture notes. We fit data to find the length of a fish as a function of its age using von Bertalanffy's equation. Then we create an allometric model for the weight of a fish as a function of its length. The questions that you answer are very similar to the ones in the Chain Rule homework problems. The fit to the von Bertalanffy equation is done using Excel's Solver very much the same way as you have done before, while the fitting to the allometric model follows the techniques learned in earlier Labs using Excel's Trendline. This is readily accessed by right clicking on the data. You may need to increase the number of significant digits given to you by Trendline. To obtain 5 significant figures, we right click on the formula, selecting Format Trendline Label, choose Number and Scientific, and finally select 5 Decimal Places. The use of Maple might be valuable for helping with the differentiations (though you can practice your skills by doing them by hand). For this problem, the first and second derivatives are zero at t = 0, because the curve is very flat there. To get the point of inflection and the value of the function and the derivative at the point of inflection you need the following commands (assuming you see the point of inflection is between t = 1 and t = 20:

> p := t -> model goes here;

> dp := diff(p(t),t); # This gives the derivative

> sdp := diff(dp,t); # This gives the second derivative

> tp := fsolve(sdp=0, t = 1..20); # This gives the point of inflection

> p(tp); evalf(%); # This gives the population at the point of inflection

> subs(t = tp, dp); # This gives the derivative at the point of inflection

Problem 2: This problem begins by fitting a quadratic function to tumor growth data. This is the classic logistic growth model, and we are attempting to find the best fitting quadratic function passing through the origin. You begin this problem by plotting the data in the second and third and fifth and sixth columns. Use Trendline as you have before to get the best quadratic through the data. Be sure to set the y-intercept equal to zero in the options and watch that you have at least 4 significant figures from the equation. Use that equation to compute equilibria (zero growth), find the vertex (maximum growth), and determine the sum of square errors between the model and the data.

The next part of the problem has you define variables a and b for the Gompertz model. The values you choose are not very sensitive, but taking a = 0.1 and b = 0.2 seems to work fine. You define the sum of square errors between the data and the Gompertz function, then use Solver as you have before to minimize this sum of square errors. You should graph the Gompertz function with about 50 evenly spaced points along the domain of your function (the p values). This is added to the graph from Part a.

The last two parts have you simulate the nonlinear discrete dynamical model starting with a variable that you define, p0. Again you find the sum of square errors between the time series data given in the table and the simulation of your model using either the logistic or Gompertz growth function. With Solver you minimize the sum of square errors by only changing p0. The rest of the questions that you need to answer are fairly routine.

Problem 3: This question examines two optimization problems from classical Calculus that involve trigonometric functions and their understanding. In each case, a diagram is provided in the lab to help you visualize the cross-sectional area of the trough. The cross-sectional area is a trapezoid. (You may want to review the area of a trapezoid in wikipedia.) For Trapezoid A, you can readily find the value of x, then the height of the trapezoid is associated with cos(q). Use your definition of cos(q) to find height as a function of x and cos(q). One base of the trapezoid is x, while the other is x plus the other legs of the triangular regions on either side of the central rectangle. The legs of these triangular regions are associated with sin(q), so use the definition of sin(q) to find the length of these legs as a function of x and sin(q). This will give you the formula for the cross-sectional area, then the volume for this trough is simply the area of the trapezoid times the length of the metal strip. Clearly, the only variable in the area function is q, so graphing the area on the given domain is relatively simple. The optimization of this volume is a problem that you can do either on paper or by taking advantage of Maple and its ability to differentiate.

For Trapezoid B, you cannot simply divide x into three equal pieces. You must start by finding both x(q) and the height h(q) (taking advantage of the definitions of the trigonometric functions). After doing this, you use the formula for the area of a trapezoid and to find the volume of the trough. In this case, the volume of Trough B is a more complicated expression of q, so you will almost certainly want to use Maple to optimize the volume.

 

Copyright © 2014 Joseph M. Mahaffy.