A method of calculating an ammonia (NH3) mass generated in a lean NOx trap (LNT) of an exhaust purification device includes sequentially calculating a NH3 mass flow at a downstream of each slice from a first slice to an n-th slice, and integrating the NH3 mass flow at the downstream of the n-th slice over a predetermined time, wherein the calculation of the NH3 mass flow at the downstream of the i-th slice comprises calculating a NH3 mass flow flowing into the i-th slice, calculating a NH3 mass flow generated at the i-th slice, and adding the NH3 mass flow generated at the i-th slice to a value obtained by subtracting the NH3 mass flow used to reduce the NOx and the O2 at the i-th slice from the NH3 mass flow flowing into the i-th slice.