“…a + be the generalized VO Riemann-Liouville integral operator of order 𝛼(t), where n − 1 < 𝛼(t) ≤ n, n ∈ IN, and g, 𝑓 ∈ C[0, T]. To simulate the VE-FIE (10), over the interval [0, T], we assume that the function g is considered such that a unique solution to the VE-FIE (10) exits on [0, T]. To describe the P-C approach, let us define a uniform grid in the interval [0, T] with N + 1 equispaced nodes t 𝑗 = 𝑗h, 𝑗 = 0, 1, • • • N, where h = T∕N and N ∈ IN.…”