/* wxMaxima 0.7.1 http://wxmaxima.sourceforge.net
Maxima 5.11.0 http://maxima.sourceforge.net
Using Lisp GNU Common Lisp (GCL) GCL 2.6.8 (aka GCL)
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
This is a development version of Maxima. The function bug_report()
provides bug reporting information.

(%i1) R:8.314;

(%i2) R;

(%i3) T:500;

(%i4) a:0.364;

(%i5) b:0.429E-6;

(%i6) algsys([R*T/(V-b) - a/V^2 = 10^6],[V]);

(%i7) B0: 104.76E-6;

(%i8) A0:507.31E-3;

(%i9) C:660;

(%i10) Beta:R*T*(B0 - (A0/R*T) -(C/T^3));

(%i11) Beta:R*T*(B0 - (A0/R*T) -(C/(T^3)));

(%i12) Beta:R*T*(B0 - (A0/(R*T)) -(C/(T^3)));

(%i13) b1:72.35E-6;

(%i14) a1:71.32E-6;

(%i15) gamma:R*T*(-B0*b1 + (A0*a1/(R*T)) - B0*C/T^3);

(%i16) gammap:(1/(R*T))*((gamma/(R*T)) - (B/(R*T))^2);

(%i17) gammap:(1/(R*T))*((gamma/(R*T)) - (Beta/(R*T))^2);

(%i18) delta:R*T*(B0*b1*C/T^3);

(%i19) deltap:(1/(R*T)^2)*((delta/(R*T)) - (3*Beta*gamma/(R*T)^2) + 2*(Beta/(R*T))^3);

(%i20) P:1013250;

(%i21) V: R*T/P + Beta/(R*T) + gammap*P + deltap*P^2;

(%i22)