Rotation diagrams measured for the CO and isotopolog lines at the source positions of IRAS 4A and 4B. The column density of each observation in each rotational level divided by the statistical weight is plotted against the excitation energy of the level. The fitted line shows the Boltzmann distribution of the rotational populations. Derived values of the rotation temperatures are presented in Table 4.