We investigate a tansient diffusion approximation of queue size distribution in $M^{X}/G^{Y}/c$ system using the diffusion process with elementary return boundary. We choose an appropriate diffusion process which approxiamtes the queue size in the system and derive the transient solution of Kolmogorov forward equation of the diffusion process. We derive an approximation formula for the transient queue size distribution and mean queue size, and then obtain the stationary solution from the transient solution. Accuracy evalution is presented by comparing approximation results for the mean queue size with the exact results or simulation results numerically.