Can anybody please help :

a_j := [-1/2 (r-q) j Delta t + 1/2 sigma^2 j^2 Delta t] / 1 + r Delta t

b_j := [1 - sigma^2 j^2 Delta t] / 1 + r Delta t

c_j := [1/2 (r-q) j Delta t + 1/2 sigma^2 j^2 Delta t] / 1 + r Delta t

with r:= 0.05, Delta t:= T/N , Delta S:=S/M, T:=0.4, q:=0, K:=1.1, S:=2, M=N:=10

Then I defined i:=N; such that f_{i,j}:=max(K - j*Delta S, 0);

And next I defined g_{k,j}:=f_{i,j}...