%function [f] = odefun(t, y, B, C)
function [f] = odefun(t, y)
B = 10;
C = 2;
    f = C*y.*(1 - y./B);
%    f = .5*(exp(t) - sin(t) - cos(t));
return
