Abstract
This paper is devoted to establish some theoretical and computational results for a coupled system of drug therapy process. The considered problem is investigated by using the concept of piecewise modeling. We apply piecewise differentiation and integration to study the proposed model. We derive some adequate results for existence and uniqueness for solution to the proposed problems by using fixed point results. Further, we apply some numerical approach based on Newton interpolation formula for equally spaced arguments to find numerical results corresponding to our proposed model. We use some real values and two sets of fractional order to present the dynamics of the considered model graphically.