could this just be `test_numerically(order*delta, (l1-l4)*(l2-l3))`