An unsteady numerical simulation method for an autorotating rotor in forward flight was developed. The flapping and rotational equations of motion of autorotation are continuously integrated for given time steps, meanwhile the induced velocity field at disc plane is obtained by the dynamic inflow theory embodying the unteadiness. The transitions from arbitrary initial states to equilibrium states were simulated. Steady autorotations as numerical solutions of equations were predicted by using two sources of blade airfoil data. The simulations using airfoil data which were obtained by a two dimensional Navier-Stokes solver in terms of angles of attack and Reynolds numbers have shown good agreements with wind tunnel experimental results.