This paper deals with the construction of numerical solution of the Black-Scholes (B-S) type equation modeling option pricing with variable yield discrete dividend payment at time $t_{d}$. Firstly the shifted delta generalized function $δ(t-t_{d})$ appearing in the B-S equation is approximated by an appropriate sequence of nice ordinary functions. Then a semidiscretization technique applied on the underlying asset is used to construct a numerical solution. The limit of this numerical solution is independent of the considered sequence of the nice type. Illustrative examples including the comparison with the exact solution recently given in [2] for the case of constant yield discrete dividend payment are presented.