A symbolic algorithm to compute immersions of polynomial systems into linear ones up to an output injection