Vibrational self-consistent field approach to anharmonic spectroscopy of molecules in solids: Application to iodine in argon matrix