Forced convection of Williamson fluid flow in porous media under constant surface heat flux conditions is investigated numerically. A model of Darcy–Forchheimer–Brinkman is used and the corresponding governing equations are expressed in dimensionless forms and solved numerically using bvp4c with MATLAB package. Boundary layer velocity, shear stress, and temperature profiles, in addition to the local Nusselt number parameter over a horizontal plate, are found. The effects of the Forchheimer parameter, Nusselt number, Darcy parameter, porous inertia, and Williamson parameter on the velocity profiles, temperature profiles, coefficient of friction, and coefficient of heat transfer are investigated. The results showed that as the Darcy parameter increases, boundary layer velocity and shear stress increase, while the temperature and Nusselt number decrease. In addition, as Williamson's parameter increases, velocity within the boundary layer, shear stress, and Nusselt number decrease while the temperature profile increases. Also, with larger values of the Forchheimer parameter, the velocity of the boundary layer, shear stress, temperature, and Nusselt number increase. Furthermore, the Nusselt number and the coefficient of friction are obtained on the surface of the horizontal plate.