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