(\frac{\partial^2 p(\theta,z)}{\partial\theta^2})_{i,j})=\frac{p_{i+1,j}-2p_{i,j}+p_{i-1,j}}{\Delta\theta^2}