用一个嵌套循环即可：

mat a = J(4,4, 9)
mat b = J(4,4, 2)
mat c = J(4,4, 0)

local col = colsof(a)
local row = rowsof(a)

forvalues i = 1/`col'{
  forvalues j = 1/`row'{
    mat c[`i',`j'] = a[`i',`j'] / b[`i',`j']
  }
}

mat list c

* ========= 结果 ===============

. mat list a

symmetric a[4,4]
    c1  c2  c3  c4
r1   9
r2   9   9
r3   9   9   9
r4   9   9   9   9

. mat list b

symmetric b[4,4]
    c1  c2  c3  c4
r1   2
r2   2   2
r3   2   2   2
r4   2   2   2   2

. mat list c

symmetric c[4,4]
     c1   c2   c3   c4
r1  4.5
r2  4.5  4.5
r3  4.5  4.5  4.5
r4  4.5  4.5  4.5  4.5

A=[aij],B=[bij],cij=aij/bij,C=[cij]