---

This study was carried out with the aim of covalent immobilization of Aspergillus oryzae beta-galactosidase and Bacillus licheniformis protease on multi-walled amino-carbon nanotubes. In this method, fractional 2k design was used to study the effect of seven continuous factors (activation pH, glutaraldehyde molarity, activation time, buffer solution pH, buffer solution molarity, MWCNT-NH3-glutaraldehyde amount and stabilization time) on the stabilization efficiency and enzyme activity. . Design-expert software was used to analyze data and draw graphs. The results showed that the aforementioned factors predict the level of enzyme activity of Bacillus licheniformis protease and Aspergillus oryzae beta-galactosidase with correlation coefficients of 0.80 and 0.92 at the rate of 77 and 88%, respectively. Also, the correlation coefficient of the covalent fixation efficiency model of Aspergillus oryzae beta-galactosidase and Bacillus licheniformis protease on multi-walled carbon nanotubes was 0.89 and 0.82, respectively, and the studied factors were able to determine the covalent fixation beta efficiency, respectively. Aspergillus oryzae galactosidase and Bacillus licheniformis protease on multi-walled amino-carbon nanotubes predict 83 and 77%, respectively.