Dynamic Engineering Systems and Mathlab assignment

Drexel University, College of Engineering 2015-2016 Academic Year

Drexel University Office of the Dean of the College of Engineering ENGR 232 – Dynamic Engineering Systems

Week 1 Laboratory Assignment

In this lab we will investigate the use of direction fields to evaluate solutions to first order differential equations.

For this lab, we will study these qualitative analysis tools via the example of population growth using the Gompertz equation as a model. This equation is used to model the growth of tumor cells, and gives an expression of the dependence of tumor size on time.

The differential equation governing the growth is:

Here is the growth rate is a reflection of the carrying capacity of the tumor.

Let ; ; and V range from 0 to 100

1. Plot the graph of vs. . Note you may have to use the vectorization operator “.” to evaluate the function of the differential equation.

2. Determine the equilibrium points (by-hand or MATLAB) and plot them on the same graph as part (1) using red circles. Use a legend to indicate the equilibrium points and the original plot.

3. Using the gradient() function, plot the derivative of versus on the same plot. That is, plot vs. .

4. Determine the stability of each equilibrium point and the region of attraction. Write these answers as comments in your code. Comment on the stability as it related to the graph.

5. From this graph, are we able to tell how long the tumor will take to grow to units if the initial size is units? Why? Write your answers as comments in the code.

We will now use the direction field tool to examine the solution.

6. Download the file from the course website in the MATLAB Resources folder. You may have to right-click and save the file to your computer. IMPORTANT! Make sure it is in the current directory you are working in.

OR

Access the dfield tools at http://math.rice.edu/~dfield/dfpp.html

7. Using the same parameters as above, plot the direction field of this differential equation. Set the time and variable limits appropriately.

8. Compare your results from part (4) to the direction field obtained to check: the stability of each equilibrium point and the region of attraction. Include these as comments in your code.

9. Estimate how long it will take for if. Hint: You can go to: WINDOW DFIELD Keyboard Input to enter the initial conditions.

10. Repeat the direction field plot with and estimate again how long it takes for if . Comment on your results as they compare to part (9).