Extendemos el método de diferencia finita no estándar de solución al estudio de modelos farmacocinéticos-farmacodinámicos. Los modelos farmacocinéticos (PK) se utilizan comúnmente para predecir las concentraciones de medicamentos que impulsan transferencias intravenosas controladas (o infusiones y transferencias orales) mientras que los modelos de interacción farmacocinética y farmacodinámica (PD) se utilizan para proporcionar predicciones de concentraciones de medicamentos que afectan la respuesta de estos medicamentos clínicos. Estructuramos un esquema de diferencia finita no estándar (NSFD) para el sistema relevante de ecuaciones que modela este proceso farmacocinético. Comparamos los resultados obtenidos con métodos estándar. El esquema es dinámicamente consistente y confiable para replicar propiedades dinámicas complejas de los modelos continuos relevantes para diferentes tamaños de paso. Este estudio proporciona asistencia para comprender el comportamiento a largo plazo del medicamento en el sistema y validar la eficiencia del esquema de diferencia finita no estándar como método de elección.