function udot=gkdv(time,u)
global h delta N alpha

udot=u;

N=length(u);

time
for i=3:N-2

            

              
            udot(i)=-1*(1/(1*h))*(1/1)*(u(i+1)+u(i)+u(i-1))*(u(i+1)-u(i-1))-delta^2*(1/(2*h^3))*(u(i+2)-2*u(i+1)+2*u(i-1)-u(i-2))+alpha*(1/(1*h))*(1/1)*(1/3)*(u(i+1)+u(i)+u(i-1)).^2*(u(i+1)-u(i-1))   ;
 
          udot(1)=0;
             udot(2)=0;
            
            udot(end)=0;
              udot(end-1)=0;
          

end