plot = plt.imshow(matrix,cmap= maps, extent= [0.0,1.0,0.0,1.0],interpolation='gaussian', aspect='auto') plt.savefig(results_dir+names+"_.png",dpi=600,bbox_inches ...
IM = int(L/dx) +1 # Grid points in x-direction alpha = 2.6*10**(-6) # diffusivity dt = 0.5*dx**2/alpha # Time step NT = round(TT/dt) +1 # Grid points for Time b ...