Abstract: An efficient method based on Singular Value Decomposition (SVD) for the computation of Lyapunov exponent is presented. This method is applicable not only to discrete systems but also to continuous systems, and it does not require frequent reorthogonalization. Numerical calculations show that the method is stable and has a higher computing efficiency, and the results are credible.