function u0 = InitialCond(y)
u0=[y;0;0;y*0+10^(-3);0;0;0;0;0;sqrt(1/(1+exp(-1)));0;0;0+0*(2.5*(y-0.5))^3+0*20*(y-0.5)^2-0*3;1];
end

