I would like to solve two nonlinear ode equations with variables coefficients with initial conditions. Below given my exact problem. This code is not giving answer.
odes_systems = {I1'[ t] == (1 - k1*t + ((k1*t)^2)/2)*(b11*I1[t] + b21*I2[t]) - (r1 + d1)*I1[t],
I2'[t] == (1 - k2*t + ((k2*t)^2)/2)*(b12*I1[t] + b22*I2[t]) - (r2 +d2)*I2[t]};
initial_condition = {I1[0] == 1, I2[0] == 1};
sol = DSolve[Join[eqn, iv], {I1, I2}, t]
where t is variables. Would be grateful if I get this solutions.