function PII=hotelling(q)
pb=100; b=1;
	

r=0.1;



P(1)=(pb-b*q(1))*q(1);

for t=2:60
	
    
    P(t)=((pb-b*q(t))*q(t))/((1+r)^(t-1)); 

end
PII=-(sum(P));
