fixed bug in computeG.m