function plotSolutionSurf(sol,y,t)
u1 = sol(:,:,1);
u2 = sol(:,:,2);
u3 = sol(:,:,3);
u4 = sol(:,:,4);
u5 = sol(:,:,5);
u6 = sol(:,:,6);
u7 = sol(:,:,7);
u8 = sol(:,:,8);
u9 = sol(:,:,9);
u10 = sol(:,:,10);
u11 = sol(:,:,11);
u12 = sol(:,:,12);
u13 = sol(:,:,13);
u14 = sol(:,:,14);

%s4=surf(y,t,(0.74.*(1+exp(-u13)).^(-1)));
s4=surf(y,t,u4);
title('alpha(y,t)')
ylabel('t')
xlabel('y')
s4.EdgeColor = 'none';
end