Numerical simulations of early star formation with “soft” equations of state