This function solves a set of delay differential equations (DDEs) that represent a target cell limited model. It takes as input the DDE model parameters and an input file containing experimental data. It plots the numerically integrated solution of the DDE and the sum of squared residuals between experimental data and model solution. It then takes upper bounds and lower bounds on some model parameters and calculates (using a local search procedure) the optimal best-fit parameters that explain the data.