
function m = momi7(f, R, parent)
    n02 = momcn(f, R, parent, 0, 2);
    n20 = momcn(f, R, parent, 2, 0);
    n11 = momcn(f, R, parent, 1, 1);
    n12 = momcn(f, R, parent, 1, 2);
    n21 = momcn(f, R, parent, 2, 1);
    n03 = momcn(f, R, parent, 0, 3);
    n30 = momcn(f, R, parent, 3, 0);
    
    m = (3 .* n21 .- n03) .* (n30 .+ n12) .* ((n30 .+ n12) .^ 2 .- 3 .* (n12 .+ n03) .^ 2) .+ (3 .* n12 .- n30) .* (n21 .+ n03) .* (3 .* (n30 .+ n12) .^ 2 .- (n21 .+ n03) .^ 2);
end
