A numerical algorithm for modelling of boson-fermion stars in dilatonic gravity