Numerical scheme for solving singular fractional partial integro‐differential equation via orthonormal Bernoulli polynomials