def lens_equation_nfw(x,ks): res = x-4.0*ks*x*aa.func_bar(x) return res
def lens_equation_nfw(x,mus): #mus = 4.0*ks res = x-mus*x*aa.func_bar(x) return res