A method and system for optimizing post prandial blood glucose excursions for persons with diabetes via computation of time-varying optimal insulin profiles are disclosed. The method builds upon a mathematical model based on the physiologic description of the human glucoregulatory system. Model parameter values are specific for each patient and meal type, and are determined from measured data by minimizing differences between both measurements and model predictions, and which are used to provide the optimal insulin profiles that are useful to diabetics in managing their chronic disease.