function heat_1d_rhs = heat_1d_rhs(tspan,u,dummy,kappa,dx,A)
heat_1d_rhs=(kappa/dx^2)*A*u;