Multiconfigurational time-dependent Hartree method for bosons: Many-body dynamics of bosonic systems