哪位大神能帮我看一下，或者在stata mata运行一下这个程序，我不造哪里写错了。。少女笨瞎

function funegarch( ut, Ezt, n, aa, bb, ccc, dd)
{
zt0=0

sigmat0=1

gzt0=(ccc+dd)*zt0-dd*Ezt

sigmat(1,1)=sqrt(exp(aa+gzt0+bb*log((sigmat0^2))))

zt(1,1)=ut(1,1)/sigmat(1,1)
for (i=1; i<=n-1; i++)
{
         if zt(i,1)<0
         {
          gzt=(ccc-dd)*zt(i,1)-dd*Ezt
                                   }
         else
          {
           gzt=(ccc+dd)*zt(i,1)-dd*Ezt
                                    }
                }
return (ut(i+1,1)/sqrt(exp(aa+gzt+bb*log((sigmat(i,1)^2)))))
                                                                            
}