Newton-Raphson equilibrium iteration scheme

Added file containing N-R iteration scheme. Currently, this does not
give the expected results. There seems to be something wrong either
with
a) The iteration itself
b) The Hessian and/or other thermodynamic calculations