function [x,u]=conditions_initialles(xmax, nPt)
dx=2*xmax/(nPt-1);
x=-xmax:dx:xmax;
u=0.5*ones(size(x));
u(x<0)=1;
u(x>0)=0;