Rust (1987) --  fixed value function iteration bug

* Use numpy.copy instead of assigning references
* transposition of u is done inside the function