drugabsorp.m

ch9例题05_ode23s口服药在人体中的吸收 matlab 程序

function DrugAbsorp    
% 口服药在人体小肠中的吸收

clear all; clc;

Dn = 0.6;  Do = 0.08;  An = 10;
[t,y] = ode23s(@ODEsfun,[0:0.0001:1],[1  1],[],Dn,Do,An);
index = find(y(:,1)<0);
y(index,1) = 0;

Fa = 1- y(end,1)^3 - y(end,2)/Do;
fprintf('\n 药物分率预计值 Fa =%6.4f\n\n',Fa)

figure(1)
comet(t(:),y(:,1));
figure(2)
comet(t(:),y(:,2));
figure(3)
plot(t,y);
title('药物颗粒半径与浓度分布变化图');
legend('药物颗粒半径变化','药物浓度分布',4);

% ------------------------------------------------------------------
function dydz = ODEsfun(t,y,Dn,Do,An)
rstar = y(1);
Cstar = y(2);
f1 = -Dn/3*(1 - Cstar)/rstar;
f2 = Dn*Do*rstar*(1 - Cstar) - 2*An*Cstar;
dydz = [f1; f2];