H = [
-6, 0, 0, 0, 0, 0, 0, 0, 0, 0;
4, -6, 0, 0, 0, 0, 0, 0, 0, 0;
4, 4, -6, 0, 0, 0, 0, 0, 0, 0;
4, 4, 4, -6, 0, 0, 0, 0, 0, 0;
0, 0, 0, 0, 2, 0, 0, 0, 0, 0;
0, 0, 0, 0, 0, 2, 0, 0, 0, 0;
0, 0, 0, 0, 0, 0, 2, 0, 0, 0;
0, 0, 0, 0, 0, 0, 0, 2, 0, 0;
0, 0, 0, 0, 0, 0, 0, 0, 2, 0;
0, 0, 0, 0, 0, 0, 0, 0, 0, 2];
lb = [0; 0; 0; 0; 2; 1; 1; 1; 1; 2];
ub = [3; 3; 3; 3; 2; 3; 3; 3; 3; 2];
A_in = [
1, -1, 0, 0, -1, 0, 0, 0, 0, 0;
-1, 1, 0, 0, -1, 0, 0, 0, 0, 0;
1, 0, -1, 0, 0, -1, 0, 0, 0, 0;
-1, 0, 1, 0, 0, -1, 0, 0, 0, 0;
0, 1, -1, 0, 0, 0, -1, 0, 0, 0;
0, -1, 1, 0, 0, 0, -1, 0, 0, 0;
1, 0, 0, -1, 0, 0, 0, -1, 0, 0;
-1, 0, 0, 1, 0, 0, 0, -1, 0, 0;
0, 1, 0, -1, 0, 0, 0, 0, -1, 0;
0, -1, 0, 1, 0, 0, 0, 0, -1, 0;
0, 0, 1, -1, 0, 0, 0, 0, 0, -1;
0, 0, -1, 1, 0, 0, 0, 0, 0, -1];
[x, obj, info, lambda] = qp(ones([10, 1]), H, zeros([10, 1]), [], [], lb, ub, [], A_in, zeros([12, 1]))