Algorithm of ARXMAX model

Hi,

I want to know algorithm of ARXMAX model when I use matlab function.

For example, like this.

sys = armax(data,[na nb nc nk])

https://kr.mathworks.com/help/ident/ref/armax.html

I want to understand the exact algorithm through the formula.

Can you give me any advise for this question?