:(function (t, u, du)
        du[1] = 0 - ((0.1 * (u[2] + u[6] + u[7]) * (u[1] + u[8] + u[9])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]) + (0.3 * (u[4] + u[6] + u[8]) * (u[1] + u[2] + u[3])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]))
        du[2] = (0.1 * (u[2] + u[6] + u[7]) * (u[1] + u[8] + u[9])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]) - (0.2 * (u[2] + u[6] + u[7]) + (0.3 * (u[4] + u[6] + u[8]) * (u[1] + u[2] + u[3])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]))
        du[3] = 0.2 * (u[2] + u[6] + u[7]) - (0.3 * (u[4] + u[6] + u[8]) * (u[1] + u[2] + u[3])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9])
        du[4] = (0.2 * (u[2] + u[6] + u[7]) + (0.3 * (u[4] + u[6] + u[8]) * (u[1] + u[2] + u[3])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9])) - 0.4 * (u[4] + u[6] + u[8])
        du[5] = (0.4 * (u[4] + u[6] + u[8]) + 0.2 * (u[2] + u[6] + u[7])) - 0
        du[6] = ((0.1 * (u[2] + u[6] + u[7]) * (u[1] + u[8] + u[9])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]) + (0.3 * (u[4] + u[6] + u[8]) * (u[1] + u[2] + u[3])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9])) - (0.2 * (u[2] + u[6] + u[7]) + 0.4 * (u[4] + u[6] + u[8]))
        du[7] = ((0.1 * (u[2] + u[6] + u[7]) * (u[1] + u[8] + u[9])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]) + 0.4 * (u[4] + u[6] + u[8])) - 0.2 * (u[2] + u[6] + u[7])
        du[8] = (0.3 * (u[4] + u[6] + u[8]) * (u[1] + u[2] + u[3])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]) - ((0.1 * (u[2] + u[6] + u[7]) * (u[1] + u[8] + u[9])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9]) + 0.4 * (u[4] + u[6] + u[8]))
        du[9] = 0.4 * (u[4] + u[6] + u[8]) - (0.1 * (u[2] + u[6] + u[7]) * (u[1] + u[8] + u[9])) / (u[1] + u[2] + u[3] + u[4] + u[5] + u[6] + u[7] + u[8] + u[9])
    end)