function truss_opt_class x0 = zeros(2,1); options = optimset('Display','iter','MaxIter',100,'MaxFunEvals',500,'LargeScale','off'); H = 2*[2.1650625e5 0; 0 8.66025e5]; f = [-1000 -500]'; A = []; b = []; Aeq = []; beq = []; lb = []; ub = []; x = quadprog(H,f,A,b,Aeq,beq,lb,ub,x0,options)