function k=standing(p,T) b=[1799 2153 2480]; Tb=[416 491 557]; for i=1:3 f(i)=b(i)*((1/Tb(i))-(1/T)); a=1.00045*p)+((15*10^-8)*p^2);...

******************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************************