I am new to Mathematica and some things I do not understand. Right now I need to write an algorithm which takes a polynomial in many variables, I denote it by $p_{i,j}$, and then I want to apply a differential operator on it. Something like $\sum_{i,j,k,l}p_{i+k-2,j+l}\partial p_{i,j}\partial p_{k,l}$ where $\partial p_{i,j}$ means differentiate the polynomial with respect to $p_{i,j}$. The first problem I am encountering right now is how to define a polynomial in these variables. I used Symbolize in Mathematica to define variables with subscript. The problem I have, how to go through this sum in the differential operator? I can't just put i and j in the subscript and then use a for loop. I want something like a list L which stores these variables and then if I want to differentiate with respect to $p_{i,j}$, I tell Mathematica to differentiate with respect to L[i,j].
I don't know how to do this and couldn't find something using Google.