In this article we consider calculation of discrete Laplacian for 1D arrays using python/numpy/scipy. We introduce the sparse matrix technique that is rather efficient for computation of certain types of diffusion PDE.

continue reading.....