For a starter indent the print ui
to see the values of ui
for all N1
. Then append your results to a list:
res = []
for n in range(0, N1):
for i in range(0,1):
u[i] = 2*Coefficient*(ui[i+1]+bi*toutside)+(1-2*Coefficient-2*bi*Coefficient)*ui[i]
for i in range(1,Nx):
u[i] = Coefficient*(ui[i+1]+ui[i-1])+(1-2*Coefficient)*ui[i]
for i in range(Nx,Nx+1):
u[i] = 2*Coefficient*(ui[i-1])+(1-2*Coefficient)*ui[i]
ui[:]= u #updates matrix for next loop
print ui
res.append(ui.copy())
print res
Produces this result:
[array([ 41.5625, 50. , 50. , 50. ]),
array([ 37.87109375, 48.41796875, 50. , 50. ]),
array([ 35.6628418 , 46.73706055, 49.70336914, 50. ]),
array([ 34.06639099, 45.21682739, 49.20280457, 49.88876343]),
array([ 32.79785633, 43.87349129, 48.58405113, 49.63152885])]