dx1/dt=-d*x1;
dx2/dt=kTL*x1-b*x2;
y1=x2;

