Scribd
Upload
53 views2 pages

Math 463: Introduction To Mathematical Biology Computer Lab Assignment #2

This document provides instructions for a MATLAB computer lab assignment on solving and plotting discrete equations. Students are asked to write MATLAB code to iterate and plot the solutions to 6 discrete models. They are also asked to simulate and discuss the effects of parameter changes on the solutions to 2 additional biological models.

Uploaded by

gs2002us
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
53 views2 pages

Math 463: Introduction To Mathematical Biology Computer Lab Assignment #2

This document provides instructions for a MATLAB computer lab assignment on solving and plotting discrete equations. Students are asked to write MATLAB code to iterate and plot the solutions to 6 discrete models. They are also asked to simulate and discuss the effects of parameter changes on the solutions to 2 additional biological models.

Uploaded by

gs2002us
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

Math 463: Introduction to Mathematical Biology

Computer Lab Assignment #2

The goal of this lab is to give you some practice solving discrete equations with
MATLAB. Your task is to write some simple code that will iterate the discrete models
below.
• Login and open MATLAB

• Click on the ’New M-file’ Icon in the top left corner of the MATLAB menu. It
looks like a blank sheet of paper.

• Type clear (return)

• Set your initial conditions, e.g. ”x(1) = 1;” Note that MATLAB does not allow
zero indices so you’ll have to shift everything up by one. The semicolon at the
end of the command tells MATLAB to perform the command, but not to print the
result to the screen.

• Define the number of iterations you want to compute: e.g. N = 15;

• Write a for-loop to compute the iterates

• Click on the ’Save’ Icon to give your file a name and save it; e.g. filename.m

• Type ’filename’ in the MATLAB command window to run your program.

For example, to solve the simple model xn+1 = −.5xn , enter the following code:
clear
x(1) = 10;
N = 20;
for n = 1:N
x(n+1) = -.5*x(n);
end
n = 1:1:N+1’;
plot(n,x,n,x,’bo’)

Use MATLAB to solve the following discrete models and plot your results to turn in.
The subplot(i,j,k) command allows you plot several graphs on the same sheet of paper.
For example:

subplot(2,2,1), plot(n,x)

breaks up your paper into 2 rows and 2 columns of graphs. The number 1 tells MATLAB
to plot the current graph in the first position.
Use the subplot command to plot problems 1-4 on the same sheet of paper. For
problems 1-4, compute the exact solution and plot it on the same axis as the numerical
solution (eg use circles for the analytical solution and a straight line for the computed
solution). For problems 5 and 6, plot the numerical solutions on the same piece of paper.
What can you say about the eigenvalues of the equations?

HINT: MATLAB cannot deal with zero as an index. This means that you will not
be able to use x0 . Instead, use the given information to start your numerical iterations
at n = 1.

1. xn+1 − 5xn + 4xn−1 = 0, x1 = 9, x2 = 33


2. xn − xn−2 = 0, x1 = 3, x2 = 5
3. xn+2 − 2xn+1 + xn = 0, x0 = 1, x1 = 3
4. xn+2 − xn+1 + xn = 0, x0 = 6, x1 = 2
5. xn+1 = xn /4 + 3yn , yn+1 = −xn /8 + yn x0 = 0, y0 = 1
6. xn+1 = xn + yn , yn+1 = 2xn x0 = 1, y0 = 0

7. Simulate the Segmental Growth Problem: an+1 − (1 + q)an − ran−1 = 0, a0 = 1 try


various values for q and r. How does changing these values effect the solution.

8. Simulate the Blood C02 problem: Cn+1 − Cn + αβCn−1 = m for relevant values of m,
α and β. HINT: Be sure to try 4αβ > and 4αβ < 1.

You might also like