This article introduces a new testing procedure for multivariate normality, employing a double estimation approach in the context of a characterizing partial differential equation (PDE). We provide a detailed theoretical framework and perform simulations to validate our method’s performance. The results indicate strong robustness and power in detecting deviations from multivariate normality across various distributions.