I have define the following ODE system

restart;
with(LinearAlgebra);
dC1 := diff(C1(t), t) = -2*a1*C1(t)+b1*C2(t);
dC2 := diff(C2(t), t) = 2*a1*C1(t)-(b1+a1)*C2(t)+2*b1*C3(t);
dC3 := diff(C3(t), t) = a1*C2(t)-(2*b1+a2)*C2(t)+b2*O1(t);
dO1 := diff(O1(t), t) = a2*C3(t)-b2*O1(t);
syst := dC1, dC2, dC3, dO1;ics := C1(0) = 1, C2(0) = 0, C3(0) = 0, O1(0) = 0;
func := {C1(t), C2(t), C3(t), O1(t)};
sol := dsolve([syst, ics], func);

After evaluating...

I have a matrix A4  

When solving a ODE system, I get "the lenght of the output exceeds the limit of 100000". This was expected so, it is OK with me. However, I would like to be able to see part of the output. Then, I can start making susbtitutions (algsubs) and shrink the output. So, how can I print just a few lines of the output?

I am having trouble getting an algebraic solution for a first order four-ODE system. The problem seems trivial to me. However, Maple is taking a very long time to solve it. The system describes a reversible "chemical" reaction that have 4 states. Basically, the system is

C1<->C2<->C3<->O1

where the rate constants are ai towards the right and bi towards the left.

Any suggestion will do.

Thanks!