A numerical variational method for calculating vibration intervals of bent triatomic molecules