Saddle-node bifurcations of multiple homoclinic solutions in ODES