Obtaining a feasible solution for underdetermined system of linear equations satisfying inequality constraints

Question

I would like to obtain a feasible solution for an under-determined system of linear equations,

$$Ax=b$$

where, $A \in \mathbb{R}^{7\times9}, \, x \in \mathbb{R}^{9\times1}\text{and } b\in\mathbb{R}^{7\times1}$.

Clearly, there is no unique solution, but I would like to obtain a solution that satisfies some non-linear inequality conditions for a few components in $x$,

\begin{align} f_1(x_1,x_2,x_3) &\le c_1 \\ f_2(x_7,x_8,x_9) &\le c_2 \\ \end{align}

To be specific, the inequalities are $|x_1 + x_2 + x_3 - k| \le c_1$ and $|x_7 + x_8 + x_9 - k| \le c_2$ , i.e. the only non-linearity here is the abs() function.

I am happy with any solution $x^*$ that satisfies the underdetermined system exactly whilst also satisfying the two constraints above.

I am currently using MATLAB as my programming environment, but pseudo-code/links to useful resources etc shall be more than enough to get started. An example MATLAB code for a dummy problem shall also be very helpful.

Answer 1

The problem was already solved in the comments section using linprog. Also CVX...

I wanted just to point out that fmincon can also be used: just set the linear equality constraints Ax=b and construct a nonlinear constraint function after the model [c,ceq,G,Geq] = nonlcon(x), where the constraints will be of the form c <= 0, ceq =0 and G and Geq are the respective gradients. This works for more general non-linear functions $f_1,f_2$ than the ones mentioned in the question. For the objective function just put a fictious one, like @(x) 0