Computing p-adic L-functions of totally real fields