Second order discretization of backward SDEs and simulation with the cubature method