A new method for solving variable coefficients fractional differential equations based on a hybrid of Bernoulli polynomials and block pulse functions. (21st March 2021)