Recent content by qazinasir

Basically in this code I am trying to assume a pressure 'Pn' to calculate the results which satisfy this relation (Error<10^-4) if this satisfy's then our assume pressure is right else if not then we use the relation P(n+1) = Pn*(Phil/Phiv) to caluclate another pressure and update the above...

% Calculate the vapor Pressure of n-heptane using SRK-EOS % P = RT / V - b - alpha*a / V(V-b) R = 8.134e-3; % MPa m3 / mole.K n = 0; n = n+1; for i=1:5 Pn = 1.6714; %(Assume Pressure)> this needs updating Tr = 0.9; Tc = 540.13; % K Pc = 2.74; % Mpa; T = Tr*Tc; w = 0.350; m =...