Formamide harmonic and anharmonic frequencies of fundamental vibrations in the gas phase and in several solvents were successfully estimated in the B3LYP Kohn-Sham complete basis set limit (KS CBS). CBS results were obtained by extrapolating a power function (two-parameter formula) to the results calculated with polarization-consistent basis sets. Anharmonic corrections using the second order perturbation treatment (PT2) and hybrid B3LYP functional combined with polarization consistent pc-n (n=0, 1, 2, 3, 4) and several Pople's basis sets were analyzed for all fundamental formamide vibrational modes in the gas phase and solution. Solvent effects were modeled within a PCM method. The anharmonic frequency of diagnostic amide vibration C=O in the gas phase and the CCl(4) solution calculated with the VPT2 method was significantly closer to experimental data than the corresponding harmonic frequency. Both harmonic and anharmonic frequencies of C=O stretching mode decreased linearly with solvent polarity, expressed by relative environment permittivity (ε) ratio (ε-1)/(2ε+1). However, an unphysical behavior of solvent dependence of some low frequency anharmonic amide modes of formamide (e.g., CN stretch, NH(2) scissoring, and NH(2) in plane bend) was observed, probably due to the presence of severe anharmonicity and Fermi resonance.