clear

q0 = zeros(1,60);       % Make a starting guess at the solution
lb = 0.0001*ones(1,60);        % Set lower bounds
ub = 200*ones(1,60);          % upper bounds

[q,NPV] = fmincon(@hotelling,q0,[],[],[],[],lb,ub,@rajhot,[])

figure(1)
plot(q)






