Equation of heat subdiffusion with neumann boundary condition