Curvas circulares

EJERCICIOS DE PROGRAMACION EN MATLAB

EJERCICIO 1

La concentración de una medicina en el cuerpo, puede ser representada con la siguiente ecuación.

Cp=DgVd*KaKa-Ke*(e-ket-e-kat)
Donde:

Dg= Dosis administrada en mg.
Vd = Volumen distribuido en lt.
Ka = Constante de absorción (h-1)
Ke = Constante de eliminación (h-1)
t = tiempo en horas desde que la medicina fue administrada

Parauna determinada medicina se dan los siguientes valores:

Dg = 150 mg.
Vd = 50 lt.
Ka = 1.6 (h-1)
Ke = 0.4 (h-1)

a) Una dosis simple es administrada en t = 0. Calcular y representar en ungráfico Cp vs t (10 horas).

b) Una dosis es administrada en t = 0. Cuatro dosis posteriores son administradas en intervalos de cuatro horas. Calcular y representar un gráfico Cp vs t (24horas).

% GRÁFICAS CP VS T
dg=input('Ingrese el valor de la dosis administrada en mg =');
vd=input('Ingrese el valor del volumen administrado en litros =');
ka=input('Ingrese el valor de la constantede absorción en h^-1 =');
ke=input('Ingrese el valor de la constante de eliminación en h^-1 =');
t1=10;
e=2.71828182846;
n=0;
for k=0:0.1:t1
cp=(dg/vd)*(ka/(ka-ke))*(e^(-ke*k)-e^(-ka*k));n=n+1;
a(1,n)=k;
b(1,n)=cp;
end
figure(1)
plot(a,b)
xlabel('tiempo ')
ylabel('cp')
m=0;
t1=0;
t2=0;
t3=0;
t4=0;
for h=0:0.1:24
if h>=4
t1=h-4;
endif h>=8
t2=h-8;
end
if h>=12
t3=h-12;
end
if h>=16
t4=h-16;
end
cp0=(dg/vd)*(ka/(ka-ke))*(e^(-ke*h)-e^(-ka*h));cp1=(dg/vd)*(ka/(ka-ke))*(e^(-ke*t1)-e^(-ka*t1));
cp2=(dg/vd)*(ka/(ka-ke))*(e^(-ke*t2)-e^(-ka*t2));
cp3=(dg/vd)*(ka/(ka-ke))*(e^(-ke*t3)-e^(-ka*t3));
cp4=(dg/vd)*(ka/(ka-ke))*(e^(-ke*t4)-e^(-ka*t4));cp=cp0+cp1+cp2+cp3+cp4;
m=m+1;
c(1,m)=h;
d(1,m)=cp;
end
figure(2)
plot(c,d)
xlabel('tiempo ')
ylabel('cp')

SOLUCIÓN

Ingrese el valor de la dosis administrada en mg =150...