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);