r= (Div(Phi u)) + ( Div(1/lambda) (Div(Phi u))) for A = A(r)