
function [x, y] = mxy(f, R, parent)
    m00 = mom(f, R, parent, 0, 0);
    m01 = mom(f, R, parent, 0, 1);
    m10 = mom(f, R, parent, 1, 0);
    
    x = m10 ./ m00;
    y = m01 ./ m00;
end
