Second update to fluid

The previous version may present some convergence issues, especially for the flux-barrier outer base points.

Now, with the proper choice of the initial guess, the solution of the nonlinear systems is more accurate.

In addition, Python code is now on par with the Matlab/Octave counterpart.