function Portfolio()
    H=[146 187 145; 187 854 104; 145 104 289]
    A=-[0.175 0.24 0.227]
    b=-0.2
    LB=zeros(3,1)
    UB=ones(3,1)
    Aeq=[1 1 1]
    beq=1

    x=quadprog(H,[0; 0; 0],A,b,Aeq,beq,LB,UB)
    
    options = optimset('LargeScale','off');
    x=quadprog(H,[0; 0; 0],A,b,Aeq,beq,LB,UB,zeros(3,1),options)
