function psi = initmom(s,d)
% d=a - parameter in the decsription of the initial momentum
% s- denote the vector of space points input for the initial momentum
psi=heaviside(s+d)-heaviside(s-d);
% the function psi is 1 between -d and d and zero otherwise