conec1=[1 2;2 3;3 4;4 5;5 6];
coord1=[-10 -5 -1 0 0 0;
	   0 0 0 0 4 10];
conec2=[1 2;2 3;3 4;4 5;5 6;6 7;7 8;8 1];
coord2=[-10 -5 -2 -2 -2 -5 -10 -10;
        3 3 3 6 9 9 9 6];
nel1=size(conec1,1);npe=size(conec1,2);
nod1=size(coord1,2);
nel2=size(conec2,1);npe=size(conec2,2);
nod2=size(coord2,2);

xhat=-0.9:0.3:0.9;nhat=size(xhat,2);
[g dg d2g] = masterp1(1,xhat);

figure 1;

for k=1:nel1
a=coord1(:,conec1(k,1));b=coord1(:,conec1(k,2));
x(1,:)=a(1)*g(1,:)+b(1)*g(2,:);
x(2,:)=a(2)*g(1,:)+b(2)*g(2,:); c=k;
scatter(x(1,:),x(2,:),20,c,"filled")
hold on;
end
for k=1:nel2
a=coord2(:,conec2(k,1));b=coord2(:,conec2(k,2));
x(1,:)=a(1)*g(1,:)+b(1)*g(2,:);
x(2,:)=a(2)*g(1,:)+b(2)*g(2,:); c=k;
scatter(x(1,:),x(2,:),20,c,"filled")
hold on;
end

hold off