Abstract
The aim of this paper is to deal with the pulsatile flow of blood in stenosed arteries using one of the known constitutive models that describe the viscoelasticity of blood witch is the generalized Oldroyd-B model with a variable-order fractional derivative. Numerical approximation for the axial velocity and wall shear stress were obtained by use of the implicit finite-difference scheme. The velocity profile is analyzed by graphical illustrations. This mathematical model gives more realistic results that will help medical practitioners and it has direct applications in the treatment of cardiovascular diseases.