Frequency dependence of susceptibility higher harmonics in YBCO samples by numerical solutions of the non linear magnetic flux diffusion equation