code: @Override
  	protected DoubleMatrix1D getFi(DoubleMatrix1D X){
  		double[] ret = new double[getMieq()];
  		for(int i=0; i<getDim(); i++){
  			ret[i] = -X.getQuick(i) + getLb().getQuick(i);
  			ret[getDim() + i] = X.getQuick(i) - getUb().getQuick(i) ;
  		}
  		return F1.make(ret);
  	}
docs: Inequality functions values at X.
  This is (-x+lb) for all bounded lb and (x-ub) for all bounded ub.