Under investigation in this paper are two local fractional partial differential systems, one is the homogeneous linear system with initial values, and other inhomogeneous non-linear boundary values. To solve these we employ variational iteration method obtain exact solutions. It shown that provides an effective mathematical tool for solving systems