Use trapezium rule to approximate area under curve (#82)