Computation of a canonical form for neutral delay-differential systems