new formula: quadprog