Using vectorial laser equation, the polarization feature of a relaxation oscillation laser is numerically analyzed. The numerical results reveal that the light field in the relaxation oscillation process is linearly polarized and the polarization direction varies with time. The single transverse-mode neodymium-doped fiber laser shows this feature.