
function u0=demo4ic(x)

%u0(1) = exp(-x^2) ;
u0(1) = exp(-0.25*x^2) * sin(5*x) ;

u0(2) = 0 ; 
