Liouvillian function