用全多重网格法求解线性方程,这是M文件functionc=MG(MK,z,g)%MK为刚度矩阵构成的向量functionc=FMG(MK,MF)