Visual Basic - como resolver este problema en visual basic

Criterio de equilibrio

Z°3 ( Z al cubo)- Z° 2 ( Z al cuadrado) (B+1) + AZ – AB =0
Al resolver esta ecuación, vamos a tener
tres valores para Z, el Z mas pequeño
es el del líquido y el Z mayor es el del vapor.

Pr= Presión Reducida
Tr= Temperatura Reducida
Pc= Presión Critica
Tc= Temperatura Critica

A= 27 Pr
__ __
64 Tr° 2

B= 1 Pr
__ ___
8 Tr

Pr= P
__
Pc

Tr = T
__
Tc
A
___
LnO =Z -1-LN(Z-B)- Z

ALGORITMO:
1. Proponer una P variamos P para cada iteración
2. Obtener las raíces de la ecuación 1. La mayor es Zv y la menor es ZL
3. Sustituir en la ec2. Determinar OL y Ov
4. Si OL = Ov OK, o por lo menos OL – Ov <=1*10 -6
5. Calcular H-Hi para el liquido ty para el vapor, con la ecuación 3
6. Si OL = Ov Proponer otra P, hasta que OL = Ov