distribution now working, need to compare to i(r)