Integrate[h^2*i*j*PDF[MultinormalDistribution[{u1,u2,u2},{{s1^2,s12,s13},{s12,s2^2,s23},{s13,s23,s3^2}}],{h,i,j}],{h,-Infinity,Infinity},{i,-Infinity,Infinity},{j,-Infinity,Infinity}]