Page 393 - Control Theory in Biomedical Engineering

Appendices for Chapter 2

function f=equations(t, m)
% t - time wector
% m - cholesterol mass vector: m(1) - cholesterol mass in I compatment,
% m(2)- cholesterol mass in II compartment
global k
global k12
global k21
global mdiet
global min
global mout
global mtis

f=zeros(2,1);
% f(1) equation that describes I compartment
f(1)=k/m(1)+ k21*m(2)-k12*m(1)+ min-mout;
% f(2) equation that describes II compartment
f(2)=-k21*m(2)+k12*m(1)-mtis+ mdiet;

% Run the two comartment model of chlesterol homeostasis
% this code supports the two-compartment model of cholesterol
homeostasis as a teaching tool
% for more: A. Wrona, J. Balbus, O. Hrydziuszko, K. Kubica,
Two-compartment model as a teaching tool for cholesterol homeostasis,
% Adv Physiol Educ 39: 372–377, 2015; doi:10.1152/advan.00141.2014
clear all
clc
global k
global k12
global k21
global mdiet
global min
global mout
global mtis
%Parameters
k=732; % the rate constant of de novo cholesterol synthesis

359