The microscale gliders are regularly affected by the local surrounding environment, like liquid rheology and physical (nearby) boundaries. This article focuses on the numerical simulations of bacterial speed over a non-Newtonian slime and its power expenditure. The flow rate generated by the swimmer, slime speed and level curves are also the point of interest. To fulfill the purpose Oldroyd-4 constant model is assumed over a rigid boundary. A complex undulating sheet is approximated as a bacterial surface. Since a slime (present below the undulating sheet) is a non-Newtonian fluid so a modeling approach of peristaltic flow problem is adopted, and dynamic equilibrium conditions are implemented for steady motion. Implicit finite difference method (FDM) is employed to calculate the numerical solution of reduced Boundary value problem (BVP). To compute the flow rate and cell speed, Broyden’s root finding algorithm is integrated with FDM. These computed values are further utilized to perceive the behavior of work done, velocity of slime and streamlines.