Newton-Raphson for VLE is working, but quite slow :(
Seems there are some performance optimizations possible by removing calls to std::vector<T>.size() since sizes will not change

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>