option = odeset('RelTol', 0.0005);
s = -2;
[x,y] = ode15s('rhs',[1 2], [1 s],option);
[m,n] = size(y);
y(m,1)