Numerical studies on many-body localization