conec=[1 2;2 3;3 4;4 5;5 6];
coord=[-10 -5 -1 0 4 10];
nel=size(conec,1);
npe=size(conec,2);
nod=size(coord,2);

xhat=-0.9:0.3:0.9;
nhat=size(xhat,2);
zhat=zeros(nhat,1);

figure 1;

[g dg d2g] = masterp1(1,xhat);

for k=1:nel

a=coord(conec(k,1));b=coord(conec(k,2));

x=a*g(1,:)+b*g(2,:);
c=k;

scatter(zhat,x,20,c,"filled")
hold on;

end

hold off