dx1/dt = -1/2500*x1 + 1/2000000*x3; dx2/dt = e2a*x7 + 1/2*x13*x11 - x5*x2 - 1/50000*x2; dx3/dt = 5*i1*x11 - 1/2*x14*x3; dx4/dt = -1/2500*x4 - 1/2000000*x3; dx5/dt = -1/5*x13*x5 - x5*x2 - x5*k,deg - x5*k2*u1*x6 - x5*k3 + x12*u1*k1 + x9*t1 + x8*t2; dx6/dt = 1/2*x4 - c5*x6; dx7/dt = -5*e2a*x7 + 1/2*x14*x3; dx8/dt = x5*x2 - x8*t2; dx9/dt = 1/5*x13*x5 - x9*t1; dx10/dt = 1/2000000*x3 - c3a*x10; dx11/dt = -1/2*x13*x11 + 1/50000*x2 - i1*x11 - x8*t2; dx12/dt = -x12*k,deg - x12*u1*k1 + k,prod; dx13/dt = -1/5*x13*x5 - x13*i1a - 1/2*x13*x11 - 1/10000*x13 + 1/2000*x14 + x10*c4a; dx14/dt = 5*x13*i1a - 1/2*x14*x3 - 1/400*x14; dx15/dt = x5*k2*u1*x6 + x5*k3 - x15*k,deg; y1 = x5 + x15 + x12; y2 = x10; y3 = x13 + x2; y4 = x3; y5 = x4; y6 = x5