4
Similar to method in example 8.2.8, page 128, reference [1]