dx1/dt=-d1*x1-d2m0*x1*x3;
dx2/dt=kTLm0*x1-b*x2;
dx3/dt=d3dm0*e0-d3*x3-d2tm0*x1*x3;
y1=x2;
x1(0)=1;
x2(0)=0;
x3(0)=1
