clear; close all; a = 3.59; b= 0.0427; R = 0.08206; T = 300; v0 = 0.05; vf = 2; v = linspace(v0,vf,1000);%m3 p = R*T./(v-b) - a./v.^2; p_ideal = R*T./v; semilogx(v,p,'b'); hold on; semilogx(v,p_ideal,'r--'); ylim([60 80]); legend('Van der Waals','Ideal gas');
全体的な違い
極値付近
さらに拡大