This paper presents a comprehensive computer simulation of hydrogenated amorphous p-i-n silicon solar cells. The physical mechanism governing solar cell operation has been modeled and solved numerically by Runge-Kutta-Gill method. Effects of gap state density, dopant impurity, diffusion length and interface recombination velocity on solar cell performance are investigated. Numerical results show that the electric field in i-region is not uniform but depends strongly on voltage and position. A rather poor fill factor may be due to the electric field variation and short diffusion length. It is found out that the life time should be improved in order to increase a fill factor and a conversion efficiency.