def _expm_multiply_interval_core_0(A, X, h, mu, m_star, s, q):
for k in range(q):<EOL><INDENT>X[k+<NUM_LIT:1>] = _expm_multiply_simple_core(A, X[k], h, mu, m_star, s)<EOL><DEDENT>return X, <NUM_LIT:0><EOL>
A helper function, for the case q <= s.