Question: how to find values on internal nodes of mesh

hi friends..!!
I am working on 1d heat equation with finite difference method and i have made rough code for it in which i considered mesh size of 4 in x direction as wel as in time space. By initial condition am known with points (0,0) (1,0) (2,0) (3,0) and (4,0). AS boundary conditions are neumman as u can see in code below i have attached. Now my problem is loops in end of program is not picking internal nodes like (1,1), (1,2), (1,3),(1,4), (2,1),...,(4,4). So please let me know what should i add so that i should have an output on all nodes. Can anyone can help me.. thankscode.mw 

This has been branched into the following page(s):
Please Wait...