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