thanks
I am solving two equations using FDM
I used nested for loops for this
[ For[iter = 1; iter < 100, iter++,
For[i = nt; i >= 1, i--,
For[j = 2; j <= nx, j++,
q[[i - 1,
j]] = Beta*Rho[[i, j]] + (1 + Alpha*(1/(nt - 1)) -
v*(nx - 1)/(nt - 1))*q[[i, j]] +
v*(nx - 1)/(nt - 1)*q[[i, j - 1]];
Rho[[i - 1,
j]] = (1 - Beta)*Rho[[i, j]] + Alpha*(1/(nt - 1))*
q[[i, j]];
]
]
]