相当不错                                                                        
                                       
                                                     
                                                     
                                             
                                                                 
                                                                                 

刚才又自己看了一些，突然有点明白了。是不是因为J（）函数本身必须在程序当中使用，即需要先定义一个函数。通过这个思路我使用了noteproblem的函数进行调试，发现可以通过，但是不能显示矩阵。刚在看struct和pointer似乎有些关联，还在继续阅读。希望大虾们能指点一下。


附有改的程序，不知是否有错求教：
mata
function noteproblem(matrix pvec)/*J函数只能作为程序里面的*/
{
pragma unset pvec
pvec  = J(1, n, NULL)
for (i=1; i<=n; i++) pvec[i] =&(J(5, 5, .))
}
end