Solving a system of coupled differential equations with dsolve_system

Question

I want to solve a system of 4 coupled differential equations:

eqs = [Eq(cP1(t).diff(t), k1*cE1(t)**3), Eq(cE1(t).diff(t), -k1 * cE1(t)**3 + k6 * cE3(t)**2), Eq(cE2(t)