function P = phasediff(T, scoeffs, vcoeffs);

P = polyval(scoeffs, T) - polyval(vcoeffs, T);