A numerical algorithm of first-order for systems driven by correlated noises