function retval = f181(x,t)
% usage: retval = f181(x,t)
% description: returns source for model parabolic equation for MOL

%retval = 6*x;
retval = 4*x.*(1-x);

end
