-->Q=eye(2,2);
-->p=[1;1];
-->[x,lagr,f]=quapro(Q,p,C,b,ci,cs,mi,x0)
f =
3.9444444
lagr =
! 0.2222222 !
! 0. !
! 1.4444444 !
x =
! 1.6666667 !
! 0.6666667 !