The present study deals with the Darcy–Brinkman–Forchheimer model for bioconvection-stratified nanofluid flow through a porous elastic surface. The mathematical modeling for MHD nanofluid flow with motile gyrotactic microorganisms is formulated under the influence of an inclined magnetic field, Brownian motion, thermophoresis, viscous dissipation, Joule heating, and stratifi-cation. In addition, the momentum equation is formulated using the Darcy–Brinkman–Forchheimer model. Using similarity transforms, governing partial differential equations are reconstructed into ordinary differential equations. The spectral relaxation method (SRM) is used to solve the nonlinear coupled differential equations. The SRM is a straightforward technique to develop, because it is based on decoupling the system of equations and then integrating the coupled system using the Chebyshev pseudo-spectral method to obtain the required results. The numerical interpretation of SRM is admirable because it establishes a system of equations that sequentially solve by providing the results of the first equation into the next equation. The numerical results of temperature, velocity, concentration, and motile microorganism density profiles are presented with graphical curves and tables for all the governing parametric quantities. A numerical comparison of the SRM with the previously investigated work is also shown in tables, which demonstrate excellent agreement.