>> xlist = [ 0.0, 7.0, 9.0, 8.0 ];
>> ylist = [ 6.0, 1.0, 3.0, 5.0 ];
>> area = quadrilateral_area ( xlist, ylist )

area =

   19.5000