Well if we account for the varying mass of the rocket, we use the equation m.dv = -v_g.dm
where m is the mass of the rocket (variable) and v_g is the velocity of the gases (500) and v is the velocity of the rocket (variable)
dividing by m and integrating you get v_f - v_i = v.ln(m_f/m_i)
throw in your numbers and you can find v_f.
Otherwise if we assume constant mass then you can just use the equation T = v_g dm/dt
where dm/dt is the rate of gas emission (100)
Equate this to F = ma, change a to (v_f - v_i)/t and solve for v_f