function rhs = heatrhs(tspan,u,dummy,A,dx)
rhs=.1*A*u;