function lj=lennard_jones_ab (a,b)
  #returns potential LJ with parameters a,b, i.e. a/r^12-b/r^6
  lj=@(x) a*power(x,-12)-b*power(x,-6);
end

