The author considers the problem defined by \(-{\underline \nabla}[a(p){\underline \nabla}p+\underline b(p)]+c(p)=f\) in \(\Omega\), \(p=-g\) on \(\partial \Omega\). He develops a mixed finite element method to approximate the solution, and proves that the solution exists uniquely. He obtains error estimates for approximations for p, \b{u} and \({\underline \nabla}.u\) where \b{u}\(=-a(p){\underline \nabla}p+\underline b(p)\). The treatment is highly abstract and there is no application to concrete examples.