function coef=planeFitViaTls(data)
% planeFitViaTls: Plane fit via total lease-squares