A neutron activation analysis method using k0-standardization has been adopted for multielement analysis using our research reactors. Gold has been used as the comparator. The input parameters, used for calculation of elemental concentration, such as subcadmium to epithermal neutron flux ratio (f), the epithermal neutron flux shape factor (alpha), and the absolute efficiency of the detector (epsilon) have been determined. The values of k0-factors and Q0 (the ratio of resonance integral to thermal neutron cross section) have been taken from the literature. The validity of applying this method in our laboratory is evaluated by analyzing the elemental concentrations with respect to the certified values in eight reference materials of different origin obtained from USGS, IAEA and NIST. The percent deviation of the measured elemental concentrations are found to be within +/- 15% to that of the certified values whereas the percent relative standard deviation ranged from 2% to 10% for most of the elements analyzed.