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