function rhs = specrhs(t,ut,dummy,k,kappa)
rhs=-kappa*k'.^2.*ut;