On the computation of Belyi functions: A numerical approach