function Vol = Vol(P,I,D,gamma)

Vol = zeros(1,length(P));
for i = 1:length(P)
  if(P(i) > I)
    Vol(i) = D-I^(1/gamma)*D/(P(i)^(1/gamma));  
  else
    Vol(i) = 0;
  end
end

end
