Post on 20-Jan-2017
PRACTICA MODELOS COMPARTIMENTALES Nombre: Germán Moreno Cuadros Fecha: 07-11-2016
Resolución de tareas.
TAREA 1.
A.
Dose =10 mg
C(t)
Vd
K= 0.28
B. R1= 0
R2= -K*Vd*C
C. = ₀ e
D. = ln₀
( )
=,
∗ = 6,4
Tarea 2.
Resolviendo con Matlab
clear all clf %PARAMETROS DEL EJERCICIO Vd = 210; %volumen de distribución [Litros] Dose = 10.*1000; % Un microgramo = 1000 miligramos Cinit= Dose./Vd; % concentración inicial del fármaco k = 0.28; %Eliminación por hora del fármaco tmax = 12; %Rango de tiempo en horas sobre el cual se va a resolver el modelo % Expresión de flujos Ri = @(t,C) 0; Re = @(t,C) k*C*Vd; % Eflux %resolviendo la ecuación diferencial sol=ode45(@Pharm,[0 tmax], Cinit,[],Vd,Ri,Re); %Graficando los resultados tstep = 0.01; tspan = 0:tstep:tmax; y = deval(sol,tspan); plot(tspan,y);grid; xlabel('tiempo (horas)') ylabel('concentracion del fármaco microgramos/litro') %Encontrando el tiempo en que la concentración en plasma va por debajo el umbral InPain = find(y<8,1,'First'); TimeForNextDose = tspan(InPain); disp(['El dolor será insoportable después de',num2str(TimeForNextDose),'hours']);
A. 10mg; 6,38 horas
B. Con 5 mg 3,9 Horas
Con 15 mg 7,82 horas
Tarea 3.
Paciente con línea intravenosa.
A. Gráfico
Ri=0 Infusión de droga constante en el tiempo
Vd
C(t)
RE= K*Vd*C
B. Código Matlab modificado.
clear all hold all %mantiene el gráfico del primer ejercicio %PARAMETROS DEL EJERCICIO Vd = 210; %volumen de distribución [Litros] Dose = 2.5*1000; % Por hora. Cinit= 0; % concentración inicial del fármaco será 0 k = 0.28; %Eliminación por hora del fármaco tmax = 12; %Rango de tiempo en horas sobre el cual se va a resolver el modelo % Expresión de flujos Ri = @(t,C) Dose;%influx Re = @(t,C) k*Vd*C; % Eflux %resolviendo la ecuación diferencial sol=ode45(@Pharm,[0 tmax], Cinit,[],Vd,Ri,Re); %Graficando los resultados tstep = 0.01; tspan = 0:tstep:tmax; y = deval(sol,tspan); plot(tspan,y);grid; xlabel('tiempo (horas)') ylabel('concentracion del fármaco microgramos/litro') %Encontrando el tiempo en que la concentración en plasma va por debajo el umbral InPain = find(y>8,1,'First'); TimeForNextDose = tspan(InPain); disp(['Dolor inaceptable despues de ',num2str(TimeForNextDose),' horas']);
TAREA 4 C₀= Dose 30mg = D
b=0,24 C2
V1 V2
K= 0,28
Ecuaciones:
22
= 1 + (− 2 2)
1 = ₀
22
= ₀ − 2 2
=
22
= − 2 2
Usando Matlab.
clear all hold all %mantiene el gráfico del primer ejercicio %PARAMETROS DEL EJERCICIO Vd = 210; %volumen de distribución [Litros]
Dose = 30*1000; % Por hora. Cinit= 0; % concentración inicial del fármaco será 0 b=0.24;% Salida de farmaco de Compartimento al compartimento 2 por hora k = 0.28; %Eliminación por hora del fármaco por hora
tmax = 12; %Rango de tiempo en horas sobre el cual se va a resolver el modelo % Expresión de flujos Ri = @(t,C) b.*Dose*exp(-b.*t);%influx Re = @(t,C) k*C*Vd; % Eflux
%resolviendo la ecuación diferencial sol=ode45(@Pharm,[0 tmax], Cinit,[],Vd,Ri,Re); % Gráficamente tstep = 0.01; tspan = 0:tstep:tmax; y = deval(sol,tspan); plot(tspan,y);grid; xlabel('tiempo (horas)') %Etiqueta del eje X ylabel('concentracion del fármaco microgramos/litro') %Etiqueta del eje Y
%Encontrando el tiempo en que la concentracion en plasma va por debajo el umbral InPain = find(y>8,1,'First'); TimeForNextDose = tspan(InPain); disp(['Dolor inaceptable despues de ',num2str(TimeForNextDose),' horas']); disp(' '); disp('fin de practica de Modelos Compartimentales');
Inyección Vía oral Vía intravenosa