We propose multiwavelength in-line digital holography with wavelength-multiplexed phase-shifted holograms and arbitrary symmetric phase shifts. We use phase-shifting interferometry selectively extracting wavelength information to reconstruct multiwavelength object waves separately from wavelength-multiplexed monochromatic images. The proposed technique obtains systems of equations for real and imaginary parts of multiwavelength object waves from the holograms by introducing arbitrary symmetric phase shifts. Then, the technique derives each complex amplitude distribution of each object wave selectively and analytically by solving the two systems of equations. We formulate the algorithm in the case of an arbitrary number of wavelengths and confirm its validity numerically and experimentally in the cases where the number of wavelengths is two and three.