Extremal Lyapunov exponents for bilinear control systems with constrained control values are computed numerically by solving discounted optimal control problems. Based on this computation a numerical algorithm to calculate stabilizing control functions is developed.