V=[
0.82 0.23 0.155 0.013 0.314
0.23 0.484 0.346 0.197 0.592
0.155 0.346 0.298 0.143 0.419
0.013 0.197 0.143 0.172 0.362
0.314 0.592 0.419 0.362 0.916
];
r = [ 1.78 0.37 0.237 0.315 0.49 ];
n=length(r);
lambda=.5;
H=V;f=zeros(n,1);
A=-eye(n);b=zeros(n,1);
Aeq=[r
     ones(1,n)];
beq=[lambda
      1];
x=quadprog(H,f,A,b,Aeq,beq)
