A new approach for obtaining cosmological constraints from Type Ia Supernovae using Approximate Bayesian Computation