On a generalized fractional boundary value problem based on the thermostat model and its numerical solutions via Bernstein polynomials