function u0 = InitialCond(y)
u0=[(-y^2+y);0;0;(-y^2+y);0;0;0;0;0;1;0;0;1;sqrt(1/(1+exp(1)))];
end

