The optimal control of a class of control unconstrained bilinear systems is considered. It is demonstrated that mere application of known necessary conditions of optimality is insufficient to solve the problem numerically. A technique overcoming this difficulty is proposed. It is demonstrated in an example involving optimal shut-down control of a nuclear reactor under xenon poisoning. The TPBVP is solved numerically using the variation of extremals method.