s_correct = fzero('findx',[-2 -1])
%
options = odeset('RelTol', 0.0005);
[x,y] = ode45('rhs',[1 2], [1 s_correct],options);
[m,n] = size(y);
xplot=x(:,1);
yplot=y(:,1);
plot(xplot,yplot)
axis([1 2 0 1]);