Numerical issues with Symbolic Derivative Operators
In order to apply the classical Van der Walls Critical Isotherm Conditions to a Carnahan-Staling based EOS, I make use of the derivative operators within a solve block. I found discrepancy between the derivative evaluated by using the symbolic derivative and the expanded derivative ( calculated symbolically ). Furthermore the solve block does not converge with the symbolic derivatives, while it does ( "better" ) with the expanded derivative.

