update

bug fix 1/r equals 0