The purpose of this paper is to develop a numerical scheme for the two-dimensional fourth-order fractional subdiffusion equation with variable coefficients and delay. Using L2?1? approximation time Caputo derivative, finite difference method second-order accuracy in temporal direction achieved. novelty introduce problem under consideration coefficients, nonlinear source term, delay constant. re...