Hello,

I'm trying to solve the following ode system:

odes := diff(a0(t), t) - diff(a1(t), t) + diff(a2(t), t) = 0, diff(a0(t), t) - diff(a2(t), t) - 16*a2(t) = 0, diff(b0(t), t) + diff(b1(t), t) + diff(b2(t), t) = 0, diff(b0(t), t) - diff(b2(t), t) - 48*b2(t) = 0, diff(a1(t), t) + 4*diff(a2(t), t) - 3*diff(b1(t), t) + 12*diff(b2(t), t) = 0, diff(a0(t), t) + diff(a1(t), t) + diff(a2(t), t) - diff(b0(t), t) + diff(b1(t), t) - diff(b2(t), t) = 0;
   
ics := a0(0) - a2(0) = 0, b0(0) - b2(0) = 1, a0(0) - a1(0) + a2(0) = 0, b0(0) + b1(0) + b2(0) = 1, a1(0) + 4*a2(0) - 3*b1(0) + 12*b2(0) = 0, a0(0) + a1(0) + a2(0) - b0(0) + b1(0) - b2(0) = 0;

dsdef := dsolve({ics, odes}, numeric);

There are 6 ordinary differential equations with 6 initial conditions. The functions are named a0(t), a1(t), a2(t), b0(t), b1(t), b2(t).

The question is that Maple solve the system analliticaly (exact dsolve) but when I run the "numerical dsolve" it returns: "Error, (in DEtools/convertsys) invalid specification of initial conditions", and I can't recognize the error/s.

Thanks.