There are physical phenomena, involving diffusion and structural vibrations, modeled by partial differential equations (PDEs) whose solution reflects their spatial distribution. Systems whose dynamics evolve on an infinite-dimensional Hilbert space, i.e., infinite-dimensional systems, are modeled by PDEs. The aim when designing a controller for infinite-dimensional systems is similar to that for finite-dimensional systems, i.e., the control system must be stable. Another common goal is to design the controller in such a way that the response of the system does not be affected by external disturbances. The controller design for finite-dimensional systems is not an easy task, so, the controller design for infinite-dimensional systems is even more challenging. The backstepping control approach is a dominant methodology for boundary feedback design. In this work, we try with the backstepping design for the boundary control of a reaction-advection-diffusion (R-A-D) equation, namely, a type parabolic PDE, but with constant coefficients and Neumann boundary conditions, with actuation in one of these latter. The heat equation with Neumann boundary conditions is considered as the target system. Dynamics of the open- and closed-loop solution of the PDE system are validated via numerical simulation. The MATLAB®-based numerical algorithm related with the implementation of the control scheme is here included.