“…10 Recently, special consideration has been paid to the application of VO fractional derivatives to model several problems, such as viscoelastic constitutive models, 11 dual phase lag bioheat problems, 12 and ultra-short pulsed laser therapy. 13 Since it is often impossible to solve VO fractional problems, analytically, varied numerical techniques have been used for such problems. Some numerical procedures for solving these class of problems are including the orthonormal Bernoulli polynomials method, 14 the Chebyshev collocation method, 15 the sixth-kind Chebyshev polynomials technique, 16 the meshless approach, 17 the Lagrange polynomials method, 18 the transcendental Bernstein expansion method, 19 and the Lagrangian piecewise interpolation technique.…”