Calculo isotermas

Enunciado: Calcule la Temperatura de ebullición del Benceno a la presión de 3 bar utilizando la construcción de Maxwell y la Ecuación de Van Der Waals.
Procedimiento: La idea principal para resolverel ejercicio es estimar una temperatura de ebullición con base en las ecuaciones de Antoine, luego con base en esta temperatura estimada y por medio de la ecuación de Van Der Waals se procederá aconstruir una línea isoterma fijando la temperatura, asignando valores de volumen y calculando la presión con la ecuación de Van der Waals.
Luego de construir esta isoterma se encontrará el punto decorte con la presión solicitada y se hallarán las áreas que se forman entre la isoterma y la constante de la presión (3bar); estas áreas deben ser muy parecidas. La resta entre estas dos áreas debe tenertendencia a cero.
Lo primero que hice fue buscar los coeficientes A, B y C para el Benceno y reemplazarlos en la ecuación de Antoine; luego los valores correspondientes son:
A =14.16
B = 2949
C =-44.56
p = 300 KPa
T = [K]
Reemplazando en:
T = B/(A-ln⁡(p))- C

T = 2949/(14.16-ln⁡(300))+ 44.56

T = 393.297 K
Ahora, desarrollando la ecuación de Van Der Waals, e igualando acero para hallar los cortes con la presión solicitada se obtiene la siguiente ecuación:

v^3- v^2*(b+(R*T)/p)+ v* a/p- (a*b)/p=0

Los valores de los coeficientes de a y b los tomé de la tabla1.6 ATKINS, Peter.

a = 18.57 (atm * 〖dm〗^6)/〖mol〗^2
b = 0.1193 ( 〖dm〗^3)/mol
R = 8.20574 * 10-2 (〖dm〗^3*atm)/(K*mol)

v^3- v^2*(0.1193+(0,0820574*393.297)/2.96077)+ v* 18.57/2.96077-(18.57*0.1193)/2.96077=0

Para esta ecuación obtengo las siguientes raíces:
V0 = 0.168491 〖dm〗^3⁄mol
V1 = 0.425276 〖dm〗^3⁄mol
V2 = 10.4424 〖dm〗^3⁄mol
Y las áreas correspondientes a estos cortes seránel área bajo la curva de la ecuación de Van Der Waals y se calculan como se muestra a continuación:

A_1= p*(v_1-v_0 )- ∫_v0^v1▒〖((R*T)/(v-b)- a/v^2 ) dv〗
A_2= ∫_v1^v2▒〖((R*T)/(v-b)- a/v^2 )...