Should get back later for whole phase diagrams

New module to realize the calculation of Px envelopes
via the continuation method.

Right now works with 2phase lines, with a starting point
of three phase lines calculation.

Added a new namelist to input files.

Now the EoS come from an external library, yaeos.

The calls are still direct to the legacy code, so an enventual
fix should occur. But first an IO must be defined on the yaeos side,
to include any kind of model from an input file.

Now the main program gets part of
it's settings from the cli library `flap`

Added the calculation of three-phase Px lines.

For now it's a beta-stage, where only a three-phase region is calculated
based on finding an intersection.
This verifies that the three-phase lines
model works properly, and later on more ways of initialization will be implemented.

Modified the tolerances to assure the same intersection is not repeated.

This probably should be revisited with more tests later.

This should be part of the `yaeos` set later

Don't write if passed max of iterations