clear
tic
  [nodes,elem,edges]=ReadMesh("./test4");
readingtime=toc
tic
  [A,b,n2d]=FEMEquation(nodes,elem,edges,1,0,-1,0,-1);
setuptime=toc
tic
  u=FEMSolveSym(nodes,A,b,n2d,0);
solvetime=toc
tic
  ShowSolution(nodes,elem,u)
graphtime=toc

np=25;
xy=[linspace(0,5,np);linspace(0,4,np)]';
[values,grad]=FEMValue(xy,nodes,elem,u)


#gset title
#gset nokey
#print("demo1.eps","-depsc");

