Fixed-point iteration

[michakraus]

We should add a simple fixed point solver. This is possible when assuming that F in the nonlinear system is of the form $F(x) = x - f(x)$. Then one can just iterate $x_{k+1} = F(x_k) - x_k$ where $x_k$ is the solution of the current iteration.