I am struggling to get sol1 in the attached code to simply complete and allow me to plot values of u[i,j]. Any suggestions on how I can proceed .
There is a mismatch in the equation A . U == b. The left-hand side is a list, the right-hand side is a list of lists. Here is a plot of u[6,1]:
A . U == b
u[6,1]
A . U == (b // Normal // Flatten); Sol1 = Solve[%, U] Plot[u[6, 1] /. Sol1, {i, 0, 10}]