_RL k123p
k123p = k12p*k3plocal
a = x3 + k1plocal*x2 + k12p*x + k123p
& 2.0*pt*k123p/a +
& 2.0*pt*k123p*da/a2 -