برای بررسی و مطالعه پیرامون جذب ، توزیع و دفع دارو در بدن نیاز به مدل سازی فارماکوکینتیکی می باشد. این مدلها بدن را به یک یا چند بیوراکتور با اختلاط کامل تشبیه می کنند. این شبیه سازی منجر به یک یا چند معادله دیفرانسیلی می گردد که با حل آنها غلظت دارو در بدن بر حسب زمان محاسبه می شود. حل تحلیلی، یکی از روشهای حل این معادلات می باشد و با استفاده از روشهای ترسیمی متداول پارامترهای فارماکوکینتیکی قابل محاسبه می باشند. از طرفی کاربرد روشهای ترسیمی مبتنی بر حل تحلیلی محدود به مدلهای ساده بودند و فاقد انعطاف پذیری لازم برای مدلهای پیچیده تر می باشد. همچنین تعیین پارامترهای مدل که با استفاده از شیب و عرض از مبداء نمودار بر حسب زمان بدست می آید خالی از خطا نمی باشد. در این پروژه، حل عددی معادلات دیفرانسیل با روش runge-kuttea و تعیین پارامترهای فارماکوکینتیکی با استفاده از یک روش بهینه سازی مبتنی بر اعداد تصادفی پیشنهاد شده است . و از داده های تجربی مطالعات بالینی موجود برای داروی افلوکساسین جهت ارزیابی روش پیشنهادی جدید استفاده شده است . نتایج بدست آمده نشان دهنده این است که روش پیشنهادی دقت و انعطاف بیشتری نسبت به روشهای متداول ترسیمی دارد.