Preliminary results in unquenched Numerical Stochastic Perturbation Theory

Introducing fermionic loops contributions in Numerical Stochastic Perturbation Theory was mainly motivated by the proposal to compute 2-3 loops for renormalization constants (and improvement coefficients). This is feasible because the computational overhead of unquenching NSPT is by far lower than for non perturbative simulations. We report on first, preliminary results for the quark propagator (basically the third loop for the critical mass) and discuss the status of the computation of other quantities.