In this paper, we try to approach chaos with numerical method. After investigating nonlinear dynamics (chaos) theory, we introduce Lyapunov exponent as chaos`s index. To look into the existence of chaos in 2-dimensional difference equation, we compute Lyapunov exponent and examine the various behaviors of solutions by bifurcation map.
