Abstract: The simultaneous differential equations on physiological pharmacokinetic model were solved by the Range-Kutta method. Concentration-time profiles in tissues or organs of rat following iv 20mg/kg ethoxybenzamine were simulated.