Dynare solves for DSGE steady-state values (nonlinear)

At present, the DSGE model and first-order conditions have been solved. There are a large number of nonlinear equations in the steady-state value part, which cannot be solved with MATLAB’s solve function. In the book, I saw that after logarithmic linearization, the BK method or undetermined coefficients are used. But after reading and observing, it seems that the steady-state value is not output, but some matrices are output. So what methods are used to solve nonlinear equations?
Can any kind-hearted person provide the code to solve the steady-state value of dsge? I will read it and learn it.
I would appreciate it very much.

See e.g. the thread