Create computeU.m