
function m = momi5(f, R, parent)
    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 = (n30 .- 3 .* n12) .* (n30 .+ n12) .* ((n30 .+ n12) .^ 2 .- 3 .* (n21 .+ n03) .^ 2) .+ (3 .* n21 .- n03) .* (n21 .+ n03) .* (3 .* (n30 .+ n12) .^ 2 .- (n21 .+ n03) .^ 2);
end
