It is worth noting that, the above ı functions covers a wider support points as compared with that used by Peskin [2] and can give the second-order approximations. This can effectively reduce the non-physical oscillations (if any). For more information about the other features of this delta function one may refer to [33].