من هم چنان منتظر جوابم دوست عزیز
یه مشکلی که وجود داره حاصل انتگرال Z/P بر حسب P (فشار) از 0 تا 26 اتمسفر هستش، این انتگرال به صورت عددی باید حل بشه، چون مرتب کردن Z برحسب P از درجه 3 میشه و باید در هر مرحله برحسب فشار، Z را حساب کنیم تا بتونیم انتگرال بگیریم و بدست آوردن Z برحسب تابعی از P تقریبا در حد توان ما نیست

در واقع تابع Z به صورت: Z^3+A*Z^2+B*Z+C=0 قابل حل هستش که به راحتی در دما و فشار مورد نظر (دقت کنید که باید در تعداد نقاط زیادی از 0 تا 26 اتمسفر انتگرال بگیریم) حاصلش محاسبه میشه و بعد با دستور سیمپسون یا تریپزودال حاصل انتگرال بدست میاد ...
حالا نمیدونم چه مقدار واضح توضیح دادم اما بعید می دونم با دست بشه حلش کرد چون همون طور که گفتم ما تابعی به صورت : Z=f(p)
z نداریم تا داخل این انتگرال قرار بدیم:
G(R)/RT=Int[(Z-1)/P.dP]h
شاید سخت به نظر بیاد اما با مطلب واقعا حلش واسم لذت بخشه ، البته بسیار interactive میشه و برای دماها و فشارهای دیگه هم جواب میده ...
در ضمن مقاله با عنوان
خواص اضافی ترمودینامیکی با پایه معادلات حالت درجه 3 وجود داره که من درخواست دادم و گرفتمش اما متاسفانه یه حل پیچیده و تجربی (تقریبی) داره
به هر حال من منتظر شمام تا اگه می خواین بنویسم برنامشو