The concentration of a drug in the body $Cp$ can be modeled by the equation:

$Cp=\frac{D_G}{V_d}\frac{k_a}{(k_a-k_e)}(e^{-k_{e}t}-e^{-k_{a}t})$

Where $D_G$ is the dosage administrated $(mg),V_d$ is the volume of distribution $(L),k_e$ is the absorption rate constant $(h^{-1}),k_a$ is the elimination rate constant $(h^{-1}),$ and $t$ is the time $($h$)$ since the drug was administered.

Create a function that bring as output the concentration of a drug in the body at time required. For a certain drug, the following quantities are given:

$D_{G}150mg,V_d=50L,k_e=0\mathrm{.}4h^{-1},k_a=1\mathrm{.}6h^{-1}$

Write as a function introduction in the comment space for help tool:

The function input variables

The function output variables

The purpose of the function

Add inside the function algorithm the command nargin

Nargin $($number of argument in$),$ is very useful to verifying if the user enter the correct quantity of argument in the function. As example if the function need $3$ argument add in the function:

if nargin $3$ errordlg$(\text{'}$Not enough input arguments'$)$ end