function idum=plotfem1(conecg,coordg,pp,conecdof,vdoft)
xh=-1:0.05:1;#drawing points
nhat=size(xh,2);
nel=size(conecg,1);npe=size(conecg,2);
nodg=size(coordg,2);
[g dg d2g] = masterlag(1,xh);
[sf dsf ddsf]=masterlag(pp,xh);
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=vdoft(conecdof(k,:))'*sf;
plot(x,f,"r","linewidth",2); hold on;
hold on;
idum=0;
end