ARMA model with fmincon

Hi,

I want to estimate an ARMA model, but not with the usual function. I need to optimize it myself. I want to use the function fmincon for this, but how can I make sure that the process is stationary?

What constraints or something do I need? And note that I have more than 1 parameter. For example:

Thanks in advance!