Add lsfit function