A non-classic optimality condition and numerical algorithm for smooth boundary controls in semi-linear first-order hyperbolic systems are presented. Additional integral control constraints are considered. The suggested approach is based on special variations of admissible continuously differentiable controls. These variations lead to a new necessary optimality condition which is a base for a numerical method. The method is applied to an optimal control problem for an age-structured population.