BENG449 Problem 8 4

BENG449 Problem 8 4 - 4/11/08 12:55 AM MATLAB Command...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon
4/11/08 12:55 AM MATLAB Command Window 1 of 3 >> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Alex Lemon % % BENG449 -- Biomedical Data Analysis % % Problem Set 8, Question 4 % % Due April 11, 2008 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Consider a model with three compartments: the absorption site (A), the % body (B) and the elimination site (E) with the following rate constants: % k0 : A --> B % km1 : B --> A % k1 : B --> E % Tabulate the eigenvalues and simulate the system. % Initialize the parameters clear A B E; k0 = 0.01; k1 = 0.035; KM1 = [0, 0.01, 0.035, 0.1]; c0 = [1; 0; 0]; N = length(KM1); S = length(c0); t = [0:250]; T = length(t); for n = 1:N % Form the transition matrix km1 = KM1(n); K = [-k0 km1 0; k0 -(k1 + km1) 0; 0 k1 0]; % Compute the eigenvalues and eigenvectors [V L] = eig(K); % Print the eigenvalues fprintf('\n%s%.6f%s\n%s%10.6f\n%s%10.6f\n%s%10.6f\n', . .. 'The eigenvalues for km1 = ', km1, ' are:', . .. ' lambda_1: ', L(1,1), .
Background image of page 1

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
Image of page 2
This is the end of the preview. Sign up to access the rest of the document.

This note was uploaded on 07/19/2008 for the course BENG 449 taught by Professor Richardcarson during the Spring '08 term at Yale.

Page1 / 6

BENG449 Problem 8 4 - 4/11/08 12:55 AM MATLAB Command...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online