In this paper, the infinite time optimal regulation problem for weakly coupled bilinear systems with quadratic performance criteria is obtained by a sequence of algebraic Lyapunov equations. This is the new approach is based on the successive approximation. In particular, the order reduction is achieved by using suitable state transformation so that the original Lyapunov equations are decomposed into the reduced-order local Lyapunov equations. The proposed algorithms not only solve optimal control problems in the weakly coupled bilinear system but also reduce the computation time. This paper also includes an example to demonstrate the procedures.