In this paper we develop a numerical method for valuing multivariate European contingent claims whose payoffs depend on more than one log-normal stochastic variables. This is achieved by means of Gauss-Hermite Integrations, applied to the principal component analysis of correlation matrix to improve our computation speed. Several examples are discussed to compare with other methods.
