Dear Mesfin, to write a code as such, you may construct a standard finite difference scheme, and enforce the condition that x-\eta and x+\gamma are nodes of the discretization, so that u(x-\eta) and u(x+\gamma) correspond to suitable node values. You only need the space step to be a submultiple of both the delays \eta and \gamma.