function y = rhs2(x)
%
%
y = 3.1 * x .* (1-x);