∫ (ln(1-q*t/m_0)+1)^2 dt