“…The use of various degrees of B-spline functions to obtain the numerical solutions of some PDEs has been shown to provide easy and simple algorithms. For instance, B-spline finite elements have been widely applied to solve elliptic equations [ 14 , 15 ], Korteweg-De Vries equation [ 16 – 18 ], Burgers’ equation [ 19 – 22 ], regularized long-wave equation [ 23 , 24 ], Fokker-Planck equation [ 25 ], advection-diffusion equation [ 26 ], and generalized equal width wave equation [ 27 ], etc., successfully. However, to the best knowledge of the authors, the B-spline FEM method is not considered for finding the approximate solution of the diffusion equation with delay.…”