Dear all, I am locking for a quadratic solver in R beeing able to solve the following sort of problem: minimize || y - Xmatrix u ||^2 subject to: Amatrix u <= Bbound and u >= lowbound While the functions pcls{mgcv} and solve.QP{quadprog} can handle the first part of the constraint it seems that they are not designed to deal with the second part (or I do not see how to specify this second constraint). Is there an existing function for this sort of problem in R? Regards, Sven Sven Lautenbach UFZ Centre for Environmental Research in the Helmholtz Association Department for Computational Landscape Ecology Permosterstr. 15 D-04318 Leipzig Germany [[alternative HTML version deleted]]