Un problema de programación cuadrática es así: minimizar $m(x)=x^tDx-2d^tx+C$ dado $Rx\geqq b$, donde $x\in\mathbb{R}^N$, $D\in\mathbb{R}^{N\times N}$, $d\in\R^N$, $R\in\R^{r\times n}$ y $b\in\R^r$.
Segundo modelo: código R call.difp <- function ( p, mx, my, nR, nC, const ) { gamma <- matrix ( data = p, nrow = nR, ncol = nC -1, byrow = T); ebeta <- exp ( gamma ) / ( 1 + apply ( exp ( gamma ), 1