function w = leibnitz(u,v)
% LEIBNITZ computes product rule for |D|
% w = |D|(uv)-u|D|v-v|D|u
w = mod_deriv(u.*v) - u.*mod_deriv(v) - v.*mod_deriv(u);
end