Approximating the leading singular triplets of a large matrix function