Somewhat working. But it is unlcear why the iteration only converges for some parameters, e.g., the number of aorsa RHSs