We consider the nonparametric estimation of the multivariate probability density function and its partial derivative with a support on nonnegative axis by dependent data. We use the class of kernel estimators with asymmetric gamma kernel functions. The gamma kernels are nonnegative; they may change their shape depending on the position on the semi-axis and possess good boundary properties for a wide class of densities. Asymptotic estimates of the multivariate density and of its partial derivatives such as biases, variances, and covariances are derived. The optimal bandwidth of both estimates is obtained as a minimum of the mean integrated squared error (MISE) by dependent data with a strong mixing. Optimal convergence rates of the MISE both for the density and its derivative are found.