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