conecg=[1 2;2 3;3 4];
coordg=[0 2 4 7];
xh=-1:0.05:1;#drawing points
nhat=size(xh,2);
nel=size(conecg,1);npe=size(conecg,2);
nodg=size(coordg,2);

pp=3;
[g dg d2g] = masterp1(1,xh);
[sf dsf ddsf]=masterlag(pp,xh);

vdofs = [0 0 2 4;
         -1 0 2 3;
          3 1 2 1];

for k=1:nel
XX(:,1)=coordg(:,conecg(k,1));
XX(:,2)=coordg(:,conecg(k,2));
x=zeros(1,nhat);
for i=1:npe
x(1,:)=x(1,:)+XX(1,i)*g(i,:);
end
f=vdofs(k,:)*sf;
plot(x,f,"linewidth",2); hold on;
hold on;
end
hold off