. clear . set obs 625 . gen x = (mod((_n-mod(_n-1,25)-1)/25,25)-12)/12 . gen y = (mod(_n-1,25)-12)/12 . gen z = normalden(x,0,2)*normalden(y,0,2) //二维乘积核 . twoway contour z x y, level(15)