R语言求二元正太分布的方差和期望的矩估计

这个代码是什么啊
好难啊
期望是二阶矩阵和协方差自定