Typo corrected:
        d(dc_i/dt)/dc_j instead of dc_i/dc_j