In this chapter, we study a problem of gradient optimal control for a bilinear reaction–diffusion equation evolving in a spatial domain Ω⊂Rn using distributed and bounded controls. Then, we minimize a functional constituted of the deviation between the desired gradient and the reached one and the energy term. We prove the existence of an optimal control solution of the minimization problem. Then this control is characterized as solution to an optimality system. Moreover, we discuss two special cases of controls: the ones are time dependent, and the others are space dependent. A numerical approach is given and successfully illustrated by simulations.