eq = a*f(x,y) + b*f(x,y).diff(x) + c*f(x,y).diff(y) - G(x,y)
    pdsolve(eq)