Summary: In this paper, the canonical correlation analysis for multivariate functional data is considered. The analysis is based on the basis functions representation of the data. The use of non-orthogonal bases is available in contrast to the approach given in the literature. The robust estimation methods of the covariance matrix are also studied in the multivariate functional canonical correlation analysis. Simulation studies and breakdown analysis suggest that the proposed methods may perform better than the classical estimator under non-normal models and in the presence of outlying observations.