This paper generalizes the two dimensional Crank Nicholson difference scheme from constant coefficient to the condition of variable coefficient with interpolation by integral.