回归方程是Y=a*X1+b*X2+c*X1*X2+d*X3+e，其中X1、X2是核心解释变量，a、b、c分别为X1、X2及两者交互项的系数，X3是其他控制变量，e是扰动项。如果X1、X2是内生变量，分别用IVX1、IVX2作为X1、X2的工具变量，那么应该怎么用stata回归？