<p>This paper addresses an optimal control problem (OCP) related to the Pennes bioheat equation, incorporating non-classical boundary conditions. To tackle this challenge, we use a discretization approach, introducing approximations for both the state and control functions. These approximations are constructed using orthonormal Bernstein basis functions and the operational matrix of differentiation corresponding to the selected basis functions. Subsequently, these approximations are incorporated into the performance index. By applying the necessary conditions for optimization, we transform the main problem into the solution of a nonlinear system of algebraic equations. Notably, our proposed scheme stands out among other direct methods for solving OCPs, as it ensures accurate satisfaction of the given initial and boundary conditions within the constraints of the problem. We discuss the convergence of the method and present various numerical examples to demonstrate the effectiveness of our approach.</p>