Hello,
you can append this code:
sNew[time_] := (A2*B[t] + A1 - rho)/eta /. {rho -> 0.96, eta -> 2} /.
s[[1]] /. t -> time;
Plot[{B[k] /. s, sNew[k]}, {k, 0, 5}]
or
Plot[{B[t] /.
s, (A2*B[t] + A1 - rho)/eta /. {rho -> 0.96, eta -> 2} /.
s[[1]]}, {t, 0, 5}]
Good feature is replace endTime=5 for
endTime = s[[1, 1, 2, 1, 1, 2]];
because you set real solved time.
M. M.