bug fixes 

1) rho function i =0 j =0
2) function k <5 instead of k <64