The present invention relates to a method for automatically determining an optimal individual dosing regimen of at least one drug for a patient suffering from a known disease, the optimal individual dosing regimen being optionally subject to at least one clinical constraint, wherein the method comprises the steps of: providing a mathematical model adapted to model a progression of said disease and an effect of the at least one drug on the progression of the disease, the model comprising individual model parameters associated with the patient; utilizing an empirical Bayesian estimation to automatically and numerically estimate the individual model parameters of the mathematical model utilizing patient data associated with the patient; automatically calculating an optimal individual dosing regimen for the mathematical model by solving an optimal control problem based on a desired progression of the disease, the estimated individual model parameters, and an initial guess for the dosing regimen; and adjusting the optimal individual dosing regimen to optionally account for at least one clinical constraint to yield the optimal individual dosing regimen subject to said at least one clinical constraint.