SED FITTING WITH MARKOV CHAIN MONTE CARLO: METHODOLOGY AND APPLICATION TO Z=3.1 LYα-EMITTING GALAXIES