The effects of thermal conductivity which depend on temperature are conversely proportional with the linear function of temperature on free convective flow where the fluid is viscous and incompressible along a heated uniform and the vertical wavy surface has been examined in this study. The boundary layer equations with the associated boundary conditions that govern the flow are converted into a nondimensional form by using an appropriate transformation. In the domain of a vertical plate that is flat, the resulting method of nonlinear PDEs is mapped and then worked out numerically by applying the implicit central finite difference technique with Newton’s quasilinearization method, and the block Thomas algorithm is well known as the Keller-box method. The outputs are obtained in the terms of the heat transferring rate, the frictional coefficient of skin, the isotherms, and streamlines. The outcomes showed that the local heat transferring rate, the local skin friction coefficient, the temperature, and the velocity all are decreasing, and both the thermal layer of boundary and velocity become narrower with the rising values of reciprocal variation of temperature-dependent thermal conductivity. On the other hand, the friction coefficient of skin, the velocity, and the temperature decrease where the friction coefficient of skin and velocity decrease by 43% and 64%, respectively, but the heat transfer rate increases by 61% approximately, and both the boundary layer thermal and velocity become thinner when the Prandtl number increases.