Question: Expanding a derivative onto an expression using Maple 2020

Hello I am trying to compute the following expression;

(x^2r+y^2+z^2) gamma^{0} partial_{x} (-xy(1-r)) gamma^{2} partial_{y} \Psi

in which r = sqrt(1-K(x^2+y^2+z^2)), K is a constant, the gamma's are the Dirac gamma matrices and \Psi is the wave function. Essentially I would like to have partial_{x} operate on the expression to the right of it and apply product rule where there will be three terms with the last one being partia_{x} partial_{y} \Psi where it evaluates on everything expect the partial and the wavefunction, combining into the term partial_{x} partial_{y} \Psi, and if the partials are the same they combine into a single second order derivative with respect to the given variable. I have done this calculation by hand but seeing that there is an extensive amount of terms I would like to automate it using Maple. I have tried using the D() calling sequence but it doesn't quite give me what I am looking for. Is there anything documentation which will give me what I am looking for?