tic
  [nodes,elem,edges]=ReadMeshTriangle("./testA.1");
readingtime=toc
tic
%  [A,b,n2d]=FEMEquation(nodes,elem,edges,'aF','fF','fF','gD','fF');
  [A,b,n2d]=FEMEquation(nodes,elem,edges,1,0,0,'gD',0);
setuptime=toc
tic
  u=FEMSolveSym(nodes,A,b,n2d,'gD');
solvetime=toc
gset nokey
tic
  ShowSolution(nodes,elem,u)
graphtime=toc

gset title
gset view 30, 30
replot
% print("hole.eps","-depsc")

uexact=gD(nodes);
maxerror=max(abs(u-uexact))

