Numerical Detection and Continuation of Saddle-node Homoclinic Bifurcations of Codimension One and Two