close all
clearvars

f=@(x,y)1;
% triangle
res=quad2D(1,f,0,0,1,0,0,1)
% res=quad2D(1,f,0,0,0,1,1,0)

% rectangle
res=quad2D(2,f,0,2,2,0)

% circle
res=quad2D(3,f,0,0,2)

% quadrilateral
res=quad2D(4,f,0,0,5,0,4,2,1,2)

% NURBS 2Dplanar
% res=quad2D(5,f,CP,W,Xi);