The distributions of X-ray diffraction (XRD) of N2+ (180 keV) implanted Si are presented. On the basis of trial and error strain function and multilayer model, according to kinematic theory of XRD, with Levenberg-Marqardt optimization method to simulate experimental curve we computed the strain distribution as the functions of depth, dose and annealing temperature. The results are analysed and discussed.