一道关于最大似然估计的题，请教stata的函数写法
变量为X1，X2，Y1，Y2，Z1，Z2
已知函数
L=ln（Φ(β1X1+β2X2-β1Y1-β2Y2)-Φ(β1X1+β2X2-β1Z1-β2Z2)）

然后用maximum likelihood 对β1，β2 进行偏微分等于0，求出β1和β2

请问各位stata高手，stata的程序应该怎么写？
谢谢各位了！