function [ro] = r0(ri, p)
ro = mod(p(1) * ri, p(2));
