function rhs1=rhs1(tspan,ut,dummy,k)
u=ifft(ut);
rhs1=-i/2*((k').^2).*ut+i*fft( (abs(u).^2).*u );