function rhs = heat_2d_rhs(tspan,u,dummy,B,N,kappa)
rhs=kappa*B*u;