drugabsorp.m

来自「《实用化工计算机模拟：MATLAB在化学工程中的应用 》这本书光盘里的程序~」· M 代码 · 共 33 行

function DrugAbsorp    
% 预测口服药在人体中的吸收
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/06/18 $

clear all
clc

Dn = 0.6;
Do = 0.08;      % Fa = 0.4869

An = 10;
[z,y] = ode23s(@Equations,[0 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
disp(' Results:')
disp('    t         r       C')
disp([z,y])
plot(z,y(:,1),'k--',z,y(:,2),'b:')
legend('r^*','C^*')

% ------------------------------------------------------------------
function dydz = Equations(z,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];