We present a new equation of state for infinite systems (symmetric, asymmetric and neutron matter) based on an extended Skyrme functional constrained by microscopic Brueckner-Bethe-Goldstone results. The resulting equation of state reproduces with very good accuracy the main features of microscopic calculations and it is compatible with recent measurements of two times Solar-mass neutron stars. We provide all necessary analytical expressions to facilitate a quick numerical implementation of quantities of astrophysical interest.