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