function fy = CNDemoFun(x)
%
% test function for Crank Nicolson
%
fy = 10*(x.*((1-x).^2));



