calculate second virail coefficient