A closed form solution is difficult to obtain due to the non linearity of the objective function (the delay function part) as well as the constraints, hence, the set of equations in 15 can be solved numerically using Newton-Raphson method [16] to obtain the optimal solution.