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

P
I 
if(P > D)
  Vol = -I^(1/gamma)*D/(P^(1/gamma))+D;  
else
  Vol = 0;
end


end