In this letter, numerical solutions of Cauchy problems are considered by multivariate Padé approximations (MPA). Multivariate Padé approximations (MPA) were applied to power series solutions of Cauchy problems that solved by using He's variational iteration method (VIM). Then, numerical results obtained by using multivariate Padé approximations were compared with the exact solutions of Cauchy problems.