NiosII is slow? Or something I am doing wrong.

I would expect this algorithm to be slow. This algorithm is performing many floating point calculations and without the floating point custom instruction all of that would end up being emulated in software using Nios II. The algorithm is also inherently slow as are most convergence algorithms and as a result you can speed it up by breaking down the algorithm and handling it in parallel (multi-core or hardware accelerators).