Area under curve trapezium method